06.12.2012 Views

Bayesian Experimental Design - Mathematical Sciences Home Pages

Bayesian Experimental Design - Mathematical Sciences Home Pages

Bayesian Experimental Design - Mathematical Sciences Home Pages

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

D-andE-optimality,buttheconditionsunderwhichtheaboveholdsdonotseemeasyto satisfy.Pilzdidnotgiveexplicitdesignsandexaminetheirpracticalimplicationsandhis<br />

estimation.Inqualitycontrolandinclinicaltrialspredictionoffutureobservationscanbe workissomewhatabstract. 2.4OtherUtilityFunctions ofspecialinterest.Inthesecasesthe<strong>Bayesian</strong>approachusespredictiveanalysiswhichcan alsobehelpfulindesigningtheexperiment.TheexpectedgaininShannoninformationon afutureobservationyn+1isusedratherthantheexpectedgainininformationonthevector Asnotedinsection1,incertainexperiments,predictioncanbemoreimportantthan<br />

p(yn+1)(priorpredictive)onyn+1istheequivalentofthequantity(3)insection2.2.The p(yn+1jy;)=Rp(yn+1j)p(jy;)d(posteriorpredictive)andthemarginaldistribution ofparameters.TheexpectedKullback-Leiblerdistancebetweenthepredictivedistribution<br />

theexpectedutility:U3()=Zlogp(yn+1jy;)p(y;yn+1j)dydyn+1: priorpredictivedistributiondoesnotdependonthedesignandthedesignthatmaximizes theexpectedgaininShannoninformationonyn+1isequivalenttothedesignthatmaximizes<br />

experiments.Inthenormallinearmodel,maximizingU3()withrespecttocorresponds ThisutilityfunctionhasbeenusedbySanMartiniandSpezzaferri(1984)foramodel selectionproblemandbyVerdinelli,PolsonandSingpurwalla(1993)foracceleratedlifetest (9)<br />

wherethenextobservationisgoingtobetakenatthepointxn+12X.Thisisequivalent tominimizingthepredictivevariance tomaximizing ?12nlog(2)+1+logh2xTn+1D()xn+1+2io;<br />

Inthespecialcaseofpredictionofyn+1ataxedpointc=xn+1,thedesignmaximizing 2n+1=2[xTn+1D()xn+1+1]: 16

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!