Bayesian Uncertainty Quantification while Leveraging Multiple Computer Model Runs

dc.contributor.authorWalsh, Stephen A.en
dc.contributor.committeechairFerreira, Marco A. R.en
dc.contributor.committeememberZick, Stephanie E.en
dc.contributor.committeememberFranck, Christopher T.en
dc.contributor.committeememberHigdon, Daviden
dc.contributor.departmentStatisticsen
dc.date.accessioned2023-06-23T08:00:54Zen
dc.date.available2023-06-23T08:00:54Zen
dc.date.issued2023-06-22en
dc.description.abstractIn the face of spatially correlated data, Gaussian process regression is a very common modeling approach. Given observational data, kriging equations will provide the best linear unbiased predictor for the mean at unobserved locations. However, when a computer model provides a complete grid of forecasted values, kriging will not apply. To develop an approach to quantify uncertainty of computer model output in this setting, we leverage information from a collection of computer model runs (e.g., historical forecast and observation pairs for tropical cyclone precipitation totals) through a Bayesian hierarchical framework. This framework allows us to combine information and account for the spatial correlation within and across computer model output. Using maximum likelihood estimates and the corresponding Hessian matrices for Gaussian process parameters, these are input to a Gibbs sampler which provides posterior distributions for parameters of interest. These samples are used to generate predictions which provide uncertainty quantification for a given computer model run (e.g., tropical cyclone precipitation forecast). We then extend this framework using deep Gaussian processes to allow for nonstationary covariance structure, applied to multiple computer model runs from a cosmology application. We also perform sensitivity analyses to understand which parameter inputs most greatly impact cosmological computer model output.en
dc.description.abstractgeneralA crucial theme when analyzing spatial data is that locations that are closer together are more likely to have similar output values (for example, daily precipitation totals). For a particular event, common modeling approach of spatial data is to observe data at numerous locations, and make predictions for locations that were unobserved. In this work, we extend this within-event modeling approach by additionally learning about the uncertainty across different events. Through this extension, we are able to quantify uncertainty for a particular computer model (which may be modeling tropical cyclone precipitation, for example) that does not provide any uncertainty on its own. This framework can be utilized to quantify uncertainty across a vast array of computer model outputs where more than one event or model run has been obtained. We also study how inputting different values into a computer model can influence the values it produces.en
dc.description.degreeDoctor of Philosophyen
dc.format.mediumETDen
dc.identifier.othervt_gsexam:38171en
dc.identifier.urihttp://hdl.handle.net/10919/115494en
dc.language.isoenen
dc.publisherVirginia Techen
dc.rightsIn Copyrighten
dc.rights.urihttp://rightsstatements.org/vocab/InC/1.0/en
dc.subjectGaussian processesen
dc.subjectuncertainty quantificationen
dc.subjectBayesian statisticsen
dc.subjectspatial statisticsen
dc.subjectnonstationary modelingen
dc.subjectmeteorologyen
dc.subjectcosmologyen
dc.titleBayesian Uncertainty Quantification while Leveraging Multiple Computer Model Runsen
dc.typeDissertationen
thesis.degree.disciplineStatisticsen
thesis.degree.grantorVirginia Polytechnic Institute and State Universityen
thesis.degree.leveldoctoralen
thesis.degree.nameDoctor of Philosophyen

Files

Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
Walsh_SA_D_2023.pdf
Size:
13.07 MB
Format:
Adobe Portable Document Format