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
   
[example-8_10_1-prior-data-conflict0][example-8_10_1-prior-data-conflict1][example-8_10_1-prior-data-conflict2]