stats.proportion.proportion_confint()

statsmodels.stats.proportion.proportion_confint

statsmodels.stats.proportion.proportion_confint(count, nobs, alpha=0.05, method='normal') [source]

confidence interval for a binomial proportion

Parameters:

count : int or array

number of successes

nobs : int

total number of trials

alpha : float in (0, 1)

significance level, default 0.05

method : string in [?normal?]

method to use for confidence interval, currently available methods :

  • normal : asymptotic normal approximation
  • agresti_coull : Agresti-Coull interval
  • beta : Clopper-Pearson interval based on Beta distribution
  • wilson : Wilson Score interval
  • jeffrey : Jeffrey?s Bayesian Interval
  • binom_test : experimental, inversion of binom_test
Returns:

ci_low, ci_upp : float

lower and upper confidence level with coverage (approximately) 1-alpha. Note: Beta has coverage coverage is only 1-alpha on average for some other methods.)

Notes

Beta, the Clopper-Pearson interval has coverage at least 1-alpha, but is in general conservative. Most of the other methods have average coverage equal to 1-alpha, but will have smaller coverage in some cases.

Method ?binom_test? directly inverts the binomial test in scipy.stats. which has discrete steps.

TODO: binom_test intervals raise an exception in small samples if one
interval bound is close to zero or one.

References

http://en.wikipedia.org/wiki/Binomial_proportion_confidence_interval

Brown, Lawrence D.; Cai, T. Tony; DasGupta, Anirban (2001). ?Interval
Estimation for a Binomial Proportion?, Statistical Science 16 (2): 101?133. doi:10.1214/ss/1009213286. TODO: Is this the correct one ?
doc_statsmodels
2017-01-18 16:19:41
Comments
Leave a Comment

Please login to continue.