
优选甘肃省常见急性传染病的预警阈值,提高预警功效。
甘肃省流行性感冒(流感)、猩红热、其他感染性腹泻、细菌性和阿米巴性痢疾(痢疾)、伤寒/副伤寒、戊型肝炎(戊肝)采用移动百分位数法,手足口病采用累积和控制图法,分别建立预警模型,通过计算灵敏度、特异度、阳性和阴性预测值、约登指数及绘制受试者工作特征曲线(receiver operating characteristic curve,ROC),选择出适合甘肃省传染病的最优预警阈值。
流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝的最优预警界值分别为P90、P80、P95、P90、P80、P90,手足口病的最优预警参数为k=1.2,H=5σ。在最优预警界值/参数下,流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝、手足口病的灵敏度分别为86.67%、100.00%、91.67%、100.00%、100.00%、100.00%、100.00%,特异度分别为86.49%、62.22%、75.00%、100.00%、97.92%、89.13%、74.47%,阳性预测值分别为72.22%、29.17%、52.38%、100.00%、80.00%、54.55%、29.41%,阴性预测值分别为94.12%、100.00%、96.77%、100.00%、100.00%、100.00%、100.00%,约登指数分别为0.73、0.62、0.67、1.00、0.98、0.89、0.74,ROC显示此预警界值/参数均为最靠近坐标图左上方的点。
甘肃省流感、其他感染性腹泻、痢疾、戊肝的预警阈值可适当调高,手足口病预警参数需调整,从而提高预警功效。
版权归中华医学会所有。
未经授权,不得转载、摘编本刊文章,不得使用本刊的版式设计。
除非特别声明,本刊刊出的所有文章不代表中华医学会和本刊编委会的观点。
2008年中国CDC基于传染病信息报告管理系统开发了传染病自动预警系统,确定了各传染病的预警阈值并于2010年12月进行统一调整后运行至今。甘肃省自传染病自动预警系统运行以来,所使用的阈值未做过调整,在实际工作中发现存在一些问题,主要是预警阳性率偏低。为提高预警效果,达到提早、及时发现传染病异常增加的目的,本研究选取甘肃省7种常见急性传染病,建立预警模型,通过分析各项指标,从而选择出适合甘肃省传染病的最优预警阈值。
2011-2016年甘肃省流行性感冒(流感)、猩红热、其他感染性腹泻、细菌性和阿米巴性痢疾(痢疾)、伤寒/副伤寒、戊型肝炎(戊肝)、手足口病7种常见传染病数据来源于中国疾病预防控制系统传染病报告信息管理系统,以周为单位整理并建立数据库。2011-2015年数据为历史基线数据,对2016年数据进行预警。
本研究采用移动百分位数法和累积和控制图法,针对不同传染病分别建立预警模型,将预警结果与历史流行水平进行对比分析,从而计算出各预警模型在不同预警界值或者不同预警参数取值下的灵敏度、特异度、阳性和阴性预测值、约登指数,并绘制受试者工作特征曲线(receiver operating characteristic curve,ROC),比较各项指标,选择出适合甘肃省常见传染病的最优预警阈值。
(1)流行的判断标准:以常见传染病2011-2016年发病率(
±2s)作为判断流行的参考标准。
(2)预警模型:对甘肃省流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝采用移动百分位数法,手足口病采用累积和控制图法,建立预警模型。①移动百分位数法:以一定置信水平的百分位数作为上控限制,即预警限,作为监视序列是否"失控"的判断依据[1]。选取2011-2015年及同期前后摆动2周的周发病率数据作为历史基线数据,对2016年的疫情数据进行预警,计算不同预警界值(P60、P65、P70、P75、P80、P90、P95),建立预警模型。②累积和控制图法:对当前值与来自历史基线的预期值之间差值累计求和,表示当前值偏离基线的程度,当合计值超过设定的预警界值,则发出预警。其基本公式:Ct=max{0,Xt-(μt+kσ)+Ct-1}[2,3]。初始值C0=0,k是该模型的重要参数,判定值设为H,用于判定疾病是否存在异常,其值的选取可影响模型的检出效果,kσ为允偏量,若均数μt偏移到μt+kσ时,希望引起警觉,当Ct≥H,发出预警信号[2,4]。根据既往文献[5,6],本研究将H和k的取值范围设定为3σ≤H≤5σ,0<k≤1.5,用移动标准差s作为σ的估计值,移动平均值作为μ的估计值。Ct通过不同的计算方式发出灵敏度由低到高的3种预警信号:C1、C2及C3,本研究选择灵敏度最高的C3,设定不同的参数,建立预警模型。
(3)评价指标:①灵敏度:预警模型根据流行标准发现真流行事件的能力;②特异度:预警模型根据流行标准判定非流行事件的能力;③阳性预测值:预警模型发出的预警信号中真流行事件所占的比例;④阴性预测值:预警模型未发出预警信号时非流行事件所占比例;⑤约登指数:反映预警模型发现流行事件与非流行事件的综合能力,其取值范围为0~1,越接近1,提示预警模型的真实性越高;⑥ROC:以灵敏度为纵坐标轴,假阳性率(1-特异度)为横坐标轴绘制,通常将曲线上最靠近坐标图左上角的点定为灵敏度和特异度最佳的临界值,其所对应的预警阈值功效最佳。
采用Excel 2010软件计算各项指标及绘制ROC。
用移动百分位数法预警甘肃省6种常见传染病发现,随着预警界值的升高,灵敏度和阴性预测值降低,特异度和阳性预测值升高;流感、猩红热、伤寒/副伤寒在P80时灵敏度达到100.00%且特异度较高,其他感染性腹泻、痢疾、戊肝在P90时灵敏度达到100.00%且特异度较高;猩红热、伤寒/副伤寒在P80时约登指数为0.62、0.98,痢疾、戊肝在P90时约登指数为1.00、0.89,其他感染性腹泻在P95时约登指数为0.67,流感在P80和P90时约登指数均为0.73,见表1。通过绘制ROC显示,流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝最靠近坐标图左上角的点分别为P90、P80、P95、P90、P80、P90,见图1。

2016年甘肃省常见传染病移动百分位数法不同预警界值预警功效比较
2016年甘肃省常见传染病移动百分位数法不同预警界值预警功效比较
| 疾病 | 评价指标 | 预警界值 | ||||||
|---|---|---|---|---|---|---|---|---|
| P60 | P65 | P70 | P75 | P80 | P90 | P95 | ||
| 流行性感冒 | 灵敏度 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 86.67 | 73.33 |
| 特异度 | 43.24 | 48.65 | 59.46 | 64.86 | 72.97 | 86.49 | 91.89 | |
| 阳性预测值 | 41.67 | 44.12 | 50.00 | 53.57 | 60.00 | 72.22 | 78.57 | |
| 阴性预测值 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 94.12 | 89.47 | |
| 约登指数 | 0.43 | 0.49 | 0.59 | 0.65 | 0.73 | 0.73 | 0.65 | |
| 猩红热 | 灵敏度 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 71.43 | 57.14 |
| 特异度 | 28.89 | 40.00 | 42.22 | 57.78 | 62.22 | 73.33 | 80.00 | |
| 阳性预测值 | 17.95 | 20.59 | 21.21 | 26.92 | 29.17 | 29.41 | 30.77 | |
| 阴性预测值 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 94.29 | 92.31 | |
| 约登指数 | 0.29 | 0.40 | 0.42 | 0.58 | 0.62 | 0.45 | 0.37 | |
| 其他感染性腹泻 | 灵敏度 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 91.67 |
| 特异度 | 27.50 | 32.50 | 37.50 | 40.00 | 45.00 | 65.00 | 75.00 | |
| 阳性预测值 | 29.27 | 30.77 | 32.43 | 33.33 | 35.29 | 46.15 | 52.38 | |
| 阴性预测值 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 96.77 | |
| 约登指数 | 0.28 | 0.33 | 0.38 | 0.40 | 0.45 | 0.65 | 0.67 | |
| 细菌性和阿米巴性痢疾 | 灵敏度 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 0.00 |
| 特异度 | 92.16 | 92.16 | 92.16 | 94.12 | 98.04 | 100.00 | 100.00 | |
| 阳性预测值 | 20.00 | 20.00 | 20.00 | 25.00 | 50.00 | 100.00 | - | |
| 阴性预测值 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 98.08 | |
| 约登指数 | 0.92 | 0.92 | 0.92 | 0.94 | 0.98 | 1.00 | 0.00 | |
| 伤寒/副伤寒 | 灵敏度 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 75.00 | 50.00 |
| 特异度 | 90.24 | 90.70 | 93.75 | 95.83 | 97.92 | 100.00 | 100.00 | |
| 阳性预测值 | 50.00 | 50.00 | 57.14 | 66.67 | 80.00 | 100.00 | 100.00 | |
| 阴性预测值 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 97.96 | 96.00 | |
| 约登指数 | 0.90 | 0.91 | 0.94 | 0.96 | 0.98 | 0.75 | 0.50 | |
| 戊型肝炎 | 灵敏度 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 83.33 |
| 特异度 | 67.39 | 69.57 | 71.74 | 78.26 | 80.43 | 89.13 | 91.30 | |
| 阳性预测值 | 28.57 | 30.00 | 31.58 | 37.50 | 40.00 | 54.55 | 55.56 | |
| 阴性预测值 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 97.67 | |
| 约登指数 | 0.67 | 0.70 | 0.72 | 0.78 | 0.80 | 0.89 | 0.75 | |
注:"-"表示细菌性和阿米巴性痢疾在P95时未发出预警信号,因此不能计算阳性预测值


累积和控制图法C3模型预警甘肃省手足口病,选取不同的预警参数计算显示,当k=1.2,H=5σ时,特异度较高、灵敏度达到100.00%且约登指数最接近1,见表2。通过绘制ROC显示,手足口病最靠近坐标轴左上角的点预警参数取值为k=1.2,H=5σ,见图2。

2016年甘肃省手足口病累积和控制图法C3模型不同参数预警功效比较
2016年甘肃省手足口病累积和控制图法C3模型不同参数预警功效比较
| 指标 | k=1.1 | k=1.2 | k=1.3 | ||||||
|---|---|---|---|---|---|---|---|---|---|
| H=3σ | H=4σ | H=5σ | H=3σ | H=4σ | H=5σ | H=3σ | H=4σ | H=5σ | |
| 灵敏度 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 40.00 | 20.00 | 20.00 |
| 特异度 | 63.83 | 65.96 | 70.21 | 65.96 | 70.21 | 74.47 | 78.72 | 87.23 | 95.65 |
| 阳性预测值 | 22.73 | 23.81 | 26.32 | 23.81 | 26.32 | 29.41 | 16.67 | 14.29 | 33.33 |
| 阴性预测值 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 | 92.50 | 91.11 | 91.67 |
| 约登指数 | 0.64 | 0.66 | 0.70 | 0.65 | 0.70 | 0.74 | 0.19 | 0.07 | 0.16 |


分别选出灵敏度、特异度、阳性预测值、阴性预测值、约登指数、ROC的最优预警阈值,综合各指标选择一个最佳平衡点,即为甘肃省常见急性传染病的最优阈值,流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝的最优预警界值分别为P90、P80、P95、P90、P80、P90,手足口病的最优预警参数为k=1.2,H=5σ,见表3。

2016年甘肃省7种常见传染病预警阈值优选结果
2016年甘肃省7种常见传染病预警阈值优选结果
| 病种 | 最优预警阈值 | 甘肃省最优阈值 | 国家设定预警阈值 | |||||
|---|---|---|---|---|---|---|---|---|
| 灵敏度 | 特异度 | 阳性预测值 | 阴性预测值 | 约登指数 | ROC | |||
| 流行性感冒 | P60~P80 | P95 | P95 | P60~P80 | P80、P90 | P90 | P90 | P80 |
| 猩红热 | P60~P80 | P95 | P95 | P60~P80 | P80 | P80 | P80 | P80 |
| 其他感染性腹泻 | P60~P90 | P95 | P95 | P60~P90 | P95 | P95 | P95 | P80 |
| 细菌性和阿米巴性痢疾 | P60~P90 | P90、P95 | P90 | P60~P90 | P90 | P90 | P90 | P80 |
| 伤寒/副伤寒 | P60~P80 | P90、P95 | P90、P95 | P60~P80 | P80 | P80 | P80 | P80 |
| 戊型肝炎 | P60~P90 | P95 | P95 | P60~P90 | P90 | P90 | P90 | P80 |
| 手足口病 | k=1.1~1.2 | k=1.3 | k=1.3 | k=1.1~1.2 | k=1.2 | k=1.2 | k=1.2 | k=1.3 |
| H=3σ~5σ | H=5σ | H=5σ | H=3σ~5σ | H=5σ | H=5σ | H=5σ | ||
注:ROC为受试者工作特征曲线
找出适宜的预警界值或预警模型参数是提高预警功效的关键。在实际工作中,根据不同传染病的发病特点以及防控要求,综合考虑灵敏度、特异度、阳性和阴性预测值、约登指数等指标,同时参照ROC,选择出不同传染病的最优阈值[7,8,9]。
移动百分位数法是一种非参数方法,不依赖数据统计学分布类型,普适性较好。本研究选取2011-2015年及同期前后摆动2周,共25个周发病率数据作为历史基线数据计算预警限,理论上对季节性进行了调整,同时前后摆动拓展了历史数据宽度,有利于模型的稳定性。累积和控制图法基于短基线数据,在探测疾病微小变化时具有较好的灵敏性和及时性。
每年通过传染病自动预警系统发出甘肃省的大量预警信号,但最终确认暴发的事件较少,预警阳性率低,甘肃省2008-2015年时间模型预警阳性率仅为0.64%[10],因此,拟合出适合甘肃省的传染病最优预警阈值意义重大。本研究用移动百分位数法预警甘肃省6种常见传染病发现,流感在P80和P90时各指标达到平衡点,但考虑到流感为丙类传染病,病例多、危害小,应适当提高特异度,减小假阳性预警信号,因此选择P90为最优阈值;猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝参考ROC,同时选择各指标的最佳平衡点为最优阈值,分别为P80、P95、P90、P80、P90。利用累积和控制图法预警手足口病,综合考虑灵敏度、特异度、阳性预测值、阴性预测值等指标,参照ROC,拟合出最优预警参数为k=1.2,H=5σ。本研究结果不同于其他省份常见传染病预警阈值研究结果[7,9],与孙乔等[11]、张洪龙等[12]对全国的传染病预警阈值的研究结果不同。
本研究存在不足。拟合出的甘肃省常见传染病的最优阈值与中国传染病自动预警系统设置的阈值不尽一致,需要对流感、其他感染性腹泻、痢疾、戊肝、手足口病的阈值进行调整,但由于中国CDC设定的传染病自动预警系统可调预警阈值范围为P60~P80,本研究优选出需要调整预警阈值的病种均不在可调范围之内。因此,如何调整甘肃省常见传染病的预警阈值,还需进一步与中国CDC进行沟通和探讨。
综上所述,传染病自动预警系统现已成为各级疾病预防控制机构早期发现传染病暴发的重要辅助工具,通过拟合优选甘肃省常见急性传染病的预警阈值,将对流感、其他感染性腹泻、痢疾、戊肝的预警阈值适当调高,手足口病预警参数调整,从而减轻基层工作者的工作量,达到提高预警功效的目的。
无





















