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