我正在尝试使用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'
,OpenBUGS在tausg
中突出显示tau
...
我做错了什么。
感谢
答案 0 :(得分:0)
您已在两个地方指定了mu - 一次为mu[i,j] <- b[i]*exp(beta),
,之后为mu<-(-1)/2*tau
。您需要更改其中一个的名称。