BUGS example 3
Gillian Raab , Napier University
June 1999



model;
{
   mean ~ dnorm( 0.0,1.0E-6)
   tau ~ dgamma(0.001,0.001)
   sd <- sqrt(1 / tau)
var<-sd*sd
   for( i in 1 : N ) {
      x[i] ~ dnorm(mean,tau)
   }
}

data;
list(N=5,x=c(2,3,3,1,0))

inits;
list(tau=1)