CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc_ Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDescO TextModels.StdModelDescTextModels.ModelDescYQTextModels.AttributesDesc1$Courier New) $Courier New "#------------------------------------------------------------------------------------------------------------------------------ Version 1 - data are loaded in log-scale #------------------------------------------------------------------------------------------------------------------------------ model{ # model's likelihood for (a in 1:LA){ for (b in 1:LB){ for (k in 1:K){ y[a,b,k] <- exp( logy[a,b,k] ) logy[a,b,k] ~ dnorm( 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 ) # # # the data are directly supplied in log-scale DATA (LIST) list( LA=3, LB=2, K=2, logy = structure( .Data=c(2.19722457733622, 2.63905732961526, 2.89037175789616, 3.36729582998647, 3.2188758248682, 3.25809653802148, 3.09104245335832, 3.2188758248682, 3.13549421592915, NA, 2.484906649788, NA), .Dim = c( 3,2,2 ) ) ) DATA (RECT FORMAT) list( LA=3, LB=2, K=2) logy[,1,1] logy[,1,2] logy[,2,1] logy[,2,2] 2.197225 2.639057 2.890372 3.367296 3.218876 3.258097 3.091042 3.218876 3.135494 NA 2.484907 NA END TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*