data {
  int<lower=0> J; // koulujen maara
  real y[J];      // keskimaaraiset hyodyt harjoittelusta jokaiselle koululle
  real<lower=0> sigma[J]; // keskimaaraisten hyotyjen keskivirheet kouluille
}

parameters {
  real theta; // theta on skalaari, eli sama harjoitusvaikutus jokaiselle koululle
}

model {
  y ~ normal(theta, sigma);
}