应用统计软件SPSS拟合生长曲线方程戴国俊1,王金玉13,杨建生2,李尚民1,原新廷1(1.扬州大学动物科学与技术学院,江苏扬州225009;2.江苏省昆山市农林局,江苏昆山215300)摘要:介绍了应用统计软件SPSS中非线性回归子程序进行生长曲线拟合的步骤,运用实例分析表明,应用SPSS软件拟合生长曲线具有交互性好,易学易用、方便快截的特点。关键词:生长曲线;统计软件SPSS;非线性回归中图分类号:S8205文献标识码:B文章编号:052925130(2006)0920028203生长曲线广泛应用于描述动植物的生长过程,尤其是家禽。生长曲线分析,实质上是把各不同时期的体重资料提供的信息综合成少数几个参数。戴国俊等[1]研究了不同环境对新扬州鸡体重、体尺生长曲线的影响。吕敏芝等[2]对不同饲养方式下马岗鹅生产曲线进行分析,采用Logistic方程描述其生长模式及特点。王志跃等[3]运用Logistic、VonBertalanffy和Gompertz三种非线性生长模型对新扬州鸡的生长曲线进行拟合,所选生长曲线模型可以较好地预测机体的生长发育情况,利用生长曲线提供的信息进行饲养管理,提高了经济效益。S1Mignon2Grasteau等[4]通过对鸡8周龄体重的正向选择,研究选择对Gompertz曲线参数的影响,结果选择提高了鸡6～16周龄的体重,同时,瞬时生长率和成熟体重也得到了提高,而且生长拐点提前了。因此,不同个体生长曲线的特征被看作是重要的综合经济指标,用于家禽品种的选育,可为选育适合不同要求的家禽品种提供依据。手工三点法进行生长曲线拟合不仅计算工作量大,而且误差也大,SPSS统计软件相对其它如SAS、STATISTICA等软件而言,具有交互性好,不需编程,易学易用,方便快捷等特点,本文介绍应用统计软件SPSS(Ver111101)拟合生长曲线方程的方法,为生长曲线拟合分析提供方便。1常用生长曲线的类型生长曲线方程大致可以分为三类,一类是表示报酬递减表现的方程,如指数函数;另一类是描述光滑S型曲线,且有一个固定拐点的方程,如Gompertz,Logistic;还有一类也是描述光滑S型曲线,但拐点是可变的方程,如VonBertalanffy,Richards等[5]。表1是3个常用描述生长过程的曲线方程及其描述生长过程的特点,3种曲线模型中待定参数A表示极限生长量(终极生长量或成熟体重),K为瞬时相对生长率,B为常数尺度。表13种常用动物生长曲线模型名称模型拐点体重(w)拐点日龄最大日增重相对生长率Gomp.Y=Ae-bexp(-kt)A/e(lnB)/kkwk(lnA-lnWi)Logi.Y=A/(1+Be-kt)A/2(lnB)/kkw/2k(1-Wi/A)Bert.Y=A(1-Be-kt)38A/27(ln3B)/k3kw/23k((A/Wi3)1/3-1)3Wi为第i日龄体重。2统计软件拟合生长曲线的方法211数据输入在SPSSVariableView界面以列方式定义2个变量,即生长时间(time)和生长量(weight)。在SPSSDataEditor的两向电子表格中,按测定的时间自上而下依次输入数据,数据输入方法如图1。收稿日期:2005212229基金项目:江苏省农业高新技术项目(BG2004316)。作者简介:戴国俊(19632),男,副教授,博士。3通讯作者。图1数据输入对话框·82·AnimalHusbandry&VeterinaryMedicine2006Vol138No19212程序选择(以Logistic模型为例)在SPSSDataEditor界面下拉菜单中点击Analyze,再选择Regression子菜单,在该子菜单中选择NonlinearRegression程序。此时出现一交互式对话框,如图2所示,接着将左侧框内的依变量Weight选入右侧Dependent框内,然后将Logistic曲线模型表达式A/(1+B3EXP(2K3time))输入到ModelExpression框中,式中A、B和K为待定参数,time为自变量,由左侧变量框内选入曲线模型相应位置,EXP是自然对数底数e的指数函数,由右下Functions框内提供的函数中选入,在输入表达式时加、减、乘、除等运算符号需严格输入。图2非线性回归对话框图3待定参数赋值对话框最后点击左侧按钮Parameters出现图3所示对话框,在该对话框中输入曲线模型中待定参数的名称(Name)以及初始值(StartingValue),每输完一个待定参数按添加按钮(Add),如有输入错误,可用Change或Remove按钮修改或去除,当所有的待定参数均输入后再按继续(Continue)。要注意的是在给模型参数初始值赋值时,最好与其真实值大致相当,因为最后模型参数是通过迭代方法计算,使平方和变化达到规定值后收敛,给出模型参数,所以如果参数所赋值相差很大,有可能使最后的计算不能收敛...
