model;
{
   theta ~ dnorm(mu,prec)I(0,1)
   y ~ dbin(theta,n)
   ypred ~ dbin(theta,npred)
}