model {

	for (j in 1:70) {
		y[j] ~ dbin(theta[j],N[j])
		theta[j] ~ dbeta(alpha,beta)
	}

	alpha ~ dexp(k)
	beta ~ dexp(l)

	k ~ dunif(0.01,1)
	l ~ dunif(0.01,1)

}