CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescx Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesch TextModels.StdModelDescTextModels.ModelDescrjTextModels.AttributesDesc1$Courier New  $Courier Newo #------------------------------------------------------------------------------------------------------------------------------ Version 1 #------------------------------------------------------------------------------------------------------------------------------ model{ # model's likelihood for (a in 1:LA){ for (b in 1:LB){ for (k in 1:K){ y[a,b,k] ~ dlnorm( mu[a,b,k], tau ) mu[a,b,k] <- mu0 + econ[a] + gender[b] + econ.gender[a,b] #mu[a,b,k]<- mu0 + econ[a] + gender[b] }}} #### 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( LA=3, LB=2, K=2, y = structure( .Data=c( 9, 14, 18, 29, 25, 26, 22, 25, 23, NA, 12, NA), .Dim = c( 3,2,2 ) ) ) DATA (RECT FORMAT) list( LA=3, LB=2, K=2) y[,1,1] y[,1,2] y[,2,1] y[,2,2] 9 14 18 29 25 26 22 25 23 NA 12 NA END TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*