1、第第7章章 参数估计参数估计总体所服从的分布类型已知总体所服从的分布类型已知/未知未知估计总体中未知的参数估计总体中未知的参数参数参数估计估计抽样抽样1 参数估计问题是利用从总体抽样得到的信息参数估计问题是利用从总体抽样得到的信息来估计总体的某些参数来估计总体的某些参数. 参数估计参数估计估计废品率估计废品率估计新生儿的体重估计新生儿的体重估计湖中鱼数估计湖中鱼数27.13点估计点估计X1,X2,Xn现从该总体抽样，得到样本现从该总体抽样，得到样本设有一个统计总体，总体的分布函数设有一个统计总体，总体的分布函数向量向量) . 为为 F(x, )，其中，其中 为未知参数为未知参数 ( 可以是可以
2、是 从样本出发构造适当的统计量从样本出发构造适当的统计量1(,)nXX作为参数作为参数 的估计量，即点估计。的估计量，即点估计。 1( ,)nxxnxx,1将将 代入估计量，得到代入估计量，得到 的估计值的估计值 4点估计点估计矩估计矩估计极大似然估计极大似然估计关键问题：如何构造统计量? 1(,)nXX5矩估计矩估计样本样本k阶原点矩阶原点矩11nkkiiAXn总体总体k阶原点矩阶原点矩kkEX 矩估计基本思想矩估计基本思想: 用样本矩估计总体矩用样本矩估计总体矩 .1lim(|()|)1nkikinXPE Xn大数定律：大数定律：K.皮尔逊皮尔逊6设总体的分布函数中含有设总体的分布函数中含
3、有k个未知参数个未知参数 k ,1(1)它的前它的前k阶原点矩都是这阶原点矩都是这k个参数的函数个参数的函数,记为：记为：1()( ,),1,iikE Xgik11212221211211( ,.,)(),1( ,.,)(),.1( ,.,)(),iinkiinkinkkkkigE XXngE XXngE XXn (2)用样本用样本i阶原点矩替换总体阶原点矩替换总体i阶原点矩阶原点矩 7(3) 解方程组，得解方程组，得 i=hi (X1, X2, Xn) (i=1,2,k);则以则以hi (X1, X2, Xn)作为作为i 的估计量的估计量 ，并，并称称hi(X1, X2, Xn)为为i 的矩
4、法估计量，而的矩法估计量，而称称hi(x1, x2, xn) 为为i 的矩法估计值。的矩法估计值。8总体期望、方差的矩估计量分别是样本均值和总体期望、方差的矩估计量分别是样本均值和样本二阶中心矩。样本二阶中心矩。例例1. 设总体设总体X的数学期望和方差分别是的数学期望和方差分别是, 2 ，求，求 , 2的矩估计量。的矩估计量。222222212222211()()()()111()nnniiiiE XE XD XEXE XuXXXnnXXXXSnnn9例例2: 已知某产品的不合格率为已知某产品的不合格率为p, 有简单随有简单随机样本机样本X1 ,X2 , Xn 求求p的矩估计量。的矩估计量。解
5、：解：E(X)=p.11niipXXn10例3：设电话总机在某段时间内接到呼唤的次数服从参数未知的泊松分布，现在收集了如下42个数据：接到呼唤次数 012345出现的频数71012832求未知参数 的矩估计。80404221x11例例4. XU(a,b),由简单随机样本由简单随机样本X1 ,X2 , Xn求求a,b的矩估计量。的矩估计量。解：解：E(X)=(a+b)/2, D(X)=(b-a)2/12.12222111()()()22111()()()1212niiE XababXXnnD XbabaMSn223(1)3(1)nnaXSbXSnn，12 矩法的优点是简单易行矩法的优点是简单易行
6、,并不需要并不需要事先知道总体是什么分布事先知道总体是什么分布 . 缺点是，当总体类型已知时，没有缺点是，当总体类型已知时，没有 充分利用分布提供的信息充分利用分布提供的信息 . 矩法特点分析：矩法特点分析：13极大似然估计极大似然估计 例：例： 设一箱中装有若干个白色和黑色的球，设一箱中装有若干个白色和黑色的球，已知两种球的数目之比为已知两种球的数目之比为3:1或或1:3，现有放回，现有放回地任取地任取3个球，有两个白球，问：白球所占的个球，有两个白球，问：白球所占的比例比例p是多少？是多少？ 14 如果只知道如果只知道0p0的泊松分布，求的泊松分布，求参数参数的极大似然估计量。的极大似然估
7、计量。21例例4. 设设X1,X2,Xn为取自总体为取自总体XU(a, b)的样的样 本本, 求求a, b的极大似然估计量的极大似然估计量. 22回顾：回顾： 设设X1,X2,Xn是取自总体是取自总体X的一个样本的一个样本其它, 010,) 1()(xxxfX求求 的矩估计量和极大似然估计量的矩估计量和极大似然估计量. 其中其中 0,101()( )(1)2E Xxf x dxxx dxX21XX1121解：解：(1)矩估计矩估计 23 (2)极大似然估计极大似然估计) 10()() 1()(1iniinxxLniixnL1ln) 1ln()(ln0)(lndLd令niixn1ln1niixn
8、12ln1247.2 点估计量的评价标准点估计量的评价标准25 评价一个估计量的好坏，不能仅仅依据一评价一个估计量的好坏，不能仅仅依据一次试验的结果，而必须由多次试验结果来衡量次试验的结果，而必须由多次试验结果来衡量 . 即确定估计量好坏必须在大量观察的基础上从即确定估计量好坏必须在大量观察的基础上从统计的意义来评价。统计的意义来评价。 常用的几条标准是：常用的几条标准是：1无偏性无偏性2有效性有效性3一致性一致性26一、一、无偏性无偏性 )(E则称则称 为为 的无偏估计的无偏估计 .),(1nXX 设设是未知参数是未知参数 的估计量，若的估计量，若 27则称则称 较较 有效有效 .1 都是参
9、数都是参数 的无偏估计量，若有的无偏估计量，若有),(122nXX 设设和和),(111nXX 212( )()DD11,()min( )nDD如果对固定的则称是 的有效估计。二、二、有效性有效性28123123123111236111333112663XXXXXXXXX123例例1：设：设X1,X2, X3是来自某总体是来自某总体X的样本，且的样本，且E(X)=,讨论讨论的以下估计量的无偏性和一致性。的以下估计量的无偏性和一致性。292,例例2：设：设X1,X2, Xn是来自某总体是来自某总体X的样本，且的样本，且 , 判断判断 的矩估计量是的矩估计量是否是无偏估计。否是无偏估计。2,DXE
10、X30三、三、一致性一致性(相合性相合性)是参数是参数 的估计量，若有的估计量，若有设设),(1nnnXX lim(|)0lim(|)1nnnnPP即则称则称是参数是参数 的一致估计量的一致估计量.n 3121,(),()nXXE XD XX设样本来自数学期望方差的总体，则 是 的一致估计量。切比雪夫大数定律切比雪夫大数定律伯努利大数定律伯努利大数定律1,npXXXp设总体为参数为 的0-1分布，为样本，则 是 的一致估计量。niinXnP11)|1(|lim327.333置信区间定义：置信区间定义：则称区间则称区间 是是 的置信度为的置信度为 的置的置信区间信区间.,21 121 和分别称为
11、置信下限和置信上限分别称为置信下限和置信上限. ),(2111nXXX ),(2122nXXX )(21 满足满足设设 是是 一个待估参数，给定一个待估参数，给定, 0 若由样本若由样本X1,X2,Xn确定的两个统计量确定的两个统计量1)(21P347.4351. 期望的区间估计期望的区间估计2已知时已知时的置信区间的置信区间2未知时未知时的置信区间的置信区间2. 求方差的区间估计求方差的区间估计 已知时已知时2的置信区间的置信区间 未知时未知时2的置信区间的置信区间单正态总体四种类型的区间估计单正态总体四种类型的区间估计36已知方差，求期望的区间估计已知方差，求期望的区间估计例例1：随机地从
12、一批服从正态分布：随机地从一批服从正态分布N(,0.022)的零的零件件16个，分别测得其长度为：个，分别测得其长度为：2.14 2.10 2.13 2.15 2.13 2.12 2.13 2.102.15 2.12 2.14 2.10 2.13 2.11 2.14 2.11估计该批零件的平均长度估计该批零件的平均长度，并求，并求的置信区间的置信区间(=0.05)37求参数求参数 的置信度为的置信度为 的置信区间的置信区间. 设设X1,Xn是取自是取自 的样本，的样本， ,2已知 ),(2 N 1,22unXunX) 10(/，利用NnX查正态分布表得查正态分布表得,2 u 1|2unXP使使
13、38例例1：随机地从一批服从正态分布：随机地从一批服从正态分布N(,0.022)的零的零件件16个，分别测得其长度为：个，分别测得其长度为：2.14 2.10 2.13 2.15 2.13 2.12 2.13 2.102.15 2.12 2.14 2.10 2.13 2.11 2.14 2.11估计该批零件的平均长度估计该批零件的平均长度，并求，并求的置信区间的置信区间(=0.05)解：解：2.14.2.112.12516X代入得查表,16,02. 0,96. 1025. 02nuu135. 2,115. 222nuXnuX的置信区间为的置信区间为(2.115,2.135).39u求置信区间的
14、步骤求置信区间的步骤(1) 构造仅与待估参数构造仅与待估参数 有关，但分布已知的有关，但分布已知的函数函数U；(2) 给定置信度给定置信度1-，得常数得常数a,b，使，使 PaUb= 1-；(3) 将将aUb变形，使得：变形，使得：),.,(),.,(212211nnXXXXXX12 ( ,)1. 区间就是 的一个置信度为的置信区间(4) 结论结论40方差未知，求期望的区间估计方差未知，求期望的区间估计例例2：随机地从一批服从正态分布：随机地从一批服从正态分布N(, 2)的零件的零件16个，分别测得其长度为：个，分别测得其长度为：2.14 2.10 2.13 2.15 2.13 2.12 2.
15、13 2.102.15 2.12 2.14 2.10 2.13 2.11 2.14 2.11估计该批零件的平均长度估计该批零件的平均长度，并求，并求的置信区间的置信区间(=0.05)41) 1(/ntnSX利用,22nStXnStX1|2tnSXP使使查查t t分布表得分布表得,2t求参数求参数 的置信度为的置信度为 的置信区间的置信区间. 设设X1,Xn是取自是取自 的样本，的样本， ,2未知),(2 N 142例例1：用一个仪表测量某物理量：用一个仪表测量某物理量9次，得到样本均次，得到样本均值为值为56.32，样本标准差为，样本标准差为0.22. 测量标准差测量标准差反映了测量仪表的精度
16、，试求反映了测量仪表的精度，试求的置信水平为的置信水平为0.95的置信区间。的置信区间。 未知未知,求方差的区间估计求方差的区间估计例例2：假设某地区：假设某地区1825岁女青年身高岁女青年身高 现抽取现抽取30名，样本均值为名，样本均值为158cm,样本方差为样本方差为(5cm)2，求，求2的置信水平为的置信水平为95%的区间估计。的区间估计。2( ,)XN 43)1()1(,)1()1(22/1222/2nSnnSn) 1() 1(222nSn利用求参数求参数 的置信度为的置信度为 的置信区间的置信区间. 设设X1,Xn是取自是取自 的样本，的样本， ),(2 N 12 1)1() 1()
17、 1(2222221nSnnP确定分位数确定分位数)1(),1(22/22/1nn使使44例例2：假设某地区：假设某地区1825岁女青年身高岁女青年身高 现抽取现抽取30名，样本均值为名，样本均值为158cm,样本方差为样本方差为(5cm)2，求，求2的置信水平为的置信水平为95%的区间估计。的区间估计。45例例3：随机地从一批服从正态分布：随机地从一批服从正态分布N(2.12, 2)的的零件零件16个，分别测得其长度为：个，分别测得其长度为：2.14 2.10 2.13 2.15 2.13 2.12 2.13 2.102.15 2.12 2.14 2.10 2.13 2.11 2.14 2.
18、11试求试求2的置信水平为的置信水平为0.95的置信区间。的置信区间。 已知已知,求方差的区间估计求方差的区间估计46求参数求参数 的置信度为的置信度为 的置信区间的置信区间. 设设X1,Xn是取自是取自 的样本，的样本， 已知已知 ),(2 N 12221(0,1)( )1,2,.,niiiXXNnin，221122/ 21/ 2),( )( )nniiiiXXnn（22211222)( )( )1niiXPnn （确定分位数确定分位数221/ 2/ 2( ),( )nn使使47例例4. 4. 对飞机的飞行速度进行对飞机的飞行速度进行1515次独立试验，测次独立试验，测得飞机的最大飞行速度得
19、飞机的最大飞行速度( (单位：单位：m/s)m/s)如下：如下：422.2 418.7 425.6 420.3 425.8 423.1 431.5 428.2438.3 434.0 411.3 417.2 413.5 441.3 423.0假设飞机最大飞行速度服从假设飞机最大飞行速度服从 求最求最大飞行速度的方差的置信度为大飞行速度的方差的置信度为0.90的置信区间。的置信区间。解：解：2的置信区间为的置信区间为( 41.407010, 142.549270 )( 41.407010, 142.549270 ).2(424.93,)N2220.0522120.952i15,424.93,1 0
20、.900.1( )(15)24.996,( )(15)7.261,(X) =1035.013504 nnn 查表代入得48n双总体双总体设总体设总体X N(1,12)，总体，总体Y N(2,22)，X1,X2,Xm来自来自X，Y1,Y2,Yn来自来自Y，且两样，且两样本相互独立。本相互独立。均值差均值差1- 2的区间估计的区间估计方差比方差比12/ 22的区间估计的区间估计49例例1.设甲乙两地区女青年身高分别服从设甲乙两地区女青年身高分别服从分布。甲地区抽取分布。甲地区抽取100名，样本均名，样本均值为值为163cm； 乙地区抽取乙地区抽取100名，样本均值为名，样本均值为160cm,求求
21、的置信水平为的置信水平为90%的置信区间。的置信区间。12(,16),(,9)NN1250u1,2已知时已知时, ,1- 2的置信区间的置信区间即得即得1- 2的置信区间的置信区间nmuYXnmuYX2221222212,122212()(0,1)XYUNmn51例例2：今抽样甲乙两地区：今抽样甲乙两地区1825岁女青年身高的数岁女青年身高的数据如下：甲地区抽取据如下：甲地区抽取50名名, 样本均值为样本均值为163cm,样本样本标准差为标准差为4cm; 乙地区抽取乙地区抽取50名名, 样本均值为样本均值为160cm,样本标准差为样本标准差为3cm。假设身高均服从正态分布并具。假设身高均服从正
22、态分布并具有公共方差，求有公共方差，求 的置信水平为的置信水平为90%的置信的置信区间。区间。1252u1,2未知未知, ,但但1=2=时时, ,1- 2的的 置信区间置信区间)2(11)(21nmtnmSYXTw2) 1() 1(2221nmSnSmSw其中53nmSnmtYXnmSnmtYXww11)2(,11)2(221- 2的置信区间：的置信区间：54u1,2未知未知, ,且且12，但容量，但容量m,n很大很大时时, , 1- 2的置信区间的置信区间221222211221)(11)(11niiniiYYnSXXmS近似代替以近似代替以nSmSuYXnSmSuYX2221222212,
23、55例例7：两位化验员：两位化验员A,B独立对某种化合物的含氯独立对某种化合物的含氯量用相同方法各作了量用相同方法各作了10次测量，测量的样本方差次测量，测量的样本方差分别为分别为S12=0.5419, S22=0.6065，设，设A,B所测量值所测量值总体为总体为X,Y，并且均服从正态分布，方差分别为，并且均服从正态分布，方差分别为22221212/0.95，求方差比的置信度为的置信区间例例8：假设人体身高服从正态分布，今抽样甲乙两：假设人体身高服从正态分布，今抽样甲乙两地区地区1825岁女青年身高的数据如下：甲地区抽取岁女青年身高的数据如下：甲地区抽取121名，样本均值为名，样本均值为16
24、4cm,样本标准差为样本标准差为4cm;乙地乙地区抽取区抽取61名，样本均值为名，样本均值为160cm,样本标准差为样本标准差为2cm，求两总体方差比的置信水平为求两总体方差比的置信水平为95%的区间估计。的区间估计。56uu1, , u2未知，方差比未知，方差比12/22的区间估计的区间估计) 1, 1(22222121nmFSSF1) 1, 1() 1, 1(22222212121nmFSSnmFP1) 1, 1(1) 1, 1(1212221222122221nmFSSnmFSSP57uu1, , u2已知，方差比已知，方差比12/22的区间估计的区间估计2211122221()( ,
25、)()miiniiXFF m nY587.5*59例例12：随机地从一批寿命服从正态分布：随机地从一批寿命服从正态分布N(,2)的的灯泡中抽取灯泡中抽取16个进行试验，分别测得样本平均个进行试验，分别测得样本平均寿命为寿命为1509.5小时小时,样本方差为样本方差为32.232. 求该批灯泡的平均寿命求该批灯泡的平均寿命的的95%置信下限。置信下限。60单侧置信区间单侧置信区间设设X分布函数为分布函数为F(x; ), 未知，给定未知，给定 (01),若由样本若由样本 X1,X2, ,Xn确定的统计确定的统计量量 ，若满足：，若满足：),.,(211nXXX12()1()1PP 或者12( ,)(,)则称随机区间或为参数 的置信度为置信度为1- 的单侧置信区间，的单侧置信区间， 分别为单分别为单侧置信下限和单侧置信上限侧置信下限和单侧置信上限12，61例例12：随机地从一批寿命服从正态分布：随机地从一批寿命服从正态分布N(,2)的的灯泡中抽取灯泡中抽取16个进行试验，分别测得样本平均个进行试验，分别测得样本平均寿命为寿命为1509.5小时小时,样本方差为样本方差为32.232.求该批灯泡的平均寿命求该批灯泡的平均寿命的的95%置信下限。置信下限。SXtn单侧置信区间（，）) 1(/ntnSX利用查正态分布表得查正态分布表得1tnSXP62
