CaliCo: a R package for Bayesian calibration

In this article, we present a recently released R package for Bayesian calibration. Many industrial fields are facing unfeasible or costly field experiments. These experiments are replaced with numerical/computer experiments which are realized by running a numerical code. Bayesian calibration intends to estimate, through a posterior distribution, input parameters of the code in order to make the code outputs close to the available experimental data. The code can be time consuming while the Bayesian calibration implies a lot of code calls which makes studies too burdensome. A discrepancy might also appear between the numerical code and the physical system when facing incompatibility between experimental data and numerical code outputs. The package CaliCo deals with these issues through four statistical models which deal with a time consuming code or not and with discrepancy or not. A guideline for users is provided in order to illustrate the main functions and their arguments. Eventually, a toy example is detailed using CaliCo. This example (based on a real physical system) is in five dimensions and uses simulated data.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[3]  T. J. Mitchell,et al.  Bayesian Prediction of Deterministic Functions, with Applications to the Design and Analysis of Computer Experiments , 1991 .

[4]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[5]  Donald R. Jones,et al.  Efficient Global Optimization of Expensive Black-Box Functions , 1998, J. Glob. Optim..

[6]  Jeong-Soo Park,et al.  A statistical method for tuning a computer code to a data base , 2001 .

[7]  A. O'Hagan,et al.  Bayesian calibration of computer models , 2001 .

[8]  Jeremy E. Oakley,et al.  Bayesian Analysis of Computer Model Outputs , 2002 .

[9]  Thomas J. Santner,et al.  Design and analysis of computer experiments , 1998 .

[10]  Dave Higdon,et al.  Combining Field Data and Computer Simulations for Calibration and Prediction , 2005, SIAM J. Sci. Comput..

[11]  Runze Li,et al.  Design and Modeling for Computer Experiments , 2005 .

[12]  James O. Berger,et al.  A Framework for Validation of Computer Models , 2007, Technometrics.

[13]  James O. Berger,et al.  Modularization in Bayesian analysis, with emphasis on analysis of computer models , 2009 .

[14]  Dirk Eddelbuettel,et al.  Rcpp: Seamless R and C++ Integration , 2011 .

[15]  Gonzalo Garcia-Donato,et al.  Bayesian Emulation, Calibration and Validation of ComputerModels , 2015 .

[16]  D. Bates,et al.  Output Analysis and Diagnostics for MCMC , 2015 .

[17]  Merlin Keller,et al.  Bayesian Model Selection for the Validation of Computer Codes , 2016, Qual. Reliab. Eng. Int..

[18]  T. Hothorn,et al.  Multivariate Normal and t Distributions , 2016 .

[19]  Merlin Keller,et al.  Adaptive Numerical Designs for the Calibration of Computer Codes , 2015, SIAM/ASA J. Uncertain. Quantification.

[20]  Long Wang,et al.  Scaled Gaussian Stochastic Process for Computer Model Calibration and Prediction , 2017, SIAM/ASA J. Uncertain. Quantification.

[21]  M. Keller,et al.  Bayesian calibration of a numerical code for prediction , 2018, 1801.01810.

[22]  James O. Berger,et al.  RobustGaSP: Robust Gaussian Stochastic Process Emulation in R , 2018, R J..

[23]  Mengyang Gu.  Jointly Robust Prior for Gaussian Stochastic Process in Emulation, Calibration and Variable Selection , 2018, Bayesian Analysis.