EXCEL在毒力回归计算中的应用

技术与方法

EXCEL在毒力回归计算中的应用

张志祥　徐汉虹

(华南农业大学昆虫毒理研究室　广州　510642)

程东美

(仲恺农业技术学院植保系　广州　510225)

CalculatingtoxicityregressionwithEXCEL.ZHANGZhi2Xiang,XUHan2Hong(LaboratoryInsectToxicity,SouthChinaAgriculturalUniversity,Guangzhou　510642,China);CHENGDong2Mei(DepartmentofPlantProtection,ZhongkaiAgrotechnicalCollege,Guangzhou　510225,China).

Abstract　Themethodofcalculatingtoxicityregression,LC50,relationcoefficient,SE,95%limiteddistanceofLC50andCTCwithEXCELwasintroduced.Keywords　EXCEL,toxicityregression摘　要　本文介绍了应用EXCEL进行毒力回归分析、计算LC50,a,b,相关系数,以及SE,LC50的95%置信区间和共毒系数的方法。关键词　EXCEL,毒力回归

-1Π2-1Π2Y-5

公式Ⅰ:Z=(2π)e

2

公式Ⅱ:w=ZΠPQ,w是权重系数,P是死亡率,Q=(1-P)。

-2-1

公式Ⅲ:SE(m)=(b×((∑nw)+(m-2

󰁡)×(∑nw(m-xx󰁡))

-1

2

))1Π,SE(m)是m的标

(

)

1　前言

EXCEL是美国微软公司研制开发的用于

个人财务分析与规划、公司营运管理与目标设定和薪水管理等的一套软件。EXCEL具有一个显著的特点,那就是它可以进行公式编辑和插入函数,并且具有连环计算的能力。作者通过使用发现,在单元格中输入公式或插入函数后,EXCEL可以将毒力回归中的所有结果计算出来,并且计算方法简单和快捷。在EXCEL中建立好计算系统后,只需输入试验浓度、供试虫数、死虫数、某一单剂在混剂中的比例,EXCEL就可将LC50,a,b,相关系数、LC50的标准误、LC50的95%置信区间和共毒系数等数据一次性

准误,m是lg(LC50),󰁡x为浓度对数平均值,n为某浓度供试总虫数。

mmm

公式Ⅳ:SE(10)=10(loge10×SE(m),SE(10)为LC50的标准误。

公式10

50

Ⅴ:LC50的95%置信限:

[1]

(lg(LC)±1196×SE(m))

公式Ⅰ～Ⅴ参考张宗炳211　计算系统的制作方法

(1988)。

计算出来。本文将详细介绍这个计算系统的制作过程。

2　计算系统的制作

当EXCEL中公式输入或函数插入完毕后,只要在显色格中输入数据,如图1所示,其它格中的数据立即就会自动显示出来,未显色格中显示出来的是数据,但实际上是所在格公式的计算结果。

现将以上计算系统的编制过程介绍如下,

收稿日期:2000212229,修回日期:2001206218

本文所采用的主要公式如下:

毒力回归方程为Y=a+bx,相关系数为

r。

　·68·昆虫知识　ENTOMOLOGICALKNOWLEDGE200239(1)

图1　用EXCEL进行毒力回归计算机上界面的表格部分下文中的“i”均表示第i行。A1～A5:Ai格中的数据为浓度。

B1～B5:表示浓度对数,Bi是Ai的对数,Bi格中的公式为“=LOG10(Ai)”。

C1～C5:Ci表示i行浓度处理的试验总虫

见公式Ⅰ。

J1～J5:此列为权重系数“w”,Ji格中的公

(1-Fi)”式为“=Ii3IiΠFiΠ,见公式Ⅱ。

K1～K5:此列为权重系数与i行浓度下供

数。

D1～D5:Di表示i行浓度处理的试验死虫数。

E1～E5:此列为死亡率,Ei格中的公式为

试虫数的乘积,即公式Ⅲ中的“nw”,Ki格中的公式为“=Ci3Ji”。K6格中的公式为“=SUM(K1:K5)”,该格为K1～K5格的总和,即公式Ⅲ中的“∑nw”。

L1～L5:此列为公式Ⅲ中的“nw(x-2

󰁡)”x,Li格中的公式为“=Ki3(Bi-AVERAGE(B1:B5))3(Bi-AVERAGE(B1:B5))”,注意在复制公式时保持“AVERAGE(B1:B5)”不变。L6格中的公式为“=SUM(L1:L5)”,该格为L1～

2

L5格的总和,即“∑nw(x-󰁡x)”。

2

M1为“(SE(m))”,其中的公式为“=1ΠD8Π

“=DiΠCi3100”。

F1～F5:此列为校正死亡率,Fi格中的公

(100-E6)”式为“=(Ei-E6)Π,E6为对照死亡

率。因计算死亡率机率值的需要,没有将Fi格的数据转化为死亡百分率,在复制此格的公式时,应注意保持公式中“E6”不变。

G1～G5:此列为校正死亡率的正态等差,Gi格中的公式为“=NORMSINV(Fi)”。“NORM2SINV()”函数本意是为了返回括号中数值(probability)的正态分布区间点,当把校正死亡率放在括号中的时候,得出的结果正好等于校正死亡率的正态等差。

H1～H5:此列为校正死亡率机率值,Hi格中的公式为“=5+Gi”。“离方差+5”等于“死亡率机率值”,所谓机率值就是正分布曲线在横

[2]

坐标上某一点的正态等差加上5。

I1～I5:Ii格中的公式为“=1ΠPOWER(PI()32,1Π2)3EXP(-1Π23(Hi-5)3(Hi-5))”,此列是为了计算权重系数公式中的参数“Z”,

D83(1ΠK6+(B9-AVERAGE(B1:B5))3(B9-(见公式Ⅲ),其中:AVERAGE(B1:B5))ΠL6)”

“D8”为毒力回归方程的b值,即斜率“;B9”为

lg(LC50)。

M2:此格为M1的平方根,即“SE(m)”,其中(见公式Ⅲ)。的公式为“=POWER(M1,1Π2)”

B7:此格为相关系数(r),格中公式为“=CORREL(B1:B5,H1:H5)”。B8:此格为毒力回归方程的截距,即a值,其中的公式为“=INTERCEPT(H1:H5,B1:B5)”。

D8:此格为毒力回归方程的斜率,即b值,其中的公式为“=SLOPE(H1:H5,B1:B5)”。

200239(1)昆虫知识　ENTOMOLOGICALKNOWLEDGE·69·　

B9:此格为lg(LC50),其中的公式为“=(5-B8)ΠD8”。

B10:此格为LC50,其中的公式为“=POW2ER(10,B9)”。

D10:此格为LC50的标准误“SE(10)”,其中的(见公式Ⅳ)。公式为“=B103LN(10)3M2”

F9和G9中的公式分别为“=B9-11963M2”和“=B9+11963M2”,即公式Ⅴ中“lg(LC50)±。1196×SE(m)”

F10和G10分别为LC50的95%置信限的上m

的公式后,单剂A、单剂B和混剂(A+B)的

LC50分别出现在B10,B20和B30处,如在H27处输入单剂A在混剂中所占的比例,在B31处

((1Π＄输入“=(1ΠB30)ΠB＄10)3H27+(1Π＄B

＄20)3(1-H27))3100”,那么共毒系数也可由计算机在一瞬间计算出来。

如果在单剂A、单剂B和混剂的显色处输入同样的数据,并在H27处输入015,那么在B31处立即就会出现“100”这个数据,这样也可以检测共毒系数公式的准确性。3　讨论限和下限,其公式分别为“=POWER(10,F9)”)。和“=POWER(10,G9)”(见公式ⅤI,J,K和L列所有的公式只是为算出M2

打下基础。

上面的公式并不是每个均要人工输入,除M列外,第1行至第5行的公式是可以相互复

制得到的,即只要将第1行的公式输好了,通过“复制”与“粘贴”键可以一下子完成第2至4列

的公式,只是需要注Fi格中的对照死亡率(E6)和Li格中的“AVERAGE(B1:B5)”始终保持不变即可。

第1～10行的公式又可通过“复制”与“粘贴”键复制到指定的地方,只需选定第1至第10行,按下“复制”键,再将光标定在第10行以

目前,可用于毒力回归的工具主要有SAS、

DPS软件和CASIO系列计算器,其中运用SAS和DPS软件进行毒力回归也非常的简单和快速。由于SAS、DPS或其它一些软件统计出来的数据需从计算机中输出,然后再在WORD中制表并输入或输入Hg等制图软件中,这些环节容易出现错误,若数据多了,校对工作量就大。因此,与这些工具相比,EXCEL在进行毒力回归上具有下列不可比拟的优点。311　只需要输入原始数据,即浓度、供试虫数、死虫数和某一单剂在混剂中的比例,EXCEL就可将毒力回归及混剂筛选的所有结果,如r,LC50,a,b,LC50的95%置信区间、LC50的标准误和共毒系数等,一步到位地全部计算出来。这是目前其它软件不具备的,就是DPS软件也不能在一次处理中将相关系数、LC50和LC50的95%置信区间同时计算出来。在计算共毒系数

下的某一格,如A11,按下“粘贴”键后,上述的公式就会瞬间出现在第11～20行,复制完毕后,第11～20行的公式不再需要人为的变动,公式随差行的变动而自动变动,如E1中公式为“=D1ΠC13100”,粘贴后,E11中的公式会自动变为“=D11ΠC113100”。

上面的公式能将毒力回归中的a,b,r,LC50的SE和LC50的95%置信区间等数据一步

到位地计算出来,同样,EXCEL也能将共毒系数自动计算出来,即只需要输入2个单剂和混剂的浓度、总虫数、死虫数和某一单剂在混剂中所占的比例,那么混剂筛选或测定中的所有结果均可求出来。

共毒系数的计算方法介绍如下:

在A11和A21处连续两次复制第1～10行

时,其它软件的不足之处就更明显了,若是单剂的原始数据输入错误,那么所有的共毒系数均需重新算过,而EXCEL就可以避免出现这种情况。　　　　　　　　　

312　EXCEL与WORD中的表格与数据可互相调用,EXCEL中的表格可直接粘贴到WORD中去,并且只需稍加改动就可形成WORD文件所需的表格,EXCEL中的数据也可直接复制到Hg或其它一些制图软件中,这样就减少了工作量,也减少了数据出错的机会。数据调到WORD后也不需再进行校对,在WORD中对数据进行

　·70·昆虫知识　ENTOMOLOGICALKNOWLEDGE200239(1)

叙述时,可通过表格转为文字功能避免人工输

入数据。由于EXCEL中的数据也可直接复制到Hg或其它一些制图软件中,因此,在绘制LD-P线时,也不需要再手工输入数据。313　由于EXCEL具有连环计算能力,只要某一格的数据有了变动,后面与之有联系的所有格的数据就会自动发生改变,因此,只要校对

EXCEL中显色格的数据,就完成了所有的较对

工作。

参　考　文　献

1　张宗炳.杀虫药剂的毒力测定.北京:科学出版社,1988.

369～372.

2　吴文君,刘惠霞,朱靖博.天然产物杀虫剂.陕西:科学

技术出版社,1998.139.

藏经楼粉螨综合防治案例　王飞生　　　　　刘敏超

(广东省韶关市贸易学校　广东　512031)　(湛江海洋大学　广东　524088)

Acaseofintegratedcontrollingofacaridsinthesutralibrary.WANGFei2Sheng　(ShaoguanTradeSchoolofGuandong,Shaoguan　512031,China);LIUMin2Chao　(ZhanjiangOceanUniversity,Zhanjian　524088,China).

Abstract　Usingelectrictemperaturerising,CaCl2drying,NaHCO2andC2H2O4emittingCO2,fumigationwithPH3andDDVP2C2H5OHapplying,theeffectofsyntheticpreventionofTyroglyphidaeinthesutralibraryisremarkable.Thismethodhasnotanyharmtoculturalrelicinthesutralibrary.

Keywords　syntheticprevention,Tyroglyphidae,protectionofculturalrellic

摘　要　应用电热器升温,CaCl2吸湿,NaHCO3和草酸晶粉释放CO2,磷化铝熏蒸,DDVP乙醇混合液综合防治某寺院藏经楼的粉螨,杀虫效果显著,且对藏经楼的文物没有危害。关键词　综合防治,粉螨,文物保护

　　广东韶关某寺院是千年古刹,寺院的藏经

楼藏有经书、古籍等上万册。藏经楼内部分书籍和建筑物,螨类危害严重,寺院方面曾使用市售的菊酯类杀虫剂进行防治,效果很不理想,加之使用这类杀虫剂不能有效杀灭螨类的卵,往往杀虫后一段时间又出现为害情况,令人十分头痛。寺院方面通过粮食部门与作者单位取得联系,要求对寺院藏经楼的螨类进行有效防治并确保文物不受损害。经组织有关技术人员进行调查研究,发现螨类危害严重的关键原因是螨类的卵没有得到有效控制,据此制定了综合防治方案,通过多种措施的防治,取得了理想的效果。具体介绍如下。

1　寺院藏经楼环境和螨类危害情况

寺院藏经楼系砖木结构,其中有多根大圆

2

木柱,面积120m,内有多排书柜和书架,藏有多种多样的经书、古籍。该寺院处于树林茂密的深山之中,每当春夏之际,当地气温一般在20～30℃,早晚云雾缭绕,空气相对湿度在75%～90%之间,十分适宜螨类生长繁殖,尤其是书虱类害虫生长繁殖快。

现场调查发现,螨类为害情况比较严重,主要发生于藏经楼的部分经书、古籍中以及木质材料(书柜、书架、柱梁)中,可见灰白色粉状的尘粒,在东南角雾气容易进入的地方危害较为

收稿日期:2001202220,修回日期:2001207202