关于Poisson分布的检验讲解

1、II承诺保证书I1引言1研究背景1研究方法及目的12Poisson分布检验的步骤和基本理论2检验步骤2检验的基本原理33关于Poisson分布检验的三个案例及实际研究7案例分析7对单位时间到来顾客数的实际研究13参考文献18英文摘要19 关于Poisson分布的检验肖秋光摘要：Poisson分布是概率论中的一种重要离散分布，在许多实际问题中都有着广泛应用本文概括了检验样本数据是否服从泊松分布的一般方法，主要是对随机数据进行图像模拟估计和利用假设检验原理对给定的临界值进行估计其中x2检验是众所周知的拟合优度检验，它能适用于任意的备择假设.另外，通过三个例子进行说明，最后用该方法对实测数据进行了分

2、析和检验，并得出了结论.关键词：Poisson分布假设检验独立变量x2统计量1引言1.1研究背景改革开放三十年来随着社会的发展、经济的增长，科学技术日新月异、人民拥有的物质日益丰富、感受到的文化也更加多元、社会的各种法规制度日臻成熟，无论是住房、保险、交通、旅游、高质量产品还是教育、饮食等.其结果是构成了大量的随机数据，而这些数据有没有什么规律可循呢？就需要我们对它进行研究.在现实生活中的许多数据经过人们大量的研究是服从泊松分布的.若通过观察记录得到了一组数据，它是否服从泊松分布，则需要我们对其进行检验.泊松分布是1837年由法国数学家泊松(PoissonS.D.17811840)首次提出的.

3、它是概率论中的一种重要的离散型随机变量的概率分布，在理论上和实践中都有广泛的应用.如110报警台24小时接到的报警次数、一定时间内发生的意外事件次数或灾害次数、布匹上的疵点数目、放射性物质放射出的粒子数目等.1.2研究方法及目的由于向110报警台的报警是一次次到来的；自然灾害是一次次发生的；放射性粒子是一个个射出的；进入商场的人是一个个到来的它们都可以看成是一种于随机时刻到来的“质点流”.要对其进行研究，首先，必须收集到有效的数据.其次，由于得到的样本数据通常是实验或统计而来，因此它不能完全的反映事物的本质.我们主要对部分数据进行抽取分析，根据部分数据对全体数据做出推断及判断.因此，研究单位时

4、间内产生的诸多随机变量有助于当事者们对各种新措施、新技术作出更为科学合理的决策.例如，商场每个时段到达的人数不一，通过调查可以确定哪个时段是人流的高峰期，可以在这个时段做一些宣传或促销产生的效益就会比其他时段高，并有效控制成本，使其用最小的投入换来最大的收益.2Poisson分布检验的步骤及基本理论2.1检验步骤2.1.1数据整理进行Poisson分布的检验时，首先要对收集到的数据进行整理.假设收集到单位时间的量为x,x,xx,然后把这些量按从小到大顺序排列起来，并查123n出其频数稍加整理制成表格如下：表1单位时间的量xi012xi频数pip0p1p2pi其中满足：x+x+x=0 xp+1x

5、p+xxp12n01ii用图像对样本数据进行模拟由于图形比较直观，而且样本数据在一定程度上能有效反映总体的分布规律，故可以用样本数据的图像模拟通过对比，对该分布进行初步判断.泊松分布的图形一般为左偏，但随九数值的增大，图形趋于对称.Fp值，则拒绝原假设H；0若a0的泊松分布，X,X,X，,X为其样本.123n考虑检验问题：H:九=九；H:九H九,现有0i=1p（x;i）=liei=e-入九刀%x!i=1i1n（x!）ii=101x）ln九i1n（x!）ii=1其中T（x,x,，x）=工x,b（九）=In九12nii=11h（x,x，x）=,c（）=en九12nn打（x!）ii=11,Tc12b

6、,T=c,j=1,2ij0,cT+b（）+b（-）0en0+（n九）qen0+（n九）c2en0j!c!0c!012-叫+j=c2+1艺加)丿j!nM=）（叫je叫+0j=0j!j=c2+1bcbcen0+-（n九en0+22（n九）c2en0c!0c!012在一般情况下上述方程不易求解，但当九不接近于零而n又不很小时，统计量0X一n九i0U=4的渐进分布为正态分布N(0,1)，则v；nX0P工Xn九u叽.=10Jo,=1i0b,b,c,c1212对一切实数u都渐近地成立(这是因为正态分布具有对称性)因此,由下式确定：a另（n九）jb（7）（n九）j2j!0c!0八0j!j=011j=c2+1

7、+b2（nX）cc02e一叫)3关于Poisson分布检验的三个案例及实际研究案例分析论反腐败与泊松分布腐败现象作为当今社会的一种非常态，它的发生、出现引起了广大群众的关注.调查显示最近几年科级腐败正在加剧，小官受贿成隐患.据悉，某检察院工作人员对某经济较落后省的320个底层官员在一年时间内的受贿金额调查纪录如下表所示.根据这些数据（金额0表示未受贿，金额1表示受贿金额大于0小于等于1，其余类同）检验受贿金额是否服从泊送分布.表21年内320个官员受贿金额（万元）统计表金额012345678910合计人数154770815225169410320来源于参考文献6用折线图像模拟数据如下：官员受贿

8、频数图系列11234567891011受贿金额图2从图形走势看，为左偏凸值分布，与泊松分布较为相似，可初步判定为泊松分布.在理论上，这里我们需要检验的是在一年的时间段内受贿官员的受贿金额是否服从泊送分布，所以可以假设H0：一年的时间内受贿官员的受贿金额服从泊送分布；H1:一年的时间内受贿官员的受贿金额不服从泊送分布x!我们知道泊送分布的概率密度函数为f（X=x）xZ,式中：:是未知参数.如果假设为真时，可以根据本数据估计九.由上表的数据可以的到在320个底层官员中，平均每一官员受贿的金额(万元)，即八0 x15+lx47+9x1+10 x0九=3.0320因此，可以用X作为九的估计值，即得到为

9、真时的概率密度函数f(Xx!根据该密度函数，就可以计算出在每一个官员的受贿金额为各个类别出现的概率，这些概率值可通过泊送分布表查得.例如，在一年内受贿金额为0万元的官员人数的概率是/(X=0)=0.498，受贿金额为1万元的概率是/(X=1)=0.1494等然后用查出的概率分别乘以样本容量n(n=320)，就可以得到各类别期望的频数例如，在320个官员中受贿金额为0万元的期望频数是0.0498x320二15.936下表列出了X2统计量的计算过程.表3X2统计量的计算过程受贿金额XI为真时的/(X二x)实际频数I期望频数e=n-/(X=x)II(ne)2IIeI00.04981515.9360.

10、055010.14944747.8080.013720.22407071.680.039430.22408171.681.211840.16805253.760.057650.10082532.2561.632260.05041616.1281.015970.021696.91280.008142.5920.981290.002710.86410万元以上0.001200.384合计1.00003203205.0068我们注意到表中，受贿金额为8，9和10万元次及以上金额的期望频数都小于5，所以将这三类归于受贿金额为7万元的合并为一类，所以合并之后的类别数k二8.这时咒2统计量为x2-丈(nie

11、i)2-5.0068ei-1i需要注意的是：根据Pearson定理，上式的咒2统计量服从自由度为k-r-1的咒2分布，其中k时类别的个数，r是估计的总体参数的个数在这里k=&r=1(只估计了一个参数九)，所以自由度为k-r-1二8-1-1二6.于是，当0.05时，查表可得x2(6)二12.592.对于样本的X2值，因为咒211ni572033835255324082731394527106来源于参考文献7在上表中的n是观测到i个粒子的时间间隔数（最后一项已经合并）若要i检验观测的数据服从泊松分布这一假设（“0.05），贝V：因为对参数为九的泊松分布是：P（X=k）=e-X,k=0,1,2k!根

12、据上表原始数据可以算得X最大似然估计X二X二3.870而X二3.870的泊松分布通过计算机计算及查表（泊松分布函数表）可得下表：表5inipinp.(nnp)2iinpi0570.20954.50.114712030.0807210.50267223830.1562407.41.461435250.2015525.50.000545320.1950508.61.076654080.1509393.50.534362730.0973253.81.452571390.0538140.30.01208450.026067.87.66739270.011229.20.165810100.005311.

13、20.12861160.00225.70.0158总和26081.0000260812.8967因此x2二12.8967,其自由度为12-1-1=10，对0.05查（咒2分布分位数X2（n）表）得X2（10）二18.307,所以我们接受H,认为观测数据服从泊松分布.p0.050卢瑟福等人经过两年时间综合多方面因素的分析，在1911年提出原子的核式模型，原子中的正电荷集中在原子中心很小的区域内，而且原子的全部质量也集中在这个区域内.原子核的半径近似为10，约为原子半径的千万分之一.卢瑟福散射实验确立了原子的核式结构，为现代物理的发展奠定了基石.这充分表明研究泊松分布具有重大意义.对印刷错误个数的

14、检验一个检验员检查了一本书的100页，并仔细记录各页中印刷错误的个数，其结果为：表6错误个数/i01234567含/个错误的页数i36401920210其频数模拟如下图：图4若要检验一页的印刷错误个数是否服从泊松分布.（取0.05）则:假设H:总体X服从泊松分布；H:总体X不服从泊松分布0136x0+40 x1+19x2+2x3+0 x4+2x5+lx6+0 x7从表中数据可得：x=1100当H成立时，九的最大似然估计为九=x=1，检验的拒绝域为：0f2X2=EL-nx2(k-r-1)npi由给出的条件可知n二100二p(X=0)=卷=03679，=p(X=1)=0.367912e-1二p(x

15、=2)=-2T=0J8397，13e-1二p(x=3)=-3T=0.0613p=p(X=4)=414e-1-4厂=0.01533，15e-1=p(X=5)=5-=0.00306616e-1p=p(X=6)=0.000511,66!p二p(X7)二1-fp二0.0000237ii=0而对于j3,有np.5，j将其合并得np=8.023，合并后k=4，查x2分布jj=3分位数x2(n)表可得：X2(4-1-1)=5.991p0.05x2362402+36.7936.79+52-100=1.44418.3978.023由1.44410.046211A203.23411.3826.525220.106

16、34-7.44130.1631611.4172.570240.18751013.1250.744050.17251012.0750.356560.132339.2614.232870.0869S6.083、80.050063.50090.025551.785100.011810.826110.00491(210.3436.52532.1112120.001910.133130.000710.049大于130.00030丿0.021丿合计1707046.5399从表中观察进入人数为0，1，8及其以上的期望频数都小于5，所以将进入人数为0，1的与进入2个的合并为一类，进入8个及以上与进入7的合并为

17、一类，那么合并后的类别数k二6，其X2统计量为：X2仝二=46.5399e1i根据Pearson定理，X2统计量服从自由度为k-r-1的x2分布，在这里k=6,r=1,所以自由度为4于是，当0.05时，查X2分布表得X2（4）=9.4880.05对于样本的X2值，有X2X2（4）落在拒绝域内所以拒绝H，接受H，即每分0.0501钟进入书城的人数不服从泊松分布，与用图像模拟得到的初步结论一致.3.4.2当取九二4.7时，有f（X=x）=”4，查泊松分布表计算得下表表9进入人数xi为真时的f（X二x）i实际频数ni期望频数(n-e)2e=nif(X二x)iei00.00915j0.637、10.0

18、42711L202.989”10.6618.180920.10054J7.03530.1574611.0182.885440.18491012.9430.669250.17381012.1660.385660.136239.5344.478070.09146、6.39880.053763.75990.028051.960100.01321l0.924I110.00561210.392f13.6643.9386120.002210.154130.000810.056大于130.00030丿0.021,合计1.00007069.98620.5377其余同上，当二0.05时，查咒2分布表得咒2（4）

19、=9.488对于样本的咒2值，仍然0.05有咒2咒2（4）落在拒绝域内，结论也一致.0.05按照众多学者研究，在一般情况下对于这种随机变量的检验应该服从Poisson分布.在这里为什么得到了相反的结论？我认为可能有两点原因：一是所采集到的样本数据可能具有特殊性，不能完全反映总体的分布，若多测几组数据进行检验结果可能会改变；二是我们假设的前提是每个人相互独立，而实际情况有许多人是有关联的，例如一个家庭四个人同时进入、两个好朋友陪同一个朋友买书等诸多情况，这个因素可能也会对结果有影响.结束语：随着社会的发展和可研究案例的增多，把自然科学领域的原理用到社会现象的研究、分析上来，对发现问题解决问题有许

20、多益处.参考文献：茆诗松等概率论与数理统计教程M.北京：高等教育出版社,2008.朱洪文应用统计M.北京：高等教育出版社,2006.王荣华等.概率论与数理统计（习题精选）.北京：北京大学出版社,2010.谢民育，吴茗，熊明.多元分布下单边备择假设的两步检验.武汉2009徐亮，丁先文等.基于经验似然的部分线性模型的统计诊断.南京，2009 HYPERLINK .en/g/20110401/10129630369.shtml .en/g/20110401/10129630369.shtmlhttp：/ HYPERLINK http:/www.phys www.phys./advlab/nuclear/10.htmPOISSONDISTRIBUTIONONTHEINSPECTIONXIAOQi