负二项分布的性质特征及在流行病学研究中的应用

1、负二项分布的性质特征及在流行病学研究中的应用         09-08-27 15:21:00     作者：韩新焕    编辑：studa20【摘要】  给出了负二项分布的分解定理，进一步研究了负二项分布的有关性质及参数 的无偏一致估计，以及在流行病学该分布的生物学意义。 【关键词】  负二项分布； 无偏一致估计； 应用负二项分布是概率论中常用的重要的离散型随机分布，它在医学中主要用于聚集性疾病及生物、微生物、寄

2、生虫分布模型等的研究。具体地说，当个体间发病概率不相等可以拟合负二项分布，如单位人数内某传染病的发病人数，某地方病、遗传病的发病人数等，这些均可通过负二项分布进行处理。本文从概率论的角度阐述负二项分布的性质及参数 的最小方差无偏估计，并且以该分布在流行病学中应用为例证讨论了其生物学意义。    1  负二项分布的概率模型    负二项分布又称帕斯卡分布（Pascal）,它有两种基本模型1：    模型：假定每次试验可能的结果只有两个：可归结为成功或失败，每次试验之间是独立，每次成功的概率均为 ，

3、直到恰好出现r（指定的一个自然数）次成功所需试验次数X,则X的概率分布为：    p(X=K)=Cr-1k-1k-1(1-)k-r=Cr-1k-1-(1-)k-r    k=r,r+1(1)    模型：假定每次试验可能的结果只有两个：可归结为成功或失败，每次试验之间是独立，每次成功的概率均为 ，试验进行到r次成功为止，记X为试验共进行的次数,则X的概率分布为3：    p(X=k)=Cr-1k+r-1k(1-)k    k=0,1,2,(2)&

4、#160;   此分布的概率是r(1-(1-)-r  的幂级数展开式的项，负二项分布由此而得名    记作 Xf(k,r,) ， 或 XNB(r,)    一个重要的特例是   r=1。   这时（2）成为    p(X=k)=(1-)k    k=0,1,2,(3)    称为几何分布。    2  性质特征  &#

5、160; 为研究负二项分布的性质，我们先给出一个重要的结论：    引理： 设XNB(r,)，则其特征函数为x(t)=r(1-(1-)eit)-r    证明： x(t)=E（eitx)=i=0Cr-1i+r-1r(1-)i eitr    =i=0Cr-1i+r-1r(1-) e)rti    =ri=0Cr-1i+r-1(1-) ert)i    =r(1-(1-)eit)-r    定理1  设: X

6、1,X2,Xr（3）的iid样本，如果    X=ri=1Xi, 则X=ri=1XiNB(r,)    证明：因为X1，X2，Xr独立同分布，又有引理知 X=ri=1Xi的特征函数为：    (t)=r(1-(1-) eit)-r    =rk=0(-r)(-r01)(-r-k+1)k! (1-) eit)k(-1)keitr    =rk=0(r+k-1)!(r-1)!k! (1-)k eit(k+1)    =k

7、=0r(1-)k eit(k+r) Cr-1r+k-1    这正是 p(X=k)=Cr-1r+k-1(1-)k 的概率分布    则 X=ri=1XiNB(r,)    定理2     设：X=X1,X2,Xn)是（1）的iid样本，则T（X）=ni=0XiNB(nr,),则有    p(T=k)=Cnr-1k-1nr(1-)k-nr    k=nr,nr+1,(4)   

8、证明：   设 的特征函数为f(t) ，那么    f(t)=x=reitxCr-1N-1N(1-)N-r =eit1-(1-)eitr    因为x是 的iid样本，所以Xi 的特征函数fi(t)=f(t),i=1,2,n    有特征函数的性质得T的特征函数为:    ni=1fi(t)eit1-(1-)eitr    由于特征函数与概率分布唯一对应，所以Tf(k,nr,) ，其概率分布便是（4）。  

9、  定理3  设：X=(X1，X2，Xn)是（1）的iid样本，则    T（X）=nr-1ni=1Xi-1， 则它是 的最小方差无偏估计。    证明： 由定理2可知    E(T(X)=k=nrnr-1k-1Cnr-1k-1nr(1-)k-nr    =k-1=nr-1 C(nr-1)-1(K-1)-1 nr-1×(1-)(k-1)-(nr-1)    =    所以T(X)是 的无

10、偏估计。     09-08-27 15:21:00     作者：韩新焕    编辑：studa20    又由于E(T(X)= ,有切贝晓夫不等式，对>0, 有    p(|T(X)-|)V(T(X)2      而    V(T(X)=k=nrnr-1k-12Cnr-1k-1nr(1-)k-nr  

11、0; =2 k=nrnr-1k-1×k-2nr-2-1×C(nr-2)-1(k-2)-1nr-2(1-)k-nr    =2 k=nr1(k-1) (k-nr)(nr-2)×C(nr-2)-1(k-2)-1nr-2(1-)k-nr < 2nr-2  k=nr C(nr-2)-1(k-2)-1nr-2(1-)(k-2)-(nr-2)    =2nr-2    所以，对>0, 都有linnp(|T(X)-|)=0 ，可见T(X)是 的一致估计。 &

12、#160;  又因为E(T(X)= ,根据 Lehmannscheff定理， 的最小方差无偏估计必存在，而T(X)=nr-1 ni=1Xi-1,只依赖T（X）= ni=1Xi , 即 T(X)=nr-1 ni=1Xi-1 是 的一致最小方差无偏估计。    3  负二项分布的最可能数和概率的最大值    如果Xk=p(X=k)=Cr-1k-1rqk-r  k=r+1,r+2,其中q=1-,则  当r-q1-q不为整数时，k0=r-q1-q时为负二项分布的唯一最可能的数，即 k=k0时，p(X=k

13、0)达到最大值。    证明：  pkpk-1=Cr-1k-1rqk-rCr-1k-2rqk-1-r=q(k-1)k-r    =>1时，kpk-1,随k增大，概率增大    <1时，k>r-q1-q, pk<PK-1,随K增大，概率减少< p>     =1时，k=r-q1-q, pk=pk-1      当r-q1-q不为整数时，则存在唯一k0=r-q1-q  满足 

14、 r-q1-q-1< k0 <R-Q1-Q< p>     使  pk0-1< pk0>pk0+1，从而k0=r-q1-q是唯一最可能的数，即 k=k0时，p(X=k0)达到最大值。    4  流行病中实例分析    假设血吸虫成虫随机地分布于人群中，即所有的个体均有同等的机会获得新感染，将致成虫在人群中呈Poisson分布。然而，由于暴露的危险性不等、易感性不一致及可能存在的获得性免疫等将导致感染的机会不等，而出现成虫集中在某一部分的人群中。一些

15、可以直接通过驱虫获得人群虫负荷分布的资料及某些尸检资料均提示蠕虫的成虫(如蛔虫、钩虫、曼氏血吸虫、鞭虫)在人群中的分布具有聚集块。聚集块内病例个体的平均数又服从 分布。即：病例的数目H服从均数为 的poisson分布,由于 是变化的，假定其概率分布可用分布表示。于是对于给定的 ，条件概率为：    p(H=h|)=he-h!,   h=0,1,2,>0    此时 的概率密度函数为：    f()=() -1e-,>0,其中>0,>0 都是参数。当和 变化时，可

16、产生一族分布曲线。可以证明，Poisson 分布就是负二项分布。    因为  ，ex= k=0 xkk!    H是离散型随机变量， 是连续型随机变量，H的边际概率函数是：    p(H=h)=JF(Z0e-hh! ()-1ed JF)    =h!() JF(Z0+h-1e(1+)d JF)    =h!(k-1)! (+h-1)!(+1)h+    =Chh+-1 (+1) (1-+1)h 

17、   =Chh+k-1 k (1-)h, (j=0,1,2)    其中，0<<+1<1, k=,   此即为负二项分布。    5  讨论    负二项分布是当poisson中参数服从 分布时所得的复合分布，分布中的参数 是不定的变化的，且其变化是有规律的。呈现的特点是病例聚集群内病例个体的密度服从 分布，病例个体间的流行病学联系与 分布有关，由于 分布的概率密度函数及图形为：    f(x)=() x-1e-x,  x>0    0,        x0    因此，负二项分布来源于poisson分布，它改进了poisson的等概条