王　蘭,　谷偉偉,　蘆凌飛

(中國(guó)礦業(yè)大學(xué) 理學(xué)院, 江蘇 徐州 221116)

?

加權(quán)p,q對(duì)稱熵?fù)p失下Poisson分布變異系數(shù)的Bayes估計(jì)

王蘭,谷偉偉,蘆凌飛

(中國(guó)礦業(yè)大學(xué) 理學(xué)院, 江蘇 徐州 221116)

為更全面研究Poisson分布的估計(jì)問(wèn)題,在加權(quán)p,q對(duì)稱熵?fù)p失函數(shù)下,討論了Poisson分布變異系數(shù)的Bayes估計(jì),并給出了Bayes估計(jì)的置信區(qū)間和多層Bayes估計(jì),得到了其具體形式。

Poisson分布; 加權(quán)p,q對(duì)稱熵?fù)p失函數(shù); Bayes估計(jì)

0　引　言

Poisson分布是1837年由法國(guó)數(shù)學(xué)家Poisson S.D.(1781-1840)首次提出的一種常用的離散分布,常與單位時(shí)間、單位產(chǎn)品或單位面積等上的計(jì)數(shù)過(guò)程相聯(lián)系,相關(guān)研究已有很多[1-4]。Poisson分布在現(xiàn)實(shí)生產(chǎn)和生活中應(yīng)用非常廣泛,是一類重要的離散分布。 其分布律為

(1)

由于隨機(jī)變量的取值有量綱,不同量綱的隨機(jī)變量用方差或標(biāo)準(zhǔn)差去比較其波動(dòng)大小不太合理,并且在取值的量綱相同的情形下,取值的大小有一個(gè)相對(duì)性問(wèn)題,取值較大的隨機(jī)變量的方差或標(biāo)準(zhǔn)差也允許大一些,所以在比較兩個(gè)隨機(jī)變量的波動(dòng)大小時(shí)常用到變異系數(shù),其具體定義如下:

定義1設(shè)隨機(jī)變量X的二階矩存在,則稱比值

為X的變異系數(shù)。

由定義1知變異系數(shù)是以其數(shù)學(xué)期望為單位去度量隨機(jī)變量取值波動(dòng)程度的特征數(shù),標(biāo)準(zhǔn)差的量綱與數(shù)學(xué)期望的量綱是一致的,所以變異系數(shù)是一個(gè)無(wú)量綱的量,是刻畫(huà)保險(xiǎn)風(fēng)險(xiǎn)最常用的指標(biāo)之一。

其中x1,x2,…,xn為隨機(jī)樣本X1,X2,…,Xn的一組實(shí)現(xiàn)值。

在統(tǒng)計(jì)決策問(wèn)題中,統(tǒng)計(jì)決策及參數(shù)估計(jì)的優(yōu)劣性在很大程度上依賴于損失函數(shù)形式的選取。常用的損失函數(shù)有平方損失函數(shù)、對(duì)稱損失函數(shù)[5]、熵?fù)p失函數(shù)[6]和LINEX損失函數(shù)[7]等。文獻(xiàn)[8-9]等已利用對(duì)稱損失函數(shù)得到相關(guān)研究的結(jié)果。文中在加權(quán)p,q對(duì)稱熵?fù)p失函數(shù)

(2)

下討論變異系數(shù)θ的估計(jì),可根據(jù)p,q的靈活選取做出最合適的估計(jì)。這里δ是待估參數(shù)θ的估計(jì)量。由損失函數(shù)的性質(zhì)易知其關(guān)于δ是嚴(yán)格凸函數(shù),并且在δ=θ處取得唯一的最小值。

1　變異系數(shù)θ的Bayes估計(jì)

下面定理給出了在損失函數(shù)(2)下變異系數(shù)θ的貝葉斯估計(jì)。

定理1記X=(X1,X2,…,Xn),在損失函數(shù)(2)下,無(wú)論先驗(yàn)分布取何種分布,變異系數(shù)θ的Bayes估計(jì)為

R(θ,δ)=E[L(θ,δ)]=

(3)

要使式(3)左端R(θ,δ)最小,則需要極小化

由于

(4)

式(4)右端關(guān)于δ求導(dǎo)并令其等于0,解得

下面定理2是在給定先驗(yàn)分布π(θ)后,結(jié)合定理1,給出了變異系數(shù)θ的精確Bayes估計(jì)。

定理2對(duì)于給定損失函數(shù)(2),若參數(shù)λ的先驗(yàn)分布為

(5)

則變異系數(shù)θ的Bayes估計(jì)為

證明易得λ的后驗(yàn)密度

h(λ|X)∝λT(X)+α-1·e-λ(n+β),

從而有

同理有

從而由定理1得

2　變異系數(shù)θ的置信區(qū)間

定理3對(duì)于Poisson分布,給定加權(quán)p,q對(duì)稱熵?fù)p失函數(shù)(2)和先驗(yàn)分布(5),對(duì)給定置信水平1-k,變異系數(shù)θ的置信區(qū)間為

由定理2,

h(λ|X)∝λT(X)+α-1·e-λ(n+β)，α,β>0，

由于2(n+β)λ～χ2(2T(X)+2α),故對(duì)給定k∈(0,1),有

1-k。

從而參數(shù)λ的置信水平1-k的置信區(qū)間為

3　變異系數(shù)θ的多層Bayes估計(jì)

定理2中的變異系數(shù)θ的Bayes估計(jì)δB(X)中的α、β仍未知,這樣的參數(shù)稱為超參數(shù)。對(duì)兩個(gè)未知參數(shù)α、β分別給定先驗(yàn)得到超先驗(yàn),這里的超先驗(yàn)和定理2中的先驗(yàn)決定了一個(gè)新的多層先驗(yàn)。θ的多層Bayes估計(jì)就是在該多層先驗(yàn)下得到。

先驗(yàn)分布(5)中的超參數(shù)α>0,β>0未知并相互獨(dú)立,利用減函數(shù)構(gòu)造法得到兩個(gè)超先驗(yàn)分別為:

h1(α)=U(0,1) ,h2(β)=U(0,c)，

(6)

其中c為常數(shù)。對(duì)于這樣的兩層先驗(yàn)分布,下面定理4給出了θ的兩層Bayes估計(jì)。

定理4對(duì)于Poisson分布(1),給定式(6)中的兩層先驗(yàn)分布和損失函數(shù)(2), 變異系數(shù)θ的兩層Bayes估計(jì)為

證明由上面討論可知變異系數(shù)θ的先驗(yàn)分布為

從而得到θ的后驗(yàn)分布為

進(jìn)而得

所以得

[1]VELLAISAMY P, SUSHMITA J. Estimating the parameter of the population selected from discrete exponential family[J]. Statistics & Probability Letters, 2008, 78(9): 1076-1087.

[2]STAMEY J D, BRATCHER T L, YOUNG D M. Parameter subset selection and multiple comparisons of Poisson rate parameters with misclassification[J]. Computational Statistics and Data Analysis, 2004, 45(3): 467-479.

[3]劉素蓉, 任海平. 加權(quán)平衡損失下泊松分布參數(shù)的Bayes估計(jì)[J]. 統(tǒng)計(jì)與決策, 2010(16): 17-19.

[4]徐凌云, 朱寧, 方愛(ài)秋, 等. 定數(shù)截尾情形指數(shù)-泊松分布參數(shù)的Bayes估計(jì)[J]. 統(tǒng)計(jì)與決策, 2010(15): 31-32.

[5]許俊美, 宋立新. 一種對(duì)稱損失下Pareto分布參數(shù)的Bayes估計(jì)[J]. 白城師范學(xué)院學(xué)報(bào), 2007, 6(3): 5-7.

[6]任海平, 宋允全. 兩個(gè)不同損失函數(shù)下一類分布族參數(shù)的Minimax估計(jì)[J]. 魯東大學(xué)學(xué)報(bào): 自然科學(xué)版, 2009, 25(3): 201-205.

[7]崔群法, 張明亮, 康會(huì)光. LINEX損失下參數(shù)未知時(shí)指數(shù)-威布爾分布的Bayes估計(jì)[J]. 河南大學(xué)學(xué)報(bào): 自然科學(xué)版, 2011, 41(4): 348-351, 393.

[8]徐寶. 壽命產(chǎn)品可靠度的貝葉斯估計(jì)[J]. 統(tǒng)計(jì)與決策, 2011(4): 153-154.

[9]王蘭. 加權(quán)p,q對(duì)稱熵?fù)p失下Reyleigh分布參數(shù)倒數(shù)的Bayes估計(jì)[J]. 貴州師范學(xué)院學(xué)報(bào), 2012, 28(3): 1-3.

[10]茆詩(shī)松, 王靜龍, 濮曉龍. 高等數(shù)理統(tǒng)計(jì)[M]. 北京: 高等教育出版社, 2006.

[11]茆詩(shī)松. 貝葉斯統(tǒng)計(jì)[M]. 北京: 中國(guó)統(tǒng)計(jì)出版社, 1999.

[12]韋博成. 參數(shù)統(tǒng)計(jì)教程[M]. 北京: 高等教育出版社, 2006.

[13]茆詩(shī)松, 程依明, 濮曉龍. 概率論與數(shù)理統(tǒng)計(jì)教程[M]. 北京: 高等教育出版社, 2004.

(編輯王冬)

Bayes estimation of Poisson coefficient of variance under weightedp,qsymmetric loss function

WANGLan,GUWeiwei,LULingfei

(College of Sciences, China University of Mining & Technology, Xuzhou 221116, China)

Aimed at a better study of the estimation of Poisson distribution, this paper discusses the Bayes estimation of Poisson coefficient of variance under the weightedp,qsymmetric loss function, gives the expression of its confidence interval and hierarchical Bayes estimation ,and offers its concrete form.

Poisson distribution; weightedp,qsymmetric loss function; Bayes estimation

2013-04-27

王蘭(1988-),女，江蘇省徐州人，碩士，研究方向：概率論與數(shù)理統(tǒng)計(jì)，E-mail:Wanglan4095@126.com。

10.3969/j.issn.1671-0118.2013.04.010

O212.8

1671-0118(2013)04-0356-04

A