秩和比法在不同灭鼠剂灭鼠前后鸟类数量变化比较中的应用

战扬孙爱峰

1.长春市二道区疾病预防控制中心,吉林长春130031；2.白城卫生职工中等专业学校，吉林白城137000

秩和比法在不同灭鼠剂灭鼠前后鸟类数量变化比较中的应用

战扬1孙爱峰2

1.长春市二道区疾病预防控制中心,吉林长春130031；2.白城卫生职工中等专业学校，吉林白城137000

目的应用秩和比法比较不同灭鼠剂灭鼠前后鸟类数量变化情况。方法在灭鼠区域设置6条样带，3年内分别采用三种毒饵灭鼠，在灭鼠前后15 d按既定方法观察样带内鸟类数量并计算鸟类密度，应用基于秩和比的可信区间法进行统计分析。结果大隆灭鼠前的RSR及其变量变换值的95%CI分别为0.9167和58.33～88.11，灭鼠后分别为0.2083和12.26～42.04；敌鼠钠盐灭鼠前分别为0.7500和45.11～74.89，灭鼠后分别为0.3056和18.67～48.45；溴敌隆灭鼠前分别为0.6667和39.85～69.63，灭鼠后分别为0.4167和25.31～55.09，除大隆灭鼠前与敌鼠钠盐灭鼠前、敌鼠钠盐灭鼠前与溴敌隆灭鼠前、溴敌隆灭鼠前与溴敌隆灭鼠后、溴敌隆灭鼠后与敌鼠钠盐灭鼠后和大隆灭鼠后、敌鼠钠盐灭鼠后与大隆灭鼠后两两组别鸟类密度之间的差异无统计学意义（P＞0.05）外，其余两两组别鸟类密度之间的差异均有统计学意义（P＜0.05）。结论单纯从鸟类密度看，应用溴敌隆在保持生态平衡，维持环境质量方面具有重要的意义，应选择溴敌隆作为灭鼠剂。

灭鼠；鸟类密度；秩和比；可信区间

吉林省西北部草原是松辽平原达乌尔黄鼠鼠疫自然疫源地，灭鼠是防制鼠疫的主要措施。但是误食毒饵可对包括鸟类在内的非靶动物带来一定的危害，影响生态系统的平衡。为评价不同毒饵对鸟类数量变化的影响，笔者应用秩和比法对三种不同灭鼠剂使用前后草原鸟类密度进行研究，探索如何选择既能降低草原鼠类危害，又能维持生态平衡的灭鼠药物。现将结果报道如下。

1 资料与方法

1.1 一般资料

2010年4月—2012年5月，在灭鼠区域设置6条样带（长3000 m，宽100 m），分别采用0.01%大隆向日葵玉米粒混合毒饵（简称大隆）、0.05%敌鼠钠盐向日葵玉米粒混合毒饵（简称敌鼠钠盐）和0.01%溴敌隆向日葵玉米粒混合毒饵（简称溴敌隆）灭鼠，于灭鼠前后15 d按照既定方法观察样带内鸟类数量，所得灭鼠前后鸟类密度观测值见表1[1]。数据来自文献，真实可靠。

表1 三种灭鼠剂灭鼠前后鸟类密度（只/km2）及其秩和比计算表

1.2 统计方法

应用基于秩和比的可信区间法[2]分析三种灭鼠剂灭鼠前后鸟类数量变化情况。为减少I类错误的发生，在计算95%CI时，按bonferroni法[3]校正检验水准。设原检验水准为0.05，则校正后的检验水准α=0.00333，故uα=2.94。

2 结果

2.1 秩和比的计算

从生态系统平衡和环境质量角度看，鸟类密度越高越好。因此，将同一样带不同灭鼠剂灭鼠前后鸟类密度观测值按从小到大顺序编秩（见表1实际观测值后括号内的数字），之后按照公式RSR=ΣR/mn（式中R表示秩次，m表示列数量，n表示样带数量）计算各列秩和比，见表1最后一列。

2.2 不同灭鼠剂灭鼠前后鸟类密度的两两比较

按照参数统计的要求，将计算得到的秩和比进行平方根反正弦代换，将计算得到的数值y按照公式y±uαSy（式中Sy=（820.7/N）0.5，N=32为所有格子数量）计算其95%CI，见表2。

表2 不同灭鼠剂灭鼠前后鸟类密度的比较

比较各组y的95%CI可知，除大隆灭鼠前与敌鼠钠盐灭鼠前、敌鼠钠盐灭鼠前与溴敌隆灭鼠前、溴敌隆灭鼠前与溴敌隆灭鼠后、溴敌隆灭鼠后与敌鼠钠盐灭鼠后和大隆灭鼠后、敌鼠钠盐灭鼠后与大隆灭鼠后两两组别鸟类密度之间的差异无统计学意义（P＞0.05）外，其余两两组别鸟类密度之间的差异均有统计学意义（P＜0.05），见表3。

3 讨论

秩和比为反映鸟类数量的一个综合信息统计量，且数值越大，鸟类密度越高。由计算得到的数值可知，鸟类密度的排序为大隆灭鼠前＞敌鼠钠盐灭鼠前＞溴敌隆灭鼠前＞溴敌隆灭鼠后＞敌鼠钠盐灭鼠后＞大隆灭鼠后；两两比较表明，除大隆灭鼠前与敌鼠钠盐灭鼠前、敌鼠钠盐灭鼠前与溴敌隆灭鼠前、溴敌隆灭鼠前与溴敌隆灭鼠后、溴敌隆灭鼠后与敌鼠钠盐灭鼠后和大隆灭鼠后、敌鼠钠盐灭鼠后与大隆灭鼠后两两组别鸟类密度之间的差异无统计学意义（P＞0.05）外，其余两两组别之间的差异均有统计学意义（P＜0.05）。

表3 不同灭鼠剂灭鼠前后鸟类密度的两两比较结果

本文所获得的信息比文献[1]更为丰富，原文只比较了三种灭鼠剂使用前后鸟类密度差异是否具有统计学意义，本文可以进行三种灭鼠剂使用前后鸟类密度两两之间的差异是否具有统计学意义，是对其所得结论的有力补充。与文献[1]的结果比较，主要差别在于溴敌隆灭鼠前后鸟类密度的差异是否具有统计学意义，原因可能与编秩时没有充分利用原始数据提供的丰富信息所致[4]，此时应利用各种编秩技巧[5]，合理分析相关资料。

以上分析可见，应用大隆和敌鼠钠盐灭鼠后鸟类密度显著下降，但是溴敌隆灭鼠后鸟类密度未改变。从灭鼠后鸟类密度的下降幅度看，大隆为62.72%，敌鼠钠盐为59.48%，溴敌隆为35.20%。单纯从鸟类密度看，从既定的灭鼠剂使用前后鸟类数量变化情况分析，应用溴敌隆在保持生态平衡，维持环境质量方面具有重要的意义，应选择溴敌隆作为灭鼠剂。

但是，在选择灭鼠剂时应综合考虑各方面的因素。如尽管溴敌隆对非靶动物毒性较低，但投药量较大，毒饵集中。再如本文所选资料为3年的监测资料，是否考虑了前面所用灭鼠剂对后续年度鸟类密度的影响，对分析结论有一定的影响；同时还要考虑灭鼠的效果，投药的方便性[6]等种种因素。

应用秩和比法分析不同灭鼠剂使用前后草原鸟类密度，可以为单纯从鸟类密度出发选择既能降低草原鼠类危害，又能维持生态平衡的灭鼠药物提供科学的依据，值得推广应用。

[1]李华宇，董建国，赵亚冰，等.草原大面积灭鼠对草原鸟类数量变化的影响[J].中国地方病防治杂志，2013，28（2）：114-115.

[2]孙爱峰.秩和比法在不同地区人体蛔虫感染现状调查结果综合评价中的应用[J].中国医药指南，2011，9（31）：226-227.

[3]方积乾主编.医学统计学与电脑实验[M].4版.上海：上海科学技术出版社，2012.

[4]史景明，孙爱峰.秩和比法在不同年度食品卫生监督质量综合评价和变化发展趋势中的应用[J].中国医药指南，2013，11（29）：594-595.

[5]赵树立，孙爱峰.秩和比法在学校食堂卫生监督质量综合评价中的应用[J].科技创新导报，2013（34）：221.

[6]王峰，任清明，刘国平，等.吉林西北部草原大面积灭鼠效果观察[J].中华卫生杀虫药械，2007，13（5）：355-356.

表3 消除碘缺乏病目标考核评估结果的因素分析

3 讨论

在本文的统计分析过程中，所以对原观测值进行相对化处理，目的在于将原始数据还原为百分制以进行混合编秩，为进行消除碘缺乏病目标考核评估结果的因素分析提供前提条件；同时，可以计算权重系数，充分利用原始数据提供的信息并体现各类指标的作用，避免经验权重对不同县市区实现消除碘缺乏病目标抽查评估结果的影响。

应用基于秩和比的可信区间法对某省实现消除碘缺乏病目标考核评估结果的因素分析表明，健康教育为主要影响因素，其次为组织管理；应用基于秩和比的分档法综合评价表明，由低到高对的县排序和等级为8、3、18、12、15和7属于差等级，17、14、6、10、13和2属于良好等级，11、1、16、9、4和5属于优秀等级，为进一步制订防治规划提供科学的依据。

原文只是以各项合计得分和最高分、最低分和全省平均分进行统计分析，本文在对不同项目县进行排序的同时，还分别利用基于秩和比的分档法和可信区间法进行聚类分析和因素分析，这有利于管理部门指导从全省宏观角度找出影响评估结果的因素，客观对项目县评估结果进行分类，获得比原文更为丰富的信息，所得结论是对原结论的补充和完善，值得推广应用。

值得注意的是，应用秩和比法所做因素分析和聚类分析与数理统计中的因子分析和聚类分析的原理[5]不同，与基于秩和比的因素分析[6]的原理亦不相同，应注意合理应用。

[参考文献]

[1]郑合明，李小烽，王传刚，等.河南省实现县级消除碘缺乏病目标考核评估结果分析[J].中国地方病防治杂志，2013，28（2）：127-128，135.

[2]孙爱峰.秩和比法在不同地区人体蛔虫感染现状调查结果综合评价中的应用[J].中国医药指南，2011，9（31）：226-227.

[3]陈雪侠，孙爱峰.秩和比法在医院医疗质量综合评价及其变化发展趋势中的应用[J].中国医药指南，2014，12（8）：249-250.

[4]方积乾.卫生统计学[M].7版.北京：人民卫生出版社，2012.

[5]方积乾.医学统计学与电脑实验[M].4版.上海：上海科学技术出版社，2012.

[6]高淑艳，孙爱峰.基于秩和比的和谐分析有关问题的探讨[J].中国卫生统计，2008，25（4）：436，442.

（收稿日期：2014-07-02）

R195

A

1672-5654（2014）09（b）-0044-02

2014-07-02）

战扬（1968-），男，吉林长白县人，本科，副主任医师，主要从事疾病预防控制和管理工作。