data {
  int<lower=0> J;
  array[J] real y;
  array[J] real<lower=0> sigma;
}
parameters {
  vector[J] eta;
  real mu;
  real<lower=0> tau;
}
transformed parameters {
  vector[J] theta;
  theta = mu + tau*eta;
}
model {
  y ~ normal(theta,sigma);
  eta ~ normal(0,1);
}