model gaussian;
# inferenza di mu e sigma, piu' distr. predittiva, 
# assumendo verosomiglianza gaussiana

const
   N = NNNN;     # number of observations
                 # NNNN viene sostituito da script run_bugs.sh
var
   X[N], mu, tau, y, sigma;

data X in "sample.dat";
inits in "mu0_tau0.in"; # questo file viene generato da set_init.R

{
   for (i in 1:N) {
      X[i] ~ dnorm(mu,tau);
   }
   mu ~ dnorm(0.0, 1.0E-4);   
   tau   ~ dgamma(1.0E-3, 1.0E-3);
   sigma <- 1.0/sqrt(tau);
   y ~ dnorm(mu,tau)
}