model {
r.obs <- 2
theta ~ dbeta(3, 27)
r ~ dbin(theta, 5) # sampling distribution
P <- step(r-r.obs-0.5) + 0.5*equals(r, r.obs) # mid-p-value
}
node mean sd MC error 2.5% median 97.5% start sample
P 0.0538 0.1759 0.001798 0.0 0.0 0.5 1001 10000
r 0.5021 0.7082 0.00735 0.0 0.0 2.0 1001 10000
theta 0.1 0.05374 5.434E-4 0.02263 0.09124 0.2269 1001 10000