model
{
# Priors
alpha.occ ~ dunif(-20, 20)
beta.occ ~ dunif(-20, 20)
alpha.p ~ dunif(-20, 20)
beta1.p ~ dunif(-20, 20)
beta2.p ~ dunif(-20, 20)
# Likelihood
for (i in 1:R) {
# Model for partially latent state
z[i] ~ dbern(psi[i]) # True occupancy z at site i
logit(psi[i]) <- alpha.occ + beta.occ * wetness[i]
for (j in 1:T) {
# Observation model for actual observations
y[i,j] ~ dbern(eff.p[i,j]) # Det.-nondet. at i and j
eff.p[i,j] <- z[i] * p[i,j]
logit(p[i,j]) <- alpha.p + beta1.p * wetness [i] + beta2.p * experience[i,j]
}
}
# Derived quantity
occ.fs <- sum(z[]) # Finite sample number of occupied sites
}