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) }