TrePreMultiple.txt

model{

for(i in 1:k){
y[i] ~ dbin(ap[i], n[i])
ap[i] <- sub.p[i]*main.Se + (1-sub.p[i])*(1-main.Sp)
#informative hyperpior for prevalence
sub.p[i] ~ dbeta(alpha,beta) T(0.001,0.999)
}

alpha <- main.ap*main.psi
beta <- main.psi*(1-main.ap)
main.ap ~ dbeta(amu, bmu) T(0.001,0.999)
main.psi ~ dgamma(apsi, bpsi) T(0.001,0.999)

#informative prior for Se and Sp
main.Se ~ dbeta(ase, bse) I(0.001,0.999)
main.Sp ~ dbeta(asp, bsp) I(0.001,0.999)

#predictions
y.pre ~ dbin(main.pstar.rep,m)

main.pstar.rep ~ dbeta(alpha,beta)
pre.pequal0 <- equals(main.ap,0)
pre.plessthan0.05 <- step(0.05-main.ap)
plessthanSetvalue <- step(perVal-main.ap)

}