馬鈴薯高產(chǎn)栽培的最優(yōu)平衡不完全區(qū)組試驗設(shè)計與分析

張 毅，方志榮，彭正松，李佩華，曲繼鵬

（西昌學(xué)院，四川 西昌 615013）

涼山州位于四川省的西南部，屬于典型的高原山地，區(qū)內(nèi)光照充足、晝夜溫差大，十分有利于馬鈴薯的生長，是四川省馬鈴薯的主產(chǎn)區(qū)，尤其是二半山（海拔1 500 m）以上的彝族聚居區(qū)，馬鈴薯不僅是農(nóng)民糧食的主要來源，也是經(jīng)濟的主要來源，在該地區(qū)發(fā)展馬鈴薯產(chǎn)業(yè)是貧困彝區(qū)農(nóng)民脫貧致富的重要渠道。2015年，涼山州馬鈴薯平均單產(chǎn)22.5 t/hm2,遠遠低于美國、比利時、新西蘭等國的46.9～48.7 t/hm2（世界糧農(nóng)組織資料）。馬鈴薯產(chǎn)量受品種[1]、栽培密度[2-3]、氮磷鉀肥配比[4-5]、微生物菌劑[6]、生長調(diào)節(jié)物質(zhì)如S-誘抗素[7-8]等的影響。為了提高涼山州馬鈴薯產(chǎn)量，有必要進行馬鈴薯高產(chǎn)栽培技術(shù)的研究。

在高產(chǎn)栽培技術(shù)的研究中，主要采用正交試驗[3,5]、二次回歸正交組合[2,4,9]等的方法。正交試驗需要進行重復(fù)試驗，要求每個區(qū)組內(nèi)的小區(qū)數(shù)要等于試驗的處理數(shù)，但由于受試實驗地的限制，導(dǎo)致區(qū)組內(nèi)的小區(qū)數(shù)小于處理數(shù)而使正交試驗受到限制。二次回歸正交組合因其處理數(shù)較多，也會因為地形、地塊大小或土壤肥力不均而使試驗受到限制。平衡不完全區(qū)組設(shè)計（Balanced Incomplete Blocks Design簡稱BIBD）作為一種析因試驗設(shè)計，可以在被試對象數(shù)目受限的條件下進行試驗設(shè)計，分析試驗因子與響應(yīng)之間的關(guān)系[10]。平衡不完全區(qū)組設(shè)計用于田間試驗主要有以下優(yōu)點：一是在一個區(qū)組內(nèi)只包含全部試驗處理的一小部分，區(qū)組內(nèi)的小區(qū)數(shù)小于試驗的處理數(shù)，因而比較容易滿足同一區(qū)組內(nèi)試驗條件相同的要求[11-12]；二是可以把要比較的幾種處理分成多個區(qū)組，只要求組內(nèi)處理的試驗條件基本相同，而允許區(qū)組間存在較大的差異[13]，從而減少了地形、地塊大小等對試驗的限制。但以上文獻雖然給出了試驗設(shè)計表，也給出了試驗結(jié)果分析的基本公式，但是由于實際操作起來比較繁瑣，對試驗設(shè)計和數(shù)據(jù)處理沒有結(jié)合具體的軟件環(huán)境給出詳細(xì)的過程，大大限制了最優(yōu)平衡不完全區(qū)組設(shè)計在田間試驗中的應(yīng)用。高化猛等[10]雖然給出了最優(yōu)平衡不完全區(qū)組設(shè)計及分析的SAS程序，但是由于研究領(lǐng)域的不同，導(dǎo)致程序使用過程中存在一定的困難。西昌學(xué)院馬鈴薯課題組在馬鈴薯高產(chǎn)栽培技術(shù)的研究過程中，進行了詳細(xì)的最優(yōu)平衡不完全區(qū)組試驗設(shè)計，在試驗數(shù)據(jù)的基礎(chǔ)上，基于SAS環(huán)境，實現(xiàn)了最優(yōu)平衡不完全區(qū)組試驗設(shè)計和方差分析，并給出了詳盡的SAS程序。本文內(nèi)容對其它作物高產(chǎn)栽培的最優(yōu)平衡不完全區(qū)組試驗的試驗設(shè)計和數(shù)據(jù)處理具有參考價值和借鑒意義。

1 試驗設(shè)計及試驗方法

1.1 選取的因素及水平

為了讓試驗設(shè)計更具有參考性，試驗選用3個馬鈴薯品種A1、A2、A3以及3種栽培措施B1、B2、B3，共9個處理（Treats），安排3個區(qū)組(Blocks)，共需27個小區(qū)，因地域所限，每個區(qū)組最多可以安排6個小區(qū)，計劃不能實施。因此，不考慮與區(qū)組有關(guān)的互作，以考察效應(yīng)Blocks，A，B，A*B為目標(biāo)，設(shè)計了一個具有3個區(qū)組，18個小區(qū)的最優(yōu)平衡不完全區(qū)組試驗，以確定影響產(chǎn)量的主要因子和最佳處理。

1.2 SAS的最優(yōu)平衡不完全區(qū)組試驗設(shè)計

在SAS中，利用factex設(shè)計試驗點，并由output語句將其輸出到指定的數(shù)據(jù)文件；optex過程完成平衡不完全區(qū)組設(shè)計及其優(yōu)化。SAS程序如下：

運行上面的SAS程序，SAS最優(yōu)平衡不完全區(qū)組試驗設(shè)計得到的結(jié)果見表1。

表1 最優(yōu)平衡不完全區(qū)組設(shè)計表

為了便于試驗的實施及結(jié)果的記錄，將表1進行整理，整理后的試驗設(shè)計見表2。

1.3 試驗實施

馬鈴薯的栽培試驗在位于普格縣五道箐鄉(xiāng)的西昌學(xué)院馬鈴薯種植基地內(nèi)完成。馬鈴薯種植基地的土壤特征為：全氮3.54 g/kg,全磷1.26 g/kg，全鉀19.83 g/kg，堿解氮317 g/kg。試驗選用了3個品種的原原種：A1（青9）、A2（大西洋）、A3（夏坡蒂）；采用了3種栽培措施：B1（對照，常規(guī)施肥）、B2（常規(guī)施肥+微生物菌劑）、B3（常規(guī)施肥+噴施s-誘抗素），試驗安排3個區(qū)組，18個小區(qū)，小區(qū)面積12.33 m2，常規(guī)大田管理，試驗在2017年4—8月完成。在收獲期，測定各小區(qū)的馬鈴薯產(chǎn)量。

2 試驗結(jié)果與分析方法

2.1 數(shù)據(jù)的獲取與整理

在馬鈴薯的收獲期，測定各小區(qū)的馬鈴薯產(chǎn)量，并將各小區(qū)產(chǎn)量的數(shù)據(jù)填入表2，以便于通過SAS程序進行后續(xù)的數(shù)據(jù)分析。

表2 最優(yōu)平衡不完全區(qū)組設(shè)計表及試驗結(jié)果

2.2 SAS的最優(yōu)平衡不完全區(qū)組試驗分析

2.2.1 應(yīng)用SAS程序進行數(shù)據(jù)處理

將表3的數(shù)據(jù)樣本創(chuàng)建為SAS數(shù)據(jù)表sasuser.malingshu01，采用glm過程對數(shù)據(jù)進行方差分析，建模語句指定因子Output的效應(yīng)為Blocks、A、B和A*B，對品種和栽培措施影響產(chǎn)量的效果進行兩兩比較。SAS程序如下：

程序輸出的方差分析結(jié)果和均值多重比較結(jié)果整理后見表3～7。

表3 馬鈴薯種植試驗因子效應(yīng)方差分析表

表4 品種的均值多重比較

表5 栽培措施的均值多重比較

2.2.2 因子效應(yīng)分析

因子效應(yīng)分析模型的決定系數(shù)R2為0.911 8，顯著性P值為0.007 7，方差分析有效。品種、品種栽培措施對馬鈴薯產(chǎn)量有顯著影響，其顯著性P值分別為0.001 9，0.003 0；區(qū)組和品種栽培互作對馬鈴薯產(chǎn)量無顯著影響（表3）。品種（因子A）LSD分析結(jié)果表明：品種對產(chǎn)量的影響從大到小排序依次為：A1、A2、A3，0.01水平上品種A1，A2對產(chǎn)量的影響顯著高于A3品種（表4）。栽培措施（因子B）LSD分析結(jié)果表明：栽培措施對產(chǎn)量影響的大小依次為：B2、B3、B1，0.01水平上栽培措施B2能明顯提高馬鈴薯的產(chǎn)量。

2.2.3 處理效應(yīng)的分析結(jié)果

處理效應(yīng)模型的決定系數(shù)R2=0.911 2，顯著性P值為0.000 7，方差分析有效。試驗處理對產(chǎn)量影響顯著，其顯著性P值為0.000 7（表6）。各處理效應(yīng)的LSD分析結(jié)果表明：試驗處理對產(chǎn)量的影響從大到小排序依次為：T2、T5、T8、T3、T4、T1、T6、T9、T7，在0.05水平上，T2處理與T5無顯著差異，顯著高于T8、T3、T4、T1、T6、T9、T7處理；在0.01水平上，T2處理與T5、T8、T3、T4處理無顯著差異，顯著高于T1、T6、T9、T7處理。

表6 玉米種植試驗處理效應(yīng)方差分析表

表7 試驗處理的均值多重比較

2.3 栽培措施的優(yōu)化

通過以上分析可以得出影響該地區(qū)馬鈴薯產(chǎn)量的主要影響因子為品種和栽培措施，品種A1、A2與B2的處理組合T2和T5均獲得最高產(chǎn)量，故確定適宜該地區(qū)土壤和氣候條件的最佳品種為A1、A2，最佳栽培方式為B2。

3 結(jié)語

最優(yōu)平衡不完全區(qū)組設(shè)計可以在被試對象數(shù)目受限的條件下進行析因試驗設(shè)計，方差分析是最優(yōu)平衡不完全區(qū)組設(shè)計的主要數(shù)據(jù)分析手段，能夠確定主要影響因子，以及各試驗因子對試驗結(jié)果影響的差異，同時能確定最佳的處理。本文圍繞上述內(nèi)容展開，通過SAS程序得到了試驗設(shè)計表，并通過方差分析法對試驗數(shù)據(jù)進行了分析，發(fā)現(xiàn)了顯著影響馬鈴薯產(chǎn)量的因子及其差異程度，并得到了最佳處理。該方法也適用于不同地域、不同品種及不同栽培措施對馬鈴薯產(chǎn)量的影響，對其它作物的高產(chǎn)栽培的最優(yōu)平衡不完全區(qū)組試驗的試驗設(shè)計和數(shù)據(jù)處理具有參考價值和借鑒意義。