CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescu Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesceTextModels.StdModelDescTextModels.ModelDescogTextModels.AttributesDesc'*H *-J # --------------- NegBin: GSS1990 DATA ------------------------------------ # --------------- SIMPLER APPROACH ------------------------------------ # ---------------- grouped data --------------------------------------------------- model{ for (i in 1:K){ n[i,1] <- males[i] n[i,2] <- females[i] } C <- 0 for(i in 1:K){ for (j in 1:2){ zeros[i,j] <- 0 zeros[i,j] ~ dpois( zeros.mean[i,j] ) zeros.mean[i,j] <- -l[i,j]*n[i,j] + C # log-likelihood for i individual lambda.star[i,j] <- (1-omega[ j ] )*lambda[ j ] + omega[ j ]*y[i] l[i,j] <- log( (1-omega[ j ])*lambda[ j ] ) + (y[i]-1)*log(lambda.star[i,j]) - loggam(y[i]+1) - lambda.star[i,j] }} beta[1] <- log(lambda[1]) beta[2] <- log(lambda[2]/lambda[1]) expbeta2<- lambda[2]/lambda[1] for (j in 1:2){ omega[j]~dbeta(0.5,0.5) lambda[j] ~ dgamma( 0.001, 0.001) di[j] <-1/pow( 1-omega[j], 2) var[j] <- lambda[j]*di[j] zeta[j] <- lambda[j]*(1-omega[j]) } } INITS list( lambda=c(1, 1), omega=c(0.5, 0.5) ) DATA list( K=27, y = c(0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, 18, 20, 22, 23, 24, 25, 27, 30, 50, 60), males = c(65, 11, 13, 14, 26, 13, 15, 7, 21, 2, 24, 6, 3, 0, 3, 3, 0, 0, 7, 0, 0, 1, 1, 0, 3, 1, 1), females = c(128, 17, 23, 16, 19, 17, 17, 3, 15, 2, 13, 10, 3, 1, 10, 1, 1, 1, 6, 1, 1, 0, 3, 1, 1, 0, 0) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ Zo * ,[ @Documents.ControllerDesc Ws,! [h$