WinBugs はdnorm
分散の代わりに精度をパラメーターとして使用することを知っています
model {
#Likelihood
for(i in 1:N1){
y1[i] ~dnorm(mu,tau)
}
sigma <- sqrt(1/tau)
#Priors
mu ~ dnorm(0,0.000001)
tau~ dgamma(taumu, taus)
}
私の質問はsigma
、次のモデルを使用するのが正しいでしょうか?
model {
#Likelihood
for(i in 1:N1){
y1[i] ~dnorm(mu,tau)
}
tau <- sqrt(1/sigma)
#Priors
mu ~ dnorm(0,0.000001)
sigma ~ dnorm(sigmamu, sigmas)
}
前もって感謝します