bayes function

Bayesian method

Calculate a binomial series lower confidence bound using Bayes' method with a Beta prior distribution.

bayes(s, n, alpha, MonteCarlo, beta.a, beta.b, ...)

Arguments

• s: Vector of successes.
• n: Vector of sample sizes.
• alpha: The significance level; to calculate a 100(1-$\alpha$)% lower confidence bound.
• MonteCarlo: Number of samples to draw from the posterior distribution for the Monte Carlo estimate.
• beta.a: Shape1 parameter for the Beta prior distribution.
• beta.b: Shape2 parameter for the Beta prior distribution.
• ...: Additional arguments to be ignored.

Returns

The 100(1-$\alpha$)% lower confidence bound.

Examples

bayes(s=c(35, 97, 59), n=c(35, 100, 60), alpha=.10, MonteCarlo=1000, beta.a=1, beta.b=1)