model
	{
	for( i in 1 : N ) {
		y[i] ~ dnegbin(theta[i],2)
		theta[i] ~ dbeta(alpha,beta)
		}
	alpha ~ dexp(1)
	beta ~ dexp(1)
	}
