數據選自Journal of the American Medical Association(http://jse.amstat.org/v4n2/datasets.shoemaker.html

)關于體溫、性別、心率的臨床數據

現對男性體溫抽樣計算下95%置信區間總體均值范圍。

1、讀取數據

import pandas as pd

import numpy as np

import seaborn as sns

import matplotlib.pyplot as plt

#讀取數據

df = pd.read_csv('http://jse.amstat.org/datasets/normtemp.dat.txt', header = None,sep = '\s+' ,names=['體溫','性別','心率'])

體溫數據.png

2、選取樣本大小，查看數據

np.random.seed(42)

#df.describe()

#樣本量為90，查看樣本數據

df_sam = df.sample(90)

df_sam.head()

樣本數據.jpeg

3、計算抽取樣本中男士體溫的均值

df3 = df_sam.loc[df_sam['性別']==1]

df3['體溫'].mean()

4、重復抽取樣本,計算其他樣本中男士體溫的均值,得到抽樣分布

boot_means = []

for _ in range(10000):

bootsample = df.sample(90, replace=True)

mean = bootsample[bootsample['性別'] == 1]['體溫'].mean()

boot_means.append(mean)

5、繪制男士體溫抽樣分布均值

男士體溫抽樣分布均值.png

6、計算抽樣分布的置信區間以估計總體均值, 置信度95%

np.percentile(boot_means, 2.5), np.percentile(boot_means, 97.5)

(97.89249519230768, 98.30741452991455)

總結

以上是生活随笔為你收集整理的python参数估计_用python求参数估计的置信区间的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。