MULTIMODEL BAYESIAN ANALYSIS OF DATA-WORTH APPLIED TO UNSATURATED FRACTURED TUFFS
|
|
- Collin Gray
- 5 years ago
- Views:
Transcription
1 MULTIMOEL BAYESIAN ANALYSIS OF ATA-WORTH APPLIE TO UNSATURATE FRATURE TUFFS M. Ye. Lu S.P. Neuman 2 and L. Xue 2 epartment of Scientific omputing Florida State University Tallahassee Florida USA s: mye@fsu.edu dl7f@fsu.edu 2 epartment of Hydrology and Water Resources University of Arizona Tucson Arizona USA s: neuman@hwr.arizona.edu xueliang983@gmail.com Abstract The management of contaminated groundwater systems requires an understanding of their response to alternative remediation strategies. Such understanding requires the collection of suitable data to help characterize the system and monitor its response to existing and future cleanup and/or containment options. It also requires incorporating such data in suitable models of water flow and contaminant transport. As the acquisition of subsurface characterization and monitoring data is costly it is imperative that the design of corresponding data collection schemes be cost-effective i.e. that the expected benefit of new information exceed its cost. A major benefit of new data is its potential to help improve one's understanding of the system in large part through a reduction in model predictive uncertainty. Traditionally value-of-information or data-worth analyses have relied on a single conceptualmathematical model of site hydrology. Yet there is a growing recognition that analyses and predictions based on a single hydrologic concept are prone to statistical bias and underestimation of uncertainty. This has led to a recent emphasis on conducting hydrologic analyses and rendering corresponding predictions by means of multiple models. Here we adopt a multimodel approach to optimum value-of-information or data-worth analyses based on model averaging within a Bayesian framewor (Neuman et al. 2). The Bayesian model averaging (BMA) approach is compatible with both deterministic and stochastic models in its maximum lielihood version (BMA) (Neuman 23) it is additionally consistent with current statistical methods of hydrologic model calibration. Implementation entails either Monte arlo simulation or lead order approximations. We implement BMA computationally on pneumatic permeability data from an unsaturated fractured tuff site near Superior Arizona USA.. INTROUTION Hydrogeologists face a daunting challenge helping to insure that contaminants in the subsurface do not pose unacceptable future riss to humans and the environment. To quantify and manage these riss one must understand their relationships to alternative pollution remediation schemes. Such understanding requires the collection of suitable data to help characterize the subsurface and monitor its response to existing and future scenarios. It also requires incorporating such data in suitable models of subsurface flow and contaminant transport. As noted by Bac (27) three strategies have traditionally been used to determine the magnitude of a data collection effort: minimizing cost for a specific level of accuracy or precision minimizing uncertainty for a given budget or responding to regulatory demands on data quantity and quality. Various combinations of these strategies have also been described such as a fitness-for-purpose approach [2]. Many today prefer a fourth approach based on value-of-information or data-worth analysis. Here the decision to collect additional data or the design of a data collection program is based on costeffectiveness. A program is considered cost-effective if the expected benefit from the new information exceeds the cost. A major benefit of new data is its potential to help improve one's understanding of the system in large part through a reduction in model predictive uncertainty. This benefit however is worth
2 the cost only if it has the potential to impact decisions concerning management of the water resource system. Value-of-information or data-worth analyses incorporating statistical decision theory have been applied to various water-related problems in the 97s and to groundwater problems in the 98s. More recent applications to groundwater resource and contamination issues have been reported among others in (Russell and Rabideau 2 ausman et al. 2). James and Freeze (993) proposed a Bayesian decision-maing framewor to evaluate the worth of data in the context of contaminated groundwater that has been widely cited in the subsequent literature. A comprehensive review focusing on health ris assessment can be found in (Yoota and Thompson 24). Additional recent publications of relevance include (Feyen and Gorelic 25 Norberg and Rosen 26). A major limitation of many existing approaches is that they rely on a single conceptual-mathematical model of the subsurface and flow/transport processes therein. Yet the subsurface is open and complex rendering its characteristics and corresponding processes prone to multiple interpretations and mathematical descriptions including parameterizations. This is true regardless of the quantity and quality of available data. Predictions and analyses of uncertainty based on a single hydrologic concept are prone to statistical bias (by committing a Type II error through reliance on an inadequate model) and underestimation of uncertainty (by committing a Type I error through under sampling of the relevant model space). Analyses of environmental data-worth which explore how different sets of conditioning data impact the predictive uncertainty of multiple models in a Bayesian context include (Freer et al. 996 Rojas et al. 2) whereas Freer et al. (996) employ Generalized Lielihood Uncertainty Estimation (GLUE) (Beven and Freer 2) Rojas et al. (2) combine GLUE with Bayesian Model Averaging (BMA) (Hoeting et al. 999). Nowa et al. (2) introduce a Bayesian approach to data worth analysis when flow and transport tae place in a random log hydraulic conductivity field. Whereas flow and transport are described by a single (linearized) model each having nown parameters other than those describing spatial variations in log hydraulic conductivity the latter is characterized by a single drift model and a continuous family of variogram models having uncertain parameters. In a similar spirit we adopt in this paper a multimodel approach to optimum value-of-information or dataworth analyses that is based on model averaging within a Bayesian framewor (Neuman et al. 2). Our approach is general in that it considers multiple models of any ind all having uncertain parameters whereas parameterizing models in the manner of (Nowa et al. 2) is elegant and computationally efficient it is unfortunately limited to a narrow range of variogram models and does not generally apply to other models such as those of flow and transport. We prefer BMA over GLUE because it (a) rests on rigorous statistical theory (b) is compatible with deterministic as well as stochastic models and (c) in its maximum lielihood () version (BMA) (Neuman 23) is consistent with current methods of hydrologic model calibration. Whereas BMA (lie the closely related approach in Nowa et al. 2) relies heavily on prior parameter statistics BMA can do without such statistics or otherwise update them on the basis of potential new data both before (during the so-called preposterior stage) and after (during the posterior stage) they are collected. We describe BMA and BMA theoretically outline the approach of Neuman et al. (2) to implement BMA computationally and illustrate the latter on pneumatic permeability data from an unsaturated fractured tuff site near Superior Arizona USA. Our methodology should be of help in designing the collection of hydrologic characterization and monitoring data in a cost-effective manner by maximizing their benefit under given cost constraints. The benefit would accrue from optimum gain in information or reduction in predictive uncertainty upon considering jointly not only traditional sources of uncertainty such as those affecting model parameters and the reliability of data but also most importantly lac of certainty about the conceptual-mathematical models that underlie the analysis and the scenarios under which the system would operate in the future. The methodology should apply to a broad range of models representing natural processes in ubiquitously open and complex earth and environmental systems.
3 2. BAKGROUN Bayesian ecision Analysis Framewor One way to cast the data-worth issue would be within a Bayesian ris-cost-benefit decision framewor such as that of Freeze et al. (99 992). Suppose without loss of generality that the data are intended to help one decide whether or not a contaminated site should be remediated. The problem can be cast in terms of a decision tree (Bac 27) containing objective functions defined for each decision alternative i 2 as i B i i Pf i i where B i is the benefit and i the investment cost ris being expressed as the product Pf i i of a ris aversion factor the probability of failure P i and the cost of failure f i. ollecting additional information generally causes the ris term to decrease due to a decrease in the probability of failure. The corresponding increase in i is the expected value (worth) of the new data. The final outcome of the analysis depends on the choice of decision rule one adopts for example maximizing would result in the largest benefit and lowest cost. i A reduction in the probability of failure comes about through a reduction in uncertainty about the expected system state present or future. The impact of hydrologic data on this expectation and the associated uncertainty are often evaluated by means of a hydrologic model. ommonly the model is considered to be certain while its parameters (and in some cases its forcing terms such as sources and boundary conditions) are treated as being uncertain due to insufficient and error-prone data. As already noted we now of only one wor that considers the impact of data on model predictive uncertainty within a Bayesian framewor by considering the model itself to be uncertain (Freer et al. 996) and one other wor that parameterizes this uncertainty (Nowa et al. 2). Below we provide bacground about Bayesian model averaging and its maximum lielihood version which we employ for this same purpose. Bayesian Model Averaging (BMA) onsider a random vector Δ the multivariate statistics of which are to be predicted with a set M of K mutually independent models M each characterized by a vector of parameters θ conditional on a discrete set of data (the case of correlated models has recently been considered in Sain and Furrer 2). In analogy to the case of a scalar (Hoeting et al. 999) we write the joint posterior (conditional) K distribution of Δ as pδ pδ MpM i.e. as the average over all models of the joint M posterior distributions pδ associated with individual models weighted by the model posterior probabilities pm. These weights are given by Bayes rule in the form K pm p MpM/ p MlpMl where l integrated lielihood of model M p M parameters pθ M the prior density of θ under model M and p M p M θ p θ M dθ is the θ being the joint lielihood of this model and its p M the prior probability of M. All probabilities are implicitly conditional on the choice of models entering into the set M. The posterior mean and covariance of Δ are given by EΔ EΔ MpM and ovδ K Δ Δ. We also consider the trace Tr ov E ov M ov E M M M scalar measure of the posterior variance of Δ. Δ which provides a
4 Maximum Lielihood Bayesian Model Averaging (BMA) BMA can be approximated through replacement of θ by an estimate θ ˆ which maximizes the lielihood p M θ. Obtaining such maximum lielihood () estimates entails calibrating each model against (conditioning on) the data using well-established statistical inverse methods. Approximating pδ M by pδ M where the subscript indicates estimation of θ was shown to be useful (for a scalar ) in the statistical literature. Neuman (23) proposed evaluating the weights pm based on a result due to Kashyap (982) according to (Ye et al. 24) K pm exp KI pm / exp KI l l p M l where KI KI min 2 2 KI ˆ 2ln 2ln ˆ N KI p M θ p θ M N ln ln F M KI being the socalled Kashyap model selection (or information) criterion for model M 2 KI min its minimum value over 2ln p M θ 2ln p θ M a negative log lielihood incorporating all candidate models and prior measurements of the parameters (if available) evaluated at θ ˆ. Here N is the dimension of θ (number of adjustable parameters associated with model M ) N is the dimension of (number of discrete data points which may include measured parameter values) and F is a normalized (by N ) observed Fisher information matrix. In the limit of large N / N KI reduces asymptotically to the 2ln ˆ BI p M θ N ln N (Ye et al. 28). Bayesian selection (or information) criterion 3. EFFET OF ATA AUGMENTATION ON UNERTAINTY We now as the question how would an augmented data base where is a potential new dataset not presently available (and thus unnown) impact the above BMA uncertainty analysis. We address the question through BMA prediction of denoted by and assessment of the corresponding predictive uncertainty. Our proposed analysis entails the following steps (for additional details see Neuman et al. 2):. Postulate a set M of K mutually independent models M with parameters θ for the desired output vector Δ 2. Obtain estimates θ ˆ of θ by calibrating each minimization of the log lielihood 2ln p M 2ln p M M against available data through θ θ then compute a corresponding parameter estimation covariance Γ and KI 3. ompute p M K exp / exp l l l 2 KI p M 2 KI p M where the subscript designates the estimation process in step 2 4. For each model M estimate EΔ M and ovδ M either through linearization or via Monte arlo simulation (we use the latter approach below): a. raw random samples (realizations) of θ from a multivariate Gaussian distribution with mean θ ˆ and covariance Γ b. Estimate EΔ M θ and ov M Δ θ for each realization of θ
5 c. Average over all realizations of ov Δ M E Δ M and θ to obtain sample estimates of 5. ompute E Δ ovδ and/or Tr ov Δ 6. Postulate a set P of I alternative geostatistical statistical or stochastic models P i with parameters π i for a potential (presently unavailable) data set the models P i may be independent of M may form extensions of M or may coincide with the latter as in the computational examples described below 7. Predict multivariate statistics of conditional on via BMA by means of the model set P using a procedure paralleling that described for Δ in steps Estimate E Δ and ovδ where the subscript designates the estimation process in step 2 with respect to an augmented data set either through linearization followed by step or via Monte arlo simulation by using the statistics of from step 7 to generate random realizations of (we use the second option below) for each realization and for each model M : a. Obtain estimates θ ˆ of θ by minimizing this negative log lielihood with respect to θ then compute the corresponding estimation covariance Γ and KI p M b. ompute c. For each model M estimate EΔ M and ovδ M via Monte arlo simulation: i. raw random samples (realizations) of θ from a multivariate Gaussian distribution with mean ˆ θ and covariance Γ E M ov Δ M θ for each ii. Estimate Δ θ and realization of θ iii. Average over all realizations of E M Δ and ov Δ M θ to obtain sample estimates of d. ompute E Δ ovδ and/or Tr ov Δ ov Δ 9. Average over all realizations of to obtain sample estimates of E Δ and/or Tr ov Δ. Repeat steps 6-9 for different sets 2 3 of potential data and select that set which maximizes the difference Tr ov E Δ Tr ov Δ Tr E ov Δ between the trace conditional on and the expected trace conditional on and. 4. APPLIATION TO AIR PERMEABILITIES OF FRATURE TUFF We apply the above approach to log air permeability data (log ) from unsaturated fractured tuff near Superior Arizona. Spatially distributed log data were obtained by Guzman et al. (996) based on a steady state interpretation of 84 pneumatic injection tests in -m-length intervals along 6 boreholes at the site (Figure ). Five of the boreholes (V2 W2A X2 Y2 Z2) are 3-m long and one (Y3) has a length of 45 m five (W2A X2 Y2 Y3 Z2) are inclined at 45 o and one (V2) is vertical. We consider two crossvalidation cases: V I in which log measured in W2a Y3 and Z2 play the role of existing data
6 Z (m) boreholes X2 and Y2 are the sites of potential new data and 2 the goal being to predict log Δ along V2 V II in which log measured in W2a X2 and Y2 play the role of boreholes V2 and Z2 are the sites of potential new data and 2 the goal being to predict log Δ along Y3 (true potential data are denoted by and their estimates by ). Given that in each case one has funds to measure log in only one borehole which among boreholes X2 and Y2 should one conduct such measurements in case V I and which among V2 and Z2 in case V II? The issue is compounded by uncertainty about the correct geostatistical model and parameters to be employed in each case. To address it we consider three alternative variogram models: exponential (Exp) spherical (Sph) and power (Pow). etails about these models their calibration and estimation of model probabilities are given in Ye et al. (24 28). Figure 2 shows how the three variogram models calibrated against and 2 compare with sample variograms in V I and V II. The corresponding posterior probabilities are listed in Table showing a preference for the exponential model in all cases. -5 W2a V2 X2 Y2 Y Z Y (m) X (m) 4 3 Figure. Spatial locations of 84 -m-scale log data at ALRS. V I (a) Exp.2 Sph Pow (c) Exp.2 Sph Pow V II.8 (b).6.4 Exp.2 Sph Pow (d).6.4 Exp.2 Sph Pow (e).4 Exp.2 Sph Pow (f).4 Exp.2 Sph Pow Figure 2. Sample and fitted variograms based on (a-b) (c-d) and (e-f). 2 Figure 3 compares measured log Δ values along borehole V2 with predicted values and their 95% confidence intervals obtained with each variogram model and with BMA based on and
7 in V I. The predictions are contained within the 95% confidence intervals which are seen to 2 narrow down slightly with data augmentation. orresponding values of Var Var and Var Var 2 are depicted in Figure 4. Augmenting the sample size from 5 () to 35 ( ) and 33 ( ) is seen to reduce predictive variance along borehole V2. As Δ Δ = 2.3 and 2 Tr ov Tr ov Tr ov Δ Tr ov Δ = 2. collecting new data along borehole X2 is seen to yield greater uncertainty reduction than doing so along borehole Y2. Table. Posterior probabilities (%) based on and. 2 V I V II Model Pow Exp Sph Pow Exp Sph p(m ) p(m ' ) p(m ' 2 ) (a) Pow - -2 (b) Exp ata index ata index - (c) Sph ata index - (d) BMA ata index Figure 3. Measured log (+) predicted (solid) and 95% confidence intervals (dashed) for (a) Pow (b) Exp (c) Sph and (d) BMA on V2 based on (blac) (blue) and (green) in V I. Since in reality the real data randomly generated data sets 2 E of 2 are unnown we predict them through the mean r. Figure 5 shows E the 2 r realizations and their 95% confidence interval in V I. E is seen to represent a smoothed version of some values of which lie outside the confidence interval suggesting a slight bias in the generated values. Variations of Var Δ E Var Δ and Var E Δ along borehole V2 are shown in Figure 6. A comparison between Figures 6 and 4 reveals that although Tr ov Δ 3.64 and Tr 3.78 ov Δ exceed Tr ov Δ by a slight amount the spatial patterns in the two figures are similar. Liewise though E VarΔ
8 2 2 in Figure 4 these too have near identical spatial patterns. One therefore expects the estimated variance reduction Var E Δ to exhibit a pattern similar to that of the true variance reduction and E VarΔ in Figure 6 tends to exceed Var 2 Δ and Var Δ 2 2 Var Δ Var Δ. The predicted variance reduction measure Tr ov E( ) Δ =.59 is seen in Table 2 to be smaller than its true counterpart ' Tr ov Tr ov ' Δ Δ = 2.3 the same applying to Tr ov E( ) Δ = Δ Δ =.. Selecting borehole X2 as the '.42 in comparison to Tr ov Tr ov ' target for collecting additional data in V I is thus seen in retrospect to have been a correct choice. Though space does not allow us to present all details of analyzing V II we nevertheless summarize the final results in Table 2. In this case the analysis indicates and validates the choice of borehole V2 as the optimum target for data collection. Variance.8.7 (a) Var( ).3 Var( ).2 Var( ) - Var( ) ata index Variance Var( ) Var( 2 ) Var( ) - Var( 2 ) ata index (b) Var Var along borehole V2 based on (a) and (b) in V I. Figure 4. Variations of Var Var and (a) ata index - -2 (b) ata index Figure 5. E (solid blue) 2 realizations r (gray) and 95% confidence interval of (red) all 2 realizations (dashed blue) in V I.
9 Variance Var( ) E Var( ) Var E( ) ata index (a) Variance Var( ) E 2 Var( 2 ) 2 Var 2 E( 2 ) ata index (b) Figure 6. Variation of Var E Var Δ and Var E Δ along predicted borehole V2 based on (a) + and (b) + 2 in V I. Table 2. Predicted and actual trace variance reductions in V I and V II. V I V II 2 2 Predicted Actual ONLUSIONS. A multimodel approach to optimum value-of-information or data-worth analyses has been proposed based on a Bayesian model averaging (BMA) framewor. We have focused on a maximum lielihood (BMA) variant of BMA that (a) is compatible with both deterministic and stochastic models (b) admits but does not require prior information about the parameters (c) is consistent with modern statistical methods of hydrologic model calibration (d) allows approximating lead predictive moments of any model by linearization and (e) updates model posterior probabilities as well as parameter estimates on the basis of potential new data both before and after such data become actually available. 2. The proposed approach should be of help to hydrogeologists in designing the collection of characterization and monitoring data at contaminated sites in a cost-effective manner by maximizing their benefit under given cost constraints. Benefits would accrue from optimum gain in information or reduction in predictive uncertainty (and ris) upon considering jointly not only traditional sources of uncertainty such as those affecting model parameters and the reliability of data but also lac of certainty about the underlying models. 6. AKNOWLEGEMENT This research was supported in part through a contract between the University of Arizona and Vanderbilt University under the onsortium for Ris Evaluation with Staeholder Participation (RESP) III funded by the U.S. epartment of Energy. The first and second authors were supported in part by NSF-EAR grant 974 and OE-ERSP grant E-S2687.
10 7. REFERENING Aaie H. (974) A new loo at statistical model identification. IEEE Transactions on Automatic ontrol A Bac P.-E. (27) A model for estimating the value of sampling programs and the optimal number of samples for contaminated soil. Environ. Geol. doi:.7/s Beven K.J. and J. Freer (2) Equifinality data assimilation and uncertainty estimation in mechanistic modelling of complex environmental systems using the GLUE methodology. Jour. Hydrology 249: 29. ausman A.M. J. oherty.. Langevin and M.. Suop (2) Quantifying data worth toward reducing predictive uncertainty. Ground Water 48(5) eutsch.v. and A.G. Journel (998) GSLIB: Geostatistical Software Library and User s Guide 2nd ed. Oxford Univ. Press New Yor. i Federico V. and S.P. Neuman (997) Scaling of random fields by means of truncated power variograms and associated spectra. Water Resour. Res. 33(5) Feyen L. and S.M. Gorelic (25) Framewor to evaluate the worth of hydraulic conductivity data for optimal groundwater resources management in ecologically sensitive areas. Water Resour. Res. 4 W39 doi:.29/23wr29. Freer J. K. Beven and B. Ambroise (996) Bayesian estimation of uncertainty in runoff predictioin and the value of data: An application of the GLUE approach. Water Resour. Res. 32(7) Freeze R.A. J. Massman L. Smith T. Sperling and B. James (99) Hydrogeological decision analysis:. A framewor. Ground Water 28(5) Freeze R.A. J. Brude J. Massman T. Sperling and L. Smith (992) Hydrogeological decision analysis: 4. The concept of data worth and its use in the development of site investigation strategies. Ground Water 3(4) Guzman A. G. A.M. Geddis M.J. Henrich.F. Lohrstorfer and S.P. Neuman (996) Summary of Air Permeability ata From Single-Hole Injection Tests in Unsaturated Fractured Tuffs at the Apache Leap Research Site: Results of Steady-State Test Interpretation. Rep. NUREG/R-636 prepared for U.S. Nuclear Regulatory ommission Washington.. Hoeting J.A.. Madigan A.E. Raftery and.t. Volinsy (999) Bayesian model averaging: A tutorial. Statist. Sci. 4(4) Hurvich.M. and.-l. Tsai (989) Regression and time series model selection in small simple. Biometria 76(2) James B.R. and R.A. Freeze (993) The worth of data in predicting aquitard continuity in hydrogeological design. Water Resour. Res. 29(7) Kashyap R.L. (982) Optimal choice of AR and MA parts in autoregressive moving average models. IEEE Transactions on Pattern Analysis and Machine Intelligence 4(2) Neuman S.P. (23) Maximum lielihood Bayesian averaging of alternative conceptual-mathematical models. Stoch. Environmental Research and Ris Assessment 7(5) OI:.7/s Neuman S.P. L.Xue M. Ye and. Lu (2) Bayesian analysis of data-worth considering model and parameter uncertainties Advances in Water Resources in press. Norberg T. and L. ROSEN (26) alculating the optimal number of contaminant samples by means of data worth analysis. Environmetrics Nowa W. F.P.J. e Barros and Y. Rubin (2) Bayesian geostatistical design: Tas-driven site investigation when the geostatistical model is uncertain. Water Resour. Res. 46 W35535 doi:.29/29wr832. Rojas R. L. Feyen O. Batelaan and A. assargues (2) On the value of conditioning data to reduce conceptual model uncertainty in groundwater modeling. Water Resour. Res. doi:.29/29wr8822 in press. Russell K.T. and A.J. Rabideau (2) ecision analysis for pump-and-treat design. Groundwater Monitoring and Remediation Summer Sain S.R. and R. Furrer (2) ombining climate model output via model correlations. Stochastic Environmental Research and Ris Assessment 24(6) OI:.7/s Thompson M. and T. Fearn (996) What exactly is fitness for purpose in analytical measurements? Analyst 2(March) Ye M. S.P. Neuman and P.. Meyer (24) Maximum Lielihood Bayesian averaging of spatial variability models in unsaturated fractured tuff. Water Resour. Res. 4 W53 doi:.29/23wr2557. Ye M. P.. Meyer and S.P. Neuman (28) On model selection criteria in multimodel analysis. Water Resour. Res. 44 W3428 doi:.29/28wr683. Yoota F. and K. Thompson (24) Value of information literature analysis: A review of applications in health ris assessment. Med. ecis. Ma
Maximum likelihood Bayesian averaging of spatial variability models in unsaturated fractured tuff
WATER RESOURCES RESEARCH, VOL. 40,, doi:10.1029/2003wr002557, 2004 Maximum likelihood Bayesian averaging of spatial variability models in unsaturated fractured tuff Ming Ye Pacific Northwest National Laboratory,
More informationPARAMETER ESTIMATION AND ORDER SELECTION FOR LINEAR REGRESSION PROBLEMS. Yngve Selén and Erik G. Larsson
PARAMETER ESTIMATION AND ORDER SELECTION FOR LINEAR REGRESSION PROBLEMS Yngve Selén and Eri G Larsson Dept of Information Technology Uppsala University, PO Box 337 SE-71 Uppsala, Sweden email: yngveselen@ituuse
More informationSensitivity analysis and assessment of prior model probabilities in MLBMA with application to unsaturated fractured tuff
WATER RESOURCES RESEARCH, VOL. 41,, doi:10.1029/2005wr004260, 2005 Sensitivity analysis and assessment of prior model probabilities in MLBMA with application to unsaturated fractured tuff Ming Ye, 1 Shlomo
More informationThe Generalized Likelihood Uncertainty Estimation methodology
CHAPTER 4 The Generalized Likelihood Uncertainty Estimation methodology Calibration and uncertainty estimation based upon a statistical framework is aimed at finding an optimal set of models, parameters
More informationEnsemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations
Ensemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations Alberto Guadagnini (1,), Marco Panzeri (1), Monica Riva (1,), Shlomo P. Neuman () (1) Department of
More informationConditional Distribution Fitting of High Dimensional Stationary Data
Conditional Distribution Fitting of High Dimensional Stationary Data Miguel Cuba and Oy Leuangthong The second order stationary assumption implies the spatial variability defined by the variogram is constant
More informationAnalysis of Regression and Bayesian Predictive Uncertainty Measures
Analysis of and Predictive Uncertainty Measures Dan Lu, Mary C. Hill, Ming Ye Florida State University, dl7f@fsu.edu, mye@fsu.edu, Tallahassee, FL, USA U.S. Geological Survey, mchill@usgs.gov, Boulder,
More informationStatistical Rock Physics
Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
More informationAccounting for conceptual model uncertainty via maximum likelihood Bayesian model averaging
Calibration ami Reliability in Groundwater Modelling: A Few Steps Closer to Reality (Proceeding of ModelCARE'2002. Prague, Czech Republic, June 2002). IAHS Publ. no. 277. 2002. 303 Accounting for conceptual
More informationBAYESIAN ESTIMATION OF UNKNOWN PARAMETERS OVER NETWORKS
BAYESIAN ESTIMATION OF UNKNOWN PARAMETERS OVER NETWORKS Petar M. Djurić Dept. of Electrical & Computer Engineering Stony Brook University Stony Brook, NY 11794, USA e-mail: petar.djuric@stonybrook.edu
More informationA first look at the performances of a Bayesian chart to monitor. the ratio of two Weibull percentiles
A first loo at the performances of a Bayesian chart to monitor the ratio of two Weibull percentiles Pasquale Erto University of Naples Federico II, Naples, Italy e-mail: ertopa@unina.it Abstract. The aim
More informationAdvanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland
EnviroInfo 2004 (Geneva) Sh@ring EnviroInfo 2004 Advanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland Mikhail Kanevski 1, Michel Maignan 1
More informationDependence of Bayesian Model Selection Criteria and Fisher Information Matrix on Sample Size
DOI 10.1007/s11004-011-9359-0 Dependence of Bayesian Model Selection Criteria and Fisher Information Matrix on Sample Size Dan Lu Ming Ye Shlomo P. Neuman Received: 11 June 2010 / Accepted: 26 August 2011
More informationAutomatic Determination of Uncertainty versus Data Density
Automatic Determination of Uncertainty versus Data Density Brandon Wilde and Clayton V. Deutsch It is useful to know how various measures of uncertainty respond to changes in data density. Calculating
More informationPRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH
PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH SURESH TRIPATHI Geostatistical Society of India Assumptions and Geostatistical Variogram
More informationLabor-Supply Shifts and Economic Fluctuations. Technical Appendix
Labor-Supply Shifts and Economic Fluctuations Technical Appendix Yongsung Chang Department of Economics University of Pennsylvania Frank Schorfheide Department of Economics University of Pennsylvania January
More informationType-curve estimation of statistical heterogeneity
WATER RESOURCES RESEARCH, VOL. 40,, doi:10.1029/2003wr002405, 2004 Type-curve estimation of statistical heterogeneity Shlomo P. Neuman Department of Hydrology and Water Resources, University of Arizona,
More informationAsymptotic analysis of three-dimensional pressure interference tests: A point source solution
WATER RESOURCES RESEARCH, VOL. 41,, doi:10.1029/2004wr003431, 2005 Asymptotic analysis of three-dimensional pressure interference tests: A point source solution Walter A. Illman 1 Department of Geoscience,
More informationNeutron inverse kinetics via Gaussian Processes
Neutron inverse kinetics via Gaussian Processes P. Picca Politecnico di Torino, Torino, Italy R. Furfaro University of Arizona, Tucson, Arizona Outline Introduction Review of inverse kinetics techniques
More informationLinear Prediction Theory
Linear Prediction Theory Joseph A. O Sullivan ESE 524 Spring 29 March 3, 29 Overview The problem of estimating a value of a random process given other values of the random process is pervasive. Many problems
More informationEstimation of Operational Risk Capital Charge under Parameter Uncertainty
Estimation of Operational Risk Capital Charge under Parameter Uncertainty Pavel V. Shevchenko Principal Research Scientist, CSIRO Mathematical and Information Sciences, Sydney, Locked Bag 17, North Ryde,
More informationA008 THE PROBABILITY PERTURBATION METHOD AN ALTERNATIVE TO A TRADITIONAL BAYESIAN APPROACH FOR SOLVING INVERSE PROBLEMS
A008 THE PROAILITY PERTURATION METHOD AN ALTERNATIVE TO A TRADITIONAL AYESIAN APPROAH FOR SOLVING INVERSE PROLEMS Jef AERS Stanford University, Petroleum Engineering, Stanford A 94305-2220 USA Abstract
More informationA MultiGaussian Approach to Assess Block Grade Uncertainty
A MultiGaussian Approach to Assess Block Grade Uncertainty Julián M. Ortiz 1, Oy Leuangthong 2, and Clayton V. Deutsch 2 1 Department of Mining Engineering, University of Chile 2 Department of Civil &
More informationMonte Carlo Simulation. CWR 6536 Stochastic Subsurface Hydrology
Monte Carlo Simulation CWR 6536 Stochastic Subsurface Hydrology Steps in Monte Carlo Simulation Create input sample space with known distribution, e.g. ensemble of all possible combinations of v, D, q,
More informationSequential Importance Sampling for Rare Event Estimation with Computer Experiments
Sequential Importance Sampling for Rare Event Estimation with Computer Experiments Brian Williams and Rick Picard LA-UR-12-22467 Statistical Sciences Group, Los Alamos National Laboratory Abstract Importance
More informationUncertainty propagation in a sequential model for flood forecasting
Predictions in Ungauged Basins: Promise and Progress (Proceedings of symposium S7 held during the Seventh IAHS Scientific Assembly at Foz do Iguaçu, Brazil, April 2005). IAHS Publ. 303, 2006. 177 Uncertainty
More informationClimate Change: the Uncertainty of Certainty
Climate Change: the Uncertainty of Certainty Reinhard Furrer, UZH JSS, Geneva Oct. 30, 2009 Collaboration with: Stephan Sain - NCAR Reto Knutti - ETHZ Claudia Tebaldi - Climate Central Ryan Ford, Doug
More informationImproved inverse modeling for flow and transport in subsurface media: Combined parameter and state estimation
GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L18408, doi:10.1029/2005gl023940, 2005 Improved inverse modeling for flow and transport in subsurface media: Combined parameter and state estimation Jasper A. Vrugt,
More informationCONSTRAINED OPTIMIZATION OVER DISCRETE SETS VIA SPSA WITH APPLICATION TO NON-SEPARABLE RESOURCE ALLOCATION
Proceedings of the 200 Winter Simulation Conference B. A. Peters, J. S. Smith, D. J. Medeiros, and M. W. Rohrer, eds. CONSTRAINED OPTIMIZATION OVER DISCRETE SETS VIA SPSA WITH APPLICATION TO NON-SEPARABLE
More informationPOPULAR CARTOGRAPHIC AREAL INTERPOLATION METHODS VIEWED FROM A GEOSTATISTICAL PERSPECTIVE
CO-282 POPULAR CARTOGRAPHIC AREAL INTERPOLATION METHODS VIEWED FROM A GEOSTATISTICAL PERSPECTIVE KYRIAKIDIS P. University of California Santa Barbara, MYTILENE, GREECE ABSTRACT Cartographic areal interpolation
More informationBayesian Approaches for Model Evaluation and Uncertainty Estimation
Bayesian Approaches for Model Evaluation and Uncertainty Estimation Tsuyoshi TADA National Defense Academy, Japan APHW 2006 Bangkok 1 Introduction In wellgauged basins Parameters can be obtained by calibrations.
More informationMachine Learning Lecture Notes
Machine Learning Lecture Notes Predrag Radivojac January 25, 205 Basic Principles of Parameter Estimation In probabilistic modeling, we are typically presented with a set of observations and the objective
More informationA characterization of consistency of model weights given partial information in normal linear models
Statistics & Probability Letters ( ) A characterization of consistency of model weights given partial information in normal linear models Hubert Wong a;, Bertrand Clare b;1 a Department of Health Care
More informationDevelopment of Stochastic Artificial Neural Networks for Hydrological Prediction
Development of Stochastic Artificial Neural Networks for Hydrological Prediction G. B. Kingston, M. F. Lambert and H. R. Maier Centre for Applied Modelling in Water Engineering, School of Civil and Environmental
More informationReservoir Uncertainty Calculation by Large Scale Modeling
Reservoir Uncertainty Calculation by Large Scale Modeling Naeem Alshehri and Clayton V. Deutsch It is important to have a good estimate of the amount of oil or gas in a reservoir. The uncertainty in reserve
More informationChecking up on the neighbors: Quantifying uncertainty in relative event location
Checking up on the neighbors: Quantifying uncertainty in relative event location The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationOn Mixture Regression Shrinkage and Selection via the MR-LASSO
On Mixture Regression Shrinage and Selection via the MR-LASSO Ronghua Luo, Hansheng Wang, and Chih-Ling Tsai Guanghua School of Management, Peing University & Graduate School of Management, University
More informationRecurrent Latent Variable Networks for Session-Based Recommendation
Recurrent Latent Variable Networks for Session-Based Recommendation Panayiotis Christodoulou Cyprus University of Technology paa.christodoulou@edu.cut.ac.cy 27/8/2017 Panayiotis Christodoulou (C.U.T.)
More informationBayesian methodology to stochastic capture zone determination: Conditioning on transmissivity measurements
WATER RESOURCES RESEARCH, VOL. 38, NO. 9, 1164, doi:10.1029/2001wr000950, 2002 Bayesian methodology to stochastic capture zone determination: Conditioning on transmissivity measurements L. Feyen, 1 P.
More informationEntropy of Gaussian Random Functions and Consequences in Geostatistics
Entropy of Gaussian Random Functions and Consequences in Geostatistics Paula Larrondo (larrondo@ualberta.ca) Department of Civil & Environmental Engineering University of Alberta Abstract Sequential Gaussian
More informationLecture : Probabilistic Machine Learning
Lecture : Probabilistic Machine Learning Riashat Islam Reasoning and Learning Lab McGill University September 11, 2018 ML : Many Methods with Many Links Modelling Views of Machine Learning Machine Learning
More informationBasic concepts in estimation
Basic concepts in estimation Random and nonrandom parameters Definitions of estimates ML Maimum Lielihood MAP Maimum A Posteriori LS Least Squares MMS Minimum Mean square rror Measures of quality of estimates
More informationDynamic Data Driven Simulations in Stochastic Environments
Computing 77, 321 333 (26) Digital Object Identifier (DOI) 1.17/s67-6-165-3 Dynamic Data Driven Simulations in Stochastic Environments C. Douglas, Lexington, Y. Efendiev, R. Ewing, College Station, V.
More informationKriging models with Gaussian processes - covariance function estimation and impact of spatial sampling
Kriging models with Gaussian processes - covariance function estimation and impact of spatial sampling François Bachoc former PhD advisor: Josselin Garnier former CEA advisor: Jean-Marc Martinez Department
More informationCROSS-VALIDATION OF CONTROLLED DYNAMIC MODELS: BAYESIAN APPROACH
CROSS-VALIDATION OF CONTROLLED DYNAMIC MODELS: BAYESIAN APPROACH Miroslav Kárný, Petr Nedoma,Václav Šmídl Institute of Information Theory and Automation AV ČR, P.O.Box 18, 180 00 Praha 8, Czech Republic
More informationCombining geological surface data and geostatistical model for Enhanced Subsurface geological model
Combining geological surface data and geostatistical model for Enhanced Subsurface geological model M. Kurniawan Alfadli, Nanda Natasia, Iyan Haryanto Faculty of Geological Engineering Jalan Raya Bandung
More informationGaussian Mixture Distance for Information Retrieval
Gaussian Mixture Distance for Information Retrieval X.Q. Li and I. King fxqli, ingg@cse.cuh.edu.h Department of omputer Science & Engineering The hinese University of Hong Kong Shatin, New Territories,
More informationarxiv: v1 [stat.ap] 27 Mar 2015
Submitted to the Annals of Applied Statistics A NOTE ON THE SPECIFIC SOURCE IDENTIFICATION PROBLEM IN FORENSIC SCIENCE IN THE PRESENCE OF UNCERTAINTY ABOUT THE BACKGROUND POPULATION By Danica M. Ommen,
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
More informationSTATS 200: Introduction to Statistical Inference. Lecture 29: Course review
STATS 200: Introduction to Statistical Inference Lecture 29: Course review Course review We started in Lecture 1 with a fundamental assumption: Data is a realization of a random process. The goal throughout
More informationTools for Parameter Estimation and Propagation of Uncertainty
Tools for Parameter Estimation and Propagation of Uncertainty Brian Borchers Department of Mathematics New Mexico Tech Socorro, NM 87801 borchers@nmt.edu Outline Models, parameters, parameter estimation,
More informationJoint inversion of geophysical and hydrological data for improved subsurface characterization
Joint inversion of geophysical and hydrological data for improved subsurface characterization Michael B. Kowalsky, Jinsong Chen and Susan S. Hubbard, Lawrence Berkeley National Lab., Berkeley, California,
More informationBayesian model selection: methodology, computation and applications
Bayesian model selection: methodology, computation and applications David Nott Department of Statistics and Applied Probability National University of Singapore Statistical Genomics Summer School Program
More informationAssessing Regime Uncertainty Through Reversible Jump McMC
Assessing Regime Uncertainty Through Reversible Jump McMC August 14, 2008 1 Introduction Background Research Question 2 The RJMcMC Method McMC RJMcMC Algorithm Dependent Proposals Independent Proposals
More informationA Stochastic Collocation based. for Data Assimilation
A Stochastic Collocation based Kalman Filter (SCKF) for Data Assimilation Lingzao Zeng and Dongxiao Zhang University of Southern California August 11, 2009 Los Angeles Outline Introduction SCKF Algorithm
More informationModeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation
Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Z. W. LI, X. L. DING Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hung
More informationF denotes cumulative density. denotes probability density function; (.)
BAYESIAN ANALYSIS: FOREWORDS Notation. System means the real thing and a model is an assumed mathematical form for the system.. he probability model class M contains the set of the all admissible models
More informationPractical Bayesian Optimization of Machine Learning. Learning Algorithms
Practical Bayesian Optimization of Machine Learning Algorithms CS 294 University of California, Berkeley Tuesday, April 20, 2016 Motivation Machine Learning Algorithms (MLA s) have hyperparameters that
More informationSimultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations
Simultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations C. Rentier,, A. Dassargues Hydrogeology Group, Departement
More informationSPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA
SPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA D. Pokrajac Center for Information Science and Technology Temple University Philadelphia, Pennsylvania A. Lazarevic Computer
More informationAUTOMATED TEMPLATE MATCHING METHOD FOR NMIS AT THE Y-12 NATIONAL SECURITY COMPLEX
AUTOMATED TEMPLATE MATCHING METHOD FOR NMIS AT THE Y-1 NATIONAL SECURITY COMPLEX J. A. Mullens, J. K. Mattingly, L. G. Chiang, R. B. Oberer, J. T. Mihalczo ABSTRACT This paper describes a template matching
More informationBasics of Uncertainty Analysis
Basics of Uncertainty Analysis Chapter Six Basics of Uncertainty Analysis 6.1 Introduction As shown in Fig. 6.1, analysis models are used to predict the performances or behaviors of a product under design.
More informationGrundlagen der Künstlichen Intelligenz
Grundlagen der Künstlichen Intelligenz Uncertainty & Probabilities & Bandits Daniel Hennes 16.11.2017 (WS 2017/18) University Stuttgart - IPVS - Machine Learning & Robotics 1 Today Uncertainty Probability
More informationA new Hierarchical Bayes approach to ensemble-variational data assimilation
A new Hierarchical Bayes approach to ensemble-variational data assimilation Michael Tsyrulnikov and Alexander Rakitko HydroMetCenter of Russia College Park, 20 Oct 2014 Michael Tsyrulnikov and Alexander
More informationGroundwater Resources Management under. Uncertainty Harald Kunstmann*, Wolfgang Kinzelbach* & Gerrit van Tender**
Groundwater Resources Management under Uncertainty Harald Kunstmann*, Wolfgang Kinzelbach* & Gerrit van Tender**, C/f &OP3 Zwnc/?, Email: kunstmann@ihw. baum. ethz. ch Abstract Groundwater resources and
More informationMustafa H. Tongarlak Bruce E. Ankenman Barry L. Nelson
Proceedings of the 0 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelspach, K. P. White, and M. Fu, eds. RELATIVE ERROR STOCHASTIC KRIGING Mustafa H. Tongarlak Bruce E. Ankenman Barry L.
More informationModule 2. Random Processes. Version 2, ECE IIT, Kharagpur
Module Random Processes Version, ECE IIT, Kharagpur Lesson 9 Introduction to Statistical Signal Processing Version, ECE IIT, Kharagpur After reading this lesson, you will learn about Hypotheses testing
More informationσ(a) = a N (x; 0, 1 2 ) dx. σ(a) = Φ(a) =
Until now we have always worked with likelihoods and prior distributions that were conjugate to each other, allowing the computation of the posterior distribution to be done in closed form. Unfortunately,
More informationValue of Information Analysis with Structural Reliability Methods
Accepted for publication in Structural Safety, special issue in the honor of Prof. Wilson Tang August 2013 Value of Information Analysis with Structural Reliability Methods Daniel Straub Engineering Risk
More informationSOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC
SOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC Alfredo Vaccaro RTSI 2015 - September 16-18, 2015, Torino, Italy RESEARCH MOTIVATIONS Power Flow (PF) and
More informationProbabilistic numerics for deep learning
Presenter: Shijia Wang Department of Engineering Science, University of Oxford rning (RLSS) Summer School, Montreal 2017 Outline 1 Introduction Probabilistic Numerics 2 Components Probabilistic modeling
More informationEmpirical Bayesian Kriging
Empirical Bayesian Kriging Implemented in ArcGIS Geostatistical Analyst By Konstantin Krivoruchko, Senior Research Associate, Software Development Team, Esri Obtaining reliable environmental measurements
More informationAsymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands
Asymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands Elizabeth C. Mannshardt-Shamseldin Advisor: Richard L. Smith Duke University Department
More informationLarge Scale Modeling by Bayesian Updating Techniques
Large Scale Modeling by Bayesian Updating Techniques Weishan Ren Centre for Computational Geostatistics Department of Civil and Environmental Engineering University of Alberta Large scale models are useful
More informationParametric Models. Dr. Shuang LIANG. School of Software Engineering TongJi University Fall, 2012
Parametric Models Dr. Shuang LIANG School of Software Engineering TongJi University Fall, 2012 Today s Topics Maximum Likelihood Estimation Bayesian Density Estimation Today s Topics Maximum Likelihood
More informationCombine Monte Carlo with Exhaustive Search: Effective Variational Inference and Policy Gradient Reinforcement Learning
Combine Monte Carlo with Exhaustive Search: Effective Variational Inference and Policy Gradient Reinforcement Learning Michalis K. Titsias Department of Informatics Athens University of Economics and Business
More informationCramér-Rao Bounds for Estimation of Linear System Noise Covariances
Journal of Mechanical Engineering and Automation (): 6- DOI: 593/jjmea Cramér-Rao Bounds for Estimation of Linear System oise Covariances Peter Matiso * Vladimír Havlena Czech echnical University in Prague
More informationACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION. <www.ie.ncsu.edu/jwilson> May 22, 2006
ACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION Slide 1 Faker Zouaoui Sabre Holdings James R. Wilson NC State University May, 006 Slide From American
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 2008 Paper 241 A Note on Risk Prediction for Case-Control Studies Sherri Rose Mark J. van der Laan Division
More informationThe Use of Copulas to Model Conditional Expectation for Multivariate Data
The Use of Copulas to Model Conditional Expectation for Multivariate Data Kääri, Meelis; Selart, Anne; Kääri, Ene University of Tartu, Institute of Mathematical Statistics J. Liivi 2 50409 Tartu, Estonia
More informationLecture 4: Probabilistic Learning. Estimation Theory. Classification with Probability Distributions
DD2431 Autumn, 2014 1 2 3 Classification with Probability Distributions Estimation Theory Classification in the last lecture we assumed we new: P(y) Prior P(x y) Lielihood x2 x features y {ω 1,..., ω K
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationHydraulic tomography: Development of a new aquifer test method
WATER RESOURCES RESEARCH, VOL. 36, NO. 8, PAGES 2095 2105, AUGUST 2000 Hydraulic tomography: Development of a new aquifer test method T.-C. Jim Yeh and Shuyun Liu Department of Hydrology and Water Resources,
More informationArtificial Intelligence & Sequential Decision Problems
Artificial Intelligence & Sequential Decision Problems (CIV6540 - Machine Learning for Civil Engineers) Professor: James-A. Goulet Département des génies civil, géologique et des mines Chapter 15 Goulet
More informationSpatial Inference of Nitrate Concentrations in Groundwater
Spatial Inference of Nitrate Concentrations in Groundwater Dawn Woodard Operations Research & Information Engineering Cornell University joint work with Robert Wolpert, Duke Univ. Dept. of Statistical
More informationUncertainty Quantification in Performance Evaluation of Manufacturing Processes
Uncertainty Quantification in Performance Evaluation of Manufacturing Processes Manufacturing Systems October 27, 2014 Saideep Nannapaneni, Sankaran Mahadevan Vanderbilt University, Nashville, TN Acknowledgement
More informationHierarchical sparse Bayesian learning for structural health monitoring. with incomplete modal data
Hierarchical sparse Bayesian learning for structural health monitoring with incomplete modal data Yong Huang and James L. Beck* Division of Engineering and Applied Science, California Institute of Technology,
More informationInfluence of parameter estimation uncertainty in Kriging: Part 2 Test and case study applications
Hydrology and Earth System Influence Sciences, of 5(), parameter 5 3 estimation (1) uncertainty EGS in Kriging: Part Test and case study applications Influence of parameter estimation uncertainty in Kriging:
More informationBayesian Learning (II)
Universität Potsdam Institut für Informatik Lehrstuhl Maschinelles Lernen Bayesian Learning (II) Niels Landwehr Overview Probabilities, expected values, variance Basic concepts of Bayesian learning MAP
More informationProbability and Information Theory. Sargur N. Srihari
Probability and Information Theory Sargur N. srihari@cedar.buffalo.edu 1 Topics in Probability and Information Theory Overview 1. Why Probability? 2. Random Variables 3. Probability Distributions 4. Marginal
More informationCoregionalization by Linear Combination of Nonorthogonal Components 1
Mathematical Geology, Vol 34, No 4, May 2002 ( C 2002) Coregionalization by Linear Combination of Nonorthogonal Components 1 J A Vargas-Guzmán, 2,3 A W Warrick, 3 and D E Myers 4 This paper applies the
More informationBayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units
Bayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units Sahar Z Zangeneh Robert W. Keener Roderick J.A. Little Abstract In Probability proportional
More informationQuality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2
Quality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2 1 Cleveland State University 2 N & R Engineering www.inl.gov ASTM Symposium on Graphite
More informationNEW ESTIMATORS FOR PARALLEL STEADY-STATE SIMULATIONS
roceedings of the 2009 Winter Simulation Conference M. D. Rossetti, R. R. Hill, B. Johansson, A. Dunkin, and R. G. Ingalls, eds. NEW ESTIMATORS FOR ARALLEL STEADY-STATE SIMULATIONS Ming-hua Hsieh Department
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 208 Petros Koumoutsakos, Jens Honore Walther (Last update: June, 208) IMPORTANT DISCLAIMERS. REFERENCES: Much of the material (ideas,
More informationA Comparison of Particle Filters for Personal Positioning
VI Hotine-Marussi Symposium of Theoretical and Computational Geodesy May 9-June 6. A Comparison of Particle Filters for Personal Positioning D. Petrovich and R. Piché Institute of Mathematics Tampere University
More informationA robust statistically based approach to estimating the probability of contamination occurring between sampling locations
A robust statistically based approach to estimating the probability of contamination occurring between sampling locations Peter Beck Principal Environmental Scientist Image placeholder Image placeholder
More information11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments
11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments M. Camporese (University of Padova), G. Cassiani* (University of Padova), R. Deiana
More informationPart 6: Multivariate Normal and Linear Models
Part 6: Multivariate Normal and Linear Models 1 Multiple measurements Up until now all of our statistical models have been univariate models models for a single measurement on each member of a sample of
More informationGaussian processes. Chuong B. Do (updated by Honglak Lee) November 22, 2008
Gaussian processes Chuong B Do (updated by Honglak Lee) November 22, 2008 Many of the classical machine learning algorithms that we talked about during the first half of this course fit the following pattern:
More information