正态总体参数的区间估计

  • 格式:doc
  • 大小:561.00 KB
  • 文档页数:7

1 第19讲 正态总体参数的区间估计

教学目的:理解区间估计的概念,掌握各种条件下对一个正态总体的均值和方差进行区间估计的方法。

教学重点:置信区间的确定。

教学难点:对置信区间的理解。

教学时数: 2学时。

教学过程:

第六章 参数估计

§6.3正态总体参数的区间估计

1. 区间估计的概念

我们已经讨论了参数的点估计,但是对于一个估计量,人们在测量或计算时,常不以得到近似值为满足,还需估计误差,即要求知道近似值的精确程度。因此,对于未知参数,除了求出它的点估计ˆ外,我们还希望估计出一个范围,并希望知道这个范围包含参数真值的可信程度。

设ˆ为未知参数的估计量,其误差小于某个正数的概率为1(01),即

ˆ{||}1P

1)ˆˆ(P

这表明,随机区间)ˆ,ˆ(包含参数真值的概率(可信程度)为1,则这个区间)ˆ,ˆ(就称为置信区间,1称为置信水平。

定义 设总体X的分布中含有一个未知参数。若对于给定的概率1(01),存在两个统计量1112(,,,)nXXX与2212(,,,)nXXX,使得

12{}1P 2 则随机区间12(,)称为参数的置信水平为1的置信区间,1称为置信下限,2称为置信上限,1称为置信水平。

注(1)置信区间的含义:若反复抽样多次(各次的样本容量相等,均为n),每一组样本值确定一个区间12(,),每个这样的区间要么包含的真值,要么不包含的真值。按伯努利大数定理,在这么多的区间中,包含真值的约占100(1)%,不包含真值的约仅占100%。例如:若0.01,反复抽样1000次,则得到的1000个区间中,不包含真值的约为10个。

(2)置信区间的长度表示估计结果的精确性,而置信水平表示估计结果的可靠性。对于置信水平为1的置信区间12(,),一方面置信水平1越大,估计的可靠性越高;另一方面区间12(,)的长度(2)越小,估计的精确性越好。但这两方面通常是矛盾的,提高可靠性通常会使精确性下降(区间长度变大),而提高精确性通常会使可靠性下降(1变小),所以要找两方面的平衡点。

在学习区间估计方法之前,我们先介绍标准正态分布的分位点概念。

设~0,1XN,若z满足条件,01PXz,则称点z为标准正态分布的分位点。例如求0.01z。按照分位点定义,我们有0.010.01PXz,则0.010.99PXz,即0.01()0.99z。查表可得0.012.327z. 又由()x图形的对称性知1zz。下面列出了几个常用的z值:

 0.001 0.005 0.01 0.025 0.05 0.10

z 3.090 2.576 2.327 1.960 1.645 1.282

2. 正态总体均值的区间估计

设已给定置信水平为1,总体2~,XN,12,,,nXXX为一个样本,2,XS分别是样本均值和样本方差。 3 (1)2已知时,的置信区间

我们知道X是的无偏估计,且有统计量/Xn~0,1N 。由标准正态分布的上分位点的定义,有

/21/XPzn

/2/21PXzXznn

这样,我们就得到了的一个置信水平为1的置信区间

/2/2,XzXznn

这样的置信区间常写成

/2Xzn

例1 从某厂生产的滚珠中随机抽取10个,测得滚珠的直径(单位:mm)如下:

14.6 15.0 14.7 15.1 14.9 14.8 15.0 15.1 15.2 14.8

若滚珠直径服从正态分布2(,)N,并且已知0.16(mm),求滚珠直径均值的置信水平为95%的置信区间。

解 计算样本均值14.92x,置信水平1=0.95,0.05,查表得/20.0251.96zz(可利用()zt查表)。由此得的置信水平为95%的置信区间为

/20.1614.921.9610Xzn

14.920.099, 14.920.099(14.821, 15.019)

注:置信水平为1的置信区间并不是唯一的。以例1来说,给定0.05,则又有

0.040.010.95/XPzzn

故 4 0.010.04,XzXznn

也是的置信水平为95%的置信区间,其区间长度为0.040.01()4.08zznn。而在对称区间0.05Xzn上,区间长度为23.92nn, 比非对称区间长度要短,较优。易知,像(0,1)N分布那样其概率密度的图形是单峰且对称的情况,当n固定时,以对称区间其长度为最短,我们选用对称区间。

(2)2未知时,的置信区间

此时不能使用/2Xzn,因为其中包含了未知参数。考虑到2S是2的无偏估计,将上述区间中的换成2SS。我们已知统计量/XSn~(1)tn,可得

/2/2(1)(1)1/XPtntnSn

/2/2(1)(1)1SSPXtnXtnnn

于是得到的一个置信水平为1的置信区间

/2(1)SXtnn

例2 在例1中,若未知,求滚珠直径均值的置信水平为95%的置信区间。

解 计算样本均值14.92x,样本标准差0.193s;置信水平1=0.95,0.05,自由度11019n,查表得/20.025(1)(9)2.26tnt。 由此得的置信水平为95%的置信区间为 5 /20.193(1)14.922.2610SXtnn

(14.92-0.138,14.92+0.138)=(14.782,15.058)

注 比较例1和例2中的置信区间,可以发现当2未知时,的置信区间区间长度要比2已知时的置信区间区间长度大,这表明当未知条件增多时,估计的精确程度变差,这也符合我们的直观感觉。

3. 正态总体方差2的区间估计

(1)已知时,2的置信区间

已知2211()niiX2~()n 但是2分布的概率密度图形不是对称的,对于已给的置信水平1,要想找到最短的置信区间是困难的。因此,习惯上仍然取对称的分位点21/2和2/2可得

2221/2/2211()()()1niiPnXn

2221122/21/2()()1()()nniiiiXXPnn

于是得到方差2的一个置信水平为1的置信区间

221122/21/2()(), ()()nniiiiXXnn

例3 在例1中,若已知14.9(mm),求滚珠直径方差2的置信水平为95%的置信区间。

解 已知14.9,置信水平1=0.95,0.05,自由度10n,查表得 6 22/20.025()(10)20.5n,221/20.975()(10)3.25n。

则方差2的置信水平为95%的置信区间为

10102222111122/21/2()()(14.9)(14.9), , ()()20.53.25nniiiiiiiiXXxxnn

0.340.34, (0.0166, 0.1046)20.53.25

(2)未知时,2的置信区间

2的无偏估计为2S,且统计量22(1)nS2~(1)n。选取分位点21/2和2/2可得

2221/2/22(1)(1)(1)1nSPnn

22222/21/2(1)(1)1(1)(1)nSnSPnn

于是得到方差2的一个置信水平为1的置信区间

2222/21/2(1)(1), (1)(1)nSnSnn

由此,我们还可以得到标准差的一个置信水平为1的置信区间

222222/21/2/21/2(1)(1)(1)(1), , (1)(1)(1)(1)nSnSnSnSnnnn

注 在实际问题中,对2做估计的时候,一般均是未知的情况。因此,我们重点掌握未知条件下求2的置信区间问题。

例4 在例1中,若未知,求滚珠直径方差2的置信水平为95%的置信区间。

解 未知,计算样本方差20.0373s,置信水平1=0.95,0.05,自由度 7 19n,查表可得22/20.025(1)(9)19.0n,221/20.975(1)(9)2.70n。

则方差2的置信水平为95%的置信区间为

2222/21/2(1)(1)90.037390.0373, , (1)(1)19.02.70nSnSnn

(0.0177,0.1243)