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)
   for( i in 1 : N ) {
llim[i]<-x[i]-0.5
ulim[i]<-x[i]+0.5
      truex[i] ~ dnorm(mean,tau) I(llim[i],ulim[i])
   }
}

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

inits;
list(truex=c(2,3,3,1,0),tau=1)