1

JAGS マニュアルを読みましたが、同じ事前分布を JAGS / R2JAGS モデルの複数のパラメーターに割り当てる方法が見つかりませんでした。

たとえば、現在、次のような多くのコードを繰り返さなければなりません。

reg.model <- function() {
  # Model structure
  for(i in 1:N){
    Y[i] ~ dnorm(mu[i], phi)
    mu[i] <- beta0 + beta1*x1[i] + beta2*x2[i]
  }
  sigma2 <- 1 / phi

  # Priors (Lots of re-typing here for beta0, beta1, beta2)
  phi   ~ dgamma(1,1)
  beta0 ~ dnorm(0, 0.01*phi)
  beta1 ~ dnorm(0, 0.01*phi)
  beta2 ~ dnorm(0, 0.01*phi)
}

このコードを乾かすには?

4

1 に答える 1