model
{
	Tau.noninformative <- 1.0E-4
	P.gamma <- 1.0E-2
	for (i in 1:N.sector) {
		# observation
		Disturbance[i] ~ dbin(p[i], N.dist)
		logit(p[i]) <- (
			a[Tr[i]]
			+ x[i]
			+ b * Distance[i]
		)
		# latent variable
		x[i] ~ dnorm(q[i], tau)
		q[i] <- (x[LeftRight[i, 1]] + x[LeftRight[i, 2]]) * 0.45
	}
	# parameters
	for (k in 1:N.transect) {
		a[k] ~ dnorm(0, Tau.noninformative)
	}
	b ~ dnorm(0, Tau.noninformative)
	tau ~ dgamma(P.gamma, P.gamma)
}