python两个总体参数的区间估计（方差比）

方差比

$\frac{s_1^2/s_2^2}{F_{a/2}}\le \frac{{\sigma }_1^2}{{\sigma }_2^2}\le \frac{s_1^2/s_2^2}{F_{1-a/2}}$

例：

为研究男女学生在生活费支出上的差异，在某大学随机抽取25名男学生和25名女学生，得到下面的结果：

男学生：${\bar{x}}_1=520,s_1^2=260$

女学生：${\bar{x}}_2=480,s_2^2=280$

试以90%的置信水平估计男女学生生活费支出方差比的置信区间

# 男学生
n1 = 25
x1_bar = 520
sigma1_2 = 260

#女学生
n2 = 25
x2_bar = 480
sigma2_2 = 280

a = 0.1

计算F分布

# 计算 自由度 

df1 = n1-1
df2 = n2-1


f_a2 = stats.f.isf(a/2,df1,df2)
f_1_a2 = stats.f.ppf(a/2,df1,df2)

计算公式

left = (sigma1_2/sigma2_2)/f_a2
right = (sigma1_2/sigma2_2)/f_1_a2

print('方差比为：{:.3f}'.format(sigma1_2/sigma2_2))
print('男女学生生活费支出方差比的90%的置信区间为（{:.3f}，{:.3f}）'.format(left,right))

方差比为：0.929
男女学生生活费支出方差比的90%的置信区间为（0.468，1.842）