謝雄澤 謝捷 尹羽豐 張清偉 白桂萍

摘要：為探明各農(nóng)藝性狀對(duì)湖北省油菜（Brassica napus L.）產(chǎn)量的影響，對(duì)2016—2017年湖北省油菜區(qū)域試驗(yàn)的參試品種（系）產(chǎn)量及各農(nóng)藝性狀進(jìn)行統(tǒng)計(jì)分析。結(jié)果表明，湖北地區(qū)油菜品種（系）在理想株型的選擇方面雖受限于株高但仍具有較大的選擇空間;通過(guò)主成分分析將油菜相關(guān)性狀進(jìn)行歸類，確定每角粒數(shù)為影響油菜產(chǎn)量的主要因素;通過(guò)多重線性回歸分析構(gòu)建了湖北省油菜產(chǎn)量預(yù)測(cè)模型，確定了產(chǎn)量構(gòu)成三要素與產(chǎn)量之間的線性關(guān)系;通過(guò)相關(guān)性分析和通徑分析發(fā)現(xiàn)，除株高和生育期外的農(nóng)藝性狀均對(duì)油菜產(chǎn)量有直接的正面影響;分枝部位對(duì)有效角果數(shù)有促進(jìn)作用，分枝數(shù)對(duì)有效角果數(shù)和每角粒數(shù)有促進(jìn)作用，但會(huì)抑制千粒重;而產(chǎn)量構(gòu)成三因素之間則存在明顯的相互抑制作用，但對(duì)產(chǎn)量貢獻(xiàn)最大的為每角粒數(shù)。

關(guān)鍵詞：油菜（Brassica napus L.）;產(chǎn)量;農(nóng)藝性狀;主成分分析;多重線性回歸分析;相關(guān)性分析;通徑分析

中圖分類號(hào)：S565.4 ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2019）24-0038-05

DOI：10.14088/j.cnki.issn0439-8114.2019.24.010 ? ? ? ? ? 開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Study on the genetic characteristics of the yield of the test variety （line）

in the rape area in Hubei province from 2016 to 2017

XIE Xiong-ze，XIE Jie，YIN Yu-feng，ZHANG Qing-wei，BAI Gui-ping

（Xiangyang Academy of Agricultural Science，Xiangyang 441057，Hubei，China）

Abstract： In order to find out the influence of various agronomic traits on the yield of rape（Brassica napus L.） in Hubei province， the statistical analysis of yield and agronomic traits of tested varieties in the rape regional trials in Hubei province from 2016 to 2017. The results show that the rape varieties （lines） in Hubei province are limited by plant height but still have a large selection space for the selection of ideal plant types. Principal component analysis was used to classify the traits of rape， and the number of horns per leaf was determined as the main factor affecting yield. Multiple linear regression analysis was used to construct a prediction model of rape yield in Hubei province， a linear relationship between the three components of production and production was determined. Correlation analysis and path analysis showed that agronomic traits except plant height and growth period have a direct positive impact on rape yield. The branch site promotes the number of effective pods， and the number of branches promotes the number of effective pods and the number of pods per pod， but inhibits thousand-grain weight. There is a significant mutual inhibition between the three factors of yield composition， but the largest contribution to yield is the number of grains per horn.

Key words： rape（Brassica napus L.）; yield; agronomic traits; principal component analysis; multiple linear regression analysis; correlation analysis; path analysis

湖北省是油菜（Brassica napus L.）種植大省，也是國(guó)內(nèi)油菜種植的主要區(qū)域之一，全省油菜產(chǎn)業(yè)的健康穩(wěn)定發(fā)展為中國(guó)植物油供給安全提供了保障[1，2]。湖北省油菜品種區(qū)域試驗(yàn)是鑒定油菜育成新品種特征特性、應(yīng)用價(jià)值及適應(yīng)區(qū)域的主要途徑，參試品種的農(nóng)藝性狀及產(chǎn)量特性代表著區(qū)域內(nèi)一定時(shí)期育種的方向與生產(chǎn)水平[3]，同時(shí)也關(guān)系到地區(qū)油菜產(chǎn)業(yè)的健康穩(wěn)定發(fā)展。高產(chǎn)、穩(wěn)產(chǎn)一直是油菜育種的主要目標(biāo)，在育種過(guò)程中對(duì)相關(guān)農(nóng)藝性狀的定向選擇是保證高產(chǎn)、穩(wěn)產(chǎn)的基礎(chǔ)[4]。油菜產(chǎn)量的形成是一個(gè)動(dòng)態(tài)的復(fù)雜過(guò)程，產(chǎn)量在形成過(guò)程中油菜的主要數(shù)量性狀依次建成，相關(guān)性狀間存在相互制約、相互補(bǔ)償?shù)年P(guān)系[5]。關(guān)于農(nóng)藝性狀對(duì)產(chǎn)量形成的相關(guān)分析已有較多報(bào)道，但由于所選定作物品種及地域等因素的不同，相關(guān)研究結(jié)論不盡一致[6-8]。研究結(jié)果表明，油菜產(chǎn)量構(gòu)成除了受有效角果數(shù)、每角粒數(shù)和千粒重3個(gè)產(chǎn)量構(gòu)成因子的影響外，其他如株高、分枝部位、有效分枝數(shù)等也是影響油菜產(chǎn)量的重要因素[9-14]。對(duì)油菜產(chǎn)量形成的相關(guān)農(nóng)藝性狀進(jìn)行分析，并理清各農(nóng)藝性狀之間客觀存在的促進(jìn)及抑制作用對(duì)于挖掘油菜產(chǎn)量潛力、提升油菜品種特性有重要的指導(dǎo)作用。但近年來(lái)有關(guān)湖北省油菜區(qū)試品種（系）產(chǎn)量成因性狀分析的研究較少，本研究以2016—2017年湖北省油菜區(qū)域試驗(yàn)油菜品種（系）的匯總數(shù)據(jù)為基礎(chǔ)，深度解析湖北省油菜品種的產(chǎn)量及相關(guān)農(nóng)藝性狀之間的關(guān)系，為高產(chǎn)油菜品種的選育及遺傳改良提供參考依據(jù)。

1 ?材料與方法

1.1 ?試驗(yàn)材料

試驗(yàn)材料源于2016—2017年湖北省油菜區(qū)域試驗(yàn)的27個(gè)品種（系），試驗(yàn)點(diǎn)包括襄陽(yáng)市、沙洋縣、武漢市、潛江市、武穴市、咸寧市咸安區(qū)。

1.2 ?研究方法

通過(guò)Excel 2016軟件對(duì)各參試油菜品種（系）的相關(guān)農(nóng)藝性狀及產(chǎn)量數(shù)據(jù)進(jìn)行匯總、整理，使用SPSS 25軟件對(duì)處理后的數(shù)據(jù)進(jìn)行方差分析、主成分分析、多重線性回歸分析、相關(guān)性分析和通徑分析，并得出研究結(jié)論。

2 ?結(jié)果與分析

2.1 ?油菜農(nóng)藝性狀及產(chǎn)量的變異性分析

作物的農(nóng)藝性狀一般是由遺傳基因和環(huán)境共同決定的，從遺傳育種角度來(lái)說(shuō)，某一性狀的變異程度越大，則該性狀的遺傳力越小，選擇的可能性越大，反之則越小[15]。由表1可以看出，與油菜產(chǎn)量相關(guān)的農(nóng)藝性狀變異系數(shù)從大到小依次為千粒重9.49%、分枝部位9.16%、每角粒數(shù)7.94%、有效角果數(shù)7.86%、分枝數(shù)6.87%、產(chǎn)量6.54%、株高4.59%和生育期0.79%。其中，油菜產(chǎn)量構(gòu)成三因素（千粒重、每角粒數(shù)、有效角果數(shù)）的變異系數(shù)均比較大，這為湖北省油菜品種（系）產(chǎn)量潛力的挖掘提供了可能;另外與油菜株型有關(guān)的農(nóng)藝性狀如分枝部位、分枝數(shù)變異系數(shù)也比較大，這為油菜理想株型的選擇提供了較大空間;但株高的變異系數(shù)相對(duì)較小，在一定程度上限制了理想株型的選擇。

2.2 ?油菜產(chǎn)量與農(nóng)藝性狀的主成分分析

主成分分析可以通過(guò)線性變換，將產(chǎn)量及與產(chǎn)量相關(guān)的農(nóng)藝性狀組合成相互獨(dú)立的能夠充分反映總體信息的指標(biāo)，便于對(duì)數(shù)據(jù)進(jìn)一步分析。在主成分分析中選取初始特征值大于1的主成分，而小于1的特征值通常被認(rèn)為該因子中得到的信息不足以證明其應(yīng)該被保留。從表2可以看出，8個(gè)變量的信息可以由3個(gè)初始特征值超過(guò)1的主成分反映77.33%，說(shuō)明提取3個(gè)主成分可較充分地反映總體變量的信息。一般來(lái)說(shuō)，累計(jì)方差百分比超過(guò)70%，即認(rèn)為比較滿意[16]。

通過(guò)各農(nóng)藝性狀的特征值及特征向量（表3）進(jìn)一步生成主成分得分系數(shù)矩陣（表4），進(jìn)而得出最終的因子得分方程，如下：

F1=0.356X1+0.404X2+0.294X7-0.060X3+0.060X4+0.038X5-0.120X6+0.106Y

F2=0.100X1+0.006X2-0.177X7+0.328X3+0.417X4-0.420X5-0.056X6+0.004Y

F3=-0.042X1-0.169X2+0.187X7+0.208X3-0.156X4+0.025X5+0.538X6+0.397Y

從上述模型可以看出，第一公因子F1基本支配了X1、X2、X7（絕對(duì)值較大的系數(shù)，下同）;第二公因子F2基本支配了X3、X4、X5;第三公因子F3基本支配了X6、Y。該模型的含義：第一公因子支配的是株高（X1）、分枝部位（X2）和生育期（X7），其中株高與分枝部位受生育期的影響，故將其定義為生育期因子;第二公因子支配的是分枝數(shù)（X3）、有效角果數(shù)（X4）和千粒重（X5），其中有效角果數(shù)和千粒重受分枝數(shù)的影響，故將其定義為分枝數(shù)因子;第三公因子支配的是每角粒數(shù)（X6）和產(chǎn)量（Y），其中每角粒數(shù)會(huì)直接影響產(chǎn)量，故將其定義為產(chǎn)量因子。該模型的啟示作用：株高和分枝部位受生育期的影響較大;分枝數(shù)對(duì)有效角果數(shù)有促進(jìn)作用，對(duì)千粒重有抑制作用;而產(chǎn)量的決定性因素則在于每角粒數(shù)。

2.3 ?油菜產(chǎn)量與農(nóng)藝性狀的多重線性回歸分析

為了解各農(nóng)藝性狀對(duì)油菜產(chǎn)量的直接影響，以株高（X1）、分枝部位（X2）、分枝數(shù)（X3）、有效角果數(shù)（X4）、千粒重（X5）、每角粒數(shù)（X6）、生育期（X7）7個(gè)性狀為自變量，以產(chǎn)量（Y）為因變量，進(jìn)行多重線性回歸分析，對(duì)回歸方程進(jìn)行顯著性檢驗(yàn)，結(jié)果表明變量X與Y之間存在極顯著線性回歸關(guān)系（F=4.329）。從表5可以看出，有效角果數(shù)、千粒重、每角粒數(shù)對(duì)應(yīng)的P值均小于0.05，達(dá)到顯著或極顯著差異，其他農(nóng)藝性狀對(duì)應(yīng)P值均大于0.05，對(duì)產(chǎn)量不存在顯著影響，最后求出回歸方程為Y=-977.243+6.534X4+293.722X5+88.355X6。分析結(jié)果的決定系數(shù)R2為0.615，模型的擬合度較好，說(shuō)明61.5%的產(chǎn)量構(gòu)成可以由有效角果數(shù)、千粒重、每角粒數(shù)來(lái)解釋。即單株有效角果數(shù)每增加1個(gè)單位（1個(gè)），單產(chǎn)增加6.534 kg/hm2，千粒重每增加一個(gè)單位（1 g），單產(chǎn)增加293.722 kg/hm2，每角粒數(shù)每增加1個(gè)單位（1粒），單產(chǎn)增加88.355 kg/hm2。

2.4 ?油菜產(chǎn)量與農(nóng)藝性狀的相關(guān)性分析

由于油菜產(chǎn)量構(gòu)成因素之間的相互影響，造成各因素對(duì)產(chǎn)量的直接作用不同程度地被間接作用所消長(zhǎng)，對(duì)經(jīng)消長(zhǎng)之后的各性狀與產(chǎn)量之間的作用進(jìn)行相關(guān)性分析（表6），進(jìn)一步了解各性狀對(duì)產(chǎn)量貢獻(xiàn)的大小及各性狀之間的相互促進(jìn)和制約關(guān)系。

從表6可以看出，油菜產(chǎn)量與各農(nóng)藝性狀之間均呈正相關(guān)。產(chǎn)量與每角粒數(shù)（r=0.444）及生育期（r=0.399）之間均為顯著正相關(guān);產(chǎn)量與株高、分枝部位、分枝數(shù)、有效角果數(shù)之間存在正相關(guān)，但相關(guān)性不顯著，與千粒重相關(guān)系數(shù)趨近于0（r=0.044），說(shuō)明千粒重對(duì)產(chǎn)量無(wú)直接影響。各農(nóng)藝性狀中，分枝部位與株高之間（r=0.846）為極顯著正相關(guān)，生育期與株高（r=0.573）及分枝部位（r=0.625）之間均為極顯著正相關(guān)，有效角果數(shù)與分枝數(shù)（r=0.406）之間為顯著正相關(guān)，千粒重與分枝數(shù)（r=-0.497）及有效角果數(shù)（r=-0.507）之間均為極顯著負(fù)相關(guān)。對(duì)產(chǎn)量構(gòu)成貢獻(xiàn)最大的每角粒數(shù)與其他農(nóng)藝性狀之間無(wú)顯著的相關(guān)性，但也可以看出株高、分枝數(shù)及生育期對(duì)每角粒數(shù)的構(gòu)成有一定的促進(jìn)作用，分枝部位、有效角果數(shù)及千粒重對(duì)其有一定的抑制作用。從分析結(jié)果看，湖北地區(qū)油菜品種（系）在當(dāng)前的栽培水平下，注重品種（系）每角粒數(shù)的提高及生育期的延長(zhǎng)能夠有效提高品種（系）的產(chǎn)量，但生育期性狀的改良難度大且會(huì)影響下季作物的茬口，不易實(shí)現(xiàn);而每角粒數(shù)與其他農(nóng)藝性狀之間存在較復(fù)雜的相互關(guān)系，在育種過(guò)程中需權(quán)衡各性狀之間的關(guān)系，謹(jǐn)慎選擇。

2.5 ?油菜產(chǎn)量與農(nóng)藝性狀的通徑分析

相關(guān)系數(shù)雖然能夠在一定程度上反映各性狀與產(chǎn)量之間的相關(guān)程度，但要理清各性狀對(duì)產(chǎn)量的作用大小還需進(jìn)行通徑分析。通徑分析可以將相關(guān)系數(shù)分解為直接作用和間接作用，從而清楚地顯示各因素的相互關(guān)系大小和相對(duì)重要性。直接作用凸顯各性狀對(duì)產(chǎn)量作用的事實(shí)，間接作用則揭示性狀間的促抑關(guān)系[17]。

由表7可以看出，分枝部位、分枝數(shù)、有效角果數(shù)、千粒重、每角粒數(shù)對(duì)產(chǎn)量構(gòu)成均存在直接的正效應(yīng)，其中每角粒數(shù)對(duì)產(chǎn)量影響效果最大，這與簡(jiǎn)單相關(guān)性分析結(jié)果一致;而株高和生育期對(duì)產(chǎn)量的直接影響為負(fù)效應(yīng)，說(shuō)明株高和生育期并不能直接促進(jìn)油菜增產(chǎn)，而是通過(guò)影響其他性狀進(jìn)而影響產(chǎn)量。從間接通徑系數(shù)可以看出，與株型相關(guān)的農(nóng)藝性狀（株高、分枝部位、分枝數(shù)）對(duì)產(chǎn)量構(gòu)成三因素的間接影響各有利弊，但影響程度均比較?。ń^對(duì)值最大為0.117、最小為0.002），而產(chǎn)量構(gòu)成三因素相互之間的間接影響則比較明顯，且均呈負(fù)值，這說(shuō)明產(chǎn)量構(gòu)成三因素在產(chǎn)量形成過(guò)程中存在明顯的相互制約關(guān)系。

與株型相關(guān)的農(nóng)藝性狀中，分枝部位對(duì)產(chǎn)量的直接影響作用較大（直接通徑系數(shù)0.535），其他性狀除有效角果數(shù)之外（間接通徑系數(shù)0.069）通過(guò)分枝部位對(duì)產(chǎn)量的影響均為負(fù)值（間接通徑系數(shù)分別是-0.122、-0.020、-0.002、-0.157和-0.034），說(shuō)明在一定程度內(nèi)提高分枝部位可以促進(jìn)有效角果的生長(zhǎng);產(chǎn)量構(gòu)成三因素通過(guò)分枝數(shù)對(duì)產(chǎn)量造成間接影響的效果較明顯，其中有效角果數(shù)（間接通徑系數(shù)0.210）與每角粒數(shù)（間接通徑系數(shù)0.200）通過(guò)分枝數(shù)對(duì)產(chǎn)量有較明顯的正面影響，而千粒重（間接通徑系數(shù)-0.274）通過(guò)分枝數(shù)對(duì)產(chǎn)量則造成較明顯的負(fù)面影響，說(shuō)明有效角果數(shù)與每角粒數(shù)的積累依賴于更多的分枝數(shù)，但分枝數(shù)過(guò)多不利于千粒重增加;而株高對(duì)產(chǎn)量的間接影響因素中，以分枝部位（間接通徑系數(shù)0.452）對(duì)產(chǎn)量的影響最為有效。

3 ?結(jié)論與討論

通過(guò)對(duì)2016—2017湖北省油菜區(qū)域試驗(yàn)匯總數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，發(fā)現(xiàn)油菜產(chǎn)量的形成與各農(nóng)藝性狀之間存在復(fù)雜的關(guān)系。油菜產(chǎn)量構(gòu)成三因素的變異系數(shù)均較大，這為湖北省油菜品種潛力挖掘提供了可能;分枝部位、分枝數(shù)變異系數(shù)較大，株高變異系數(shù)相對(duì)較小，說(shuō)明在油菜理想株型的選擇方面雖然受限于株高但仍具有較大的空間。通過(guò)主成分分析將油菜生產(chǎn)的相關(guān)性狀信息簡(jiǎn)化為生育期因子、分枝數(shù)因子和產(chǎn)量因子，以便于對(duì)各性狀進(jìn)行歸類，并將油菜產(chǎn)量的主要影響因素鎖定為每角粒數(shù)。在此基礎(chǔ)上通過(guò)多重線性回歸分析構(gòu)建湖北省油菜產(chǎn)量預(yù)測(cè)模型Y=-977.243+6.534X4+293.722X5+88.355X6（Y為產(chǎn)量，X4為有效角果數(shù)、X5為千粒重、X6為每角粒數(shù)），確定了產(chǎn)量構(gòu)成三要素與產(chǎn)量之間的線性關(guān)系。隨后通過(guò)相關(guān)性分析和通徑分析對(duì)影響油菜產(chǎn)量各因素之間的促進(jìn)和抑制關(guān)系進(jìn)行解析，發(fā)現(xiàn)各農(nóng)藝性狀除株高和生育期外均對(duì)油菜產(chǎn)量有直接的正面影響;另外與株型有關(guān)的農(nóng)藝性狀中，分枝部位對(duì)有效角果數(shù)有促進(jìn)作用，分枝數(shù)對(duì)有效角果數(shù)和每角粒數(shù)有促進(jìn)作用，但會(huì)抑制千粒重;而產(chǎn)量構(gòu)成三因素之間則存在明顯的相互抑制作用，但對(duì)產(chǎn)量貢獻(xiàn)最大的性狀為每角粒數(shù)。

綜合上述分析結(jié)果，湖北省在油菜高產(chǎn)品種選育過(guò)程中應(yīng)當(dāng)注意兼顧產(chǎn)量構(gòu)成三因素的協(xié)調(diào)發(fā)展，可側(cè)重于選育每角粒數(shù)偏多的品種以獲得高產(chǎn)。而在理想株型的選擇方面，可以在合適生育期和株高的前提下，選擇分枝部位相對(duì)較高、分枝數(shù)偏多的油菜品種，以提高角果面積指數(shù)，充分發(fā)揮光能利用效率，提高產(chǎn)量。但此株型在生長(zhǎng)后期存在一定的倒伏風(fēng)險(xiǎn)，在育種時(shí)需注重品種的抗倒伏特性。

參考文獻(xiàn)：

[1] 陳 ?靜.湖北省油菜產(chǎn)業(yè)發(fā)展問(wèn)題研究[D].武漢：華中農(nóng)業(yè)大學(xué)，2014.

[2] 劉小林.基于因子分析的湖北省油菜生產(chǎn)影響因素研究[J].農(nóng)村經(jīng)濟(jì)與科技，2015，26（12）：15-16.

[3] 李秀萍，杜德志，張冬曉，等.國(guó)家春油菜新品種區(qū)域試驗(yàn)10年研究進(jìn)展[J].中國(guó)種業(yè)，2012（12）：55-56.

[4] 張 ?毅.我國(guó)冬油菜區(qū)域試驗(yàn)品種的高產(chǎn)穩(wěn)產(chǎn)和適應(yīng)性分析[J].中國(guó)油料作物學(xué)報(bào)，2018，40（3）：54-61.

[5] 胡海珍，盧華平.甘藍(lán)型油菜數(shù)量性狀與產(chǎn)量性狀的相關(guān)分析與通徑分析[J].安徽農(nóng)業(yè)科學(xué)，2014，42（19）：6195-6197.

[6] 張培竹，彭少丹，馬忠琴，等.云南小油菜主要農(nóng)藝性狀的相關(guān)及通徑分析[J].甘肅農(nóng)業(yè)科技，2009（8）：12-15.

[7] 張振蘭，李永紅，鄭 ?磊，等.黃淮區(qū)試試驗(yàn)點(diǎn)油菜主要農(nóng)藝性狀分析[J].中國(guó)農(nóng)學(xué)通報(bào)，2018，34（27）：51-55.

[8] 李宏軍，楊 ?鴻，朱傳霞，等.湖南省油菜品種農(nóng)藝性狀分析[J].作物雜志，2015（3）：41-44.

[9] 胡學(xué)聰，劉道敏.冬油菜區(qū)試品種性狀分析[J].安徽農(nóng)業(yè)科學(xué)，2018，46（22）：28-30.

[10] 吳桂仙，楊臘梅，曾 ?林，等.不同品種和密度對(duì)山地油菜產(chǎn)量及其構(gòu)成因素的影響[J].安徽農(nóng)業(yè)科學(xué)，2018，46（10）：56-59.

[11] 宋 ?稀，劉鳳蘭，鄭普英，等.高密度種植專用油菜重要農(nóng)藝性狀與產(chǎn)量的關(guān)系分析[J].中國(guó)農(nóng)業(yè)科學(xué)，2010，43（9）：1800-1806.

[12] 張 ?芳，趙永國(guó)，谷鐵城，等.2001-2010年國(guó)家審定冬油菜品種的產(chǎn)量與主要性狀分析[J].中國(guó)油料作物學(xué)報(bào)，2012，34（3）：239-244.

[13] SONG X，LIU F L，ZHENG P Y，et al. Correlation analysis between agronomic traits and yield of rapeseed（Brassica napus L.） for high-density planting[J].Scientia agricultura sinica，2010， 43（9）：1800-1806.

[14] 郭修遠(yuǎn)，賴玉嫦，段 ?寧，等.不同甘藍(lán)型雜交油菜產(chǎn)量構(gòu)成因素與產(chǎn)量的相關(guān)分析[J].廣東農(nóng)業(yè)科學(xué)，2014，41（8）：39-42.

[15] 崔嘉成.雜交油菜優(yōu)異親本重要農(nóng)藝性狀的遺傳分析[D].北京：中國(guó)農(nóng)業(yè)科學(xué)院，2013.

[16] 趙 ?倩，姜鴻明，孫美芝，等.山東省區(qū)試小麥產(chǎn)量與產(chǎn)量構(gòu)成因素的相關(guān)和通徑分析[J].中國(guó)農(nóng)學(xué)通報(bào)，2011，27（7）：42-45.

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