CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc1$Courier NewR $Courier Newb #------------------------------------------------------------------------------------------------------------------------------ Chapter 7 Example 1 Poisson data Aircraft damage dataset (Montogomery etal 2006) #------------------------------------------------------------------------------------------------------------------------------ model{ # Poisson model likelihood for (i in 1:30){ damage[i] ~ dpois( lambda[i] ) log(lambda[i]) <- beta[1] + beta[2] * type[i] + beta[3] * bombload[i] + beta[4] * airexp[i] } # # prior for (j in 1:4){ beta[j]~dnorm( 0.0, 0.001 ) B[j] <- exp( beta[j] ) } # # profiles # values for bombload profiles[1,1] <- ranked( bombload[], 1 ) # minimum of bombload profiles[2,1] <- mean(bombload[]) # mean of bombload profiles[3,1] <- 0.5*( ranked( bombload[], 15 )+ranked( bombload[], 16 )) #median profiles[4,1] <- ranked( bombload[], 30 ) #max # values for airexp profiles[1,2] <- ranked( airexp[], 30 ) #max experience profiles[2,2] <- mean(airexp[]) #mean profiles[3,2] <- 0.5*( ranked( airexp[], 15 )+ranked( airexp[], 16 )) #median profiles[4,2] <- ranked( airexp[], 1 ) #min experience for (k in 1:4){ a4.profile[k] <- exp( beta[1] + beta[3]*profiles[k,1] + beta[4]*profiles[k,2] ) a6.profile[k] <- a4.profile[k]*B[2] # this is equivalent to setting exp( beta[1] + beta[3]*profile[k,1] + beta[4]*profile[k,2] ) } } INITS list( beta=c(0,0,0,0) ) DATA (LIST) list( 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) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*