python 信賴區間

import numpy as np

import scipy.stats as stats

# 估計 95% 信心水準下，老年人每星期看電視平均時間的信賴區間？

# alpha = 0.05

# 假設樣本平均時間的抽樣分配服從常態分佈: xbar ~ N(mu_xbar, sigma_xbar)

# 抽樣人數 n = 100 ; 樣本平均時間 xbar = 21.2 小時

# 母體看電視時間的標準差 sigma = 8 小時

def Interval_estimation(n, xbar, sigma, alpha):

    # standard error of the mean:

    se = sigma / np.sqrt(n)

    # margin of error:

    # Hint: Z_(0.05/2) for Two-tails

    E = stats.norm.ppf(1 - alpha/2, loc=0, scale=1) * se

    # (1 - alpha)% Confidence Interval of Population mean with Known Variance:

    Interval = [xbar - E, xbar + E]

    return Interval

Interval_estimation(n=100, xbar=21.2, sigma=8, alpha=0.05)