扩展功能
文章信息
- 谢忠杭, 詹美蓉, 欧剑鸣, 洪荣涛, 陈武
- XIE Zhong-hang, ZHAN Mei-rong, OU Jian-ming, HONG Rong-tao, CHEN Wu
- 2010-2015年福建省戊型肝炎空间聚集性及变化规律研究
- Spatial clustering and changing of hepatitis E in Fujian province, 2010-2015
- 疾病监测, 2017, 32(5): 377-381
- Disease Surveillance, 2017, 32(5): 377-381
- 10.3784/j.issn.1003-9961.2017.05.007
-
文章历史
- 收稿日期:2016-09-02
戊型肝炎(戊肝)是由戊肝病毒引起的一种人兽共患病,是急性肝炎的常见病因,每年急性戊肝的发病人数超过300万,死亡人数约7万[1]。戊肝临床症状类似甲型肝炎,但在孕妇、老年人等感染者中有较高( > 20%)的病死率[2];该病经肠道传播为主,易引起大规模暴发或流行,1986-1988年在我国南疆地区因戊肝病毒污染饮用水发生了世界上最大的一次戊肝流行[3]。传染病的发生、传播、流行与地理环境关系密切,传染病疫情数据间往往存在空间自相关性,研究传染病空间分布特征时采用空间分析方法要优于传统数理统计方法,因此近年空间分析方法被广泛运用于戊肝等传染性疾病的研究[4-10]。本研究对福建省2010-2015年戊肝监测资料进行空间自相关分析,旨在揭示福建省戊肝疫情在县域层面上的空间分布规律,以识别防控工作的重点区域,为聚集区域的进一步研究提供依据,同时为制定防控策略提供参考。
1 资料与方法 1.1 资料来源福建省2010-2015年戊肝监测资料来源于中国疾病预防控制信息系统中福建省2010-2015年戊肝报告卡(以现住址、发病日期统计,仅包括临床诊断和确诊病例的已终审卡,但不包括已删除卡)。人口资料来源于国家统计局,研究区域包括福建省10个区市中的88个县(市、区),福建省矢量化县界地图由中国疾病预防控制中心提供。
1.2 空间自相关分析理论空间分析方法是研究疾病流行空间聚集性的有效方法,通过空间分析模型,可以清晰地了解疾病发生的热点地区,探索疾病聚集的位置和范围[11]。其中,空间自相关分析是研究某空间单元的属性值是否相似于其邻近空间单元的属性值,包括全局空间和局部空间自相关分析,属于探索性空间数据分析技术。
1.2.1 全局空间自相关分析全局空间自相关分析可以从整体上反映某空间单元属性值的空间分布情况,判断其是否具有聚集性,常用的统计量有Moran′s I、Geary′s C、Getis和Join Count等。本研究选用全局Moran′s I统计量来评估福建省戊肝发病率在县域空间层面上所表达的模式是聚类模式、离散模式还是随机模式,其计算公式[12]:
式中n是空间单元总数,在本研究中为福建省88个县(市、区);wij是空间单元i和j的空间权重,在本研究中采用基于Queen原则的二进制空间邻接矩阵,即若两个空间单元有公共边界或公共点,则wij=1,否则wij=0;xi、xj分别代表空间单元i和j的属性值,在本研究中为某年某县(市、区)的戊肝发病率。Moran′s I的取值范围在-1~1之间,以标准化统计量Z检验其统计学意义。|Z|≥1.96,则认为其空间自相关性有统计学意义。I > 0,则提示存在正的空间自相关,研究对象的空间分布呈聚类模式,|I|值越大,空间聚集性也越强[13];I < 0,提示存在负的空间自相关,研究对象的空间分布呈离散模式,|I|值越大,空间差异性也越大;I=0,提示不存在空间自相关,研究对象的空间分布呈随机模式,无空间聚集性。
1.2.2 局部空间自相关分析局部空间自相关分析可以探测出具体的自相关区域,常用的方法有局部Moran′s I、Moran散点图和局部Getis指数(Gi*)等。本研究使用Gi*统计量识别具有统计显著性的热点和冷点,其计算公式[13]:
式中n是空间单元总数,xj是空间单元j的属性值,wij是空间单元i和j的空间权重。Gi*值为正,提示该区域为一个高值的空间聚集,称为热点(Hot Spot);Gi*值为负,提示该区域为一个低值的空间聚集,称为冷点。|Gi*|=3的要素反映置信度为99%的统计显著性,|Gi*|=2的要素反映置信度为95%的统计显著性,|Gi*|=1的要素反映置信度为90%的统计显著性,Gi*=0时要素的聚类没有统计学意义。
1.3 统计分析使用SAS 9.2软件建立福建省戊肝病例数据库,计算分年度、分区(县)的发病率,比较戊肝疫情特征在热点区域与非热点区域间的差异性;用χ2检验进行率/构成比的比较;用Cochran-Armitage Test for Trend进行趋势检验。使用ArcGIS 10.2软件绘制县域发病率地图、进行空间自相关分析。P < 0.05为差异有统计学意义。
2 结果 2.1 发病概况福建省2010-2015年共报告戊肝6 367例,各年度发病率依次为2.42/10万、3.26/10万、3.36/10万、3.24/10万、2.34/10万和2.45/10万,发病率总体呈下降趋势(Z=-3.827,P < 0.001),年均为2.85/10万,无聚集性疫情报告。
2.2 全局空间自相关分析由表 1可见,各年度戊肝发病率的空间分布均呈聚类模式:2010-2013年呈中等程度空间聚集(Moran′s I值介于0.329~0.498之间),2013-2015年空间聚集程度呈减弱态势(2015年Moran′s I值为0.124)。
年份 | Moran′s I值 | E(I)值 | Var(I)值 | Z值 | P值 |
2010 | 0.357 2 | -0.011 5 | 0.003 7 | 6.082 7 | < 0.000 1 |
2011 | 0.410 2 | -0.011 5 | 0.003 5 | 7.079 7 | < 0.000 1 |
2012 | 0.329 4 | -0.011 5 | 0.003 3 | 5.935 4 | < 0.000 1 |
2013 | 0.498 2 | -0.011 5 | 0.004 0 | 8.014 4 | < 0.000 1 |
2014 | 0.216 3 | -0.011 5 | 0.003 9 | 3.653 9 | 0.000 3 |
2015 | 0.124 0 | -0.011 5 | 0.003 5 | 2.304 3 | 0.021 2 |
由图 1、表 2可见,2010-2015年分别有8、8、6、11、8、9个热点区域,分别有0、0、0、1、1、1个冷点区域;2010 2015年间的热点区域,90.00%(45/50) 为三明市的县(区);三明市的县(区)中,在2010-2015年内成为热点区域的县(区)数占三明市所辖县(区)总数的83.33%(10/12)。
年份 | 热点 | 冷点 | |||
Gi*=3 | Gi*=3 | Gi*=1 | Gi*=-1 | ||
2010 | 梅列1,明溪1, 沙县1, 永安1, 三元1 | 无 | 宁化1, 将乐1, 建宁1 | 无 | |
2011 | 梅列1, 明溪1, 沙县1, 永安1, 三元1, 宁化1 | 将乐1 | 建宁1 | 无 | |
2012 | 梅列1, 明溪1, 沙县1, 永安1, 三元1 | 宁化1 | 无 | 无 | |
2013 | 梅列1, 明溪1, 沙县1, 永安1, 三元1, 将乐1, 泰宁1 | 宁化1, 建宁1 | 罗源2, 大田1 | 长泰3 | |
2014 | 梅列1, 明溪1, 沙县1, 永安1 | 大田1, 福安4 | 三元1, 霞浦4 | 龙海3 | |
2015 | 梅列1, 明溪1, 沙县1, 永安1, 大田1 | 建宁1, 新罗5 | 将乐1, 华安3 | 长泰3 | |
注:“1”表示该县(区)归三明市所辖;“2”表示该县(区)归福州市所辖;“3”表示该县(区)归漳州市所辖;“4”表示该县(区)归宁德市所辖;“5”表示该县(区)归龙岩市所辖。 |
热点区域(Gi*≥1) 在6年间累计发病853例、发病率为7.40/10万,非热点区域(Gi* < 1) 累计发病5 514例、发病率为2.60/10万,发病率在热点区域与非热点区域间的差异有统计学意义(χ2=884.785,P < 0.001)。确诊病例所占构成比,在热点区域(为94.37%,805/853) 与非热点区域(为92.05%,5 861/6 367) 间差异有统计学意义(χ2=5.715,P=0.017)。热点区域发病月份分布总体呈以2月为主峰、7月为次峰的双峰分布,非热点区域总体呈1-4月为高发期、5-12月为低发期的双平台分布。在热点区域中城市病例占56.04%(478/853),与非热点区域(30.67%,1 691/5 514) 间的差异有统计学意义(χ2=211.687,P < 0.001)。热点区域中男女性别比为2.29 : 1,非热点区域为2.50 : 1,两者间差异无统计学意义(χ2=1.120,P=0.290)。热点区域发病的年龄以45~岁组(占15.36%,131/853) 为峰值的单峰分布,非热点区域以55~岁组(占12.62%,696/5 514) 为峰值的单峰分布,且其峰值相对较低、波峰相对平缓。热点区域前5位高发职业依次为农民(占35.99%,307/853)、离退人员(占15.71%,134/853)、工人(占14.30%,122/853)、家务及待业(占10.20%,87/853) 和干部职员(占7.27%,62/853),非热点区域中这5个职业的发病构成比分别为37.25%(2 054/5 514)、7.02%(387/5 514)、3.59%(198/5 514)、19.91%(1 098/5 514)、3.72%(205/5 514),相比热点区域离退人员、工人和干部职员的发病构成比显著较高(χ2分别为74.263、177.570、23.181,均P < 0.001),家务及待业的发病构成比显著较低(χ2=46.015,P < 0.001),农民的发病构成比差异无统计学意义(χ2=0.503,P=0.478)。
3 讨论本研究揭示了福建省戊肝疫情在县域层面上的空间分布规律:(1) 存在空间聚集性。(2) 发病热点区域具有集中性和连续性。2010-2015年,90.00%(45/50) 的热点区域为三明市的县(区);三明市的梅列区、明溪县、沙县、永安市连续6年成为热点区域,三元区连续5年、宁化县连续4年、大田县连续3年成为热点区域;将乐县、建宁县有4年成为热点区域。(3) 发病热点区域呈现出扩散趋势。2010-2012年,热点区域集中在三明市的8个县(区);2013-2015年,热点区域涵盖了三明市的10个县(区),且包括了其他4个设区市的5个县(区)。
福建省戊肝空间分布的聚集性,与湖南省相比有明显差异[4]:空间集聚程度福建省呈减弱态势,而湖南省呈增强态势,热点区域分布福建省呈现出集中性、连续性和扩散性,而湖南省呈现出交替性和阶段性,这除了受地区间报告质量的差异影响外,还与下述因素有关。(1) 区域间人群的戊肝检测比例不同。我国2009年颁布的《食品安全法实施条例》要求食品从业人员每年均需进行戊肝检测,而三明市近年小型餐饮业发展迅速,辖区居民餐饮类从业人员比例相对较高,相应的戊肝检测比例亦较高,这可能与热点区域的形成有关,以致热点区域的戊肝发病呈现出相应的季节性(食品行业人员通常在每年2-3月或7-8月进行从业健康体检)和高发人群相对年轻等特征。(2) 区域间疑似病例的确诊比例不同。目前在我国的临床实践中,戊肝诊断率与实际感染率存在不符[14],漏、误诊现象不可避免,而“三明医改模式”可能有助于提高当地医生确诊意识和疑似病例的确诊比例,以致热点区域的戊肝确诊病例构成比显著高于非热点区域。(3) 其他影响因素不同。热点区域中离退人员、干部职员以及工人的占比显著高于非热点区域;热点区域中戊肝发病以城市病例为主(占56.04%),显著高于非热点区域的30.67%,提示热点区域的形成,难以简单用饮用水污染[15]和/或农作物污染[16]增加戊肝传染风险来解读,提示尚存在其他影响因素。
综上所述,福建省戊肝发病存在空间聚集性,且空间聚集现象具有集中性、连续性和扩散趋势,三明市是防控的重点区域,应进一步深入调研以加强其防控。
志谢: 本研究得到福建省各级医疗卫生机构及中国疾病预防控制中心的大力支持,谨致谢意。作者贡献:
谢忠杭 ORCID:0000-0002-6432-1386
谢忠杭:负责数据整理、分析与论文撰写
詹美蓉:负责数据收集、文章的构思与修改
欧剑鸣、洪荣涛、陈武:研究设计与指导,文章修改
[1] | Rein DB, Stevens GA, Theaker J, et al. The global burden of hepatitis E virus genotypes 1 and 2 in 2005[J]. Hepatology, 2012, 55: 988–997. DOI:10.1002/hep.25505 |
[2] | Vivek R, Nihal L, Illiayaraja J, et al. Investigation of an epidemic of hepatitis E in nellore in south India[J]. Trop Med Int Health, 2010, 15(11): 1333–1339. DOI:10.1111/tmi.2010.15.issue-11 |
[3] | Meng XJ. Recent advances in hepatitis E virus[J]. J Viral Hepat, 2010, 17(3): 153–161. DOI:10.1111/jvh.2010.17.issue-3 |
[4] | Liu Y, Liang WJ, Li JH, et al. Characteristic of spatial-temporal distribution of hepatitis E in Hunan province, 2006-2014[J]. Chinese Journal of Epidemiology, 2016, 37(4): 543–547. (in Chinese) 刘意, 梁维君, 李俊华, 等. 湖南省2006-2014年戊型肝炎时空分布特征[J]. 中华流行病学杂志, 2016, 37(4): 543–547. |
[5] | Yu YN, Lyu HJ. The spatial autocorrelation of hepatitis E[J]. Chinese Journal of Health Statistics, 2014, 31(1): 113–115. (in Chinese) 俞以内, 吕惠进. 戊型病毒性肝炎的空间自相关分析[J]. 中国卫生统计, 2014, 31(1): 113–115. |
[6] | Zheng YH, Liu XF, Ren XW, et al. Spatial distribution of hepatitis A in Gansu, 2013[J]. Disease Surveillance, 2015, 30(8): 624–628. (in Chinese) 郑芸鹤, 刘新凤, 任晓卫, 等. 2013年甘肃省甲型病毒性肝炎发病的空间聚集性分析[J]. 疾病监测, 2015, 30(8): 624–628. DOI:10.3784/j.issn.1003-9961.2015.08.005 |
[7] | Tang L, Zhang Y, Xing DG, et al. Spatial data based study of distribution of hepatitis C in Chongqing[J]. Chinese Journal of Epidemiology, 2016, 37(1): 80–84. (in Chinese) 唐路, 张燕, 幸奠国, 等. 基于空间数据分析技术的重庆市丙型肝炎发病研究[J]. 中华流行病学杂志, 2016, 37(1): 80–84. |
[8] | Sun HQ, Xiao GX, Guo Y, et al. Spatial clustering and changing trend of hepatitis C prevalence during 2008-2012 in China[J]. Chinese Journal of Public Health, 2014, 30(3): 286–289. (in Chinese) 孙海泉, 肖革新, 郭莹, 等. 中国2008-2012年丙肝流行规律及空间聚集性分析[J]. 中国公共卫生, 2014, 30(3): 286–289. DOI:10.11847/zgggws2014-30-03-10 |
[9] | Xiao GX, Hu YH, Ma JQ, et al. Spatial clustering and changing trend of hand-foot-mouth disease during 2008-2011 in China[J]. Chinese Journal of Epidemiology, 2012, 33(8): 808–812. (in Chinese) 肖革新, 胡跃华, 马家奇, 等. 中国2008-2011年手足口病空间聚集性及变化趋势研究[J]. 中华流行病学杂志, 2012, 33(8): 808–812. |
[10] | Gou FX, Liu XF, Liu DP, et al. Spatial-temporal distribution of hepatitis B in Gansu province, 2009-2014[J]. Chinese Journal of Epidemiology, 2016, 37(1): 85–89. (in Chinese) 苟发香, 刘新凤, 刘东鹏, 等. 甘肃省2009-2014年乙型肝炎报告发病率的时空聚集性研究[J]. 中华流行病学杂志, 2016, 37(1): 85–89. |
[11] | Wang T, Xue FZ, Chen YJ, et al. The spatial epidemiology of tuberculosis in Linyi city, China, 2005-2010[J]. BMC Public Health, 2012, 12: 885. DOI:10.1186/1471-2458-12-885 |
[12] | Liu G, Bi RT, Wang SJ, et al. The use of spatial autocorrelation analysis to identify PAHs pollution hot spots at an industrially contaminated site[J]. Environ Monit Assess, 2013, 185(11): 9549–9558. DOI:10.1007/s10661-013-3272-6 |
[13] | Amaral PV, Anselin L. Finite sample properties of Moran's I test for spatial autocorrelation in tobit models[J]. Papers Reg Sci, 2014, 93(4): 773–781. DOI:10.1111/pirs.12034 |
[14] | Zhang Q, Tang GP, Liu HH, et al. Investigation on the causes of increased hepatitis E cases reported in Guizhou province[J]. Chinese Journal of Epidemiology, 2015, 36(3): 228–231. (in Chinese) 张倩, 唐光鹏, 刘慧慧, 等. 贵州省戊型肝炎报告病例增加原因的调查[J]. 中华流行病学杂志, 2015, 36(3): 228–231. |
[15] | Teshale EH, Denniston MM, Drobeniuc J, et al. Decline in hepatitis E virus antibody prevalence in the united states from 1988-1994 to 2009-2010[J]. J Infect Dis, 2015, 211(3): 366–373. DOI:10.1093/infdis/jiu466 |
[16] | Pavio N, Meng XJ, Doceul V. Zoonotic origin of hepatitis E[J]. Curr Opin Virol, 2015, 10: 34–41. DOI:10.1016/j.coviro.2014.12.006 |