data {
  int<lower=0> J;
  real y[J];
  real<lower=0> sigma[J];
}

parameters {
  real mu;
  real<lower=0> tau;
  real eta[J];    // korvataan theta parametrilla eta
}

// lisataan 'transformed parameters'-blokki, jonne maaritellaan theta etan funktiona
transformed parameters {
  real theta[J];
  for(j in 1:J) 
    theta[j] = mu + eta[j];
}

model {
  eta ~ normal(0, tau);  // eta on koulun harjoitusvaikutuksen poikkeama populaatiojakauman keskiarvosta mu
  y ~ normal(theta, sigma);
}