model { ############################################################################################### ##-----------------------------------------## ##---------- PARAMETERS & PRIORS ----------## ##-----------------------------------------## #----- SPATIAL PROCESS -----# for(i in 1:n.individuals){ sxy[i,1] ~ dunif(0, x.max) sxy[i,2] ~ dunif(0, y.max) pOK[i] <- habitat.mx[trunc(sxy[i,2])+1, trunc(sxy[i,1])+1] # CHECK IF INSIDE HABITAT MASK. OK[i] ~ dbern(pOK[i]) # ONES TRICK }#i #----- DEMOGRAPHIC PROCESS -----# psi ~ dunif(0,1) #----- DETECTION PROCESS ------# sigma ~ dunif(0,25) alpha <- -1/(2*sigma*sigma) p0 ~ dunif(0,1) for (i in 1:n.individuals){ for(j in 1:n.detectors){ d2[i,j] <- pow(sxy[i,1]-detector.xy[j,1],2) + pow(sxy[i,2]-detector.xy[j,2],2) p[i,j] <- p0*exp(alpha * d2[i,j]) }#j }#i ############################################################################################## ##--------------------------------## ##---------- LIKELIHOOD ----------## ##--------------------------------## for (i in 1:n.individuals){ z[i] ~ dbern(psi) for (j in 1:n.detectors){ y[i,j] ~ dbin(p[i,j] * z[i], n.trials[j]) # REALIZATION OF DETECTION OF INDIVIDUAL i AT DETECTOR j }#j }#i ############################################################################################## ##----------------------------------------## ##---------- DERIVED PARAMETERS ----------## ##----------------------------------------## N <- sum(z) ############################################################################################## }