第6章 概率统计方法模型上.doc

第6章 在各种传染病的流行过程中，无论健康人还是病人，任何两个人之间接触的机会都是随机的，而且当健康人与病人接触时，健康人是否被传染也是一个随机的事件。我们通过建立传染病随机模型来分析这些随机规律。 假设人群总的规模为N，在总人群中，病人的数量为m，健康人的数量为s，即满足N=m+s。在人们的日常生活中，任意两人之间（包括健康人和病人）接触的概率相同，每人平均与k个人接触。当健康人和病人接触时，被传染的概率为p。在以上假设的参数中，m和s通常是已知的，k和p可以通过专家的经验和统计数据获得。我们分析的目的是寻找健康人群中每天平均被感染的人数与已知参数之间的关系，以及初始参数对传染病的扩散速度和流行趋势的影响。 我们首先以每一名健康人为研究对象，探讨其每天被感染的概率，而每一名健康人被一名指定病人接触并传染的概率等于每名健康人与指定传染者接触的概率乘以接触时感染的概率。记人群中任意两人接触的概率为q，则对每一名健康人来说，其每天接触的人数服从二项分布，分布函数为 , （6.1.1） 这个分布的期望为k，即，进而。这样，一名健康人被一名指定病人接触并感染的概率为 . 进一步，对人群中的每一名健康人来说，其每天不被感染的概率为，被感染的概率为 . （6.1.2） 所以，对人群中的所有健康人来说，每天被感染的人数服从二项分布，分布函数为 , （6.1.3） 每天被感染的人数期望为，标准差为 为了得到简明的结果，对进行近似计算，由于通常人群的总数，且根据Talyor展开，得 ， 因此， . （6.1.4） 通过式（6.1.4）可以看出平均每天被感染的人数与s、m、p和k之间的关系。进而可以度量平均每天被感染人数的相对误差即 （6.1.5） 由式（6.1.4）可以看出，对于健康人群来说，每天平均被感染的人数与人群中每人每天平均接触的人数k，健康人与病人接触时被感染的概率p成正比。当n，p，k都确定的情况下，时，也就是在整个人群中，病人和健康人的数量各占一半时，每天被感染的人数达到最大。 为了对传染病的传染过程有一个直观的了解，假设一个人口总量n=10000的人群，在日常生活中，平均每人每天接触的人数k=18，健康人与病人接触时被感染的概率p=10%，对于不同的m，平均每天被感染人数与相对误差的变化趋势如图6.1.1和图6.1.2所示。可见被感染人数随着病人数量的增大而增大，直到病人数量占总人群数量的一半时达到最大，随后呈下降趋势。随着病人人口的增加每天被感染人数的相对误差一直呈减少趋势，尤为明显的是病人数量增长的前期，相对误差急剧减少。 图6.1.1 平均每天被感染人数的趋势图 图6.1.2 平均每天被感染人数的相对误差趋势 R编程如下： crb - function(m, n=10000, p=0.1, k=18) #函数 { u-(m*(n-m)*p*k)/(n-1);u } m-1:10000 plot(1:10000,crb(m), xlab=m, ylab=平均每天被的传染人数,type=l, col=blue) crb1 - function(s, n=10000, p=0.1, k=18) #相对误差函数 { miugama-((n-1-m*p*k)/((n-m)*m*p*k))^0.5;miugama } m-1:6000 plot(1:6000,crb1(m), xlab=m, ylab=相对误差,type=l, col=red) （2）企鹅繁殖模型 企鹅的繁殖过程是一个典型的随机不确定模型。首先，每只母企鹅下蛋的数量是随机的，服从泊松分布，其次，每个企鹅蛋是否可以成功孵化也是不确定的。针对这一问题，我们在合理假设的基础上，建立概率模型，求企鹅后代个数的期望值。 根据人们的统计，企鹅生蛋的个数是服从参数为的泊松分布，即 （6.1.6） 而每个生蛋能发育成企鹅的概率为p，且每个生蛋能否发育成企鹅是彼此独立的随机事件。令代表企鹅后代的个数，且有。可见取非负的整数值0, 1, 2, …，对于的概率，我们利用全概率公式 (6.1.7) 注意到每个生蛋发育成小企鹅是相互独立的，且发育成小企鹅的概率为p，因此，实际上反映了有k个生蛋，每个生蛋独立发育，恰好发育成 k 个企鹅的概率。显然，它是一个伯努利试验，因而