# File name "burn2.jags"
model{
  # likelihood
  for(i in 1:nSite) {
    logit(psi[i]) <- b0[fire[i]] + bElev * elev[i]
    z[i] ~ dbern(psi[i])
    y[i] ~ dbin(p * z[i], nOcc[i])
  }
  # priors
  for(i in 1:nCat) {
    psi0[i] ~ dbeta(1,1)
    b0[i] <- logit(psi0[i])
  }
  bElev ~ dunif(-5, 5)
  p ~ dbeta(1,1)
}