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)
}
このコードを乾かすには?