程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
您现在的位置: 程式師世界 >> 編程語言 >  >> 更多編程語言 >> Python

用Python 計算t分布的置信區間

編輯:Python

設置隨機種子

確保數組的隨機

np.random.seed(1)

生成正態分布

norm_dist = stats.norm(loc, scale)

其中, loc = 數組均值, scale = 數組方差

產生隨機數組

sample = norm_dist.rvs(size=10)
size 為樣本數量

計算t分布的置信區間

當樣本數量比較小<30, 使用t分布

stats.t.interval(0.95,df,mu,se)
df= len(sample)-1 #自由度
mu= np.mean(sample) # 均值
std = np.std(sample,ddof=1) # 樣本標准方差
se = std/ np.sqrt(len(sample) # 標准誤差
0.95 = 置信水平

用到的庫
np = import numpy as np
stats = from scipy import stats
np. 分布名稱.rvs:產生服從指定分布的隨機數

完整代碼

#引入用於計算的庫
import numpy as np
import pandas as pd
import scipy as sp
from scipy import stats
#設置隨機種子
np.random.seed(1)
#設置正態分布
norm_dist=stats.norm(loc=5,scale=0.09)
#采集樣本
sample = norm_dist.rvs(size=20)
#自由度
df=len(sample)-1
#均值
mu=np.mean(sample)
#標准方差, 計算樣本用ddof=1, 總體用ddof=0
std = np.std(sample,ddof=1)
#標准誤差
se = std/np.sqrt(len(sample))
#置信區間
interval = stats.t.interval(0.95,df,mu,se)

參考《用python動手學統計學》


  1. 上一篇文章:
  2. 下一篇文章:
Copyright © 程式師世界 All Rights Reserved