spss進(jìn)行籽粒灌漿擬合的詳細(xì)步驟

SPSS進(jìn)行灌漿擬合的詳細(xì)過程一、取樣：進(jìn)行灌漿曲線擬合，取樣非常重要，務(wù)必準(zhǔn)確標(biāo)記抽穗開花時間，同時持續(xù)較長時間的取樣，一般花后間隔3-5天取樣，并持續(xù)至充分的黃熟收獲取樣時間。二、灌漿擬合例如：水稻花后每5天取樣，直至成熟，樣品80℃烘箱處理48h，質(zhì)量衡定后脫去穎殼稱質(zhì)量，并換算成千粒質(zhì)量，數(shù)據(jù)如下取樣時間（天）干物重（g/1000粒）50.921107.5811514.3422016.3922516.0263016.292（一）參數(shù)初始值的設(shè)定籽粒灌漿擬合實(shí)質(zhì)是曲線回歸。曲線回歸擬合是否成功，擬合的好壞與參數(shù)初始值的定義密切相關(guān)。那么灌漿的參數(shù)初始值應(yīng)該定義成多少呢？這個難道了很多人。灌漿擬合采用的是Richards方程，為W=A（1+Be-kt）-1\N。當(dāng)N=1時，Richards方程就等于Logistics方程：Y=A/（1+Be-kt）。因此在進(jìn)行Richards方程擬合之前可進(jìn)行Logistics方程的擬合，以獲得其他參數(shù)B、K、A的初始值。1、打開spss軟件并輸入數(shù)據(jù)，如圖點(diǎn)擊【分析】、【回歸】、【曲線估計(jì)】。2、打開曲線估計(jì)對話框，將x選入自變量，y選入因變量，同時在模型處選中l(wèi)ogistic，在【下限】框內(nèi)，輸入比最大y值略大一些的數(shù)值，本處可輸入16.5。然后勾選上顯示anova表格，如下圖所示3、點(diǎn)擊確定輸出以下結(jié)果從系數(shù)表中可知b0=1.366，b1=0.756從觀測值數(shù)據(jù)中選取3對觀察值，x為3個等間距值，如5、10、15。3對觀測值可為（x1，y1）=（5,0.921），（x2，y2）=（10,7.581），（x3，y3）=（15,14.342）.按照以下公式計(jì)算出A值：A=(y22*(y1+y3)-2*y1*y2*y3)/(y22-y1*y3)=15.293然后B=b0*A=20.89，K=-lnb1=-ln0.756=0.279。因此得到初始值：A=15.293，B=20.89，K=0.279，N=1.（二）Richards方程擬合1、如圖所示選擇【非線性】2打開非線性對話框，將y選入因變量，同時將Richard方程按照spss的語法規(guī)范輸入模型。如下圖所示然后在非參數(shù)對話框內(nèi)，點(diǎn)擊【參數(shù)】，打開參數(shù)對話框，并輸入之前算出的參數(shù)初始值3、點(diǎn)擊確定，輸出結(jié)果A=16.255，B=33，K=0.407，N=0.583最后得到Richards方程為：y=16.255（1+33*e-0.407*x）-1/0.583對上述Richards方程求一階導(dǎo)數(shù)，可得下式灌漿速率(GR，即單位時間千粒重的增長量，gd-1)依時間x的灌漿速率方程GR=dy/dx=AkBe-kx/[N（1+Be-kx）