Dear the developer of JAGS, I have a problem when I am using dbern and pnorm. My problem is that a JAGS code runs successfully when For (I in N) { d[i] <- [some code] p[i] <- pnorm(d[i],0,1) ip[i] <- ifelse(p[i]>0.5, 1-(1e-10), 1e-10) y[i] ~ dbern(ip[i]) } But, the JAGS code did not run when [1] For (I in N) { d[i] <- [some code] p[i] <- pnorm(d[i],0,1) ip[i] <- ifelse(p[i]>0.5, 1, 0) y[i] ~ dbern(ip[i]) } Or, when [2] For (I in N) { d[i] <- [some code] p[i] <- pnorm(d[i],0,1) y[i] ~ dbern(p[i])...