#  8.6 Maple: Confidence profile examples:  Chapter 27

model;
{
   theta.s <- (theta.sm * r) * 100000
   e.d <- theta.s - theta.n
   r ~ dbeta( 0.5, 0.5)
   phi.s ~ dbeta( 0.5, 0.5)
   phi.n ~ dbeta( 0.5, 0.5)
   theta.em ~ dbeta( 0.5, 0.5)
   theta.lm ~ dbeta( 0.5, 0.5)
   theta.sm <- phi.s * theta.em + 1 - phi.s
   theta.nm <- phi.n * theta.em + (1 - phi.n) * theta.lm
   theta.n <- (theta.nm * r) * 100000
   r.r ~ dbin(r,n.r)
   r.s ~ dbin(phi.s,n.s)
   r.n ~ dbin(phi.n,n.n)
   r.em ~ dbin(theta.em,n.em)
   r.lm ~ dbin(theta.lm,n.lm)
}