Computer model validation with functional output

被引:211
作者
Bayarri, M. J.
Berger, J. O.
Cafeo, J.
Garcia-Donato, G.
Liu, F.
Palomo, J.
Parthasarathy, R. J.
Paulo, R.
Sacks, J.
Walsh, D.
机构
[1] Univ Valencia, Dept Stat & Operat Res, E-46100 Valencia, Spain
[2] Duke Univ, Dept Stat Sci, Durham, NC 27708 USA
[3] GM Corp, Res & Dev, Warren, MI 48090 USA
[4] Univ Castilla La Mancha, Dept Econ, Albacete 02071, Spain
[5] Rey Juan Carlos Univ, Stat & Operat Res Dept, Madrid 28945, Spain
[6] Univ Tecn Lisboa, Dept Matemat, ISEG, P-1200781 Lisbon, Portugal
[7] Natl Inst Stat Sci, Res Triangle Pk, NC 27709 USA
[8] Massey Univ, Inst Informat & Math Sci, Auckland, Albany, New Zealand
关键词
computer models; validation; functional data; bias; Bayesian analysis; uncertain inputs;
D O I
10.1214/009053607000000163
中图分类号
O21 [概率论与数理统计]; C8 [统计学];
学科分类号
020208 ; 070103 ; 0714 ;
摘要
A key question in evaluation of computer models is Does the computer model adequately represent reality? A six-step process for computer model validation is set out in Bayarri et al. [Technometrics 49 (2007) 138-154] (and briefly summarized below), based on comparison of computer model runs with field data of the process being modeled. The methodology is particularly suited to treating the major issues associated with the validation process: quantifying multiple sources of error and uncertainty in computer models; combining multiple sources of information; and being able to adapt to different, but related scenarios. Two complications that frequently arise in practice are the need to deal with highly irregular functional data and the need to acknowledge and incorporate uncertainty in the inputs. We develop methodology to deal with both complications. A key part of the approach utilizes a wavelet representation of the functional data, applies a hierarchical version of the scalar validation methodology to the wavelet coefficients, and transforms back, to ultimately compare computer model output with field output. The generality of the methodology is only limited by the capability of a combination of computational tools and the appropriateness of decompositions of the sort (wavelets) employed here. The methods and analyses we present are illustrated with a test bed dynamic stress analysis for a particular engineering system.
引用
收藏
页码:1874 / 1906
页数:33
相关论文
共 28 条
  • [1] Circuit optimization via sequential computer experiments: design of an output buffer
    Aslett, R
    Buck, RJ
    Duvall, SG
    Sacks, J
    Welch, WJ
    [J]. JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES C-APPLIED STATISTICS, 1998, 47 : 31 - 48
  • [2] BAYARRI M, 2006, 165 NAT I STAT SCI
  • [3] A framework for validation of computer models
    Bayarri, Maria J.
    Berger, James O.
    Paulo, Rui
    Sacks, Jerry
    Cafeo, John A.
    Cavendish, James
    Lin, Chin-Hsu
    Tu, Jian
    [J]. TECHNOMETRICS, 2007, 49 (02) : 138 - 154
  • [4] BAYARRI MJ, 2005, 163 NAT I STAT SCI
  • [5] CAFEO J, 2001, FRAMEWORK VERIFICATI
  • [6] Craig P. S., 1997, Case Studies in Bayesian Statistics, V3, P36
  • [7] Bayesian forecasting for complex systems using computer simulators
    Craig, PS
    Goldstein, M
    Rougier, JC
    Seheult, AH
    [J]. JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION, 2001, 96 (454) : 717 - 729
  • [8] BAYESIAN PREDICTION OF DETERMINISTIC FUNCTIONS, WITH APPLICATIONS TO THE DESIGN AND ANALYSIS OF COMPUTER EXPERIMENTS
    CURRIN, C
    MITCHELL, T
    MORRIS, M
    YLVISAKER, D
    [J]. JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION, 1991, 86 (416) : 953 - 963
  • [9] EASTERLING RG, 2001, SAND20010243
  • [10] Fuentes M, 2003, INT STAT REV, V71, P201