model{
  theta ~ dbeta(1,1)
  y ~ dbin(theta,n)
}