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 mu;            // koulukohtaisten harjoitusvaikutusten keskiarvo
  real<lower=0> tau;  // koulukohtaisten harjoitusvaikutusten keskihajonta, huomaa maarittelyvali (0,\infty),
                      // tama PITAA antaa parametreille
  real theta[J];      // koulukohtaiset harjoitusvaikutkset
}

model {
  // kaytetaan priorina tasajakaumaa p(mu,tau) \propto 1, joten prioria ei tarvitse maarittaa erikseen
  theta ~ normal(mu, tau);  // harjoitusvaikutuksen theta jakauma
  y ~ normal(theta, sigma); // keskimaaraisen hyodyn jakaumat kouluittain
}