4. 區(qū)間估計

還以為你被上節(jié)課的內(nèi)容唬住了~終于等到你，還好沒放棄！

本節(jié)我們將說明兩個問題：總體均值 的區(qū)間估計和總體比例 的區(qū)間估計。

區(qū)間估計經(jīng)常用于質(zhì)量控制領(lǐng)域來檢測生產(chǎn)過程是否正常運(yùn)行或者在“控制之中” ，也可以用來監(jiān)控互聯(lián)網(wǎng)領(lǐng)域各類數(shù)據(jù)指標(biāo)是否在正常區(qū)間。

一個總體均值的區(qū)間估計

• 大樣本的情況下

  • 已知，

  • 未知，

• 小樣本的情況下

  • 已知，
  • 未知，

另外補(bǔ)充一個公式，樣本量 這個了解就好，大部分情況下是不缺數(shù)據(jù)的，盡可能選數(shù)據(jù)量稍大些的數(shù)據(jù)。

把以上過程編寫成Python的自定義函數(shù)：

import numpy as np
import scipy.stats
from scipy import stats as sts


def mean_interval(mean=None, sigma=None,std=None,n=None,confidence_coef=0.95):
    """
    mean:樣本均值
    sigma: 總體標(biāo)準(zhǔn)差
    std: 樣本標(biāo)準(zhǔn)差
    n:   樣本量
    confidence_coefficient：置信系數(shù)
    confidence_level:置信水平 置信度
    alpha:顯著性水平
    功能：構(gòu)建總體均值的置信區(qū)間
    """
    alpha = 1 - confidence_coef
    z_score = scipy.stats.norm.isf(alpha / 2)            # z分布臨界值
    t_score = scipy.stats.t.isf(alpha / 2, df = (n-1) )  # t分布臨界值
   
    if n >= 30: 
        if sigma != None:
            me = z_score * sigma / np.sqrt(n)
            print("大樣本，總體 sigma 已知：z_score:",z_score)
        elif sigma == None:
            me = z_score * std / np.sqrt(n)
            print("大樣本，總體 sigma 未知 z_score",z_score)
        lower_limit = mean - me
        upper_limit = mean + me
    if n < 30 :
        if sigma != None:
            me = z_score * sigma / np.sqrt(n)
            print("小樣本，總體 sigma 已知 z_score * sigma / np.sqrt(n) n z_score = ",z_score)
        elif sigma == None:
            me = t_score * std / np.sqrt(n)
            print("小樣本，總體 sigma 未知 t_score * std / np.sqrt(n) n t_score = ",t_score)
            
        print("t_score:",t_score)
        lower_limit = mean - me
        upper_limit = mean + me
    
    return (round(lower_limit, 1), round(upper_limit, 1))

應(yīng)用：網(wǎng)站流量UV區(qū)間估計:

某網(wǎng)站流量UV數(shù)據(jù)如下[52,44,55,44,45,59,50,54,62,46,54,42,60,62,43,42,48,55,57,56]，我們研究一下該網(wǎng)站的總體流量uv均值，我們先把數(shù)據(jù)放進(jìn)來

import numpy as np
data = np.array([52,44,55,44,45,59,50,54,62,46,54,42,60,62,43,42,48,55,57,56])

計算一下均值為：

x_bar = data.mean()
x_bar
# 51.5

樣本標(biāo)準(zhǔn)差為：

x_std = sts.tstd(data,ddof = 1) #  ddof=1時,分母為n-1;ddof=0時,分母為n
x_std
# 6.840283158189472

進(jìn)行區(qū)間估計：

mean_interval(mean=x_bar, sigma=None,std= x_std,  n=n, confidence_coef=0.95)

輸出結(jié)果：

小樣本，總體 sigma 未知 t_score * std / np.sqrt(n) 
t_score =  2.093024054408263
(48.3, 54.7)

于是我們有95%的把握，該網(wǎng)站的流量uv介于 [48, 55]之間。

值得一提的是，上面這個案例的數(shù)據(jù)是實際上是公眾號山有木兮水有魚 的按天統(tǒng)計閱讀量……有人可能要說了，你這數(shù)據(jù)也太慘了，而且舉個案例都是小樣本。我想說，小樣本的原因是這新號一共發(fā)了也沒幾天，至于數(shù)量低，你幫忙動動小手轉(zhuǎn)發(fā)轉(zhuǎn)發(fā)，這數(shù)據(jù)也就高了~希望下次舉例的時候這個能變成大樣本，均值怎么著也得個千兒八百的，感謝感謝！

一個總體比例的區(qū)間估計

其中樣本量

def proportion_interval(p=None, n=None, confidence_coef =0.95):
    """
    p: 樣本比例
    n: 樣本量
    confidence_coef: 置信系數(shù)
    功能：構(gòu)建總體比例的置信區(qū)間
    """
    alpha = 1 - confidence_coef
    z_score = scipy.stats.norm.isf(alpha / 2)  # z分布臨界值
    
    me = z_score * np.sqrt((p * (1 - p)) / n) 
    lower_limit = p - me
    upper_limit = p + me
    
    return (round(lower_limit, 3), round(upper_limit, 3))

下期將為大家?guī)?a href="http://www.3lll3.cn/view/205064.html" style="text-decoration: none; color: #1e6bb8; word-wrap: break-word; font-weight: bold; border-bottom: 1px solid #1e6bb8;">《Python統(tǒng)計學(xué)極簡入門》之假設(shè)檢驗

這里分享一個你一定用得到的小程序——CDA數(shù)據(jù)分析師考試小程序。 它是專為CDA數(shù)據(jù)分析認(rèn)證考試報考打造的一款小程序?？梢詭湍憧焖賵竺荚?、查成績、查證書、查積分，通過該小程序，考生可以享受更便捷的服務(wù)。 掃碼加入CDA小程序，與圈內(nèi)考生一同學(xué)習(xí)、交流、進(jìn)步！