不知道什么原因,个人感觉已经都赋予初始值了,为什么在加载完后 还会出现 ”this chain contains uninitialized values“提示呢?望大神指导,不胜感激
代码如下:
model {
for (i in1:n) {
for (j in 1:K){
y[i,j] ~ dnorm( mu[i,j], tau )
mu[i,j] <- m + a
}
a~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, a=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0))