背景与目的麻疹是由麻疹病毒引起的急性感染性疾病,其传染源为处在传染期的麻疹病人,主要传播途径为飞沫传播,是严重威胁公共卫生安全的重要传染性疾病。中国向世界卫生组织承诺2012年消除麻疹,但由于各种原因,全国麻疹仍然有升高趋势。山东省人口密度较高,流动性大,因此根据麻疹发病特征,探索近年来麻疹疫情特征仍具有重要的公共卫生意义。常见的流行病学研究方法包含观察性研究和实验性研究,而理论流行病学方法研究较少。常规的常微分方程传染病动力学模型,具有建模简单、计算较快的特点,但该方法估计的参数对应马尔科夫链的均值,未考虑随机因素对发病数的影响,同时常微分方程模型输出曲线平滑,对发病数据较少的地市并不适用。卡尔曼滤波针对实际发生病例数和真实发病数之间的区别加入随机因素,是一种经典的研究方法,但要求发病数据、噪声和参数满足正态性假设,因此该方法具有一定局限性。状态空间模型(Space State Model)则假设发病数据、噪声和参数满足任意分布,是理论流行病学研究中的新方法。近年来缺少针对疫苗保护率进行的专项研究,缺乏疫苗保护率等数据,因此针对重要的不确定参数应当进行多种假设。理论流行病学以基础再生数(R0)来衡量麻疹的传染性,该参数是单个病例平均能够传染的二代病例数,该值的大小说明麻疹传染性的高低。本研究利用基础再生数表明不同地市麻疹暴发的风险,并为山东省麻疹消除计算所需要的疫苗保护率,旨在为山东省麻疹消除的政策制定提供科学指导。资料与方法本研究资料来源于传染病网络直报系统,提取了山东省2014～2015年麻疹个案数据(包含发病时间和地点等信息)。并根据麻疹发病时间和地点,建立合适的状态空间模型,利用迭代滤波(Iterated Filter)对参数进行估计,置信区间采用轮廓似然函数估计。根据国家麻疹消除标准,麻疹消除目标定义为发病率低于百万分之一。通过统计模拟估计在不同疫苗保护率假设下,麻疹消除的概率。结果1.山东省2014～2015年麻疹基础再生数在麻疹发病数量保持不变的情况下,当疫苗保护率分别为50%、60%、70%和80%时,相应的基础再生数分别为2.49(95%CI:2.45～2.55)、2.61(95%CI:2.55～2.69)、3.34(95%CI:3.29～3.48)和 5.07(95%CI:4.96～5.24)。2.山东省2014～2015年17地市麻疹基础再生数山东省17地市麻疹基础再生数从高到低依次为4.20(济宁)、4.12(泰安)、3.93(菏泽)、3.92(德州)、3.81(潍坊)、3.81(莱芜)、3.47(烟台)、3.45(日照)、3.42(济南)、3.33(东营)、3.27(枣庄)、3.21(淄博)、2.98(滨州)、2.91(威海)、2.87(青岛)、2.39(临沂)和2.35(聊城)。全省呈现西高东低,南高北低的趋势。3.消除麻疹目标疫苗保护率根据模拟计算结果,假设真实疫苗保护率为60%,若通过疫苗免疫使疫苗保护率提高至75%,麻疹消除概率为30.32%;疫苗保护率提高至78.13%,麻疹消除概率为95.15%;疫苗保护率提高至80.00%,麻疹的消除概率为99.80%。假设真实疫苗保护率为70%,若通过免疫措施使疫苗保护率提高至81.55%时,则有95.30%的概率实现麻疹消除,当目标免疫率达到82.50%时,麻疹的消除概率为99.22%。结论1.不同城市之间麻疹传播能力不同,全省呈现西高东低,南高北低的趋势。2.如果真实世界疫苗保护率较低,则意味着人群中更多的人需要进行免疫保护,而如果人群疫苗免疫保护率高,则疫苗保护率稍加提高,即可达到消除麻疹的目的。3.状态空间模型由状态模型和观测模型组成,该模型考虑到观测数据的误差,能够有效地模拟仿真山东省麻疹实际疫情数据。