监测
甘肃省2016年7种常见急性传染病预警阈值研究
中华流行病学杂志, 2018,39(3) : 352-356. DOI: 10.3760/cma.j.issn.0254-6450.2018.03.020
摘要
目的

优选甘肃省常见急性传染病的预警阈值,提高预警功效。

方法

甘肃省流行性感冒(流感)、猩红热、其他感染性腹泻、细菌性和阿米巴性痢疾(痢疾)、伤寒/副伤寒、戊型肝炎(戊肝)采用移动百分位数法,手足口病采用累积和控制图法,分别建立预警模型,通过计算灵敏度、特异度、阳性和阴性预测值、约登指数及绘制受试者工作特征曲线(receiver operating characteristic curve,ROC),选择出适合甘肃省传染病的最优预警阈值。

结果

流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝的最优预警界值分别为P90P80P95P90P80P90,手足口病的最优预警参数为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显示此预警界值/参数均为最靠近坐标图左上方的点。

结论

甘肃省流感、其他感染性腹泻、痢疾、戊肝的预警阈值可适当调高,手足口病预警参数需调整,从而提高预警功效。

引用本文: 成瑶, 刘新凤, 孟蕾, 等.  甘肃省2016年7种常见急性传染病预警阈值研究 [J] . 中华流行病学杂志, 2018, 39(3) : 352-356. DOI: 10.3760/cma.j.issn.0254-6450.2018.03.020.
参考文献导出:   Endnote    NoteExpress    RefWorks    NoteFirst    医学文献王
扫  描  看  全  文

正文
作者信息
基金 0  关键词  0
English Abstract
评论
阅读 0  评论  0
相关资源
引用 | 论文 | 视频

版权归中华医学会所有。

未经授权,不得转载、摘编本刊文章,不得使用本刊的版式设计。

除非特别声明,本刊刊出的所有文章不代表中华医学会和本刊编委会的观点。

2008年中国CDC基于传染病信息报告管理系统开发了传染病自动预警系统,确定了各传染病的预警阈值并于2010年12月进行统一调整后运行至今。甘肃省自传染病自动预警系统运行以来,所使用的阈值未做过调整,在实际工作中发现存在一些问题,主要是预警阳性率偏低。为提高预警效果,达到提早、及时发现传染病异常增加的目的,本研究选取甘肃省7种常见急性传染病,建立预警模型,通过分析各项指标,从而选择出适合甘肃省传染病的最优预警阈值。

资料与方法
1.资料来源:

2011-2016年甘肃省流行性感冒(流感)、猩红热、其他感染性腹泻、细菌性和阿米巴性痢疾(痢疾)、伤寒/副伤寒、戊型肝炎(戊肝)、手足口病7种常见传染病数据来源于中国疾病预防控制系统传染病报告信息管理系统,以周为单位整理并建立数据库。2011-2015年数据为历史基线数据,对2016年数据进行预警。

2.研究方法:

本研究采用移动百分位数法和累积和控制图法,针对不同传染病分别建立预警模型,将预警结果与历史流行水平进行对比分析,从而计算出各预警模型在不同预警界值或者不同预警参数取值下的灵敏度、特异度、阳性和阴性预测值、约登指数,并绘制受试者工作特征曲线(receiver operating characteristic curve,ROC),比较各项指标,选择出适合甘肃省常见传染病的最优预警阈值。

(1)流行的判断标准:以常见传染病2011-2016年发病率(±2s)作为判断流行的参考标准。

(2)预警模型:对甘肃省流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝采用移动百分位数法,手足口病采用累积和控制图法,建立预警模型。①移动百分位数法:以一定置信水平的百分位数作为上控限制,即预警限,作为监视序列是否"失控"的判断依据[1]。选取2011-2015年及同期前后摆动2周的周发病率数据作为历史基线数据,对2016年的疫情数据进行预警,计算不同预警界值(P60P65P70P75P80P90P95),建立预警模型。②累积和控制图法:对当前值与来自历史基线的预期值之间差值累计求和,表示当前值偏离基线的程度,当合计值超过设定的预警界值,则发出预警。其基本公式:Ct=max{0,Xt-(μt+)+Ct-1}[2,3]。初始值C0=0,k是该模型的重要参数,判定值设为H,用于判定疾病是否存在异常,其值的选取可影响模型的检出效果,为允偏量,若均数μt偏移到μt+时,希望引起警觉,当CtH,发出预警信号[2,4]。根据既往文献[5,6],本研究将Hk的取值范围设定为3σH≤5σ,0<k≤1.5,用移动标准差s作为σ的估计值,移动平均值作为μ的估计值。Ct通过不同的计算方式发出灵敏度由低到高的3种预警信号:C1、C2及C3,本研究选择灵敏度最高的C3,设定不同的参数,建立预警模型。

(3)评价指标:①灵敏度:预警模型根据流行标准发现真流行事件的能力;②特异度:预警模型根据流行标准判定非流行事件的能力;③阳性预测值:预警模型发出的预警信号中真流行事件所占的比例;④阴性预测值:预警模型未发出预警信号时非流行事件所占比例;⑤约登指数:反映预警模型发现流行事件与非流行事件的综合能力,其取值范围为0~1,越接近1,提示预警模型的真实性越高;⑥ROC:以灵敏度为纵坐标轴,假阳性率(1-特异度)为横坐标轴绘制,通常将曲线上最靠近坐标图左上角的点定为灵敏度和特异度最佳的临界值,其所对应的预警阈值功效最佳。

3.统计学分析:

采用Excel 2010软件计算各项指标及绘制ROC。

结果
1.不同预警界值预警功效比较:

用移动百分位数法预警甘肃省6种常见传染病发现,随着预警界值的升高,灵敏度和阴性预测值降低,特异度和阳性预测值升高;流感、猩红热、伤寒/副伤寒在P80时灵敏度达到100.00%且特异度较高,其他感染性腹泻、痢疾、戊肝在P90时灵敏度达到100.00%且特异度较高;猩红热、伤寒/副伤寒在P80时约登指数为0.62、0.98,痢疾、戊肝在P90时约登指数为1.00、0.89,其他感染性腹泻在P95时约登指数为0.67,流感在P80P90时约登指数均为0.73,见表1。通过绘制ROC显示,流感、猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝最靠近坐标图左上角的点分别为P90P80P95P90P80P90,见图1

点击查看表格
表1

2016年甘肃省常见传染病移动百分位数法不同预警界值预警功效比较

表1

2016年甘肃省常见传染病移动百分位数法不同预警界值预警功效比较

疾病评价指标预警界值
P60P65P70P75P80P90P95
流行性感冒灵敏度100.00100.00100.00100.00100.0086.6773.33
 特异度43.2448.6559.4664.8672.9786.4991.89
 阳性预测值41.6744.1250.0053.5760.0072.2278.57
 阴性预测值100.00100.00100.00100.00100.0094.1289.47
 约登指数0.430.490.590.650.730.730.65
猩红热灵敏度100.00100.00100.00100.00100.0071.4357.14
 特异度28.8940.0042.2257.7862.2273.3380.00
 阳性预测值17.9520.5921.2126.9229.1729.4130.77
 阴性预测值100.00100.00100.00100.00100.0094.2992.31
 约登指数0.290.400.420.580.620.450.37
其他感染性腹泻灵敏度100.00100.00100.00100.00100.00100.0091.67
 特异度27.5032.5037.5040.0045.0065.0075.00
 阳性预测值29.2730.7732.4333.3335.2946.1552.38
 阴性预测值100.00100.00100.00100.00100.00100.0096.77
 约登指数0.280.330.380.400.450.650.67
细菌性和阿米巴性痢疾灵敏度100.00100.00100.00100.00100.00100.000.00
 特异度92.1692.1692.1694.1298.04100.00100.00
 阳性预测值20.0020.0020.0025.0050.00100.00-
 阴性预测值100.00100.00100.00100.00100.00100.0098.08
 约登指数0.920.920.920.940.981.000.00
伤寒/副伤寒灵敏度100.00100.00100.00100.00100.0075.0050.00
 特异度90.2490.7093.7595.8397.92100.00100.00
 阳性预测值50.0050.0057.1466.6780.00100.00100.00
 阴性预测值100.00100.00100.00100.00100.0097.9696.00
 约登指数0.900.910.940.960.980.750.50
戊型肝炎灵敏度100.00100.00100.00100.00100.00100.0083.33
 特异度67.3969.5771.7478.2680.4389.1391.30
 阳性预测值28.5730.0031.5837.5040.0054.5555.56
 阴性预测值100.00100.00100.00100.00100.00100.0097.67
 约登指数0.670.700.720.780.800.890.75

注:"-"表示细菌性和阿米巴性痢疾在P95时未发出预警信号,因此不能计算阳性预测值

点击查看大图
图1
2016年甘肃省6种常见传染病受试者工作特征曲线
点击查看大图
图1
2016年甘肃省6种常见传染病受试者工作特征曲线
2.不同预警参数预警功效比较:

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

点击查看表格
表2

2016年甘肃省手足口病累积和控制图法C3模型不同参数预警功效比较

表2

2016年甘肃省手足口病累积和控制图法C3模型不同参数预警功效比较

指标k=1.1k=1.2k=1.3
H=3σH=4σH=5σH=3σH=4σH=5σH=3σH=4σH=5σ
灵敏度100.00100.00100.00100.00100.00100.0040.0020.0020.00
特异度63.8365.9670.2165.9670.2174.4778.7287.2395.65
阳性预测值22.7323.8126.3223.8126.3229.4116.6714.2933.33
阴性预测值100.00100.00100.00100.00100.00100.0092.5091.1191.67
约登指数0.640.660.700.650.700.740.190.070.16
点击查看大图
图2
2016年甘肃省手足口病的受试者工作特征曲线
点击查看大图
图2
2016年甘肃省手足口病的受试者工作特征曲线
3.甘肃省常见急性传染病预警阈值优选:

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

点击查看表格
表3

2016年甘肃省7种常见传染病预警阈值优选结果

表3

2016年甘肃省7种常见传染病预警阈值优选结果

病种最优预警阈值甘肃省最优阈值国家设定预警阈值
灵敏度特异度阳性预测值阴性预测值约登指数ROC
流行性感冒P60P80P95P95P60P80P80P90P90P90P80
猩红热P60P80P95P95P60P80P80P80P80P80
其他感染性腹泻P60P90P95P95P60P90P95P95P95P80
细菌性和阿米巴性痢疾P60P90P90P95P90P60P90P90P90P90P80
伤寒/副伤寒P60P80P90P95P90P95P60P80P80P80P80P80
戊型肝炎P60P90P95P95P60P90P90P90P90P80
手足口病k=1.1~1.2k=1.3k=1.3k=1.1~1.2k=1.2k=1.2k=1.2k=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种常见传染病发现,流感在P80P90时各指标达到平衡点,但考虑到流感为丙类传染病,病例多、危害小,应适当提高特异度,减小假阳性预警信号,因此选择P90为最优阈值;猩红热、其他感染性腹泻、痢疾、伤寒/副伤寒、戊肝参考ROC,同时选择各指标的最佳平衡点为最优阈值,分别为P80P95P90P80P90。利用累积和控制图法预警手足口病,综合考虑灵敏度、特异度、阳性预测值、阴性预测值等指标,参照ROC,拟合出最优预警参数为k=1.2,H=5σ。本研究结果不同于其他省份常见传染病预警阈值研究结果[7,9],与孙乔等[11]、张洪龙等[12]对全国的传染病预警阈值的研究结果不同。

本研究存在不足。拟合出的甘肃省常见传染病的最优阈值与中国传染病自动预警系统设置的阈值不尽一致,需要对流感、其他感染性腹泻、痢疾、戊肝、手足口病的阈值进行调整,但由于中国CDC设定的传染病自动预警系统可调预警阈值范围为P60~P80,本研究优选出需要调整预警阈值的病种均不在可调范围之内。因此,如何调整甘肃省常见传染病的预警阈值,还需进一步与中国CDC进行沟通和探讨。

综上所述,传染病自动预警系统现已成为各级疾病预防控制机构早期发现传染病暴发的重要辅助工具,通过拟合优选甘肃省常见急性传染病的预警阈值,将对流感、其他感染性腹泻、痢疾、戊肝的预警阈值适当调高,手足口病预警参数调整,从而减轻基层工作者的工作量,达到提高预警功效的目的。

利益冲突
利益冲突

参考文献
[1]
杨维中.传染病预警理论与实践[M].北京人民卫生出版社2012.
YangWZ. Early warning of infectious disease theory and practice[M]. BeijingPeople’s Medical Publishing House, 2012.
[2]
王全意.疾病监测信息报告管理系统数据分析手册[M].北京中国协和医科大学出版社2010.
WangQY. Data analysis manual for the disease surveillance information reporting and management system[M]. BeijingPeking Union Medical College Press, 2010.
[3]
杨筱婷刘新凤孟蕾.甘肃省流感流行预警方法探研[J].中华流行病学杂志201637(3):430-433. DOI:10.3760/cma.j.issn.0254-6450.2016.03.029.
YangXT, LiuXF, MengL, et al. Study on early warning method for influenza epidemic in Gansu province[J]. Chin J Epidemiol, 2016, 37(3):430-433. DOI:10.3760/cma.j.issn.0254-6450.2016.03.029.
[4]
LiZJ, LaiSJ, BuckeridgeDL, et al. Adjusting outbreak detection algorithms for surveillance during epidemic and non-epidemic periods[J]. J Am Med Inform Associat, 2012, 19(el):e51-53. DOI:10.1136/amiajnl-2011-000126.
[5]
CarpenterTE. Evaluation and extension of the Cusum technique with an application to Salmonella surveillance[J]. J Vet Diagn Invest, 2002, 14(3):211-218. DOI:10.1177/104063870201400304.
[6]
李静杨鹏吴双胜.应用累积和模型探测北京市猩红热流行起始时间的预警分析[J].中华流行病学杂志201334(5):526-530. DOI:10.3760/cma.j.issn.0254-6450.2013.05.025
LiJ, YangP, WuSS, et al. Early detection on the onset of scarlet fever epidemics in Beijing, using the Cumulative Sum[J]. Chin J Epidemiol2013, 34(5):526-530. DOI:10.3760/cma.j.issn.0254-6450.2013.05.025
[7]
王棠明小燕李十月.湖北省手足口病控制图法预警效果评价[J].中国公共卫生201531(10):1324-1326. DOI:10.11847/zgggws2015-31-10-25.
WangT, MingXY, LiSY. Efficiency of control graph method in early warning for hand-foot-mouth disease in Hubei province, China[J]. Chin J Public Health, 2015, 31(10):1324-1326. DOI:10.11847/zgggws2015-31-10-25.
[8]
王瑞平春雅丽郭晓芹.上海市松江区手足口病控制图法预警界值优选研究[J].中国卫生统计201128(3):264-265. DOI:10.3969/j.issn.1002-3674.2011.03.013.
WangRP, ChunYL, GuoXQ, et al. A study of early detection for hand-foot-mouth disease in Songjiang district of Shanghai[J]. Chin J Health Stat, 2011, 28(3):264-265. DOI:10.3969/j.issn.1002-3674.2011.03.013.
[9]
王瑞平陈立凌刘成.苏州市重点传染病移动百分位数法预警界值优选分析[J].中国卫生统计201431(4):598-601,605.
WangRP, ChenLL, LiuC, et al. Study of early detection for principal communicable diseases by moving percentile method in Suzhou City[J]. Chin J Health Stat, 2014, 31(4):598-601, 605.
[10]
刘新凤成瑶刘东鹏.甘肃省2008-2013年传染病自动预警系统时间模型运行结果分析[J].中国公共卫生201733(4):634-637. DOI:10.11847/zgggws2017-33-04-30.
LiuXF, ChengY, LiuDP, et al. Performance efficiency of temporal model in China infectious diseases automated-alert and response system in Gansu province, 2008-2013[J]. Chin J Public Health, 2017, 33(4):634-637. DOI:10.11847/zgggws2017-33-04-30.
[11]
孙乔赖圣杰李中杰.移动百分位数法不同阈值在传染病暴发预警中的效果比较[J].中华流行病学杂志201132(5):450-453. DOI:10.3760/cma.j.issn.0254-6450.2011.05.006.
SunQ, LaiSJ, LiZJ, et al. Comparison on the different thresholds on the 'moving percentile method’ for outbreak detection[J]. Chin J Epidemiol, 2011, 32(5):450-453. DOI:10.3760/cma.j.issn.0254-6450.2011.05.006.
[12]
张洪龙孙乔赖圣杰.移动百分位数法分地区设定预警阈值对传染病预警效果的影响分析[J].中华预防医学杂志201448(4):265-269. DOI:10.3760/cma.j.issn.0253-9624.2014.04.006.
ZhangHL, SunQ, LaiSJ, et al. Analysis of effect on infectious diseases outbreak detection performance by classifying provinces for moving percentile method[J]. Chin J Prev Med, 2014, 48(4):265-269. DOI:10.3760/cma.j.issn.0253-9624.2014.04.006.
 
 
展开/关闭提纲
查看图表详情
回到顶部
放大字体
缩小字体
标签
关键词