空间扫描统计量在手足口病空间聚集性探测中的应用

空间扫描统计量在手足口病空间聚集性探测中的应用

手足疾病（hmfd）是由各种肠道病毒（iv）引起的常见传染病，通常发生在5岁以下的儿童。在大多数情况下，症状轻微，主要表现为发热、手、脚、口腔等部位的皮疹和皮疹。一些儿童可能会导致心肌炎、肺水肿和病毒性脑膜脑炎等并发症。个别重症患儿如果病情发展快,会导致死亡。自1974年Schmidt等首次报道以来,全世界很多国家和地区先后报道了EV71的感染和流行情况。2008年5月2日,我国将手足口病纳入丙类法定报告传染病管理。北京市顺义区自2006年初开始有手足口病例报告,其后发病率呈逐年上升趋势。疫情资料表明,手足口病具有一定的时间和空间的聚集倾向或趋势。本文采用Tango等提出的Flexible空间扫描统计量对顺义区手足口病的空间聚集性进行研究,探讨手足口病发病的空间聚集性特征,为进一步探索影响手足口病发病的因素及制定有效的防控措施和卫生决策提供理论依据。数据和方法1.调查单位的信息清理和分析2009年1月至2010年12月北京市顺义区手足口病疫情资料来源于“中国疾病预防控制信息系统”,参照网络直报工作要求对资料进行初步清理,并以“发病日期”计算手足口病发病时间以及病例的发病年龄,现住详细地址以及报告地区编码进行空间聚集性分析。人口资料来源于“中国疾病预防控制基本信息系统”的数据(统计局人口数);顺义区MapInfo基础地理信息系统数据,包括区界及乡镇级以上界线等来自北京市疾病预防控制中心传地所。2.fluser-p值的计算和分析本研究采用Flexible空间扫描统计量分析软件FleXScanV3进行病例的空间聚集性探测。Flexible空间扫描统计量是2005年Tango等在Kulldorff圆形扫描统计量基础上提出的,主要是为了探测圆形扫描窗口探测不到的不规则形状聚集区域。对每一个扫描窗口,均根据Poisson分布计算出理论发病数,然后依据实际发病数和理论发病数构造检验统计量对数似然比(loglikelihoodratio,LLR),用LLR来评价扫描窗口内发病数的异常程度。通常LLR值越大,则该窗口发病数异常程度越高。通常将LLR最大的窗口定为发病数异常程度最高的窗口,然后对该窗口的统计学意义进行评价,即计算P值。Flexible空间扫描统计量计算的具体过程如下。将空间扫描统计量S定义为所有可能的扫描窗口Z中最大的似然比:其中L(z)是扫描窗口Z的似然函数值,L0是基于无效假设得到的似然函数值。nz为扫描窗口Z中的实际发病数,μ(Z)为根据无效假设得到的扫描窗口Z中预期发病数,N是整个地区的总发病数,μ(G)为根据无效假设得到的整个地区的预期发病数,且有μ(G)=N。则Poisson模型的似然比为:由于扫描统计量的概率分布极为复杂,本研究采用Kulldorff等提出的蒙特卡罗法(MonteCarlomethod)计算检验统计量的P值,其迭代次数选择999。此外,除了对LLR最大的窗口进行评价外,扫描统计量还可对LLR较大的其他窗口进行统计学意义的评价,从而尽可能找出所有发病数异常的区域。Flexible空间扫描统计量方法中扫描窗口采用预先设定最大相邻区域数K的方法,从而避免不合理的区域出现。参考以往文献可知,聚集最多覆盖地理区域个数不宜超过总地理区域个数的10%~15%,顺义区辖区共25个乡镇和街道,因此本研究在分析时使用cluster最多覆盖地理区域数K=3。由于顺义区手足口病发病情况在不同时间和不同地区差异较大,因此本研究对手足口病报告发病例数分乡镇分月份进行统计,利用人口数和报告发病数计算理论发病人数,以期获得更为全面和准确的聚集性探测结果。将数据导入FleXScanV3软件进行运算,并用MapInfo软件将探测结果进行可视化呈现。结果1.病例纳入分析2009年1月至2010年12月顺义区共报告手足口病3376例,其中25例属于疑似病例,2例为地址不详病例,将其删除,余下3349例病例纳入分析,其中临床诊断病例3244例,实验室诊断病例105例。将报告病例数按月份做图可知,顺义区手足口病发病高峰主要在4~9月份,以夏秋季发病最高,呈明显的季节性(图1)。2.足口病的聚集期在季节和季节内聚集的区域分布根据P值计算NOR(nulloccurrencerate),计算公式:NOR=1/P天。NOR是指按假设不存在聚集性情况下,出现虚假聚集信号的时间间隔。例如P=0.001,表示每1000天可能出现这样一个虚假聚集信号。由结果可知,不同月份顺义区手足口病空间聚集程度和空间聚集区域明显不同。手足口病高发月份空间聚集程度较低发月份强,空间聚集主要集中在4~9月份,这与手足口病的季节性发病特点一致。在手足口病发病率较低的月份,手足口病病例在地区间成散在的分布,无明显空间聚集性发生。随着发病率的上升,手足口病发生空间聚集cluster的个数逐渐上升,聚集区域也逐渐增多。另外,手足口病发生聚集的区域主要在人口密集,外来人口聚集的乡镇,有的地区甚至连续几个月为聚集性区域。这基本符合手足口病城乡结合部相邻的街道或外来人口多的地区发病较高的特点。例如,2010年5~7月份,高丽营、天竺、后沙裕地区一直为手足口病病例聚集地区,其中高丽营镇4~8月份,后沙峪镇5~10月份均连续5个月为手足口病病例聚集地区。此类地区人群生活、居住的环境卫生条件较差,托幼机构管理水平低,人口稠密,居民手足口病等传染病防控知识缺乏,这些都为手足口病的传播和流行提供了条件。3.带模型的7d临床资料分析采用MapInfo软件可以直观地显示Flexible空间扫描统计量分析结果,图中阴影所示的区域即为研究所确定的可能的手足口病发病的空间聚集区域(以2010年3、4月份分析结果为例,见图2~3)。flusgerf空间聚集性探索目前,传染病空间聚集性的研究是全球公共卫生领域的热点问题之一。Kulldorf等于1997年提出圆形扫描统计量分析,该方法采用大小和位置动态变化的圆形扫描窗口进行聚集性探测。该方法探测到的聚集性地理区域形状主要是圆形或者近似圆形。假如某传染病沿河流或交通路线传播,则聚集性的形态是狭长的条状,这种情况下,采用圆形扫描窗口就不能很好地对聚集性进行探测。2005年Tango等在Kulldorff圆形扫描统计量基础上提出Flexible空间扫描统计量方法,主要是为了在空间维度上探测圆形扫描窗口探测不到的不规则形状的聚集区域,通过对疾病未知的空间聚集性进行探索以达到早期预警的目的。与Kulldorff等提出的圆形扫描统计量不同,Flexible空间扫描统计量预先不设定扫描窗口的形状,而是以动态变化的不规则型空间扫描窗口对在地理边界上有连接的不同区域集进行扫描,同时该方法将扫描区域限制在每个起始区域周边的一个较小的邻域内。有资料显示,在模拟数据的聚集性探测中,对圆形或近似圆形的聚集区域,Kulldorff圆形扫描方法有稍高的检验效能;而在非圆形的聚集区域时,Flexible扫描统计量的检验效能优于Kulldorff圆形扫描统计量,并且探测出的聚集区域更符合实际的情况。目前,扫描统计量方法已广泛应用于许多医学研究领域,用于疾病空间聚集性的探测以及早期预警。同时该方法还可以实现对聚集性的定位。因该方法采用的概率模型对不同区域间非均匀的人口密度进行了校正以消除其潜在的影响,且采用大小可变的扫描窗口等特点又使其检验效能较其他分析方法显著提高,成为当前疾病时间聚集性或空间聚集性分析研究中的热点。而且,扫描统计量方法在分析之前对聚集性的大