model {
Y ~ dbin(0.5, 8)
P2 <- step(2.5 - Y) # does Y = 2, 1 or 0?
}

   node   mean   sd   MC error   2.5%   median   97.5%   start   sample
   P2   0.1448   0.3519   0.003317   0.0   0.0   1.0   1   10000
   Y   4.004   1.417   0.01291   1.0   4.0   7.0   1   10000

[example-2_1_2-coins0][example-2_1_2-coins1]