radon_no_pool_chr.stan

data {
  int<lower=1> N;
  int<lower=1> J; // number of counties
  array[N] int<lower=1, upper=J> county;
  vector[N] x;
  vector[N] y;
}
transformed data {
  matrix[N, 1] cov = [x']';
}
parameters {
  vector[1] beta;
  real mu_a;
  real<lower=0> sigma_a;
  real<lower=0> sigma_y;
  vector<offset=mu_a, multiplier=sigma_a>[J] alpha;
}
model {
  beta ~ std_normal();
  mu_a ~ std_normal();
  alpha ~ normal(mu_a, sigma_a);
  sigma_a ~ cauchy(0, 2.5);
  sigma_y ~ cauchy(0, 2.5);
  
  y ~ normal_id_glm(cov, alpha[county], beta, sigma_y);
}