Learn R Programming

statpsych (version 1.7.0)

ci.bayes.prop: Bayesian credible interval for a proportion

Description

Computes a Bayesian credible interval for a population proportion using the mean and standard deviation of a prior Beta distribution along with sample information. The mean and standard deviation of the posterior Beta distribution are also reported. For a noninformative prior, set the prior mean to .5 and the prior standard deviation to 1/sqrt(12) (which corresponds to a Beta(1,1) distribution). The prior variance must be less than m(1 - m) where m is the prior mean.

Usage

ci.bayes.prop(alpha, prior_mean, prior_sd, f, n)

Value

Returns a 1-row matrix. The columns are:

  • Posterior mean - posterior mean of Beta distribution

  • Posterior SD - posterior standard deviation of Beta distribution

  • LL - lower limit of the credible interval

  • UL - upper limit of the credible interval

Arguments

alpha

alpha level for 1-alpha credibility interval

prior_mean

mean of prior Beta distribution

prior_sd

standard deviation of prior Beta distribution

f

number of participants who have the attribute

n

sample size

References

Gelman2004statpsych

Examples

Run this code
ci.bayes.prop(.05, .4, .1, 12, 100)

# Should return:
# Posterior mean Posterior SD       LL        UL
#      0.1723577   0.03419454 0.1111747 0.2436185


Run the code above in your browser using DataLab