《spss进行籽粒灌浆拟合的详细步骤》由会员分享，可在线阅读，更多相关《spss进行籽粒灌浆拟合的详细步骤（10页珍藏版）》请在金锄头文库上搜索。
1、SPSS 进行灌浆拟合的详细过程一、取样： 进行灌浆曲线拟合，取样非常重要，务必准确标记抽穗开花时间， 同时持续较长时间的取样，一般花后间隔3-5 天取样，并持续至充分 的黄熟收获取样时间。例如：水稻花后每5天取样，直至成熟，样品80C烘箱处理48 h,质量衡定后脱去颖壳称质量，并换算成千粒质量，数据如下取样时间（天）干物重（g/1000粒）50.921107.5811514.3422016.3922516.0263016.292（一）参数初始值的设定籽粒灌浆拟合实质是曲线回归。曲线回归拟合是否成功，拟合的 好坏与参数初始值的定义密切相关。那么灌浆的参数初始值应该定义 成多少呢？这个难道了很多
2、人。灌浆拟合采用的是Richards方程，为W=A （1+B e-kt） -in。当N=1 时，Richards 方程就等于 Logistics 方程：Y=A/ (1+ B e-k t)。因 此在进行 Richards 方程拟合之前可进行 Logistics 方程的拟合，以获 得其他参数B、K、A的初始值。1、打开spss软件并输入数据，如图点击【分析】、【回归】、【曲线估 计】。y土未标遁1-IBM SPSS StoTistics 額堀藹蛭器KPFO US.El flSiV teISQi UA iJflriAl 直轉過EH& 2实用ISffiJi 匿口朗留助H 国r鮭*垂 CT)b斓均回叩*
3、一 G如供粽空 *S3IVrii克量-*-* 115 DO翌210 OO7 53广又线性掺型*混台嗣凶*315 0014 34Am OO16 19相天ICIEo3S.OO16 03:手汪It:hfn 巧V：217 u：、四目埶戕性翅損砂 丽址* L：.由瓏计. 翔部甘母=-.30 001S 29丁 r= LagiBbQ 吗 # - LogiElk- 力石F fflPftJbriu1516毎总蛆|型因磁ifl労折芒H；l刃弋凰|. fflr.i l -iwi.V7i? 1- -、工亍壬圧1、El Roc 曲即 *QA1HEG：.192021!?2、打开曲线估计对话框，将 x 选入自变量， y 选
4、入因变量，同时在模型处选中logistic，在【下限】框内，输入比最大y值略大一些的数值，本处可输入16.5。然后勾选上显示anova表格，如下图所示3、点击确定输出以下结果LogisticRR方调鏗尺方估计值的标淮933.871.0391.120自变邑为xANOVA平方和df均方FSigI；怕34 389134.3S927.049.0075.08541 271总计39 4755x-未标淮f匕勲数1Sig.E掃准谣X(戦.7561.366.0411.434.39310.551.953.000.395因型量为 ln1 IX- 1 M6.5DDi -从系数表中可知 b0=1.366，b1=0.75
5、6从观测值数据中选取3对观察值，x为3个等间距值，如5、10、15。3 对观测值可为 %, y1)= (5, 0.921), (x2，y2)= (10,7.581), (x3， y3) =(15,14.342) .按照以下公式计算出 A 值：A=(y22*(y1+y3)-2*y1*y2*y3) / (y22-y1*y3)=15.293然后 B=b0*A=20.89，K=-ln b=-ln0.756=0.279。因此得到初始值：A=15293, B=2089, K=0279, N=1.（二） Richards方程拟合1、如图所示选择【非线性】2 打开非线性对话框，将 y 选入因变量，同时将 Ri
6、chard 方程按照 spss的语法规范输入模型。如下图所示I國数和特跌变星旧:函數组:全部算术CDF与菲中心CDF转换当前日期耐间曰期运算日明创淫抽取日期SHHQ0OI3OS口EIEEEE然后在非参数对话框内，点击【参数】，打开参数对话框，并输入之前算出的参数初始值非线性回归：蓼数使用上一券析的起始值型）帮助B(2a K(a.279JN.1)(1 4 B c EXPt-K *)*mE鏈到IJ三保存蟲也空卬琳总旦lF)C|确定 融且h 重置迟) 取莎二| m兰邛CBF与圭中心GDF当前曰期耐间曰卷走豆曰肚咙他虫1_1軼L*J和根芟去过m M3、点击确定，输出结果鑫皴估计憶feif95% I嚨A
7、10.255.13515.67316.837B33.00027.930-07.171153.172K.407041.231.584N.583.214-.3381.505寥纹他计值冏相尖注AEKNA1.000-.448-.540-.370E-.4481.000.956.968K-.5409561.000.857N-.370.968.8571.000ANOVAa二方和df均方l-lf=l1048.3734262.093痩差.0832041索更正的总计1048.4566i.更正的鼠计199.8905因娈量：a. R A. = 1 -(裁平方和已更正的平方 .fili = 1.000 -A=16.25
8、5，B=33，K=0.407， N=0.583最后得到 Richards 方程为：y=16255 (1+33*e-0.407*x) -1/0.5834 2 0J 1 J00101520253036花后天数对上述Richards方程求一阶导数，可得下式灌浆速率（GR,即 单位时间千粒重的增长量，g d-i）依时间x的灌浆速率方程GR=dy / dx= A k B e-k x / N（1+B e-kx）（1+N）/N，根据参数 A、B、K、N 可计算以下参数： 生长速率G为最大的H期与最大生长速率G：对方程（1）求二阶导数，可得G随时间f而改变的速率，令其为0,可得- lii.v小、1 iimx*
9、G =（彳）KK将 tmax G 作为 x 代入 Richard 方程求得的 y 值即为 Gmaxmax Gmax 平均牛长速率乙与活跃牛长期Q:对方程（2）积分得到:(4)(5)_1 ，- w = ANKG = 6t =A w 0 d/2( V + 2)活跃生长期D为生长终值上除以石，即A 2(% + 2)D 6K 划仆灌浆过稈的前、中、后期，牛长速率方稈 G具有两个拐点”求其对的二阶导数，并令为零， 可得两个拐点在t坐标上的值0和心为2 | 2 卜 3v - nW 十 5/1 = -In I）/K（6）2BN1十3N -X!N十行十5q = ln（） / K（ 7 ）-IS假定达99%A时为实际灌浆终值期心，依（1）得LOO N I /（）-1
