model gaussian; # Semplicemente una singola misura, con sigma nota const N = 1, # number of observations sigma = 1.5; # sigma fissa var X[N], mu, tau; data X in "gauss_1m.dat"; inits in "gauss_1m.in"; # l'inizializzazione non e' fondamentale # in questo semplice caso (provare a commentare) { tau <- 1/(sigma*sigma); for (i in 1:N) { X[i] ~ dnorm(mu,tau); } mu ~ dnorm(0.0, 1.0E-4); }