CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc5 - TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc'* *&R#------------------------------------------------------------------------------------------------------------------------------ Version 1: hierarchical #------------------------------------------------------------------------------------------------------------------------------ model { for (i in 1:n) { for (j in 1:K){ y[i,j] ~ dnorm( mu[i,j], tau ) mu[i,j] <- m + a[i] } a[i]~dnorm( 0, tau.a) } # prior distributions m ~ dnorm( 0.0, 0.001) tau~dgamma(0.001, 0.001) tau.a~dgamma(0.001, 0.001) s2 <- 1/tau s2.a <- 1/tau.a ts2 <- s2+s2.a cor <- s2.a/ts2 s <- sqrt(s2) s.a <- sqrt(s2.a) } DATA list(n=20, K=2, y=structure(.Data=c(108, NA, 91, NA, NA, NA, 104, NA, 99, 97, 95, 98, 93, 97, 99, 96, 90, 100, 92, 95, 101, 89, 97, 97, 97, 100, 96, 95, 106, 100, 100, 98, 90, 99, 88, 98, 92, 92, 100, 101), .Dim = c(20, 2) ) ) INITS list(m=0.0, tau=1.0, tau.a=1.0) #------------------------------------------------------------------------------------------------------------------------------ #Version 3: simple no randome effects #------------------------------------------------------------------------------------------------------------------------------ model { for (i in 1:n) { for (j in 1:K){ y[i,j] ~ dnorm( m, tau ) } } # prior distributions m ~ dnorm( 0.0, 0.001) tau~dgamma(0.001, 0.001) s2 <- 1/tau s <- sqrt(s2) } INITS list(m=0.0, tau=1.0) #------------------------------------------------------------------------------------------------------------------------------ #Version 4: fixed effects #------------------------------------------------------------------------------------------------------------------------------ model { for (i in 1:n) { for (j in 1:K){ y[i,j] ~ dnorm( mu[i,j], tau ) mu[i,j] <- m + a[i] } a[i]~dnorm( 0, 0.001) } # prior distributions m ~ dnorm( 0.0, 0.001) tau~dgamma(0.001, 0.001) s2 <- 1/tau s <- sqrt(s2) } INITS list(m=0.0, tau=1.0) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*