CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescF Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesc6 TextModels.StdModelDescTextModels.ModelDesc@8TextModels.AttributesDesc1$Courier NewR $Courier New3 Q#------------------------------------------------------------------------------------------------------------------------------ Chapter 8 Example 1 Poisson data Aircraft damage dataset (Montogomery etal 2006) #------------------------------------------------------------------------------------------------------------------------------ model{ # vectors used for the calculation of DIC for (i in 1:n){ y1[i] <- damage[i] y2[i] <- damage[i] y3[i] <- damage[i] y4[i] <- damage[i] y5[i] <- damage[i] y6[i] <- damage[i] y7[i] <- damage[i] y8[i] <- damage[i] y1[i]~dpois( lambda[i,1] ) y2[i]~dpois( lambda[i,2] ) y3[i]~dpois( lambda[i,3] ) y4[i]~dpois( lambda[i,4] ) y5[i]~dpois( lambda[i,5] ) y6[i]~dpois( lambda[i,6] ) y7[i]~dpois( lambda[i,7] ) y8[i]~dpois( lambda[i,8] ) } # Poisson model likelihood for (k in 1:8){ for (i in 1:n){ log(lambda[i,k]) <- beta[1,k] + gamma1[k]*beta[2,k] * type[i] + gamma2[k]*beta[3,k] * bombload[i] + gamma3[k]*beta[4,k] * airexp[i] }} # # prior for (k in 1:8){ for (j in 1:4){ beta[j,k]~dnorm( 0.0, 0.001 ) } } # } INITS list( beta = structure( .Data= c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), .Dim = c(4, 8) ) ) DATA (LIST) list( n=30, damage = c(0, 1, 0, 0, 0, 0, 1, 0, 0, 2, 1, 1, 1, 1, 2, 3, 1, 1, 1, 2, 0, 1, 1, 2, 5, 1, 1, 5, 5, 7), type = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), bombload = c(4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 7, 7, 7, 10, 10, 10, 12, 12, 12, 8, 8, 8, 14, 14, 14), airexp = c(91.5, 84, 76.5, 69, 61.5, 80, 72.5, 65, 57.5, 50, 103, 95.5, 88, 80.5, 73, 116.1, 100.6, 85, 69.4, 53.9, 112.3, 96.7, 81.1, 65.6, 50, 120, 104.4, 88.9, 73.7, 57.8), gamma1=c(0,1,0,1,0,1,0,1), gamma2=c(0,0,1,1,0,0,1,1), gamma3=c(0,0,0,0,1,1,1,1) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*