Байесовская подборка моделей с рьягами

Question

Изображение модели

Я пытаюсь выполнить выборку Гиббса с моделью в файле изображения, используя rjags.Я написал модель следующим образом.

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)
}

Тем не менее, все гаммы равны нулю в образце, который я получил.Есть ли способ исправить мою модель, чтобы получить другие результаты?