model
{
	Tau.noninformative <- 1.0E-4
	for (i in 1:N.sector) {
		Disturbance[i] ~ dbin(p[i], N.dist)
		logit(p[i]) <- (
			(x[i] - N.size * 0.5)
			+ b[1, Tr[i]] * Distance[i]
		)
		x[i] ~ dbin(q[i], N.size)
		logit(q[i]) <- (
			b[2, Tr[i]]
			+ b[3, Tr[i]] * (x[LeftRight[i, 1]] + x[LeftRight[i, 2]])
		)
	}
	for (j in 1:N.b) {
		for (k in 1:N.transect) {
			b[j, k] ~ dnorm(0, Tau.noninformative)
		}
	}
}