winbugs:应为逗号管理员

Question

全!

我使用winbugs来做简单的线性回归。但是，系统总是给出错误消息，需要一个逗号。

以下是我的模型声明：

model {
   for (i in 1:I)
    {
      Z[i] ~ dnorm(beta0 + beta1 * X[i], tau)
     }
   tau <- 1/(sigma*sigma)
   sigma ~ dunif(0, 100)
   beta0 ~ dnorm(0, 1E-6)
   beta1 ~ dnorm(0, 1E-6)
 }

它有什么问题？谢谢。

Answer 1

与在dnorm(beta0 + beta1 * X[i], tau)中不同，WinBugs不支持参数分布的表达式。

你的问题的解决方案是

model {
  for (i in 1:I)
  {
      Z[i] ~ dnorm(mu[i], tau)
      mu[i]<- beta0 + beta1 * X[i]
  }
  tau <- 1/(sigma*sigma)
  sigma ~ dunif(0, 100)
  beta0 ~ dnorm(0, 1.0E-6)
  beta1 ~ dnorm(0, 1.0E-6)
 }