我在用winbugs标定随机系数模型，但是在模型compile时不停的报multiple definitions of node Y.不知道是什么原因啊。我的模型如下
model
{
for(p in 1 : N) {
Y[p] ~ dnorm(mu[p], tau[p])
mu[p] <- (u[1]+alpha[station[p], 1] )+ (u[2]+alpha[station[p],2]) *XB[p]+(u[3]+alpha[station[p], 3] ) *XC[p]+(u[4]+alpha[station[p], 4]) *XD[p]+(u[5]+alpha[station[p], 5] ) *XE[p]+(u[6]+alpha[station[p], 6]) *XF[p]+(u[7]+alpha[station[p], 7] ) *XG[p]+(u[8]+alpha[station[p], 8]) *XH[p]+(u[9]+alpha[station[p], 9] )*XI[p]
sigma2[p] <- 1 / tau[p]
tau[p] ~ dgamma(1,1)
}
# Priors for random coefficients:
for (j in 1 : M) {
alpha[j, 1:9 ] ~ dmnorm(gamma[1:9 ], T[1:9 ,1:9 ]);
alpha1[j] <- alpha[j,1]
# Priors for fixed coefficients:
for (k in 1 : 9) { u[k] ~ dnorm(0.0, 0.0001) }
}
# Hyper-priors:
gamma[1:9] ~ dmnorm(mn[1:9 ], prec[1:9 ,1:9 ]);
T[1:9 ,1:9 ] ~ dwish(R[1:9 ,1:9 ], 9)
Sigma.alpha[1 : 9,1 : 9] <- inverse(T[1:9,1:9 ])
}
另外，我在学习winbugs时觉得很痛苦，因为没有什么参考资料，每次报错也不知道错误的原因是什么，只能上google海搜，现在两个月了这个模型还没有标定出来，觉得学习的很没有效率~不知道大家是怎么学习winbugs的，使用了什么参考教材，大家讨论一下呗~
