Rjags ile tahminler nasıl oluşturulur?

JJS dilinde belirtilen bir modelde MCMC çalıştırmak için rjags kullandım. Bu modeli çıkartmanın ve onunla tahminler yapmanın iyi bir yolu var mı (parametrelerimin posterior dağılımlarını kullanarak)? Modeli R'de yeniden belirtebilir ve parametre posterlerimin modlarını takabilirim; Sadece bunu yapmanın daha az gereksiz bir yolu olup olmadığını merak ediyorum.

Http://sourceforge.net/p/mcmc-jags/discussion/610037/thread/0ecab41c aynı soruyu sorduğuna inanıyorum .

Genellikle tahminleri JAGS'de yapabilirsiniz. Aşağıda bağımlı değişken olarak FEV (akciğer kapasitesi ile ilgili bir şey) ve prediktör olarak yaş ve sigara içimi göstergesi olan bir regresyon örneği verilmiştir.

FEV20'ler ve FEV20ns, 20 yaşında sigara içen ve 20 yaşında sigara içmeyen için öngörülen FEV değerleridir.

model
{
for(i in 1:n){
    FEV[i] ~ dnorm(mu[i],tau)
    mu[i] <- beta[1] + beta[2]*Age[i] + beta[3]*Smoke[i]  + beta[4]*Age[i]*Smoke[i]
}

#priors
beta[1] ~ dnorm(0,0.001)
beta[2] ~ dnorm(0,0.001)
beta[3] ~ dnorm(0,0.001)
beta[4] ~ dnorm(0,0.001)
tau ~ dgamma(0.001,0.001)
sigma<-1/sqrt(tau) 

## Predict the FEV for a 20 year old smoker and for a 20 year old nonsmoker
mu20s <-  beta[1] + (beta[2]+beta[4])*20 + beta[3]
mu20ns <-  beta[1] + beta[2]*20 
FEV20s ~ dnorm(mu20s,tau)
FEV20ns ~ dnorm(mu20ns,tau)
}

Örnek: Bayesci Fikirler ve Veri Analizi