论文部分内容阅读
目的了解阿坝州棘球蚴病流行特征,阐明棘球蚴病流行程度和范围,探索棘球蚴病流行影响因素以及有效的棘球蚴病插值预测方法。方法：①采用分层整群等比例随机抽样方法,以阿坝州牧区、半农半牧、农区和城镇各层人口数为基础做等比例抽样,采用直接标准化法计算全州各县标准化患病率；运用PearsonX2检验、Cochran-Mantel-Haensze检验、Fisher’Exact Test以及Cochran-Armitage趋势检验等方法对人群棘球蚴病患病率以及性别、年龄、民族、文化程度、职业、宗教信仰、居住方式等不同属性分组的患病率做差异性检验,分析阿坝州人体棘球蚴病分布特征。②应用ArcGIS9.3软件进行专题制图,以直观显示人体棘球蚴患病率、家畜患病率以及犬的感染率的地区分布图。③应用直方图法、正态QQ Plot图法进行数据分布检验；应用Voronoi图来查找离群值；三维趋势分析识别数据的全局趋势；变异函数分析探讨其是否具有空间自相关性以及适合拟合何种变异函数模型,进而选择适宜Kriging插值方法。④对阿坝州人群棘球蚴病患病率进行Indicator Kriging插值预测分析；采用交叉验证法对插值结果的精度进行评价；运用全域Moran’s I、Getis法进行空间聚集性分析。⑤运用Spearman等级相关法和Cochran-Mantel-Haensze检验分析各因素与人群棘球蚴病患病率之间的相关性。⑥建立负二项回归模型进行棘球蚴病流行影响因素分析,运用Pearsonx2检验、信息准则(AIC)以及Schwarz贝叶斯准则(BIC)确定模型的拟合优度。结果①地区分布：若尔盖县的标准化患病率最高为2.786%(1159/41602),其次壤塘、马尔康、红原、阿坝、金川等5县标准化患病率依次为1.923%(800/41602)、1.832%(762/41602)、1.543%(642/41602)、1.219%(507/41602)和1.017%(423/41602),汶川的标化患病率最低为0.029%(12/41602)。②人群分布：阿坝州棘球蚴病患病率为0.885%(368/41602)；男、女患病率分别为0.69%(130/18948)和1.05%(238/22654),女性高于男性(产15.64,P<0.001)；棘球蚴病患病率随年龄增长呈升高趋势(Z=9.4583,P<0.001),随教育程度提高相应下降(Z=11.5902,P<0.001)；藏族、回族、汉族、羌族患病率分别为1.227%(352/28692)、0.259%(1/386)、0.177%(12/6786)、0.052%(3/5735),四个民族患病率差异有统计学意义(χ2=115.407,P<0.001)；牧民与半农半牧人群患病率分别为1.987%(145/7298)和1.364%(63/4618),高于其他职业人群(χ2155.3097,P<0.001)；信仰佛教人群患病率1.394%(304/21803)高于无宗教信仰人群(0.329%)(χ2=132.2495,P<0.001)。③Indicator Kriging插值：棘球蚴病患病率高概率区主要分布于阿坝州西部及西北部地区,九寨沟至金川之间狭长地带为过渡地带,过渡地带东部地区是棘球蚴病低概率区。过渡地带随着阈值水平提高向西偏移。④交叉验证：4个阈值患病率Indicator Kriging预测插值估计偏性较小,估计偏差变异度较小,估计值与实测值的一致性较好。⑤空间聚集性：阿坝州全境内人群棘球蚴病患病率存在高患病率聚集区(Moran’s I和Getis系数分别为0.36和0.06,二者P<0.001)。⑥影响因素：以海拔高度、生产类型为自变量建立的负二项回归模型的拟合优度较好(PearsonX2=1.0670, AIC=690.1650, BIC=710.1902);海拔高度的回归系数为0.364,P<0.001,生产类型以农区为参考,牧区和半农半牧的回归系数分别为0.9398、0.7377,P<0.001；海拔高度的RR值为1.43914,95%CI为(1.20581,1.73007),牧区和半农半牧的RR值为2.55953和2.09121,95%CI分别为(1.46626,4.45492)和(1.29517,3.38347)。结论①阿坝州西部及西北部流行程度较东部及东南部为重；女性、藏族、牧民、信奉佛教以及受教育水平较低人群患病率较高。②负二项回归模型适用于棘球蚴病流行影响因素分析,海拔高度、生产类型是棘球蚴病流行的重要影响因素。③Indicator Kriging插值方法能较好地预测棘球蚴病流行范围和流行程度,可为其它同类地区流行现状评估提供参考。