[R] package R2WinBUGS question; was: Multilevel
Uwe Ligges
ligges at statistik.tu-dortmund.de
Sat Jan 29 18:35:47 CET 2011
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)
> }
>
>
>
More information about the R-help
mailing list