model   {
for(i in 1:N){
PregStatus[i] ~ dbin(regression_prob[i], 1)
logit(regression_prob[i]) <- intercept[Site[i]] + b1 * x1[i] + GROUP_effect[GROUP[i]]


for(Site_iterator in 1:25){
    intercept[Site_iterator] ~ dnorm(mu_int, tau_int)
}

#imperfect detection of x1
tmp[i] <- dnorm(N,1,1)
p1[i]<-tmp[i]/(1+tmp[i])

#modeling disease status d1
d1[i]~dbern(p1[i])

#probability of getting a positive diagnostic test result p2
p2[i]<-ifelse(d1[i]==1,SN,1-SP)

#modeling of diagnostic test result x1
x1[i]~dbern(p2[i])