CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesck c TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDesc  TextModels.AttributesDesc1$Courier New $Courier New #------------------------------------------------------------------------------------------------------------------------------ Version 3- three way simple #------------------------------------------------------------------------------------------------------------------------------ model{ # CR dummy variables for (i in 1:n){ D.gender[i]<-equals(g[i],2) D.econ2[i]<-equals(e[i],2) D.econ3[i]<-equals(e[i],3) D.city[i]<-equals(ci[i],2) } # STZ dummy variables #for (i in 1:n){ # Dstz.gender[i]<-equals(g[i],2)-equals(g[i],1) # Dstz.econ2[i]<-equals(e[i],2)-equals(e[i],1) # Dstz.econ3[i]<-equals(e[i],3)-equals(e[i],1) # Dstz.city[i]<-equals(ci[i],2)-equals(ci[i],1) #} # # model's likelihood for (i in 1:n){ y[i] ~ dnorm( mu[i], tau ) mu[i] <- mu0 + gender * D.gender[i] + econ[2]*D.econ2[i] + econ[3]*D.econ3[i] + city*D.city[i] + gender.econ[2]*D.gender[i]*D.econ2[i] + gender.econ[3]*D.gender[i]*D.econ3[i] + gender.city*D.gender[i]*D.city[i] + econ.city[2]*D.econ2[i]*D.city[i] + econ.city[3]*D.econ3[i]*D.city[i] + gender.econ.city[2]*D.gender[i]*D.econ2[i]*D.city[i] + gender.econ.city[3]*D.gender[i]*D.econ3[i]*D.city[i] } #### set zero all nuisance parameters econ[1] <- 0.0 econ.city[1]<-0.0 gender.econ[1]<-0.0 gender.econ.city[1]<-0.0 # priors mu0~dnorm( 0.0, 1.0E-04) gender ~ dnorm( 0.0, 1.0E-04) city ~ dnorm( 0.0, 1.0E-04) gender.city ~ dnorm( 0.0, 1.0E-04) for (k in 2:3){ econ[k] ~ dnorm( 0.0, 1.0E-04) gender.econ[k] ~ dnorm( 0.0, 1.0E-04) econ.city[k] ~ dnorm( 0.0, 1.0E-04) gender.econ.city[k] ~ dnorm( 0.0, 1.0E-04) } tau ~dgamma( 0.01, 0.01) s <- sqrt(1/tau) # precision } INITS list( mu0=1.0, econ=c(NA, 0,0), gender=0, city=0, gender.city=0, gender.econ=c(NA,0,0), econ.city=c(NA,0,0), gender.econ.city=c(NA,0,0), tau=1.0) DATA (LIST) list(n=10) g[] e[] ci[] y[] 1 1 1 9 1 2 1 25 1 3 1 23 2 1 1 18 2 2 1 22 2 3 1 12 1 1 2 14 1 2 2 26 2 1 2 29 2 2 2 25 END TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*