B.9 - JAGS model to perform the same Monte Carlo evaluation done in Appendix B.8

(Only the model is provided - steering R commands are left as exercise.)

model {
  n.I.1 ~ dbin(p1, ns1)
  n.NI.1 <- ns1 - n.I.1 
  nP.I.1 ~ dbin(pi1, n.I.1)
  nP.NI.1 ~ dbin(pi2, n.NI.1)
  nP.1 ~ sum(nP.I.1, nP.NI.1)
  fP.1 <- nP.1 / ns1

  n.I.2 ~ dbin(p2, ns2)
  n.NI.2 <- ns2 - n.I.2 
  nP.I.2 ~ dbin(pi1, n.I.2)
  nP.NI.2 ~ dbin(pi2, n.NI.2)
  nP.2 ~ sum(nP.I.2, nP.NI.2)
  fP.2 <- nP.2 / ns2

  D.fP <- fP.2 - fP.1

  pi1 ~ dbeta(r1, s1)
  pi2 ~ dbeta(r2, s2) 
}