Calibration of imperfect mathematical models by multiple sources of data with measurement bias

10/27/2018
by   Mengyang Gu, et al.
0

Model calibration involves using experimental or field data to estimate the unknown parameters of a mathematical model. This task is complicated by discrepancy between the model and reality, and by possible bias in the data. We consider model calibration in the presence of both model discrepancy and measurement bias using multiple sources of data. Model discrepancy is often estimated using a Gaussian stochastic process (GaSP), but it has been observed in many studies that the calibrated mathematical model can be far from the reality. Here we show that modeling the discrepancy function via a GaSP often leads to an inconsistent estimation of the calibration parameters even if one has an infinite number of repeated experiments and infinite number of observations in a fixed input domain in each experiment. We introduce the scaled Gaussian stochastic process (S-GaSP) to model the discrepancy function. Unlike the GaSP, the S-GaSP utilizes a non-increasing scaling function which assigns more probability mass on the smaller L_2 loss between the mathematical model and reality, preventing the calibrated mathematical model from deviating too much from reality. We apply our technique to the calibration of a geophysical model of Kīlauea Volcano, Hawai`i, using multiple radar satellite interferograms. We compare the use of models calibrated using multiple data sets simultaneously with results obtained using stacks (averages). We derive distributions for the maximum likelihood estimator and Bayesian inference, both implemented in the "RobustCalibration" package available on CRAN. Analysis of both simulated and real data confirm that our approach can identify the measurement bias and model discrepancy using multiple sources of data, and provide better estimates of model parameters.

READ FULL TEXT

page 5

page 24

page 26

page 27

research
07/10/2018

A theoretical framework of the scaled Gaussian stochastic process in prediction and calibration

The Gaussian stochastic process (GaSP) is a useful technique for predict...
research
01/05/2022

RobustCalibration: Robust Calibration of Computer Models in R

Two fundamental research tasks in science and engineering are forward pr...
research
02/02/2018

Stochastic Kriging for Inadequate Simulation Models

Stochastic kriging is a popular metamodeling technique for representing ...
research
02/17/2023

Bayesian calibration with summary statistics for the prediction of xenon diffusion in UO2 nuclear fuel

The evolution and release of fission gas impacts the performance of UO2 ...
research
07/29/2022

HUG model: an interaction point process for Bayesian detection of multiple sources in groundwaters from hydrochemical data

This paper presents a new interaction point process that integrates geol...
research
06/08/2023

Computational methods for fast Bayesian model assessment via calibrated posterior p-values

Posterior predictive p-values (ppps) have become popular tools for Bayes...
research
04/13/2020

Embedded model discrepancy: A case study of Zika modeling

Mathematical models of epidemiological systems enable investigation of a...

Please sign up or login with your details

Forgot password? Click here to reset