CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDescTextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc1$Courier New  $Courier NewI #------------------------------------------------------------------------------------------------------------------------------ Version 2 #------------------------------------------------------------------------------------------------------------------------------ model{ # model's likelihood for (i in 1:n){ y[i] ~ dnorm( mu[i], tau ) mu[i] <- mu0 + econ[ e[i] ] + gender[ g[i] ] + econ.gender[e[i],g[i]] # mu[a,b,k] <- mu0 + econ[ e[i] ] + gender[ g[i] ] } #### CR Constraints econ[1] <- 0.0 gender[1] <- 0.0 econ.gender[1,1] <-0.0 for (a in 2:LA){econ.gender[a,1]<-0.0} for (b in 2:LB){econ.gender[1,b]<-0.0} # priors mu0~dnorm( 0.0, 1.0E-04) for (a in 2:LA){econ[a]~dnorm( 0.0, 1.0E-04)} for (b in 2:LB){gender[b]~dnorm( 0.0, 1.0E-04)} for (a in 2:LA){ for (b in 2:LB){ econ.gender[a,b]~dnorm( 0.0, 1.0E-04) }} tau ~dgamma( 0.01, 0.01) s <- sqrt(1/tau) # precision for (a in 1:LA){ for (b in 1:LB){ mean.spq[a,b] <- mu0 + econ[a] + gender[b] + econ.gender[a,b] }} } INITS list( mu0=1.0, econ=c(NA, 0,0), gender=c(NA, 0), econ.gender=structure(.Data=c(NA, NA, NA, 0, NA, 0), .Dim=c(3,2)), tau=1.0 ) DATA (LIST) list(n=10, LA=3, LB=2) g[] e[] y[] 1 1 9 1 2 25 1 3 23 2 1 18 2 2 22 2 3 12 1 1 14 1 2 26 2 1 29 2 2 25 END TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*