贝叶斯模型与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)
}

但是，在我获得的样本中，伽玛都为零。 有没有一种方法可以修正我的模型以获得不同的结果？