package R2WinBUGS question; was: Multilevel
Can you please 1. reply to me personally, not only to the list (which I won't read all the time) and 2. cite what you wrote before. I tend to delete messages that are not of importance for me and rely on the fact that old information is cited. Now I do not have the data nor the code anymore and I really do not start to search in the web or in the list archives. Best, Uwe Ligges
On 29.01.2011 18:20, martanair wrote:
I have R2WinBUGS
and this is my radon1.bug
# Bugs code for multilevel model for radon
# with bsmt as an individual predictor
# varying-intercept model
model {
for (i in 1:n){
y[i] ~ dnorm (y.hat[i], tau.y)
y.hat[i]<- a[county[i]] + b*x[i]
}
b ~ dnorm (0, .0001)
tau.y<- pow(sigma.y, -2)
sigma.y ~ dunif (0, 100)
for (j in 1:J){
a[j] ~ dnorm (mu.a, tau.a)
}
mu.a ~ dnorm (0, .0001)
tau.a<- pow(sigma.a, -2)
sigma.a ~ dunif (0, 100)
}