data { int J; // number of trials real y[J]; // estimated log relative risk real sigma[J]; // se of log relative risk } parameters { real mu; real tau; real eta[J]; } transformed parameters { real theta[J]; for (j in 1:J) theta[j] = mu + tau * eta[j]; } model { eta ~ normal(0, 1); y ~ normal(theta, sigma); }