data = c(1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
posterior <- prop_model(data)
median(posterior)

0.19

quantile(posterior, c(0.05, 0.95))

5%    95% 
0.06  0.39

sum(posterior > 0.07) / length(posterior)

0.93