statsmodels を使って二項比率の信頼区間を求めるデモプログラム。

binomial_proportion_demo.py

# Requirements:
#   pip install statsmodel pandas numpy scipy

from statsmodels.stats import proportion as prop


def demo(num_clicks, num_impressions):
    alpha = 0.05

    methods = [
        'normal',
        'beta',
        'wilson',
        'agresti_coull'
    ]

    print "---"
    for method in methods:
        print '%s: %s' % (method, prop.proportion_confint(num_clicks, num_impressions, alpha, method))


if __name__ == '__main__':
    demo(5, 1000)
    demo(0, 1000)