TimeT 发表于 2014-3-4 22:33 
首先,我声明我没读你的模型(抱歉,太多东西了)。不过对你的问题invalid or unexpected token sceanned倒 ...
非常感谢你啊,我按照你说的方法改了 ,现在显示expected comma,这又是什么意思啊?模型现在修改如下:
model
{
for( i in 1 : N ) {
for( j in 1 : T ) {
Y[i , j]<-mu0+r0[ i]+(mu1+r1
) * x1[ i,j]+(mu2+r2[i ]) * x2[i,j]+(mu3+r3[ i]) *x3[i,j]+(mu4+r4)*x4[i,j]+(mu5+r5)*x5[i,j]+(mu6+r6)*x6[i,j]+(mu7+r7)*x7[i,j]+(mu8+r8)*TIME[i,j]+e
r[r0,r1,r2,r3,r4,r5,r6,r7,r8] ~ dmnorm(0.0 ,T[ , ])
e~dnorm(mu.e,prec.e)
}
T[ , ] ~dwish (R[ , ],v)
}
mu.e<-0.0
prec.e~dgamma(0.001,0.001)
sigma.e <- 1 / sqrt(prec.e)
mu0 ~ dnorm(0,1.0E-6)
mu1 ~ dnorm(0,1.0E-6)
mu2 ~ dnorm(0,1.0E-6)
mu3 ~ dnorm(0,1.0E-6)
mu4 ~ dnorm(0,1.0E-6)
mu5 ~ dnorm(0,1.0E-6)
mu6 ~ dnorm(0,1.0E-6)
mu7 ~ dnorm(0,1.0E-6)
mu8 ~ dnorm(0,1.0E-6)
}
我这是一个随机系数模型。另外,我还有一些问题想请教您。如果是多个变量服从多元正态分布,那么这个多元正太分布的均值是一个标量,还是向量啊?还有我这行代码r[r0,r1,r2,r3,r4,r5,r6,r7,r8] ~ dmnorm(0.0 ,T[ , ]),想表示的是r0到r8这9个变量服从多元正态分布,为什么总提示expected comma?
另外我想问问有没有什么比较好的参考书啊,我现在看着全英文的感觉比较吃力,老师又催的很紧。