§7.3正态总体均值与方差的置信区间
一、单个正态总体均值的置信区间
设(X 1 ,X 2 ,⋯,X n )是来自正态总体N(μ,σ 2 )的样本.
1.设σ 2 已知,求μ的置信度为1−α的置信度区间.
(X ¯ ¯ ¯ −σn √ μ α2 ,X ¯ ¯ ¯ +σn √ μ α2 )
2.设σ 2 未知,求μ的置信度为1−α的置信区间.
由于S 2 是σ 2 的无偏估计,因此用S 2 代替σ 2 ,有t=X ¯ ¯ ¯ −μSn √ ∼t(n−1)
由附表3查得t α2 (n−1),有
P⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ∣ ∣ ∣ ∣ ∣ ∣ X ¯ ¯ ¯ −μSn √ ∣ ∣ ∣ ∣ ∣ ∣ <t α2 (n−1)⎫ ⎭ ⎬ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ =1−α,
即P{X ¯ ¯ ¯ −Sn √ t α2 (n−1)<μ<X ¯ ¯ ¯ +Sn √ t α2 (n−1)}=1−α,
于是得到μ的1−α置信区间为
(X ¯ ¯ ¯ −Sn √ t α2 (n−1),X ¯ ¯ ¯ +Sn √ t α2 (n−1))
例1.某车间生产的螺杆直径服从正态分布N(μ,σ 2 ),今随机地从中抽取5只测得直径值为22.3,21.5,22.0,21.8,21.4.
(1)已知σ=0.3,求μ的0.95置信区间;
(2)如果σ未知,求μ的0.95置信区间.
解:n=5,x ¯ =15 ∑ i=1 5 x i =21.8.
(1)已知σ=0.3,α=0.05,查表得μ α2 =1.96,因此μ的0.95置信区间为
(x ¯ −σn √ μ α2 ,x ¯ +σn √ μ α2 )=(21.8−0.35 √ ×1.96,21.8+0.35 √ ×1.96)=(21.537,22.063)
(2)σ未知,s 2 =15−1 ∑ i=1 5 (x i −x ¯ ) 2 =0.135,s=0.367.查表得t α2 (5−1)=t 0.025 (4)=2.7764,因此μ的0.95置信区间为(x ¯ −sn √ t α2 (n−1),x ¯ +sn √ t α2 (n−1))=(21.8−0.3675 √ ×2.7764,21.8+0.3675 √ ×2.7764)=(21.345,22.255)
二、单个正态总体方差的置信区间
1.设μ已知,求σ 2 的置信度为1−α的置信区间
由于
χ 2 =1σ 2 ∑ i=1 n (X i −μ) 2 ∼χ 2 (n)
对于给定置信度1−α,查表可得χ 2 1−α2 (n)及χ 2 α2 (n),使
P{χ 2 1−α2 (n)<1σ 2 ∑ i=1 n (X i −μ) 2 <χ 2 α2 (n)}=1−α
P⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ∑ i=1 n (X i −μ) 2 χ 2 α2 (n) <σ 2 <∑ i=1 n (X i −μ) 2 χ 2 1−α2 (n) ⎫ ⎭ ⎬ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ =1−α,
因此,σ 2 的1−α的置信区间为
⎛ ⎝ ⎜ ⎜ ⎜ ∑ i=1 n (X i −μ) 2 χ 2 α2 (n) ,∑ i=1 n (X i −μ) 2 χ 2 1−α2 (n) ⎞ ⎠ ⎟ ⎟ ⎟
2.设μ未知,求σ 2 的置信度为1−α的置信区间
由于
(n−1)S 2 σ 2 ∼χ 2 (n−1)
对于给定置信度1−α,查表可得χ 2 1−α2 (n−1)及χ 2 α2 (n−1),使得
P{χ 2 1−α2 (n−1)<(n−1)S 2 σ 2 <χ 2 α2 (n−1)}=1−α
即P⎧ ⎩ ⎨ (n−1)S 2 χ 2 α2 (n−1) <σ 2 <(n−1)S 2 χ 2 1−α2 (n−1) ⎫ ⎭ ⎬ =1−α
因此,方差σ 2 的置信度为1−α的置信区间为
⎛ ⎝ (n−1)S 2 χ 2 α2 (n−1) ,(n−1)S 2 χ 2 1−α2 (n−1) ⎞ ⎠
而标准差σ的1−α的置信区间为
⎛ ⎝ ⎜ n−1 − − − − − √ Sχ 2 α2 (n−1) − − − − − − − − √ ,n−1 − − − − − √ Sχ 2 1−α2 (n−1) − − − − − − − − − − √ ⎞ ⎠ ⎟
例2.从正态总体N(μ,σ 2 )中抽取容量为5的样本,其观测值为1.86,3.22,1.46,4.01,2.64(1)已知μ=3,求σ 2 的0.95置信区间;(2)如果μ未知,求σ 2 的0.95置信区间.
解:n=5,α=0.05
(1)已知μ=3,查表得χ 2 1−α2 (n)=χ 2 0.975 (5)=0.831,χ 2 α2 (n)=χ 2 0.025 (5)=12.833.由已知数据算得∑ i=1 5 (x i −μ) 2 =4.8693,因此σ 2 的0.95置信区间为
⎛ ⎝ ⎜ ⎜ ⎜ ∑ i=1 n (X i −μ) 2 χ 2 α2 (n) ,∑ i=1 n (X i −μ) 2 χ 2 1−α2 (n) ⎞ ⎠ ⎟ ⎟ ⎟
=(4.869312.833 ,4.86930.831 )=(0.379,5.860)
(2)μ未知,查表得χ 2 1−α2 (n−1)=χ 2 0.975 (4)=0.484,χ 2 α2 (n−1)=χ 2 0.025 (4)=11.143,由已知数据算得x ¯ =2.64,(n−1)s 2 =4.2141,因此σ 2 的0.95置信区间为
⎛ ⎝ ⎜ n−1 − − − − − √ Sχ 2 α2 (n−1) − − − − − − − − √ ,n−1 − − − − − √ Sχ 2 1−α2 (n−1) − − − − − − − − − − √ ⎞ ⎠ ⎟
=(4.214111.143 ,4.21410.484 )=(0.378,8.707)
三、两个正态总体均值差的置信区间
设总体X∼N(μ 1 ,σ 2 1 ),总体Y∼N(μ 2 ,σ 2 2 ),X与Y独立,(X 1 ,X 2 ,⋯,X m )与(Y 1 ,Y 2 ,⋯,Y n )分别来自X与Y的相互独立的样本,并设它们的样本均值分别为X ¯ ¯ ¯ ,Y ¯ ¯ ¯ ,样本方差分别为S 2 1 ,S 2 2 .
1.设σ 2 1 和σ 2 2 都已知,求μ 1 −μ 2 的置信度为1−α的置信区间.
由于X ¯ ¯ ¯ 与Y ¯ ¯ ¯ 相互独立,且X ¯ ¯ ¯ ∼N(μ 1 ,σ 2 1 m ),Y ¯ ¯ ¯ ∼N(μ 2 ,σ 2 2 n ),于是可知X ¯ ¯ ¯ −Y ¯ ¯ ¯ ∼N(μ 1 −μ 2 ,σ 2 1 m +σ 2 2 n ),从而μ=X ¯ ¯ ¯ −Y ¯ ¯ ¯ −(μ 1 −μ 2 )σ 2 1 m +σ 2 2 n − − − − − − − − √ ∼N(0,1)
对于给定的置信度1−α,查表得μ α2 ,使
P⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ |X ¯ ¯ ¯ −Y ¯ ¯ ¯ −(μ 1 −μ 2 )|σ 2 1 m +σ 2 2 n − − − − − − − − √ <μ α2 ⎫ ⎭ ⎬ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ =1−α
即P⎧ ⎩ ⎨ X ¯ ¯ ¯ −Y ¯ ¯ ¯ −σ 2 1 m +σ 2 2 n − − − − − − − − √ μ α2 <μ 1 −μ 2 <X ¯ ¯ ¯ −Y ¯ ¯ ¯ +σ 2 1 m +σ 2 2 n − − − − − − − − √ μ α2 ⎫ ⎭ ⎬ =1−α
从而得到μ 1 −μ 2 的置信度为1−α的置信区间为
⎛ ⎝ X ¯ ¯ ¯ −Y ¯ ¯ ¯ −σ 2 1 m +σ 2 2 n − − − − − − − − √ μ α2 ,X ¯ ¯ ¯ −Y ¯ ¯ ¯ +σ 2 1 m +σ 2 2 n − − − − − − − − √ μ α2 ⎞ ⎠
2.设σ 2 1 =σ 2 2 =σ 2 为未知,求μ 1 −μ 2 的置信度为1−α的置信区间.
(X ¯ ¯ ¯ −Y ¯ ¯ ¯ −t α2 (m+n−2)S w 1m +1n − − − − − − √ ,X ¯ ¯ ¯ −Y ¯ ¯ ¯ +t α2 (m+n−2)S w 1m +1n − − − − − − √ )
其中S 2 w =(m−1)S 2 1 +(n−1)S 2 2 m+n−2
例3.设总体X∼N(μ 1 ,4),总体Y∼N(μ 2 ,6),分别独立地从这两个总体中抽取样本,样本容量依次为16和24,样本均值依次为16.9和15.3,求两个总体均值差μ 1 −μ 2 的置信度为0.95的置信区间.
解:由题设可知m=16,n=24,x ¯ =16.9,y ¯ =15.3,σ 2 1 =4,σ 2 2 =6,1−α=0.95,α=0.05,查附表1得μ α2 =μ 0.025 =1.96.从而可得μ 1 −μ 2 的置信度为0.95的置信区间为
(16.9−15.3−1.96×416 +624 − − − − − − − − √ ,16.9−15.3+1.96×416 +624 − − − − − − − − √ )=(0.214,2.986)
例4.为了估计磷肥对某种农作物增产的作用,选20块条件大致相同的地块进行对比试验.其中10块地施磷肥,另外10块地不施磷肥,得到单位面积的产量(单位:kg)如下:施磷肥:620,570,650,600,630,580,570,600,600,580不施磷肥:560,590,560,570,580,570,600,550,570,550设施磷肥的地块单位面积产量X∼N(μ 1 ,σ 2 ),不施磷肥的地块单位面积产量Y∼N(μ 2 ,σ 2 ).求μ 1 −μ 2 的置信度为0.95的置信区间.
解:由题设,两个正态总体的方差相等,但σ 2 未知,m=10,n=10,α=0.05,1−α=0.95,x ¯ =600,y ¯ =570,s 2 1 =64009 ,s 2 2 =24009 ,s w =(m−1)s 2 1 +(n−1)s 2 2 m+n−2 − − − − − − − − − − − − − − − − − − √ =22,查表得t α2 (m+n−2)=t 0.025 (18)=2.1009,因此μ 1 −μ 2 的置信度为0.95的置信区间为
(x ¯ −y ¯ −s w 1m +1n − − − − − − √ t α2 (m+n−2),x ¯ −y ¯ +s w 1m +1n t α2 (m+n−2) − − − − − − − − − − − − − − − − − − √ )
=(600−570−22×110 +110 − − − − − − − − √ ×2.1009,600−570+22×110 +110 − − − − − − − − √ ×2.1009)
=(9,51)
四、两个正态总体方差比的置信区间
1.设μ 1 ,μ 2 均已知,求σ 2 1 /σ 2 2 的置信水平为1−α的置信区间.
由于1σ 2 1 ∑ i=1 m (X i −μ 1 ) 2 与1σ 2 2 ∑ j=1 n (Y j −μ 2 ) 2 相互独立,且
1σ 2 1 ∑ i=1 m (X i −μ 1 ) 2 ∼χ 2 (m),
1σ 2 2 ∑ j=1 n (X j −μ 2 ) 2 ∼χ 2 (n),
从而可知
F=∑ i=1 m (X i −μ 1 ) 2 ∑ j=1 n (Y j −μ 2 ) 2 ⋅nσ 2 2 mσ 2 1 ∼F(m,n)
对于给定的置信度1−α,查表得F 1−α2 (m,n)及F α2 (m,n),使
P⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ F 1−α2 (m,n)<∑ i=1 m (X i −μ 1 ) 2 ∑ j=1 n (Y j −μ 2 ) 2 ⋅nσ 2 2 mσ 2 1 <F α2 (m,n)⎫ ⎭ ⎬ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ =1−α
即
因此,σ 2 1 σ 2 2 的置信度为1−α的置信区间为
⎛ ⎝ ⎜ ⎜ ⎜ ⎜ n∑ i=1 m (X i −μ 1 ) 2 m∑ j=1 n (Y j −μ 2 ) 2 1F α2 (m,n) ,n∑ i=1 m (X i −μ 1 ) 2 m∑ j=1 n (Y j −μ 2 ) 2 1F 1−α2 (m,n) ⎞ ⎠ ⎟ ⎟ ⎟ ⎟
2.设μ 1 与μ 2 都未知,求σ 2 1 /σ 2 2 的置信水平为1−α的置信区间
⎛ ⎝ S 2 1 S 2 2 ⋅1F α2 (m−1,n−1) ,S 2 1 S 2 2 F α2 (n−1,m−1)⎞ ⎠
例5.设总体X∼N(24,σ 2 1 ),总体Y∼N(20,σ 2 2 ).从总体X和Y中独立地抽的样本值如下:
总体X:23,22,26,24,22,25;
总体Y:22,18,19,23,17
求σ 2 1 /σ 2 2 的置信度为0.95的置信区间.
解:已知μ 1 =24,m=6;μ 2 =20,n=5.由已知数据可算得
∑ i=1 6 (x i −24) 2 =14,∑ j=1 5 (y j −20) 2 =27.
因1−α=0.95,故α=0.05.查附表5,可得F 0.025 (6,5)=6.98,F 0.025 (5,6)=5.99
从而可得σ 2 1 /σ 2 2 的置信度为0.95的置信区间为
(5×146×27×6.98 ,5×14×5.996×27 )=(0.06,2.59)
例6.从参数μ 1 ,μ 2 ,σ 2 1 ,σ 2 2 都未知的两正态总体N(μ 1 ,σ 2 1 ),N(μ 2 ,σ 2 2 )中分别独立地抽取样本,它们的样本容量分别为m=10,n=8,样本方差分别为s 2 1 =3.6,s 2 2 =2.8,求二总体方差比σ 2 1 /σ 2 2 的置信度为0.95的置信区间.
解:这里1−α=0.95,α=0.05,查F分布表得:F α2 (m−1,n−1)=F 0.025 (9,7)=4.82,F α2 (n−1,m−1)=F 0.025 (7,9)=4.20,σ 2 1 /σ 2 2 的置信度为0.95的置信区间为
⎛ ⎝ s 2 1 s 2 2 ⋅1F α2 (m−1,n−1) ,s 2 1 s 2 2 ⋅F α2 (n−1,m−1)⎞ ⎠
=(3.62.8 ×14.82 ,3.62.8 ×4.20)
=(0.27,5.42)