DAVE HIGDON*, MARC KENNEDYâ , JAMES C. CAVENDISH*, JOHN A. CAFEO§, AND. ROBERT D. RYNEЦ. Abstract. We develop a statistical approach for ...
!#"$%&'( !"$)*+ ',-& *$ .& %/ 0 #& 1+%+ 2(3547698;:=7@7A BDC37E(FHG76I@;@;6I2KJML BDN 3OC6KPQFOR;FI35476S@O2O:=P+87T/BDN >78O@U3RDFI3OVW6K>MX B(YSZI[ E(>7\I6IE^]_2)RDE5J?@D6K` acb^dfehghikjle m)n*opqosr outwvyx%z-{=|}zh|~w{|~wuzltWzlxqxqv zl}vyD}
zlzls|of~wf~wqcq/uof=|}zl~w |#~wx/ofpq~ws|~wvy/{O|qzh| zloczlp/o-;~ |.|qozl~wp.vl;zfvyxq/|of){~w qtzh|~wvyv pqoft R] /xq~wuzltwt B |qocuvyx//|of){~w /tzh|~wvy*fv+pqo v pqoftw{(z xq /{~wzlt{+{|of,zlqpof+/~wof{Dz){os|DvlW~wqxq+|{Kc{vyo;+/vDzlqp{xkoff~
ofp B vl|qof{Dqq qvD R 3 tw~w~ |ofpzlvyq |cvl7
oftwp$p
zh|}z%vy|qo|/oxq /{~wzlt;{+{|ofzlozurlzl~wtzlqtwo|v~w/vy!zlkvy/|c|qo qq qvu;~wqxq+|{(zlqpzltw{v?|v ~w/vyzlkvy/|K|qo;qqfof|zl~w |c|
zh|(~w{Kzl{{v+f~zh|ofpD~ |czM{~w)qtzh|~wvy
zl{ofp xqoup/~wf|~wvy R] qoMzlxqxqv zl}y~ r ofqofoMzltwtwvD{Kvyu qzlfqffvyofq| zl|~w~w q|%vyou zlqqpqf~wofq-|zl~w v |pqoft pqq~w/ocxq|/v#|{tw~w¡~wo ~ R |uzhzl|~wtwvy~wqq}{)zh|vy~wvyk|¢}q£ o q kof)vl;{~w)qtzh|~wvyq{|
zh| zl.¤o uzl~woupvy/|u£ pqzlq~w{puqofx
/zluqof=s|}zl~wk os||K~wof.of|q|oqoMvy{/~w{ of/rltzhzh||~w~wvyvycfxqv+p/v+o)fofzl{/{p|
|zhq| o/zl~wfof|tw
pqzl{t5|x/q o+{zl~wsu|zlqt5zl{=tD/
{=|ofoftwp£ p
zh|z#vy|/o|qo xq /{~wuzlt5{+{|of R ] qoDof{qt |~wqzl
zlt +{~w{S+~woftwpq{Ixqofpq~ws|~wvyq{Izlqp|qof~wKzl{{v f~zh|ofpcqqfof|zl~w |~wof{K;/~wtwo;zlffvyq |~w/ vy qt ¥ |~wxqtwo7{vy/fof{;vl^/qfou=|}zl~w | R n*o?q{oMz \ zu of{~zlvy /tzh|~wvyczlqpoft cvy < zlq{{~zlxqv+fof{{(v pqoftw{I|v v pqoftSqq qvu;*/qs|~wvyq{vl(|/ov+pqoft~wqxq/|{ R] qoof{|~wczh|~wvy*~w{Muzl~wofp*vy/|q{~w/ C zlyvur%}
zl~w C vy |o F zltwv R)] q~w{?os|qv pqvytwvyl%~w{Mzlxqxqtw~wofp#|v|Kvos¦/zlxqtwof{7z}qzlyofp.x
zl=|~wftwo zlufoftwof}zh|vy£§zlqpz {xkvl|Koftwpq~w/xqv fou{{ R ¨© ª «)¬ gl^dym G os Kvypq{uMuzltw~wq}zh|~wvy £§fvyxq/|ouMos¦+x¤of~wou |{£Wxqoup/~wf|zlq~wtw~ |¤£Wqqfof|zl~w |.+qzl/¥ |~
uzh|~wvy £ < zl/{{~zlxqv fof{{£/v+p/out¤ryzltw~wp
zh|~wvy £
{~w qtzh|vyK{f~wofqfo ac®°¯ dh±5b ² © jlecj ³´ikdudhµw¶§jyiqeyµ ¬ · dymO¸y¹ ï3O@O2E5J?@;6 ñ Î×WÌQÎhÙ¤ß,Ö Ùà.ߧÊ^ÎÌ+ËÙkØÔÎhÑܧÓqФâ ó ×WÌ ÓkßWߧÏÔÕÖ ÓkÎÔÙ¤Í9ÙkØÑuÎhÓkÎÔÕÑfÎhÔÖ+ÓkÏàÌ Îh×WÙ^ÜWÑ.ÎÙӤͧÓkÏÐ^ÑuÌ/Ñ æ ÔÍ Þ Ù¤Ï Þ ÔÍWÛÖ ÙàßWÊ^ÎÌ+ËcÑuÔà*ʧÏÓkÎÔ٤ͧÑ-ÒSÌ Û ÓkÍQÜ^Ì+Ö+Ó¤Ü^Ì/Ñ-ÏÓkÎÌ Ë/â ó Ù*ÜWÓ
Îh̤åà*ʧÖl×Ù¤ØSÎh×WÌÑuÎhÓkÎÔÕÑfÎhÔÖ+ÓkϧàÌ Î×WÙ^Ü^ÙÏÙÛ¤Ð*ËhÌ ÏÌ Þ Ó¤ÍÎ?ÎhÙ#ÎקÌÓkͧӤÏÐ^ÑÔÑMÙkØKÖ Ù¤àßWÊWÎÌ Ë ÑÔàÊWÏÓkÎÔ٤ͧÑÖ+ÓkÍ1ÒÌ#Ëh٤ʧۤ×WÏÐÒWËhÙ¤ò¤Ì+Í$ÔÍ ÎÙ.ÎקÌØÚÙ¤ÏÏÙ
èÔÍWÛÍWÙ¤Í§Ì ðWÖ ÏʧÑÔ Þ ÌÖ+Ó
ÎÌ+Û¤ÙËÔÌ+Ñ+â @ Ì ð^ßSÌ ËhÔàÌ+ÍÎlÓkÏÜWÌ+ÑÔÛÍ9ñ Ü^Ì ÎhÌ Ëhà.ÔÍWÔÍWÛÔÍWßWÊ^Î1ÑuÌ ÎuÎÔÍWÛ Ñ Ó
ÎèקÔÖl×,ÎhÙäÖ Ó¤ËËhÐÙ¤ÊWÎ Ó ÑÌ A ÊWÌ ÍÖ Ì%Ù¤ØÑuÔà*ÊWÏÕÓ
ÎhÔÙÍÜ^Ì/ÑuÔۤͧÑö",!.-%*,6L@J3",9K=Rq>*,=1D&YvV",!J8 * 9K=^A!?"%@E1;9)Rl\ ouyof B 3R ¾ R
2 zuD~wp B/2)Rq8 of}yofczl B+3)R+VKRqC.RqP ~ | B zl/p C.R n*of{| B ofpq{ RwBq> ¦+vyp $ q~ r of{~ |c¾of{{ B
ý U
8D:=7@BkG76S@O@D6I2KJ BqFK354O6I@;2O:=P/8B
FI3OVW6K>ï3O@O2E5J?@;6 posterior model uncertainty
3
3
2
2 y(x), η(x,θ)
y(x), η(x,θ)
prior uncertainty
1 0 −1
0 −1
−2 −3
1
0
0
.5
0.2
−2
1
θ
0.4
0.6
0.8
−3
1
0
0
.5
1
θ 0.2
0.4
x
3
3
2
2
1 0
−2
−2 0.4
0.6 x
0.8
1
0 −1
0.2
1
1
−1
0
0.8
calibrated prediction
y(x), ζ(x)
δ(x)
posterior model discrepancy
−3
0.6 x
0.8
1
−3
0
0.2
0.4
0.6 x
Å%ÆÈÇM É ÉÈË R=ÌOÍ/Î2Ï ÍP×Ñ!ØÖÍ/Õ&)1K3*E:,G7*? 7$'C8 1;9K6",!J94",1D&L",=",8 Q,G 9KG}BkN vyqn1~wtwos 8æP vyq{ B ¼y¹y¹y¼ R ð ¼ Ë» ñc É cÉ ±þ ªú ö%môHÿHýbNÉ ô"úM#ô^B q>*,=1D&YvV",!J8 *OG 1D",1;9KG 1;94@ ",8Z$L&E1435*):,G}BP x/~wqyof¥ 4 oftzl : q B^º ::: R ð ¼y¸ ñmø7fÉ òmZÉ ±þ öÆ ÿ#÷PÉ RôHÿýHú ö>Æô¯ôHö%y÷¬ô¯ÿHýQÉ ;ú¼ÿ Höô^B 7G Q,$BA 1;*%1
