CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescDocuments.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesc}TextModels.StdModelDescTextModels.ModelDesc TextModels.AttributesDesc'*+ '* '* '*PS .*PSxStdFolds.FoldDesc  o6 *PS/Click on one of the arrows to display the modelx0 x  n6 *PS.Click on one of the arrows to display the datax0 xx r T6 *PSlist(tau=1,b=c(0,0))8xx: x  o6 *PS/Click on one of the arrows to display the model j*PS    Nx0 xx  6 *PSXlist(x=c(15,26,10,9,15,20,18,11,8,20,7,9,10,11,11,10,12,42,17,11,10), y=c(95,71,83,91,102,87,93,100,104,94,113,96,83,84,102,100,105,57,121,86,100), CPO=c(0.035063,0.020899,0.010772,0.024637,0.024444,0.035249,0.03367, 0.03453,0.033512,0.028571,0.019175,0.033113,0.010772,0.014934, 0.032342,0.035137,0.025132,0.014975,0.000138,0.019113,0.035137)).xx0 xx s U6 *PS list(tau=1,b=c(0,0))8xx  3**uTTextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDesczTextRulers.AttributesDescD ZGo Z%$.6?Homework 5 Stat 544 - Spring 2005 Problem 1 Part 1 Modelmodel { for (i in 1:21) {y[i] ~ dnorm(m[i],tau) m[i] <- b[1]+b[2]*x[i] r[i] <- y[i]-m[i] f[i] <- sqrt(tau/6.28)*exp(-0.5*tau*pow(r[i],2)) # Using R take inverses of average h to get estimate of CPO k[i] <- 1/f[i]} for (j in 1:2) {b[j] ~ dnorm(0,0.0000001)} tau ~ dgamma(1,0.001) }  Data list(x=c(15,26,10,9,15,20,18,11,8,20,7,9,10,11,11,10,12,42,17,11,10), y=c(95,71,83,91,102,87,93,100,104,94,113,96,83,84,102,100,105,57,121,86,100))  Inits Click on one of the arrows to display the initial values Part 2-4 Model model { for (i in 1:21) {y[i] ~ dnorm(m[i],tau) m[i] <- b[1]+b[2]*x[i] r[i] <- y[i]-m[i] f[i] <- sqrt(tau/6.28)*exp(-0.5*tau*pow(r[i],2)) a[i] <- CPO[i]/f[i] # divergence criteria g1[i] <- a[i]*log(a[i]) # Kullback 1 g2[i] <- -log(a[i]) # Kullback 2 g3[i] <- 0.5*abs(a[i]-1) # L1 norm g4[i] <- pow(a[i]-1,2)} # chi2 norm # priors for (j in 1:2) {b[j] ~ dnorm(0,0.0000001)} tau ~ dgamma(1,0.001) s2 <- 1/tau} Data Click on one of the arrows to display the data Inits Click on one of the arrows to display the initial values TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc aY?$ ZGo * ,[ @Documents.ControllerDesc @[s ^h