OpenBUGS 代码在生存中给出错误“预期逗号”

Question

我正在尝试使用 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 在tau 中突出显示tausg... 我做错了什么。 谢谢

Answer 1

您在两个地方指定了 mu - 一次是 mu[i,j] <- b[i]*exp(beta),，然后是 mu<-(-1)/2*tau。您需要更改其中一个的名称。