	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
	}
