黃金超
(滁州職業(yè)技術學院 基礎部,安徽 滁州 239000)
指數(shù)分布族參數(shù)的經(jīng)驗Bayes檢驗問題已有很多研究,如文獻[1-5]對其做了不同程度的統(tǒng)計推斷研究,并獲得一些重要有意義的結果。彭家龍等[6]研究了Cox模型參數(shù)的經(jīng)驗Bayes檢驗,在適當?shù)臈l件下,收斂速度的階可任意接近但絕大多數(shù)研究經(jīng)驗Bayes檢驗問題的文獻,都是利用密度函數(shù)的通常核估計來構造核函數(shù)。本文將采用密度函數(shù)的遞歸核估計和EB檢驗函數(shù)的單調(diào)性,重新構造一類廣義指數(shù)分布族參數(shù)的檢驗函數(shù)。在一定的條件下獲得了收斂速度的階可任意接近O(n-1)。改進了文獻[6]的EB檢驗相應結果,推廣了現(xiàn)有文獻的相應結果。
考慮如下廣義指數(shù)分布族[6]:設:
這里 q(x)和 Q(x)為連續(xù)q(x)=Q′(x)>0, Q(x)>0,且參數(shù)空間為
考慮分布族式(1)中參數(shù)θ的如下EB檢驗問題
其中,θ0>0為已知常數(shù)。
對檢驗函數(shù)式(2),設損失函數(shù)為“線性損失”。
設k>0且為常數(shù),D={d0,d1}是行動空間,d0表示接受H0,d1表示否定H0,I[A]表示A的示性函數(shù)。
設G(θ)為θ未知先驗分布,且:
式(4)為隨機化判別函數(shù),則δ(x)的風險函數(shù)為:
此處:
其中:
為r.v.X的邊緣分布,而:
由式(6)至式(8)得α(x)的另一表達式:
其中,f(1)(x)、q(1)(x)分別為f(x)、q(x)一階導數(shù),u(x)
由Cauchy-Schwarz不等式和式(8)可得:
其中,φ(1)(x)為φ(x)一階導數(shù),由式(10)可知φ(x)是單調(diào)遞減且連續(xù)。
本文假定:
在假定式(11)成立下,先驗分布G(θ)是非退化的,由介值定理可知,一定存在點aG∈(0,∞),使得φ(aG)=θ0。又由式(9)可知:
因此,由式(5)可知Bayes判決函數(shù)為:
其Bayes風險為:
在式(13)中,當先驗分布G(θ)已知,且δ(x)=δG(x)時,R(G)可以達到的,但此處G(θ)未知,因此δG(x)無使用價值,于是引入EB方法。
設X1,X2,…,Xn和X是iid樣本,令密度函數(shù)為f(x)如式(7),iid樣本作如下假定:
假定Cs,α表示R1中一族密度函數(shù),其s階導數(shù)存在,的正整數(shù)。
令Kr(x)(r=0,1,…,s-1)是有界的Borel可測,在區(qū)間(0,1)之外為0,且滿足(B):
(B2)Kr(x)在R1上是可微的,且
本文假定先驗分布G(θ)非退化,且屬于下列先驗分布類:
這里A1、A2為給定的常數(shù),0<A1<A2<∞,通常取A1為充分小A2充分大。
記f(0)(x)=f(x),f(r)(x)為f(x)的第r階導數(shù),r=0,1,…,s。定義f(r)(x)的遞歸核估計[7]為:
其中,hn↓0且hn>0,Kr(x)是滿足條件(B)的核函數(shù),此估計具有一種遞歸性質(zhì),即:
定義α(x)的估計量:
由假設(C),結合式(12),定義EB檢驗函數(shù)為:
EB檢驗的構造方法由文獻[8]最先提出來的。
令En表示對 r.v.X1,X2,…,Xn的聯(lián)合分布求均值,則δn(x)的全面Bayes風險為:
假定c0,c1,c2,…表示與n無關的正常數(shù):
引理 1:設fn(r)(x) 由式(14)定義,其中X1,X2,…,Xn為 iid樣本序列,假定條件(A)和(B)成立,hn↓0,當時,則有:
證明:由Cr不等式可知,對r=0,1有:
由式(14)和核函數(shù)的性質(zhì)可知:
由Taylor展開得:
將式(20)代入式(19)可得:
再由f(x)∈Cs,α, 及|Kr(t)| ≤C,X1,X2,…,Xn單調(diào)遞減可知:
由式(22)可得:
故有:
將式(24)和式(25)代入式(18),結論成立。
注1:當 0<λ≤1時可任意接近Rn-R(G)。
定理1:設X1,X2,…,Xn為來自分布族(1)iid樣本序列。R(G)、Rn分別由式(13)和式(17)給出,且假設(A)—(C)成立,當時,有:
這里s≥3的正整數(shù)。
證明:由式(13)和式(17)可知:
由式(9)和式(15),Cr不等式與引理1可知:
當x∈(0,A1)時,由式(12)和式(16)可知,δn(x)=0,δG(x)=0,x∈(A2,∞)時,δn(x)=1,δG(x)=1。 故En(δn(x))-δG(x)=0,因此:
當x∈(A1,aG)時,由式(12)和式(16)可知,δG(x)=0,,利用Markov不等式和式(27)有:
當12<λ<1時,且是 第 2 類 瑕 積 分, 瑕 點 為x=aG,由式(8)、式(10)和比較判別法則:
同理可證:
因此,將式(28)至式(31)代入式(26),得:
注2:文獻[6]在舍入數(shù)據(jù)下構造的經(jīng)驗Bayes檢驗,得到了收斂速度階為,其中 0<λ<1,s≥4,當λ→1,s→∞時,可任意接近。本文利用遞歸核估計和EB檢驗的單調(diào)性,重新建立EB檢驗,得到收斂速度的階為。其中 0<λ<1,s≥3,當λ→1,s→∞ 時可任意接近于O(n-1)。文獻[6]在舍入數(shù)據(jù)下的條件,構造EB檢驗比本文復雜;文獻[6]定理的條件要求s≥4比本文s≥3要強,本文引理的證明比文獻[6]引理1要簡單些。由于本文利用EB檢驗的單調(diào)性,在一定的條件下,獲得了收斂速度的階的結果比文獻[6]快了約1倍,改進文獻[6]中的相應結果。
在式(1)中,設Q(x)=x2則r.v.X為Rayleigh分布,即f(x|θ)=2xθe-θx2I(x>0), 設θ的先驗分布族為:
β和r為常數(shù)且β>0,r>0,有
由式(33)可知f(x)關于x任意階可導函數(shù)且有界,即f(x)∈Cs,α,條件(A)成立,在假定(B)成立,故驗證(C)成立即可。
即假定(C)也成立,故定理1成立。