理论流行病学课件

1、11、概述、概述2、用途、用途3、建立、建立4、实例、实例4、应用、应用2 定义：定义：用数学公式明确地和定量地表达病因、宿主和环用数学公式明确地和定量地表达病因、宿主和环境之间构成的疾病流行规律，同时从理论上探讨不同防境之间构成的疾病流行规律，同时从理论上探讨不同防治措施的效应。治措施的效应。 发展史：发展史： harmerharmer（19601960）：决定流行过程动态规律为易感者和接触率）：决定流行过程动态规律为易感者和接触率 rossross（19111911）：确定性模型）：确定性模型 reed & frostreed & frost（19281928）：）：ree

2、d-frostreed-frost模型模型 armitage & dollarmitage & doll（1960s)1960s)：肿瘤形成模型（非传染病）：肿瘤形成模型（非传染病） 1970s1970s：计算机的应用，提出多状态、时间序列、时空聚集性模：计算机的应用，提出多状态、时间序列、时空聚集性模型型 1980s1980s：混沌论；协同论；灰色模型：混沌论；协同论；灰色模型3第二节第二节 流行病数学模型的用途流行病数学模型的用途 1、研究疾病分布规律、研究疾病分布规律 （1）研究疾病地区聚集性：二项分布、）研究疾病地区聚集性：二项分布、poisson分布分布 （2）研究疾

3、病时间分布特征）研究疾病时间分布特征 （3）研究疾病人群分布）研究疾病人群分布 2、研究疾病流行过程、研究疾病流行过程空间模型、多等级模型空间模型、多等级模型4提出假设选择模型结构模型参数估计建立模型现场观察资料模型性质不满意模型性质满意拟合度理想模型性质满意拟合度欠佳拟合度比较5（一）模型所适用的条件（一）模型所适用的条件1 1、封闭人群、封闭人群2 2、经空气飞沫传播的急性呼吸道传染病、经空气飞沫传播的急性呼吸道传染病3 3、固定的有效接触率、固定的有效接触率（每个个体在单位时间（每个个体在单位时间内与其他个体相互交往发生有效接触的概率）内与其他个体相互交往发生有效接触的概率）4 4、易感

4、者转归：、易感者转归：有效接触后获感染并传染给有效接触后获感染并传染给其他易感者，病后获免疫力。其他易感者，病后获免疫力。5 5、上述各条件在流行过程中保持不变、上述各条件在流行过程中保持不变6（二）明确流行病学等级及其相互转移流程（二）明确流行病学等级及其相互转移流程s(t) ：第：第t 代易感者；代易感者； s(t+1)：第：第(t+1)代的易感者代的易感者c(t)：第：第t 代的病例及传染者；代的病例及传染者； c(t+1)：第：第(t+1)代的病例代的病例i(t)：第：第t 代的免疫者；代的免疫者； i(t+1)：第：第(t+1)代的免疫者代的免疫者s(t)c(t)i(t)s(t+1)

5、s(t+2)c(t+1)c(t+2)i(t+1)i(t+2)7（三）确定模型的数学表达式及其参数（三）确定模型的数学表达式及其参数1、确定型模型：、确定型模型：初值一经确定，参数不再变动初值一经确定，参数不再变动（1）确定模型：下一代将发生的病例数）确定模型：下一代将发生的病例数c(t+1) 是有效接触率是有效接触率p0与与t 代代病例数和病例数和t 代易感人数的乘积代易感人数的乘积 c(t+1)=p0ct st（2） 确定关键变量确定关键变量(参数参数)：有效接触率：有效接触率p0（有效接触指能引起传染（有效接触指能引起传染的接触）。以的接触）。以q 表示该人群中单位时间内任何两个体不发生有

6、效表示该人群中单位时间内任何两个体不发生有效接触的概率，则接触的概率，则q=1-p。如该人群在。如该人群在t 时间有时间有ct 个病例，则个病例，则qct 为一为一名易感者和个名易感者和个ct 病例不发生感染的概率，而病例不发生感染的概率，而1- qct 则为一名易感者则为一名易感者至少和一例患者发生有效接触的概率。即传染率，而某代出现的至少和一例患者发生有效接触的概率。即传染率，而某代出现的病例数为：病例数为： 某代病例数某代病例数=传染率传染率上代余下的易感者数。由此，上代余下的易感者数。由此，计算公式为：计算公式为：（3）数学表达式：）数学表达式： c(t+1)=st(1-qct) s(

7、t+1)=st-c(t+1) i(t+1)=it+ct 82、随机性模型：、随机性模型：确定性模型仅在大人群、传染者接触人多的情确定性模型仅在大人群、传染者接触人多的情况下可作为随机性模型的估计值，对新一代病例数作出点估计。况下可作为随机性模型的估计值，对新一代病例数作出点估计。其其假定有效接触率不变假定有效接触率不变，这与实际不符。，这与实际不符。如以可变的参数计算新如以可变的参数计算新一代病例数的概率区间，先按给定的初值计算出一代病例数的概率区间，先按给定的初值计算出t+1代不同病例的代不同病例的概率，再选择较大概率的事件作为新出始值计算下一代出现不同概率，再选择较大概率的事件作为新出始值

8、计算下一代出现不同病例数的概率，即为病例数的概率，即为。（1）模型：）模型： 式中：式中：pc(t+1)|st,ct为在为在t 时间内存在时间内存在st 个易感者和个易感者和ct 个病个病例的条件下例的条件下,，t+1时间出现时间出现c(t+1) 个病人的概率。也用个病人的概率。也用r表示，表示， 是是0至至st 之间的任何整数。之间的任何整数。)(1111)()1 ()!( !,|ttcsctrcttttttttqqcscscscp9例：例：在在5个易感者中发生了个易感者中发生了1个病例，假定个病例，假定p=0.2，下一，下一代发生代发生0、1、2、3、4病人数的概率是：病人数的概率是：0例

9、例:1例例:2例例:3例例:4例例:0016. 012 . 01)8 . 0()8 . 01 ()!44( ! 4! 40256. 08 . 02 . 04)8 . 0()8 . 01 ()!34( ! 3! 41536. 08 . 02 . 06)8 . 0()8 . 01 ()!24( ! 2! 44096. 08 . 02 . 01)8 . 0()8 . 01 ()!14( ! 1! 44096. 08 . 011)8 . 0()8 . 01 ()!14( ! 0! 44)44(14113)34(13122)24(1213)14(1114)04(10110（四）参数估计及模型拟合（四）参

10、数估计及模型拟合 根据经验和实际资料假定参数值根据经验和实际资料假定参数值 分别将数个假定参数值代入分别将数个假定参数值代入 比较各假定参数获得的曲线和实际比较各假定参数获得的曲线和实际流行曲线的拟合度，取最佳拟合者流行曲线的拟合度，取最佳拟合者 （最小卡方值法、列线图法、最大（最小卡方值法、列线图法、最大似然法）似然法）11水痘流行期间儿童总数水痘流行期间儿童总数n=196过去患过此次未感染者：过去患过此次未感染者：40无患病史，此次未感染者：无患病史，此次未感染者：60无患病史，此次感染者：无患病史，此次感染者：96全部流行期间：全部流行期间：79天天病例呈代出现病例呈代出现其它条件符合要

11、求其它条件符合要求 代代 例数例数 1 1 2 2 3 14 4 38 5 34121、假定有效接触率、假定有效接触率p： 0.01、0.02、0.03 则则q：0.99、0.98、0.972、利用公式计算各代期望、利用公式计算各代期望病例数：结果见表病例数：结果见表（p=0.02)3、以期望病例数为理论数，、以期望病例数为理论数，与各代实际发病例数比与各代实际发病例数比较：卡方检验较：卡方检验 tto22)(13143个个p值中最小的卡方值为值中最小的卡方值为p=0.02，选此进一步拟合。，选此进一步拟合。154、进一步拟合，结果见表、进一步拟合，结果见表结果：以结果：以p=0.0231的卡

12、方值为最小，但的卡方值为最小，但p0.05，说明模型结构，说明模型结构有问题有问题16三、修正模型后拟合（三、修正模型后拟合（reed-frostreed-frost模型衍化式）模型衍化式）1 1、免疫屏障的作用、免疫屏障的作用（1 1）假定）假定2 2名免疫者保护一名易感者（阈值为名免疫者保护一名易感者（阈值为33%33%），拟合优度颇），拟合优度颇佳佳)1(20)1(cttttqisc17（2 2）假定）假定1 1名免疫者保护名免疫者保护1 1名易感者（阈值为名易感者（阈值为50%50%）)1 (0)1(cttttqisc（3）将公式改写为：）将公式改写为：ttttiscnkc01)(1以

13、以k3.6或或p=k/155=0.023226拟合，则效果颇佳拟合，则效果颇佳182 2、隐性感染作用：、隐性感染作用： 式中：式中：b b为流行过程中隐性感染与显性为流行过程中隐性感染与显性感染的比例常数；感染的比例常数； 为第代为第代t t 已经积已经积累的感染者。累的感染者。)1)(0) 1(cttittqcscbtic0b19四、模型抽象研究采用公式：)1 (0)1(cttttqiscttttiscnkc01)(120（一）有效接触率和隔离的变动对流行过程的影响（一）有效接触率和隔离的变动对流行过程的影响21（二）隔离对流行过程的影响以以p=0.05，每代隔离，每代隔离15新病例为例，则：新病例为例，则：tcsi第一代第一代1例例11499 0第二代第二代24.95例例225474 1第二代第二代25例，隔离例，隔离5例例3303 171 26第三代病例数第三代病例数303.434145 26329第三代第三代303例，隔离例，隔离61例，例，242例仍传播例仍传播第四代第四代144.9例，免疫人数超过易感者数例，免