data { int J; real y[J]; real sigma[J]; } parameters { real mu; real tau; vector[J] eta; } transformed parameters { vector[J] theta; theta = mu + tau*eta; } model { eta ~ normal(0.0, 1.0); y ~ normal(theta, sigma); }