第24卷第2期许昌学院学报Vol24No.2加05年3月JOURNALOFXUCHANGUNIVERSITYMar..2005文章墒号:167l一9824(2005)02~0018一‘)4
一致最小方差无偏估计的判定及其求法
王芬玲,樊明智,
(许昌学院数学系,河南许昌461000)
摘要:给出了一致最小方差无偏估计的判定方法及存在的充要条件,并给出了一致最小方
差无偏估计的几种求法,最后通过实例说明了这些方法的应用.
关键词:一致最小方差无偏估计;无偏估计;充分统计量;C—R下界
中图分类号:0212.1文献标识码:A
O引言
当样本容量n>1时,可估参数(统计中,一般将存在无偏估计的参数函数称为可估参数)的无偏估计不惟一设g(∞为可估参数,我们把g(∞的所有无偏估计组成的类记为岵,如何从坛中选取一个较好的估计是我们关心的问题.根据均方误差原则,一个无偏估计其方差越小越好,对固定大小的样本来说,估计餐的方差不能任意小.那么,在所有无偏估计中,估计量的方差下界是什么?有无达到这个下界的估计量?如存在,如何把它找出来?为此本文将给出一致最小方差无偏估计的判定方法及其求法.
1主要结果
定义设g(∞为可估参数,如果r(x)是g(口)的无偏估计,且对略中任一个估讹妒(x),有
-缸b(r(x))≤协目(驴(x))V口∈@
则称为r(x)为g(∞的一致最小方差无偏估计(uniforrnlyMinif肌mVadanceunbiasedEstimate),简记为UMVUEflI.
引理l设s(x)是分布族{,(x,口),臼∈@}的充分统计量,妒(盖)是g(疗)的无偏估计t令r(x)=E(妒(x)ls(x)),则r(x)也是g(目)的无偏估计,且
%埽(r(x))≤忆印(妒(卫))V疗∈@
引理2设s(丑)是分布族{,(丑,口),日∈@}的完备充分统计量,g(口)为可估参数,则g(目)的uMVuE存在,它是s(Ⅳ)的函数且在几乎处处意义下是惟一的.
定理1设丁(.jf)是g(∞的无偏估计,%啊(丁(x))(∞,则r(引为g(日)的uMVuE的充要条件是对任~O的无偏估计P(膏),若忆r(P(x))<m,则有c锄(f(x),r(趵)=O
证明令
“={妒(置):岛(妒(石))=g(∞,仇御(x)<+∞,V曰∈@f;
‰={妒(盖):局(P(x))=0,%啊驴(盖)<+∞,V口∈@},
设r(x)是g(鲫的uMVuE,任取驴(x)∈‰及^∈R,则
妒’(x)=^妒(置)+r(x)∈U,且有沌印(r(x))≤ya唯(A驴(x)+r(x))
即^2-缸妇(垆(x))+2^co蜥(p(X),r(_]|『))≥0.
由^∈月的任意性知:c鲫日(妒(x),r(x))=o,V臼∈@、
收稿日期:2004一11—19作者简介:王芬玲(1968一),士,河南都陵人,讲师,研究方向:基础数学
第24卷第2期王芬玲,等:一致最小方差无偏估计的判定及其求法19反之,设妒(Ⅳ)∈%都有‰口(妒(剐,丁(x))=oV口∈@要证r(x)是g(∞的uMVuE,若妒’(并)∈u,则有r(x)一妒’(x)∈【,o,由假设条件得
c。%(r(x)一∥(卫),r(x))=0
易((r(x)一97(x))r(x))=0
岛(r(x))一E。(r(x))妒’(x))=0
由许瓦兹不等式得:一
(Fd(产(Ⅳ)))2=(E口(r(x)驴’(x)))2≤岛(r(x))岛((妒’(x))2)
从而
E口(严(x)))≤E口((妒’(x))2)
又
¨E。(丁(x))=F日(节7(x))=g(8),.‘.V廿日(T(x))≤I缸7。(于’(x))
由(妒’(x))∈u之任意性可知,r(膏)是g(∞的uMVuE.
定理2设总体分布族"(x,日),口∈@}满足如下条件:
(ii)分布旌有共同支撑^,且对任意x∈一,v睁∈@,掣存在;(i)@为开区间;
(iii)岛(掣)=』掣出-0'v日∈@;(Iv),(口):岛(i掣)2:』(曼掣)了(石,口)如,ov日∈@
x。,x2,…,x。i.i.d.是抽自总体的一个样本,i(x)是g(目)的无偏估计,且满足
刍』...』;(x。,・・・x。)(rl,(x。,口))dx,-・dx。:』・・・』;(x,,…z。)品(J●,(z。,目))dz。・-・dx。
则有
’
‰印(i(J))≥上i嘉;字v日∈@(1)
把(1)式称为cramer-R帅不等式,简称c.R不等式,[g7(口)]2/n,(口)称为c-RF界.
定理2的证明见文献[2]或[3].
2uMvuE的几种求法
由r述讨论,只要完备充分统}|量存在,可估参数uMvuE一定存在.因为充分统计量有一个重要作用是降低无偏估计的方差.
引理l、2提供了两种求uMvuE的方法.
方法l若s(x)为完备充分统计量,寻求函数^(s(x))为g(∞的无偏估计,则^Cs(x))为g(∞的UMVUE
例1殴样本x∥・x.是来自U(0,p)的一个样本,求疗的uMvuE
解由凶子分解定理知x(。)=rrm(盖∥一盖。)是克分统计量,其分布密度函数为
p(f,∞=m“一1/∥,0<f<口
可验证该分布族是完备的,因而x㈨也是8的完备充分统计量.
因为删(。)=J:芳也=了鲁目,故生}x(。)是口的uMVuE.
方法2首先找一个无偏估计p(x),然后将之对完备充分统计昔求条件期望,则r(x)=E(P(x)Is(x))是譬(∞的uMVuE.
例2设某种产品的寿命服从带有未知参数A的指数分布,对固定时问to,产品在t。前失效的概率肌(x≤fo)=J—e“~兮g(^),求g(^)的uMVuE.
许昌学院学报2005年3月
解先寻找一个完箭允分统计量.
抽取,t个产品进行试验,得样本盖∥・・x。,由指数族的性质知:
s(x)=∑x。为^的一个完备充分统计量.
令学(x)=,(J.;%)足g(^)的无偏估计,因为E(妒(盖))=肌(xl≤to)=g(^),则在s(置)=s的条件下求条件期望:
,,
当oc詈≤1时,有
E(P(鼻)ls(x)=s?=,(x。≤州s(丑)=s)=p(志≤导)
=fcn一,,c-一z,“一2at=-一(・一导)“一1
当导≥1时,E(妒(x)|5(x)=s)=1,
故g(^)的uMVuE为:
咖):『1_(・一高)”1瑚果s(如“.
1l,如果s(x)<“
方法3利用定理1,uMvuE存在的充要条件来寻求uMvuE.
例3没样本x∥‘’x。,i.id是来自N(p,。2)的一个样本,证明i=吉蚤x,及s2=_与善(x.一x)2分别为肚,a2的uMVuF
证明令参数守问@={(严,口2)l卢∈尺,口2>0},V臼∈@,氏∈u,有:
印。=(2—2)一{J…j8。(x∥‘’,%)exp【一刍蚤(置一p)2】出r…如。=o(2)
两边关于Ⅳ求导,得:。
(z一2)一;卜・』蹦m…^)妄窨(≈一p)“一【一去砉(≈一严)2]岫…出。=。(3)
因此
●
(2—2)一;』…』岛(z∥..,z。)iexp[_刍骞(毛一岸)2]出,…如。=o、(4)
即毋(a。i):o,又困E孑=p,由定理1知i:÷∑x。及为Ⅳ的uMvuE.
(3)式两边关于“求导,得:
(2—2)一j』¨・』蹦m…^)[骞(戡一一)】2e砷【一刍耋(‰一一)2】也,…出。=。(5)
(2)式两边关于一2求导,得:
(2—2)一;j.-fd。(x∥‘‘,“)骞(甄一一)2ex一[一刍薯(z。一一)2】出-…如。=。(6)
(6)式减去上倍的(5)式,得:
(z一2)一;卜-』蹦钆…,剞耋(她一i)2e印【一去骞(*。一一)2】出,…也。=。
即岛(%52):0,V臼∈@,又因E庐2=口2,由定理1知52及为d2的uMVuE.
方法4利用定理2,c.R不等式提供了一种验证某些估计是uMvuE的方法.在定理2的条件成立时,由c—R不等式知.无偏估计的方差不低于c—R下界.因此,若某个无偏估计的方差达到了c.R下界,则它必
第24卷第2期王芬玲,等:一致最小方差无偏估计的判定及其求法
是一个UMvUE.
倒4设样本J∥‘‘咒,i.i.d是来自口(1,P),求p的uMVuE.
解设总体分布为,(g,p)=矿(1一p)1~,*=o,l,o<p<1
容易验证;,(x,p)lp∈p(o,1)}满足定理2的条件
因为:
故P的无偏估计的c—R下界为ib=掣.‰)=耳[毒・帕川卜乓【者与]2:南%(并-p)2=万尚
i作为p的无偏估汁,有:
%砸)=去耋地。(置)丢删刊=世≯
i达到了c—R下界,故i为p的uMvuE.
方法5通过解方程来求uMVuE.若s(工)为完备充分统计量,可估两数F(口)的uMvuE.是方程E(^(S(x)))=g(日)的唯一解.
如例1也可采用方法5来解决,即先找出^(5(x))为日的uMVuE,解方程:E(^(s(x)))=口,再求出^(5(x)).
因为盖(。)也是目的完备充分统计量,记_:r㈨取值为‘,设^(f)是口的元偏估计(即为目的uMvuE)
r日
从而有E(h(£))=l矗(£)-nP。1/扫”出=臼,v口>o
JO
r日r日
即:l^(t)n,。1出:俨¨=I(n+1)‘q£,v口>o
对上式两边求导得:
^(£)・n£8—1:(n+1)r8日U^(£):!!型f,V口>0n
-
因而^(5(x)):型x㈨是口的uMvuE.
涟:方法1比较简单,但适用范匍有一定的局限性;方法2虽计算量较大,但适用范围广,是一种普遍使用的方法;方法3需要构造。的无偏估计,方法灵活但需要一定的技巧;方法4要求密度函数满足一定的条件(即正则条件1—4)方可使用;方法5足一种不常用的方法,但在完备充分统计量s(丑)容身确定时,利用参数估计的无偏性建立等式还是比较简单的.
参考文献
[1]茆诗松,王静龙,濮晓龙高等数理统计[M].北京:高等教育出版社,1998
[2]中山大学教学力学系概率论与数理统计(下)[MJ.北京:高等教育出版社.1980
[3]华东师范大学数学系概率论与数理统计教程[M]北京:高等教育出版社,1983
责任编校:周伦
JudgementandSolutionMethodsofUMVUE
WANGFen。ling,FANMing‘zhi
(上卸Ⅱ玎舢眦0,肘珊^删n£溉,x∽矗。昭‰洫倦毋,x眦^o昭461000,现iM)
Abstract:Thi8pape。pmvidesmethodsofjudgement0fUMVUE,necessaryandsufficiemconditionofexistence,severals0】u石ons抽UMVUE,ande。p】ajns£b。pr靴出丑l印p】ica石on。f出eseme出ods呐aafewe髓“1ples.
Keywords:UMVUE;Unbiasedestimate;su珩c;entstaIist记;C—R10werbound
第24卷第2期许昌学院学报Vol24No.2加05年3月JOURNALOFXUCHANGUNIVERSITYMar..2005文章墒号:167l一9824(2005)02~0018一‘)4
一致最小方差无偏估计的判定及其求法
王芬玲,樊明智,
(许昌学院数学系,河南许昌461000)
摘要:给出了一致最小方差无偏估计的判定方法及存在的充要条件,并给出了一致最小方
差无偏估计的几种求法,最后通过实例说明了这些方法的应用.
关键词:一致最小方差无偏估计;无偏估计;充分统计量;C—R下界
中图分类号:0212.1文献标识码:A
O引言
当样本容量n>1时,可估参数(统计中,一般将存在无偏估计的参数函数称为可估参数)的无偏估计不惟一设g(∞为可估参数,我们把g(∞的所有无偏估计组成的类记为岵,如何从坛中选取一个较好的估计是我们关心的问题.根据均方误差原则,一个无偏估计其方差越小越好,对固定大小的样本来说,估计餐的方差不能任意小.那么,在所有无偏估计中,估计量的方差下界是什么?有无达到这个下界的估计量?如存在,如何把它找出来?为此本文将给出一致最小方差无偏估计的判定方法及其求法.
1主要结果
定义设g(∞为可估参数,如果r(x)是g(口)的无偏估计,且对略中任一个估讹妒(x),有
-缸b(r(x))≤协目(驴(x))V口∈@
则称为r(x)为g(∞的一致最小方差无偏估计(uniforrnlyMinif肌mVadanceunbiasedEstimate),简记为UMVUEflI.
引理l设s(x)是分布族{,(x,口),臼∈@}的充分统计量,妒(盖)是g(疗)的无偏估计t令r(x)=E(妒(x)ls(x)),则r(x)也是g(目)的无偏估计,且
%埽(r(x))≤忆印(妒(卫))V疗∈@
引理2设s(丑)是分布族{,(丑,口),日∈@}的完备充分统计量,g(口)为可估参数,则g(目)的uMVuE存在,它是s(Ⅳ)的函数且在几乎处处意义下是惟一的.
定理1设丁(.jf)是g(∞的无偏估计,%啊(丁(x))(∞,则r(引为g(日)的uMVuE的充要条件是对任~O的无偏估计P(膏),若忆r(P(x))<m,则有c锄(f(x),r(趵)=O
证明令
“={妒(置):岛(妒(石))=g(∞,仇御(x)<+∞,V曰∈@f;
‰={妒(盖):局(P(x))=0,%啊驴(盖)<+∞,V口∈@},
设r(x)是g(鲫的uMVuE,任取驴(x)∈‰及^∈R,则
妒’(x)=^妒(置)+r(x)∈U,且有沌印(r(x))≤ya唯(A驴(x)+r(x))
即^2-缸妇(垆(x))+2^co蜥(p(X),r(_]|『))≥0.
由^∈月的任意性知:c鲫日(妒(x),r(x))=o,V臼∈@、
收稿日期:2004一11—19作者简介:王芬玲(1968一),士,河南都陵人,讲师,研究方向:基础数学
第24卷第2期王芬玲,等:一致最小方差无偏估计的判定及其求法19反之,设妒(Ⅳ)∈%都有‰口(妒(剐,丁(x))=oV口∈@要证r(x)是g(∞的uMVuE,若妒’(并)∈u,则有r(x)一妒’(x)∈【,o,由假设条件得
c。%(r(x)一∥(卫),r(x))=0
易((r(x)一97(x))r(x))=0
岛(r(x))一E。(r(x))妒’(x))=0
由许瓦兹不等式得:一
(Fd(产(Ⅳ)))2=(E口(r(x)驴’(x)))2≤岛(r(x))岛((妒’(x))2)
从而
E口(严(x)))≤E口((妒’(x))2)
又
¨E。(丁(x))=F日(节7(x))=g(8),.‘.V廿日(T(x))≤I缸7。(于’(x))
由(妒’(x))∈u之任意性可知,r(膏)是g(∞的uMVuE.
定理2设总体分布族"(x,日),口∈@}满足如下条件:
(ii)分布旌有共同支撑^,且对任意x∈一,v睁∈@,掣存在;(i)@为开区间;
(iii)岛(掣)=』掣出-0'v日∈@;(Iv),(口):岛(i掣)2:』(曼掣)了(石,口)如,ov日∈@
x。,x2,…,x。i.i.d.是抽自总体的一个样本,i(x)是g(目)的无偏估计,且满足
刍』...』;(x。,・・・x。)(rl,(x。,口))dx,-・dx。:』・・・』;(x,,…z。)品(J●,(z。,目))dz。・-・dx。
则有
’
‰印(i(J))≥上i嘉;字v日∈@(1)
把(1)式称为cramer-R帅不等式,简称c.R不等式,[g7(口)]2/n,(口)称为c-RF界.
定理2的证明见文献[2]或[3].
2uMvuE的几种求法
由r述讨论,只要完备充分统}|量存在,可估参数uMvuE一定存在.因为充分统计量有一个重要作用是降低无偏估计的方差.
引理l、2提供了两种求uMvuE的方法.
方法l若s(x)为完备充分统计量,寻求函数^(s(x))为g(∞的无偏估计,则^Cs(x))为g(∞的UMVUE
例1殴样本x∥・x.是来自U(0,p)的一个样本,求疗的uMvuE
解由凶子分解定理知x(。)=rrm(盖∥一盖。)是克分统计量,其分布密度函数为
p(f,∞=m“一1/∥,0<f<口
可验证该分布族是完备的,因而x㈨也是8的完备充分统计量.
因为删(。)=J:芳也=了鲁目,故生}x(。)是口的uMVuE.
方法2首先找一个无偏估计p(x),然后将之对完备充分统计昔求条件期望,则r(x)=E(P(x)Is(x))是譬(∞的uMVuE.
例2设某种产品的寿命服从带有未知参数A的指数分布,对固定时问to,产品在t。前失效的概率肌(x≤fo)=J—e“~兮g(^),求g(^)的uMVuE.
许昌学院学报2005年3月
解先寻找一个完箭允分统计量.
抽取,t个产品进行试验,得样本盖∥・・x。,由指数族的性质知:
s(x)=∑x。为^的一个完备充分统计量.
令学(x)=,(J.;%)足g(^)的无偏估计,因为E(妒(盖))=肌(xl≤to)=g(^),则在s(置)=s的条件下求条件期望:
,,
当oc詈≤1时,有
E(P(鼻)ls(x)=s?=,(x。≤州s(丑)=s)=p(志≤导)
=fcn一,,c-一z,“一2at=-一(・一导)“一1
当导≥1时,E(妒(x)|5(x)=s)=1,
故g(^)的uMVuE为:
咖):『1_(・一高)”1瑚果s(如“.
1l,如果s(x)<“
方法3利用定理1,uMvuE存在的充要条件来寻求uMvuE.
例3没样本x∥‘’x。,i.id是来自N(p,。2)的一个样本,证明i=吉蚤x,及s2=_与善(x.一x)2分别为肚,a2的uMVuF
证明令参数守问@={(严,口2)l卢∈尺,口2>0},V臼∈@,氏∈u,有:
印。=(2—2)一{J…j8。(x∥‘’,%)exp【一刍蚤(置一p)2】出r…如。=o(2)
两边关于Ⅳ求导,得:。
(z一2)一;卜・』蹦m…^)妄窨(≈一p)“一【一去砉(≈一严)2]岫…出。=。(3)
因此
●
(2—2)一;』…』岛(z∥..,z。)iexp[_刍骞(毛一岸)2]出,…如。=o、(4)
即毋(a。i):o,又困E孑=p,由定理1知i:÷∑x。及为Ⅳ的uMvuE.
(3)式两边关于“求导,得:
(2—2)一j』¨・』蹦m…^)[骞(戡一一)】2e砷【一刍耋(‰一一)2】也,…出。=。(5)
(2)式两边关于一2求导,得:
(2—2)一;j.-fd。(x∥‘‘,“)骞(甄一一)2ex一[一刍薯(z。一一)2】出-…如。=。(6)
(6)式减去上倍的(5)式,得:
(z一2)一;卜-』蹦钆…,剞耋(她一i)2e印【一去骞(*。一一)2】出,…也。=。
即岛(%52):0,V臼∈@,又因E庐2=口2,由定理1知52及为d2的uMVuE.
方法4利用定理2,c.R不等式提供了一种验证某些估计是uMvuE的方法.在定理2的条件成立时,由c—R不等式知.无偏估计的方差不低于c—R下界.因此,若某个无偏估计的方差达到了c.R下界,则它必
第24卷第2期王芬玲,等:一致最小方差无偏估计的判定及其求法
是一个UMvUE.
倒4设样本J∥‘‘咒,i.i.d是来自口(1,P),求p的uMVuE.
解设总体分布为,(g,p)=矿(1一p)1~,*=o,l,o<p<1
容易验证;,(x,p)lp∈p(o,1)}满足定理2的条件
因为:
故P的无偏估计的c—R下界为ib=掣.‰)=耳[毒・帕川卜乓【者与]2:南%(并-p)2=万尚
i作为p的无偏估汁,有:
%砸)=去耋地。(置)丢删刊=世≯
i达到了c—R下界,故i为p的uMvuE.
方法5通过解方程来求uMVuE.若s(工)为完备充分统计量,可估两数F(口)的uMvuE.是方程E(^(S(x)))=g(日)的唯一解.
如例1也可采用方法5来解决,即先找出^(5(x))为日的uMVuE,解方程:E(^(s(x)))=口,再求出^(5(x)).
因为盖(。)也是目的完备充分统计量,记_:r㈨取值为‘,设^(f)是口的元偏估计(即为目的uMvuE)
r日
从而有E(h(£))=l矗(£)-nP。1/扫”出=臼,v口>o
JO
r日r日
即:l^(t)n,。1出:俨¨=I(n+1)‘q£,v口>o
对上式两边求导得:
^(£)・n£8—1:(n+1)r8日U^(£):!!型f,V口>0n
-
因而^(5(x)):型x㈨是口的uMvuE.
涟:方法1比较简单,但适用范匍有一定的局限性;方法2虽计算量较大,但适用范围广,是一种普遍使用的方法;方法3需要构造。的无偏估计,方法灵活但需要一定的技巧;方法4要求密度函数满足一定的条件(即正则条件1—4)方可使用;方法5足一种不常用的方法,但在完备充分统计量s(丑)容身确定时,利用参数估计的无偏性建立等式还是比较简单的.
参考文献
[1]茆诗松,王静龙,濮晓龙高等数理统计[M].北京:高等教育出版社,1998
[2]中山大学教学力学系概率论与数理统计(下)[MJ.北京:高等教育出版社.1980
[3]华东师范大学数学系概率论与数理统计教程[M]北京:高等教育出版社,1983
责任编校:周伦
JudgementandSolutionMethodsofUMVUE
WANGFen。ling,FANMing‘zhi
(上卸Ⅱ玎舢眦0,肘珊^删n£溉,x∽矗。昭‰洫倦毋,x眦^o昭461000,现iM)
Abstract:Thi8pape。pmvidesmethodsofjudgement0fUMVUE,necessaryandsufficiemconditionofexistence,severals0】u石ons抽UMVUE,ande。p】ajns£b。pr靴出丑l印p】ica石on。f出eseme出ods呐aafewe髓“1ples.
Keywords:UMVUE;Unbiasedestimate;su珩c;entstaIist记;C—R10werbound