概率统计Python计算:双正态总体已知总体方差总体均值差双侧假设的Z检验

在这里插入图片描述
X X X Y Y Y相互独立且 X X X~ N ( μ 1 , σ 1 2 ) N(\mu_1,\sigma_1^2) N(μ1,σ12) Y Y Y~ N ( μ 2 , σ 2 2 ) N(\mu_2,\sigma_2^2) N(μ2,σ22),其中 σ 1 2 \sigma_1^2 σ12 σ 2 2 \sigma_2^2 σ22是已知的。来自 X X X Y Y Y的容量分别为 n 1 n_1 n1 n 2 n_2 n2的样本均值为 X ‾ \overline{X} X Y ‾ \overline{Y} Y。对显著水平 α \alpha α,检验双侧假设 H 0 : μ 1 − μ 2 = δ , H 1 : μ 1 − μ 2 ≠ δ H_0:\mu_1-\mu_2=\delta,H_1:\mu_1-\mu_2\not=\delta H0:μ1μ2=δ,H1:μ1μ2=δ。由于检验统计量 X ‾ − Y ‾ − δ σ 1 2 / n 1 + σ 2 2 / n 2 \frac{\overline{X}-\overline{Y}-\delta}{\sqrt{\sigma_1^2/n_1+\sigma_2^2/n_2}} σ12/n1+σ22/n2 XYδ~ N ( 0 , 1 ) N(0,1) N(0,1),与利用ztestL或ztestR计算已知总体方差计算两总体均值差单侧假设的Z检验相仿,我们可以用函数ztest2计算双侧假设 H 0 : μ 1 − μ 2 = δ H_0:\mu_1-\mu_2=\delta H0:μ1μ2=δ的检验。
例1设甲、乙两厂生产同型号的灯泡,其寿命 X X X Y Y Y分别服从正态分布 N ( μ 1 , σ 1 2 ) , N ( μ 2 , σ 2 2 ) N(\mu_1,\sigma_1^2),\\N(\mu_2, \sigma_2^2) N(μ1,σ12),N(μ2,σ22),已知它们寿命的标准差分别为84h和96h,现从两厂生产的灯泡中各取60只,测得灯泡的平均寿命甲厂为1295h,乙厂为1230h。在显著水平 α = 0.05 \alpha=0.05 α=0.05下,问能否认为两厂生产的灯泡寿命有无显著差别?
解: 按题意,需对假设
H 0 : μ 1 − μ 2 = 0 , H 1 : μ 1 − μ 2 ≠ 0 H_0:\mu_1-\mu_2=0, H_1:\mu_1-\mu_2\not=0 H0:μ1μ2=0,H1:μ1μ2=0
作显著水平为 α = 0.05 \alpha=0.05 α=0.05的双侧检验。下列代码完成本例计算。

import numpy as np								#导入numpy
xmean=1295										#样本均值
ymean=1230										#样本均值
xsigma2=84**2									#总体方差
ysigma2=96**2									#总体方差
n1=60											#样本容量
n2=60											#样本容量
alpha=0.05										#显著水平
Z=(xmean-ymean)/np.sqrt(xsigma2/n1+ysigma2/n2)	#检验统计量值
accept=ztest2(Z, alpha)							#双侧假设Z检验
print('mu1-mu2=0 is %s.'%accept)

第2~8行按题面设置各项数据,第9行计算检验统计量值 X ‾ − Y ‾ − δ σ 1 2 / n 1 + σ 2 2 / n 2 \frac{\overline{X}-\overline{Y}-\delta}{\sqrt{\sigma_1^2/n_1+\sigma_2^2/n_2}} σ12/n1+σ22/n2 XYδ为Z(注意,此时 δ = 0 \delta=0 δ=0),第10行调用函数ztest2检验双侧假设 H 0 : μ 1 − μ 2 = 0 H_0:\mu_1-\mu_2=0 H0:μ1μ2=0。运行程序,输出

mu1-mu2=0 is False.

表示拒绝假设 H 0 : μ 1 − μ 2 = 0 H_0:\mu_1-\mu_2=0 H0:μ1μ2=0,即两厂生产的灯泡寿命有显著差别。
例2 全市高三学生进行数学毕业会考,随机抽取10名男生和8名女生的会考成绩如下:
男生: 65 , 72 , 89 , 56 , 79 , 63 , 92 , 48 , 75 , 81 女生: 78 , 69 , 65 , 61 , 54 , 87 , 51 , 67 \text{男生:}65,72,89,56,79,63,92,48,75,81\\ \text{女生:}78,69,65,61,54,87,51,67 男生:65,72,89,56,79,63,92,48,75,81女生:78,69,65,61,54,87,51,67
若男生的会考成绩 X X X~ N ( μ 1 , 1 0 2 ) N(\mu_1, 10^2) N(μ1,102),女生的会考成绩 Y Y Y~ N ( μ 2 , 9. 5 2 ) N(\mu_2,9.5^2) N(μ2,9.52),试问男生和女生的平均成绩是否相同( α = 0.05 \alpha=0.05 α=0.05)?
解: 按题意,需对双侧假设
H 0 : μ 1 − μ 2 = 0 , H 1 : μ 1 − μ 2 ≠ 0 H_0:\mu_1-\mu_2=0, H_1:\mu_1-\mu_2\not=0 H0:μ1μ2=0,H1:μ1μ2=0
作Z检验。下列代码完成本例计算。

import numpy as np#导入numpy
x=np.array([65, 72, 89, 56, 79, 63, 92, 48, 75, 81])#样本数据
y=np.array([78, 69, 65, 61, 54, 87, 51, 67])		#样本数据
xmean=x.mean()										#样本均值
ymean=y.mean()										#样本均值
xsigma2=10**2										#总体方差
ysigma2=9.5**2										#总体方差
n1=x.size											#样本容量
n2=y.size											#样本容量
alpha=0.05											#显著水平
Z=(xmean-ymean)/np.sqrt(xsigma2/n1+ysigma2/n2)		#检验统计量值
accept=ztest2(Z, alpha)								#双侧假设Z检验
print('mu1-mu2=0 is %s.'%accept)

第11行计算检验统计量值 X ‾ − Y ‾ − δ σ 1 2 / n 1 + σ 2 2 / n 2 \frac{\overline{X}-\overline{Y}-\delta}{\sqrt{\sigma_1^2/n_1+\sigma_2^2/n_2}} σ12/n1+σ22/n2 XYδ为Z(注意,此时 δ = 0 \delta=0 δ=0),第12行调用函数ztest2检验双侧假设 H 0 : μ 1 − μ 2 = 0 H_0:\mu_1-\mu_2=0 H0:μ1μ2=0。运行程序,输出

mu1-mu2=0 is True.

表示接受假设 H 0 : μ 1 − μ 2 = 0 H_0:\mu_1-\mu_2=0 H0:μ1μ2=0,即男女生的平均成相同。
写博不易,敬请支持:
如果阅读本文于您有所获,敬请点赞、评论、收藏,谢谢大家的支持!
代码诚可贵,原理价更高。若为AI学,读正版书好
返回《导引》

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值