model;
{
 y ~ dbin(pi,n)
 pi ~ dbeta(newalpha,newbeta)I(0,1)
 newalpha <- y + alpha
 newbeta <- n - y + beta
}

data;
list(n=16,alpha=2,beta=4)

inits;
list(pi=0.5,y=8)