model{ b.eff~dunif(-100, 100) b.spec~dunif(-100, 100) psi~dunif(0,1) for (i in 1:2){ alpha[i]~dunif(-100, 100) ##monitor p.out p.out.l[i]<-alpha[leech_type[i]]+ b.eff*mean_eff + b.spec*mean_spec p.out[i]<-exp(p.out.l[i])/(exp(p.out.l[i])+1) } for (j in 1:J){ z[j]~dbern(psi) for (k in 1:maxocc){ l.psp[j,k]<-alpha[leech_type[k]] + b.eff*efftot[j,k] + b.spec*spec[j,k] p[j,k]<-exp(l.psp[j,k])/(1+exp(l.psp[j,k])) p.eff[j,k]<-z[j]*p[j,k]*effabs[j,k] y[j,k]~dbern(p.eff[j,k]) } } }