模擬結(jié)果圖

按照兩個(gè)樣本方差比分布的理論推導(dǎo)可以得到

這里

為樣本方差，

為總體方差。

這里將對(duì)這一結(jié)論進(jìn)行基于Python的計(jì)算機(jī)模擬。代碼如下。

1 生成兩個(gè)總體，總體1服從正態(tài)分布

，總體2服從正態(tài)分布

。

import matplotlib.pyplot as plt #加載可視化工具

import seaborn as sns #加載可視化工具

import numpy as np #加載numpy工具

import pandas as pd

np.random.seed(1234)# 設(shè)置隨機(jī)數(shù)種子，保證每次的隨機(jī)數(shù)相同

pdata1=4+3*np.random.randn(50000) #生成服從標(biāo)準(zhǔn)正態(tài)分布的50000個(gè)數(shù)據(jù)

np.random.seed(2345)

pdata2=7+2*np.random.randn(50000) #生成服從標(biāo)準(zhǔn)正態(tài)分布的50000個(gè)數(shù)據(jù)

sns.distplot(pdata1,bins=50,kde=True,label="$N(4,3^2)$") #繪制直方圖，50個(gè)分箱(直方)，kde表示需要密度曲線

sns.distplot(pdata2,bins=50,kde=True,label="$N(7,2^2)$") #繪制直方圖，50個(gè)分箱(直方)，kde表示需要密度曲線

plt.title("Population")

plt.legend()

plt.show()

總體分布

2 隨機(jī)抽樣5000次，樣本量分別為50和40，得到兩個(gè)樣本方差

的分布，以及

的分布

from tqdm import tqdm

import statistics

loop_var=pd.DataFrame(columns=('vars1', 'vars2','vars_ratio'))

loop=5000

n_size=[50,40]

for loop in tqdm(range(loop)):

np.random.seed(loop+1)

sample1 = np.random.choice(pdata1, size=n_size[0])

vars1 = statistics.variance(sample1)

np.random.seed(loop+1+5000)

sample2 = np.random.choice(pdata2, size=n_size[1])

vars2 = statistics.variance(sample2)

vars_ratio=(vars1/9)/(vars2/4)

loop_var=loop_var.append(pd.DataFrame({'vars1':[vars1],'vars2':[vars2],'vars_ratio':[vars_ratio]}),ignore_index=True)

loop_var.head()

3 畫(huà)出兩個(gè)樣本方差

的分布圖，以及

的分布圖

import matplotlib.pyplot as plt #加載可視化工具

import seaborn as sns #加載可視化工具

plt.figure()

fig = plt.figure(figsize=(8,5))

sns1=sns.distplot(loop_var["vars1"],bins=30,kde=True,label="Distribution of sample variance 1") #繪制直方圖，30個(gè)分箱(直方)，kde表示需要密度曲線

sns2=sns.distplot(loop_var["vars2"],bins=30,kde=True,label="Distribution of sample variance 2")

sns3=sns.distplot(loop_var["vars_ratio"],bins=30,kde=True,label="Distribution of sample variance ratio")

#plt.xlabel("X")

plt.legend()

plt.show()

抽樣分布

動(dòng)動(dòng)手，讓統(tǒng)計(jì)學(xué)變得更有趣！！

與50位技術(shù)專(zhuān)家面對(duì)面20年技術(shù)見(jiàn)證，附贈(zèng)技術(shù)全景圖

總結(jié)

以上是生活随笔為你收集整理的python方差选择样本_两个样本方差比的分布（来自Python的实验）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。