rjags - 使用 rjags 进行贝叶斯模型选择

Question

模型图像

我正在尝试使用 rjags 对图像文件中的模型进行 Gibbs 采样。我写的模型如下。

model {
  for (i in 1:N) {
    y[i] ~ dnorm(mu[i], lambda)
    mu[i] <- ((c(1, gamma)*beta)%*%X[i,])
  }
  for (i in 1:(p+1)) {
    beta[i] ~ dnorm(bhat[i], 0.00001)
  }
  for (i in 1:p) {
    gamma[i] ~ dbern(eta)
  }
  eta ~ dbeta(1,1)
  lambda ~ dgamma(v0/2, v0*s0/2)
}

但是，在我获得的样本中，伽玛值都为零。有没有办法修复我的模型以获得不同的结果？

标签： rjags