model
		{
		   for( i in 1 : N ) {
		      p[i] ~ dbeta(alpha, beta)
			r[i] ~ dbin(p[i], n[i])
		   }
		alpha ~ dexp(1.0)
		beta ~ dexp(1.0)
		}
