CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesce ] TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDesc TextModels.AttributesDesc1$Courier NewB $Courier New0 #------------------------------------------------------------------------------------------------------------------------------ Version 4- three way using matrices and vectors model gender*econ #------------------------------------------------------------------------------------------------------------------------------ model{ # # Creating the design matrix X for (i in 1:n){ X[i,1]<- 1.0 # constant term # CR dummy variables X[i,2]<-equals(g[i],2) # gender X[i,3]<-equals(e[i],2) # econ2 X[i,4]<-equals(e[i],3) # econ3 X[i,5]<-equals(ci[i],2) # city # STZ dummy variables # X[i,2]<-equals(g[i],2)-equals(g[i],1) # X[i,3]<-equals(e[i],2)-equals(e[i],1) # X[i,4]<-equals(e[i],3)-equals(e[i],1) # X[i,5]<-equals(ci[i],2)-equals(ci[i],1) # specification of interaction terms X[i,6]<-X[i,2]*X[i,3] # gender*econ2 X[i,7]<-X[i,2]*X[i,4] # gender*econ3 #X[i,8]<-X[i,2]*X[i,5] # gender*city #X[i,9]<-X[i,3]*X[i,5] # econ.city2 #X[i,10]<-X[i,4]*X[i,3] # econ.city3 #X[i,11]<-X[i,2]*X[i,3]*X[i,5] # gender.econ.city2 #X[i,12]<-X[i,2]*X[i,4]*X[i,5] # gender.econ.city3 } # # model's likelihood for (i in 1:n){ y[i] ~ dnorm( mu[i], tau ) mu[i] <- inprod( X[i,], beta[] ) } # priors for (j in 1:4){ beta[j]~dnorm( 0.0, 1.0E-04) } beta[5]<-0.0 for (j in 6:P){ beta[k]~dnorm( 0.0, 1.0E-04) } tau ~dgamma( 0.01, 0.01) s <- sqrt(1/tau) # precision } INITS list( beta=c(0,0,0,0,0,0,0,0,0,0,0,0), tau=1.0) DATA (LIST) list(n=10, P=7) 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*