我正在尝试使用 OpenBUGS 拟合分层模型,代码如下:
model
{
for (i in 1:N) {
for (j in 1:M){
# Survival times bounded below by censoring times:
t[i,j] ~ dweib(r,mu[i,j]) I(cen[i,j],);
mu[i,j] <- b[i]*exp(beta);
}
# Random effects:
b[i] ~dlnorm(mu, tau)
}
# Priors:
beta ~ dnorm(0, 0.001);
tau ~ dgamma(1.0E-3, 1.0E-3);
r~ dgamma(1, 1.0E-3);
sigma2<-1/tau
mu<-(-1)/2*tau
}
但是,我得到了
错误 : 'expected a comma'
,并且 OpenBUGStausg
在tau
...中突出显示了我做错了什么。谢谢