秩和比法在不同滅鼠劑滅鼠前后鳥類數量變化比較中的應用

戰揚孫愛峰

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-），男，吉林長白縣人，本科，副主任醫師，主要從事疾病預防控制和管理工作。