"Imagine, a grown man too ignorant to believe in the Devil!" — R.A. Lafferty, Past Master
R Tutorial: A Bayesian Estimate of Proportion
This is an old chestnut in Bayesian statistics, using the conjugate beta prior to find a beta posterior distribution for a proportion. If you’re unfamiliar with the calculation of the posterior distribution, there’s a link in the tutorial.