CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDescb Z TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc'*Y *#------------------------------------------------------------------------------------------------------------------------------ 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) # # calculation of residuals for (i in 1:n) { for (j in 1:K){ res[i,j]<- y[i,j]-mu[i,j] }} # Calculation of R2 R2 <- 1 - pow( sd(res[1:n,1:K])/sd(y[1:n,1:K]), 2 ) } DATA list(n=20, K=2, y=structure(.Data=c(108, 98, 91, 94, 93, 96, 104, 99, 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 2: multivariate #------------------------------------------------------------------------------------------------------------------------------ model { for (i in 1:n) { y[i, 1:K] ~ dmnorm( mu[i, 1:K], T[1:K, 1:K] ) for (j in 1:K){ mu[i,j] <- m } } T[1:K,1:K] <- inverse( S[1:K, 1:K] ) ts2 <- s2+s2.a S[1,1] <- ts2 S[2,2] <- ts2 S[1,2] <- s2.a S[2,1] <- s2.a s2 <- 1/tau s2.a <- 1/tau.a # prior distributions m ~ dnorm( 0.0, 0.001) tau~dgamma(0.001, 0.001) tau.a~dgamma(0.001, 0.001) cor <- s2.a/ts2 s <- sqrt(s2) s.a <- sqrt(s2.a) } TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*