李必富

摘要：對(duì)紛繁復(fù)雜的數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析篩選出目標(biāo)品種，是育種家必須面對(duì)的問(wèn)題。采用3種不同高產(chǎn)穩(wěn)產(chǎn)統(tǒng)計(jì)分析方法對(duì)實(shí)際試驗(yàn)數(shù)據(jù)進(jìn)行處理，并比較分析結(jié)果，表明應(yīng)用高穩(wěn)系數(shù)法對(duì)玉米進(jìn)行高產(chǎn)穩(wěn)產(chǎn)統(tǒng)計(jì)分析最為簡(jiǎn)便、可靠。

關(guān)鍵詞：玉米；高產(chǎn)穩(wěn)產(chǎn)；統(tǒng)計(jì)分析；常規(guī)方法；Eberhart-Russell模型；高穩(wěn)系數(shù)法

中圖分類號(hào)：S513 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1674-1161（2017）04-0064-03

選育高產(chǎn)穩(wěn)產(chǎn)品種是玉米育種的最終目的，而環(huán)境與基因互作是玉米品種最終產(chǎn)量形成的根本，因此，對(duì)選育的品種在投放市場(chǎng)前需經(jīng)過(guò)多年的種植鑒定分析。初期育種規(guī)模小，僅用肉眼經(jīng)驗(yàn)觀察即可；而現(xiàn)在育種規(guī)模逐年擴(kuò)大，新組合鑒定數(shù)量和鑒定點(diǎn)次成倍增加，特別是品種審定制度的推行，均要求“靠數(shù)據(jù)說(shuō)話”。

在玉米育種中，系統(tǒng)客觀地評(píng)價(jià)一個(gè)品種的高產(chǎn)穩(wěn)產(chǎn)性有多種方法。傳統(tǒng)的方法采用平均數(shù)值（x）較對(duì)照增產(chǎn)百分比來(lái)估算品種的產(chǎn)量性能，以SSR測(cè)驗(yàn)法估算品種間產(chǎn)量差異大小（高產(chǎn)性），以標(biāo)準(zhǔn)差（S）、變異系數(shù)（CV）估測(cè)產(chǎn)量的穩(wěn)定性。1966年，Eberhart和Russell提出以回歸模型分析品種穩(wěn)定性的方法，通過(guò)計(jì)算第一穩(wěn)定系數(shù)（bi）和第二穩(wěn)定系數(shù)（S2di）評(píng)價(jià)品種的穩(wěn)定性。1994年，溫振明提出高穩(wěn)系數(shù)法（HSC），可以在一個(gè)數(shù)列里面表現(xiàn)出品種的豐產(chǎn)性和穩(wěn)定性。這3種方法原理各不相同，計(jì)算結(jié)果體現(xiàn)出的結(jié)論也不盡相同。本研究通過(guò)試驗(yàn)數(shù)據(jù)分析，比較三者在實(shí)際生產(chǎn)生活中的應(yīng)用效果。

1 材料與方法

1.1 試驗(yàn)資料

數(shù)據(jù)資料來(lái)源于2016年鐵嶺市農(nóng)業(yè)科學(xué)院在遼寧省內(nèi)異地鑒定區(qū)域試驗(yàn)。試驗(yàn)共6個(gè)品種（含對(duì)照），5個(gè)點(diǎn)次；種植采用壟作，5行區(qū)，3次重復(fù)，密度為60 000株/hm2。

1.2 數(shù)據(jù)統(tǒng)計(jì)方法

1.2.1 常規(guī)方法 計(jì)算平均產(chǎn)量■i=ΣXij/n，其中i=1，…，v，j=1，…，t，與對(duì)照的增產(chǎn)比可反映品種的豐產(chǎn)性。標(biāo)準(zhǔn)差Si=sqrt[ΣX2ij-（ΣXij）2/n]/（n-1），變異系數(shù)CVi=Si/■i×100，可反映品種的穩(wěn)定性。品種間產(chǎn)量差異采用SSR測(cè)驗(yàn)法計(jì)算，反映品種間差異的顯著性。

1.2.2 Eberhart-Russell模型 Eberhart-Russell模型主要通過(guò)計(jì)算兩個(gè)穩(wěn)定參數(shù)來(lái)考量品種的高產(chǎn)穩(wěn)產(chǎn)性能，將環(huán)境變量以指數(shù)形式體現(xiàn)。以每一個(gè)試點(diǎn)1 a中參試品種的平均產(chǎn)量減去所有試點(diǎn)產(chǎn)量的平均值，作為環(huán)境指數(shù)I，Ij=■.j-■..，計(jì)算每一個(gè)品種的平均產(chǎn)量相對(duì)于環(huán)境指數(shù)的回歸系數(shù)b，稱為第一穩(wěn)定參數(shù)，bi=Σ（Xij·Ij）/ΣIj2。因?yàn)閎i的平均數(shù)等于1，因此可以用測(cè)驗(yàn)bi與1的差異的t測(cè)驗(yàn)值反映差異的顯著性。將直線回歸離差S2di稱為第二穩(wěn)定系數(shù)，S2di=Σδ2ij/（t-2）-S2e/r，其中δ2ij為第i個(gè)品種的離回歸平方和，S2e為合并試驗(yàn)誤差。定義bi=1，S2di=0的品種為穩(wěn)定性品種。S2di值小，說(shuō)明這一品種的產(chǎn)量是隨環(huán)境指數(shù)變化而變化的，對(duì)這一品種在各個(gè)環(huán)境下的產(chǎn)量變動(dòng)可以較可靠地估測(cè)；S2di值大，說(shuō)明這一品種的產(chǎn)量變動(dòng)不能可靠地估測(cè)，品種的產(chǎn)量穩(wěn)定性差。

總平方和SS=ΣX2ij-（ΣXij）2/t，回歸平方和SS=（ΣXij·Ij）2/ΣIj2，離回歸平方和Σδ2ij=總平方和SS-回歸平方和SS，回歸離差S2di=Σδ2ij/（t-2）-S2e/r。最后用F測(cè)驗(yàn)法檢驗(yàn)穩(wěn)定性差異的顯著性，F(xiàn)=ΣΣδ2ij/（t-2）/（S2e/r）

1.2.3 高穩(wěn)系數(shù)法（HSC） 溫振明提出的高穩(wěn)系數(shù)計(jì)算公式為HSCi=（Ga-Gi）/Ga×100，以百分?jǐn)?shù)的形式顯示，其中Ga為目標(biāo)品種的穩(wěn)定產(chǎn)量，Gi為參試品種的穩(wěn)定產(chǎn)量。育種中常簡(jiǎn)單計(jì)，以比對(duì)照增產(chǎn)10%作為目標(biāo)品種的穩(wěn)定產(chǎn)量（即Ga=1.1×■ck），以參試品種的多點(diǎn)平均值減去標(biāo)準(zhǔn)差作為參試品種的穩(wěn)定產(chǎn)量（即Gi=■i-Si），因此，可以HSCi=（1.1×■ck-■i+Si）/（1.1×

■ck）×100來(lái)計(jì)算品種的高穩(wěn)系數(shù)。HSC法同時(shí)考量了品種的豐產(chǎn)性和穩(wěn)產(chǎn)性，HSCi值越小，說(shuō)明品種的高產(chǎn)穩(wěn)產(chǎn)性越好。

1.2.4 相關(guān)性計(jì)算 利用Excel相關(guān)性計(jì)算公式計(jì)算。

2 結(jié)果與分析

2.1 供試品種的產(chǎn)量計(jì)算和穩(wěn)定性分析

5個(gè)供試品種及對(duì)照在5個(gè)試驗(yàn)點(diǎn)的平均產(chǎn)量結(jié)果及各項(xiàng)穩(wěn)定性參數(shù)計(jì)算結(jié)果見(jiàn)表1和表2。

由表1和表2可知：平均產(chǎn)量超過(guò)對(duì)照的有v1，v4，v3這3個(gè)品種，但均沒(méi)有超過(guò)10%，互相之間產(chǎn)量的差異也沒(méi)有達(dá)到顯著水平，但v1與對(duì)照產(chǎn)量差異達(dá)到極顯著水平，v4和v3則與對(duì)照差異不顯著。

根據(jù)計(jì)算原理，bi=1，S2di=0的品種穩(wěn)定，因此可以通過(guò)測(cè)定bi與1的差異及S2di與0的差異來(lái)反映品種的穩(wěn)定性差異。僅有v2的t測(cè)驗(yàn)達(dá)到極顯著，說(shuō)明該品種的穩(wěn)定性極差。S2di與0的F測(cè)驗(yàn)值均不顯著，表明本試驗(yàn)通過(guò)S2di值計(jì)算不能體現(xiàn)品種間的穩(wěn)定性差異，只能通過(guò)S2di值大小簡(jiǎn)單評(píng)判品種的穩(wěn)定性大小。通過(guò)結(jié)合產(chǎn)量及標(biāo)準(zhǔn)差變異系數(shù)分析，v1，v4，v3均為高產(chǎn)穩(wěn)產(chǎn)品種，其中v4品種的標(biāo)準(zhǔn)差和變異系數(shù)最小，穩(wěn)定性最好；通過(guò)結(jié)合產(chǎn)量及穩(wěn)定性參數(shù)bi分析，v1，v4，v3為高產(chǎn)穩(wěn)產(chǎn)品種，其中v1品種的bi值最接近1，穩(wěn)定性最好；通過(guò)結(jié)合產(chǎn)量和穩(wěn)定性參數(shù)S2di分析，v1，v4，v3為高產(chǎn)穩(wěn)產(chǎn)品種，其中v1品種的S2di值最接近0，穩(wěn)定性最好。HSCi值綜合考量了品種的產(chǎn)量和變異情況，根據(jù)數(shù)據(jù)可知高產(chǎn)穩(wěn)產(chǎn)品種為v1，v4，v3，其中v1品種的HSCi值最小，高產(chǎn)穩(wěn)產(chǎn)性最好。

2.2 高產(chǎn)穩(wěn)產(chǎn)相關(guān)數(shù)據(jù)排序及相關(guān)性分析

以品種產(chǎn)量平均數(shù)大小及穩(wěn)定性參數(shù)表示穩(wěn)定性高低進(jìn)行排序，結(jié)果見(jiàn)表3。

對(duì)表3中產(chǎn)量與穩(wěn)定性參數(shù)進(jìn)行相關(guān)性計(jì)算，結(jié)果見(jiàn)表4。將較對(duì)照增產(chǎn)的品種單獨(dú)進(jìn)行各穩(wěn)定性參數(shù)的相關(guān)性分析，結(jié)果見(jiàn)表5。將增產(chǎn)品種和對(duì)照一起進(jìn)行相關(guān)穩(wěn)定性參數(shù)的相關(guān)性分析，結(jié)果見(jiàn)表6。

由表4、表5和表6可知：僅高穩(wěn)系數(shù)與平均產(chǎn)量成正相關(guān)，其余指數(shù)與平均數(shù)相關(guān)性變化不穩(wěn)定，不能單獨(dú)體現(xiàn)出品種的高產(chǎn)豐產(chǎn)性，只能結(jié)合兩個(gè)數(shù)據(jù)一起簡(jiǎn)單比較分析。本試驗(yàn)平均數(shù)排序和HSCi值排序結(jié)果相同，則通過(guò)表4中數(shù)據(jù)可得高穩(wěn)系數(shù)與標(biāo)準(zhǔn)差和變異系數(shù)的相關(guān)性中等。

3 結(jié)論與討論

結(jié)合試驗(yàn)數(shù)據(jù)，采用3種不同計(jì)算分析方法得出的品種高產(chǎn)穩(wěn)產(chǎn)結(jié)論基本相同，只是在穩(wěn)定性大小方面有所差異。從計(jì)算過(guò)程來(lái)看，常規(guī)方法簡(jiǎn)單，筆算或者使用計(jì)算器很容易計(jì)算；Eberhart-Russell模型考量的變量較多，計(jì)算較為復(fù)雜，使用計(jì)算機(jī)軟件也需要多個(gè)步驟才能計(jì)算；而高穩(wěn)系數(shù)法既考量了豐產(chǎn)性又考量了穩(wěn)產(chǎn)性，計(jì)算過(guò)程也非常簡(jiǎn)單。因此，在玉米高產(chǎn)穩(wěn)產(chǎn)鑒定試驗(yàn)統(tǒng)計(jì)分析時(shí)，建議采用高穩(wěn)系數(shù)分析法。這與劉海燕、李恒嶺等人的結(jié)論一致。

由于試驗(yàn)容量的限制，試驗(yàn)數(shù)據(jù)偏少，誤差較突出，各穩(wěn)定性參數(shù)的相關(guān)性分析結(jié)論不夠嚴(yán)謹(jǐn)，有待于進(jìn)一步研究。

參考文獻(xiàn)

[1] 蓋鈞鎰.試驗(yàn)統(tǒng)計(jì)方法[M].北京：中國(guó)農(nóng)業(yè)出版社，2013.

[2] 胡秉民，耿旭.作物穩(wěn)定性分析法[M].北京：科學(xué)出版社，1993.

[3] 溫振明，張永科.用高穩(wěn)系數(shù)法估算玉米雜交種高產(chǎn)穩(wěn)產(chǎn)性的探討[J].作物學(xué)報(bào)，1994，20（4）：508-512.

[4] 劉海燕.利用高穩(wěn)系數(shù)法綜合分析玉米新品種的研究[J].玉米科學(xué)，2003，11（1）：37-38.

[5] 李恒嶺.對(duì)2012—2013年遼寧省玉米區(qū)域試驗(yàn)中熟組參試品種的高穩(wěn)系數(shù)法評(píng)價(jià)[J].遼寧農(nóng)業(yè)科學(xué)，2016（3）：77-79.