# 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 (a[county[i]], tau.y)
  }
  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)
}