基于農(nóng)藥毒力計算的6個統(tǒng)計分析軟件比較

常菊花　何月平

摘要 選用了能夠用于農(nóng)藥毒力數(shù)據(jù)機率值分析的通用軟件（SAS、SPSS和DPS）和專門程序（Polo、BA和EPA機率值分析程序），以二化螟對三唑磷的室內(nèi)毒力數(shù)據(jù)為例，比較分析不同統(tǒng)計分析軟件的計算結(jié)果。結(jié)果表明，當(dāng)對照組死亡率為0時，6個軟件程序計算得到的LD值（LD5、LD50和LD95）、截距和斜率等基本相等，但是LD值的95%置信限有差異。當(dāng)對照組死亡率不為0時，發(fā)現(xiàn)6個程序計算得到的毒力資料數(shù)值都有差異，其中LD（LD5、LD50和LD95）、截距和斜率的數(shù)值相差較小，但是LD值的95%置信限差異較大。

關(guān)鍵詞 SAS；SPSS；Polo；EPA；DPS；BA

中圖分類號 S127 文獻標(biāo)識碼 A 文章編號 0517-6611（2014）03-00746-03

Abstract Six statistical analysis software （SAS， SPSS， DPS， Polo， BA and EPA） were used to calculate the pesticide toxicology data. When the mortality in control group was 0， the Lethal dose（LD5， LD50 and LD95）， intercept and slope were almost equal， but the 95% confidence limits calculated by the six statistical analysis software were different. When the mortality in control group was not 0， the pesticide toxicology data calculated by six statistical analysis software （SAS， SPSS， DPS， Polo， BA and EPA） were all different. There were slight differences between the Lethal dose （LD5， LD50 and LD95）， intercept and slope， but there were large differences between the 95% confidence limits calculated by the six statistical analysis software.

Key words SAS； SPSS； Polo； EPA； DPS； BA

在農(nóng)藥生物測定中，致死中量（LD50）或有效中量（ED50）和回歸線的斜率（b值）是表示藥劑對生物效力的代表性數(shù)值。為了要測定LD50或b值，通常采用的統(tǒng)計分析方法是機率值分析法（Probit Analysis）[1]。機率值分析法和計算致死中量的方法很多，但傳統(tǒng)的手工計算方法十分復(fù)雜，不僅要花費大量時間，而且很容易出錯。隨著電子計算機技術(shù)的發(fā)展，出現(xiàn)了不少軟件或程序應(yīng)用于生物測定數(shù)據(jù)的機率值分析和致死中量的計算，在很大程度上簡化了運算程序[2-6]。筆者采用機率值分析軟件（Polo、BA、EPA）和生物統(tǒng)計軟件（SPSS、SAS、DPS）計算農(nóng)藥毒力數(shù)據(jù)，比較分析不同程序的計算結(jié)果，擬推薦出可供廣泛使用且結(jié)果精確的程序或軟件用于生物測定數(shù)據(jù)的統(tǒng)計分析。

1 材料與方法

1.1 程序和測試數(shù)據(jù)

以室內(nèi)采用毛細管點滴法測定三唑磷對水稻二化螟幼蟲的室內(nèi)毒力數(shù)據(jù)為例（表1），采用SAS 8.1、SPSS 19.0、DPS 7.05版、PoloPlus （Probit and Logit Analysis， LeOra Software）、EPA （EPA probit analysis program used for calculation LC/EC values version 1.5）、BA（Bioassay Data Processing and Management System）6個統(tǒng)計分析軟件，計算毒力回歸方程、LD5、LD50、LD95 及其95%置信限，以及卡平方值和相關(guān)系數(shù)等，比較分析不同程序的計算結(jié)果。

1.2 軟件使用和結(jié)果輸出

各軟件用于機率值分析的操作步驟如下，在表2～3中列出6個軟件計算所得的LD5、LD50和LD95及其95%置信限、回歸方程的截距、斜率及其標(biāo)準(zhǔn)誤、卡平方值或相關(guān)系數(shù)等，所有數(shù)值保留小數(shù)點后3位有效數(shù)字。

1.2.1 Polo軟件（Probit and logit Analysis）。Polo軟件是LeOra軟件公司研發(fā)的專門用于生物測定資料的統(tǒng)計分析軟件，其開發(fā)是基于BLISS算法，并包含parallelism和equality假設(shè)測驗[6]。PoloPlus版本是與微軟windows兼容的最新版Polo程序。在PoloPlus程序中打開保存的毒力數(shù)據(jù)文件，在參數(shù)選擇“Choose options”窗口中，選擇劑量的log轉(zhuǎn)換“Convert doses to logarithms （standard）”，在數(shù)學(xué)模型“Mathematical Model”中選擇機率值“Probit”，輸入所要獲得的致死濃度（如LD5、LD50、LD95）。完成數(shù)據(jù)格式檢查“Check data”和計算操作后，PoloPlus會在“Display results”和“Display summary”界面中顯示數(shù)據(jù)分析結(jié)果。這個程序能夠計算顯示出毒力回歸資料的斜率和截距及其標(biāo)準(zhǔn)誤、致死濃度（LD1～LD99）及其95%置信限和卡方檢驗結(jié)果。毒力數(shù)據(jù)在TXT或DOC文本文件中按照如下格式輸入，濃度為表1中劑量（mg/ml）乘以點滴量0.04 μl/頭。

2.2 對照死亡率不為0時的程序計算結(jié)果比較

當(dāng)對照死亡率不為0時，使用Polo、DPS、EPA和BA等程序時，直接將對照死亡個數(shù)輸入到程序中。采用SAS程序編碼進行分析之前，先將死亡率進行校正后再輸入到軟件中進行計算，經(jīng)校正后的死亡蟲數(shù)保留小數(shù)點后2位有效數(shù)值，如0.4 μg/頭的劑量下的校正死亡蟲數(shù)為11.38；使用SPSS軟件時，表3中的SPSS-1計算處理組是先將校正死亡蟲數(shù)輸入到SPSS中（校正后的死亡蟲數(shù)保留小數(shù)點后2位有效數(shù)值），再按照對照死亡率為0來進行計算分析；SPSS-2計算處理組為采用軟件帶有的“Natural response rate”選擇用于設(shè)定對照組的自然死亡率（如表1中的對照死亡率為0.033）。

6個統(tǒng)計軟件計算得到的結(jié)果見表3。SPSS-1和SPSS-2兩組得到的結(jié)果有差異，SPSS-1和SAS計算結(jié)果基本一致，SPSS-2和Polo計算結(jié)果基本一致。另外，①LD值（LD5、LD50和LD95）：BA、SAS和SPSS-1計算得到的LD值相等或相近，Polo、DPS、SPSS-2和EPA得到的LD值相近，如BA、SAS和SPSS-1計算得到的LD50值（0.806）要高于Polo、DPS、SPSS-2和EPA的LD50值（0.797～0.801）。②LD值的95%置信限：BA、DPS和EPA計算得到的LD值得95%置信限跨度要窄于POLO、SAS、SPSS-1和SPSS-2得到的置信限。③斜率和截距：SAS、BA和SPSS-1得到的斜率（1.939～1.940）要略高于Polo、DPS和SPSS-2得到的斜率（1.926），EPA算得斜率為1.923。SAS和SPSS算得的截距數(shù)值加上數(shù)值5后就與其他程序得到的結(jié)果一致。④卡平方值或相關(guān)系數(shù)：EPA算得相關(guān)系數(shù)值最?。?.881）。BA和DPS得到的相關(guān)系數(shù)分別為0.959和0.968。