尽管我的 OpenBugs 代码没有任何问题,但我不断收到相同的错误消息。
包含先验的模型如下所示:
model{
for (i in 1:n) {
y[i] ~ dnorm(mu[i], tau)
mu[i] <- z1
}
## Priors
z1 ~ dnorm(200, 0.0025)
beta ~ dnorm(5, 0.0001)
p0 ~ dunif(0, 100)
Limit[1] ~ dnorm(0, 0.0001)
Limit[2] ~ dnorm(0, 0.0001)
ranked_limits <- sort(Limit)
pi_1<- ranked_limits[1]
pi_2 <- ranked_limits[2]
z_hat ~ dnorm(133, 0.0004)
tau ~ dgamma(0.001, 0.001)
sigma2 <- 1/tau
}
每次我尝试检查模型时,它都会在第 16 行返回错误“Expected Right Parenthesis”。
该行的代码行是:
pi_1 <- ranked_limits[1]
我收到此错误是否有任何特殊原因。
通过进行一些测试,当我注释掉该行时,下面的行返回相同的错误。如果我随后也注释掉该行,则下面的行将返回错误。
我在这里缺少什么吗?