	model
	{
		for (i in 1 : K) {
			r0[i]  ~ dbin(p0[i], n0[i])
			r1[i] ~ dbin(p1[i], n1[i])
			logit(p0[i]) <- mu[i]
			logit(p1[i]) <- mu[i] + logPsi[i]
			logPsi[i]    <- alpha + beta1 * year[i] + beta2 * (year[i] * year[i] - 22) + b[i]
			b[i] ~ dnorm(0, tau)
			mu[i]  ~ dnorm(0.0, 1.0E-6)
		}
		alpha  ~ dnorm(0.0, 1.0E-6)
		beta1  ~ dnorm(0.0, 1.0E-6)
		beta2  ~ dnorm(0.0, 1.0E-6)
		tau    ~ dgamma(1.0E-3, 1.0E-3)
		sigma <- 1 / sqrt(tau)
	}
