CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesce Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDescU TextModels.StdModelDescTextModels.ModelDesc  TextModels.AttributesDesc1$Courier New ;$Courier Newb .**uTTextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescTextRulers.AttributesDescL Zo Z%$.6?HLQ*uTgL Zo Z%$.6?HLQ =*:~*uTgL Zo Z%$.6?HLQ#------------------------------------------------------------------------------------------------------------------------------ # Political Party data # Data takes from Agrest (2002), page 303, Table 7.15, problem 7.3 #------------------------------------------------------------------------------------------------------------------------------ # Model 1: USING MULTI(P,N) #------------------------------------------------------------------------------------------------------------------------------ model{ # model's likelihood for (i in 1:n){ for( k in 1:K ) { # linear predictors eta[i,k] <- beta[1,k] + beta[2,k]*gender[i] + beta[3,k]*race[i] expeta[i,k]<-exp(eta[i,k]) # probabilities (link function) p[i,k] <- expeta[i,k]/sum(expeta[i,1:K]) } N[i] <- sum( y[i,1:K] ) # stochastic part y[i,1:K] ~ dmulti( p[i,1:K], N[i] ) } # for (j in 1:P){ # coefficients for the baseline category are constrained to zero beta[j,1] <- 0.0 # independent normal low information priors for (k in 2:K){ beta[j,k] ~ dnorm( 0.0, 0.001) } } } INITS list( beta=structure(.Data=c(NA, 0, 0, NA, 0, 0, NA, 0, 0), .Dim = c(3, 3))) DATA (LIST) list( n=4, K=3, P=3, gender = c(1, 1, 0, 0), race = c(1, 0, 1, 0), y = structure( .Data=c(132, 42, 172, 56, 176, 6, 129, 4, 127, 12, 130, 15), .Dim = c(4, 3) ) )   node mean sd MC error 2.5% median 97.5% start sample beta[1,2] 1.361 0.1758 0.006274 1.025 1.36 1.703 5001 5000 beta[1,3] -1.158 0.2577 0.01392 -1.687 -1.147 -0.6838 5001 5000 beta[2,2] 0.249 0.1981 0.005818 -0.1405 0.2482 0.6359 5001 5000 beta[2,3] -0.07157 0.1597 0.004306 -0.3789 -0.0712 0.2468 5001 5000 beta[3,2] -3.26 0.209 0.006158 -3.676 -3.259 -2.857 5001 5000 beta[3,3] 1.333 0.2553 0.01367 0.8502 1.325 1.85 5001 5000  TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc aY?$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*