2016—2021年厦门市流行性腮腺炎流行特征分析与趋势预测

林艺兰，张怡盾，郭志南

厦门市疾病预防控制中心，福建 厦门 361021

流行性腮腺炎（以下简称流腮）是因腮腺炎病毒感染引发的急性呼吸系统传染病，呈全球流行，我国的发病率也一直居高不下[1]。作为丙类传染病，流腮传染性强，常在幼托机构、学校等集体场所暴发流行，极大干扰正常教学工作，严重影响我国儿童青少年身心健康[2]。了解流腮的流行特征和准确预测其变化趋势，有利于政府部门风险研判并合理配置卫生资源。GM（1，1）模型较好地填补概率统计的缺陷，可有效地降低原始数列随机波动性的影响，具备运算方便、所需信息较少等优点[3]。本研究收集分析厦门市2016—2021年流腮的流性特征和变化趋势，并用GM（1，1）模型预测未来三年厦门市流行性腮腺炎发病率，以期为制定综合防控措施提供参考依据。

1 资料与方法

1.1 资料来源 流腮发病资料和厦门市常住人口资料来自国家传染病监测系统，选择发病日期为2016年1月1日—2021年12月31日、常住地址为厦门市、诊断结果为流腮的发病数据，剔除重复报告病例。不同地区、年份和性别人群发病率比较时，以2010年全国人口普查数据为标准人口计算标化率，以消除年龄结构差异的影响。

1.2 趋势分析和预测方法 用率的年度变化百分比（annual mean percentage change,APC）评价发病率的时间变化趋势。用线性回归法计算APC，公式：APC（%）=(ea-1)×100%，y=ax+b，其中y表示各年发病率的自然对数值，x为年份，斜率的检验采用t检验。用拟和GM（1，1）模型预测发病率。建模采用累加法生成数据，先进行事前检验，阈值为，其中n为用于建模的发病率个数，若事前检验未通过，尝试对发病率进行平方根转换；用后验差比值（，其中S e表示残差标准差，S x表示发病率标准差）评价模型拟合效果，C≥0.65为不合格模型，不能进行预测，尝试减少建模的发病率个数重新拟合模型，C＜0.35模型为优秀。用平均绝对百分比误差（mean absolute percentage error,MAPE）评价预测精度，其中表示实际值，xt表示预测值。

1.3 统计分析 用Excel 2007进行数据整理和汇总，用SAS 9.3进行统计分析和建模。计量资料符合正态分布使用描述其集中趋势，不符合正态分布使用M（Q1，Q3）描述其集中趋势，检验水准取α=0.05。

2 结果

2.1 流腮报告病例情况 2016—2021年，厦门市共报告流腮2 332例，年发病率在3.99/10万～13.60/10万间波动，平均发病率9.20/10万，标化发病率10.30/10万，其中男性平均发病率10.98/10万，标化发病率11.86/10万；女性平均发病率7.27/10万，标化发病率8.43/10万，男性发病率高于女性。

从年龄分布看，发病年龄中位数为7（5，11）岁。1～14岁儿童共报告发病1 976例（占报告发病总数的84.73%），平均发病率57.77/10万，以学生（1 043例）和幼托机构儿童（684例）为主。发病高峰出现在5～9岁组，平均发病率71.30/10万。见图1。

图1 2016—2021年厦门市流腮发病率年龄分布

从月份分布看，各月均有发病病例，呈现双高峰分布，发病高峰在5月（286例）和11月（295例），2月和8月发病数相对较少。见图2。

图2 2016—2021年厦门市流腮发病例数月份分布

2.2 流腮发病率变化趋势 2016—2021年各年厦门市流腮的发病情况见表1。全市发病率APC为-21.61%（t=-3.49，P＜0.05），提示全市发病率年均下降21.61%，其中男性发病率年均下降23.47%（t=-3.95，P＜0.05），女性发病率呈波动变化（t=-2.74，P＞0.05）。

表1 2016—2021年厦门市流腮发病情况

2.3 流腮发病率预测结果 2020年厦门市流腮发病率下降明显，数据经平方根转换后，事前检验值为0.656 7，超过阈值，无法拟合模型。尝试去除部分年份后，最终男、女及总体人群均用2018—2021年发病率数据建模。后验差检验显示，C在0.15～0.30之间，三个模型均为优秀，可用于预测，建模参数和预测结果见表2。全市发病率预测值与实际值的MAPE为11.14%，预测2022—2024年厦门市流腮发病率分别为1.73/10万、0.93/10万和0.50/10万。

表2 GM（1，1）模型预测厦门市流腮发病率结果

3 讨论

流腮传染性强，在空间狭小、人员密集的场所易引起暴发流行，处理不当还可能并发脑膜炎、胰腺炎、睾丸炎等疾病，甚至导致死亡，危害不容忽视[4]。2016年全国流腮报告发病率12.84/10万[5]，本研究结果显示2016年厦门市流腮报告发病率13.68/10万，略高于全国平均水平，远高于福建省2016年发病率（7.7/10万）[6]，说明厦门市流腮发病率处于全省中上水平。男性发病率高于女性的流行特点与多数地区报道一致[7-8]，可能与女性卫生习惯较好有关。发病人群主要集中于1～14岁，以学生和托幼机构儿童为主，这与学校和托幼机构人群密切接触频率高、易形成利于流腮传播的环境有关。发病率最低的月份是2月和8月，正是学生寒暑假期间，进一步提示流腮防控应重点关注幼托儿童和中小学生，以托幼机构、学校等集体单位为重点防控场所，严格落实入学预防接种查验和每日晨检制度。

流腮的防控措施包括行为干预和疫苗接种。沿海省份外来人口多，常以疫苗接种为主要抓手[9]。福建省2008年将麻腮风联合疫苗纳入一类疫苗免疫规划项目，腮腺炎的发病率上升趋势逐渐得到遏制，全省发病率从2011年的43.2/10万下降至2017年的7.5/10万[6]。本研究结果也显示，厦门市流腮发病率存在下降趋势，年均下降21.61%，与浙江、江西等省份变化趋势一致[10-11]，下降原因除与疫苗接种有关外，也与新冠疫情发生以来居民佩戴口罩、避免去人口密集场所等行为有关。提示良好的个人卫生行为对流腮感染的发生有很好的预防作用，今后应进一步加强该方面的宣传。男性发病率下降明显，而女性发病率存在波动，可能与男性此前发病率较高，降幅相对更加明显有关，也可能是新冠疫情期间行为习惯的改变对男性发病率影响更大。

目前用于流行性腮腺炎发病预测的数学模型主要有动力学模型、神经网络模型、自回归平均移动模型、灰色模型等[5]。动力学模型参数设定易受主观和经验因素影响，较难解释。神经网络模型尤适用于处理模糊不清或不精确的信息，但其选取神经元及网络层数过程可能出现较大误差进而影响预测准确性[12]。自回归平均移动（autoregressive integrated moving average,ARIMA）模型短期预测精确度高且实用性强，长期预测需不断加入新值，拟合修正模型[13]。灰色模型中的GM（1，1）模型对建模原始数据要求较低，建模相对简单，被广泛应用于疾病发病率预测[14]。本研究用GM（1，1）模型拟合厦门市流腮发病率，预测值与实际值误差（MAPE）为11.14%，优于彭颖等用ARIMA模型预测武汉市流腮发病率的MAPE（23.64%）[15]，说明GM（1，1）模型也可用于流腮发病率预测。受新冠肺炎疫情防控措施的影响，2020年厦门市流腮发病率出现较大幅度下降，数据波动过大，导致本研究建模只能用4年发病率数据拟合模型，建模数据偏少，有待于进一步收集数据以不断修正模型，提高模型稳定性和预测准确性。根据预测结果，在维持现有防控措施下，2022—2024年厦门市流腮发病率存在继续下降趋势。2022年12月开始，厦门市也逐步放开新冠肺炎管控措施，短期内人们仍将保持戴口罩习惯，注意减少人群聚集，有利于控制流腮发病；长期上看，人们的防控意识可能逐渐淡化，有必要加强宣传教育，提高人们自我防护意识，防止流腮发病出现反弹或暴发。

综上所述，厦门市近年来流腮发病率存在下降趋势，但仍处在福建省中上水平，应进一步加强疫苗接种和行为干预。