馬尾松球果各標(biāo)志值與出種粒數(shù)回歸分析補報

徐浩峰，韋蓉靜，莫周衛(wèi)

（1.都勻市馬寨馬尾松良種基地；2.黔南州林業(yè)科學(xué)研究所，貴州 都勻 558000）

通過多年的生產(chǎn)實踐，馬尾松種子出種率在2%～2.8%之間浮動。影響出種率的因素較多，如樹木的營養(yǎng)狀況、干旱或多雨、球果處理程度等方面，但很多因素不易測定和控制，我們想從較為直觀易測的數(shù)據(jù)如球果的果重、果長、果徑、種鱗數(shù)幾個標(biāo)志值來研究與出種粒數(shù)之間是否有線性關(guān)系。

1 試驗材料、球果處理以及標(biāo)志值的測定方法

1.1 試驗材料

供試的馬尾松球果在2011年采摘于貴州省都勻市馬尾松良種基地1.5代種子園，分別從48個家系中選擇無病蟲害、無機械損傷、大小長短在該家系內(nèi)適中的球果各10個，共選擇球果480個做為總體。

1.2 球果處理以及標(biāo)志值的測定方法

用電子天平對各球果重量稱量，用游標(biāo)卡尺測定果長、果徑，人工計數(shù)種鱗數(shù)，并分別裝入網(wǎng)袋按家系和序列編號。球果經(jīng)石灰水浸泡48h后晾至不滴水，用烘箱燒烤至球果完全開放，種子自然脫落時取出去翅除雜質(zhì)。本次測定的種子要求種粒飽滿、無變質(zhì)、無雜質(zhì)，用電子天平稱重并人工計數(shù)。

2 球果各標(biāo)志值及數(shù)據(jù)統(tǒng)計

2.1 球果各標(biāo)志值

從480組數(shù)據(jù)中隨機抽取96組數(shù)據(jù)做為樣本，抽取率為20%，抽取數(shù)據(jù)如下：

表1 馬尾松球果與種子樣本標(biāo)志值測定

續(xù)表1

續(xù)表1

2.2 數(shù)據(jù)統(tǒng)計

表2 樣本標(biāo)志值數(shù)據(jù)統(tǒng)計

通過表1的數(shù)據(jù)統(tǒng)計出表2，從表2得出樣本各標(biāo)志值的標(biāo)準(zhǔn)差分別為22.530＞16.431＞7.358＞0.902＞0.320；種粒數(shù)＞種鱗數(shù)＞果重＞果長＞果徑，變異系數(shù)比較34.3%＞31.2%＞19.1%＞14.9%＞10.1%；果重＞種粒數(shù)＞種鱗數(shù)＞果長＞果徑，從標(biāo)準(zhǔn)差和變異系數(shù)中得出數(shù)據(jù)分布分散程度最大的是種粒數(shù)、果重、種鱗數(shù)三項標(biāo)志值，種粒數(shù)最大值103、最小值7、極差110；果重最大值44.3、最小值9.8、極差34.5；種鱗數(shù)最大值130、最小值48、極差130，但上述各標(biāo)志值還不能得出與種粒數(shù)有無相關(guān)性的結(jié)論，還要通過以下回歸分析來討論。

3 回歸分析與檢驗

通過把種粒數(shù)做為因變量，其它標(biāo)志值做為自變量處理，組合成果重－種粒數(shù)、果長－種粒數(shù)、果徑－種粒數(shù)、種鱗數(shù)－種粒數(shù)共四組數(shù)據(jù)進行一元線性回歸與檢驗得出以下分析：

表3 一元回歸統(tǒng)計表

表3中四組數(shù)據(jù)一元回歸的相關(guān)系數(shù)與表值Ra比較得0.530＞0.379＞R0.01＝0.262＞0.227＞R0.05＝0.201＞0.147；即果徑－種粒數(shù)R值＞果重－種粒數(shù)R值＞種鱗數(shù)－種粒數(shù)R值＞表值R0.05＞果長－種粒數(shù)R值，果徑、果重與種粒數(shù)的直線回歸相關(guān)系數(shù)極顯著，種鱗數(shù)－種粒數(shù)的直線回歸相關(guān)系數(shù)在表值R0.01和表值R0.05之間相關(guān)系數(shù)顯著，它們所配合的樣本直線回歸方程有意義，果長－種粒數(shù)的直線回歸相關(guān)系數(shù)小于表值R0.05，表明回歸相關(guān)系數(shù)不顯著，樣本直線回歸方 程無線性意義。

表4 方差分析表

表4方差分析表中通過方差分析得到F檢驗值，所得結(jié)論與回歸統(tǒng)計表中的相關(guān)系數(shù)檢驗結(jié)果一致，直線回歸關(guān)系果徑－種粒數(shù)、果重－種粒數(shù)呈極顯著水平，種鱗數(shù)－種粒數(shù)呈顯著水平，直線回歸方程有意義，果長－種粒數(shù)直線回歸關(guān)系不顯著，直線回歸方程無意義。

表5 一元回歸參數(shù)表

由表5一元回歸參數(shù)表得出，三組數(shù)據(jù)回歸方程的回歸常數(shù)和回歸系數(shù)，因表3、表4中果長－種粒數(shù)回歸關(guān)系不顯著，也就不討論該項的回歸參數(shù)，通過表5得回歸方程如下：果重－種粒數(shù)回歸方程：＝47.2975＋1.1634xi；果徑－種粒數(shù)回歸方程＝－45.9152＋37.3219xi；種鱗數(shù) －種粒數(shù)回歸方程＝45.5129＋0.3116xi

該三組回歸方程對估計種粒數(shù)都有意義，前兩個方程的置信度更高，能構(gòu)成一元回歸線性關(guān)系說明果重、果徑、種鱗數(shù)對種粒數(shù)的正相關(guān)，三組數(shù)據(jù)越大種粒數(shù)也越大。但因變量（種粒數(shù)）與自變量（果重、果徑、種鱗數(shù)）的是否有相互關(guān)系通過一元直線回歸表現(xiàn)出一定的相關(guān)性，取得兩變量之間相互關(guān)系的直線近似表達式，但只是一一對應(yīng)的變量之間的關(guān)系式，較為片面不能反映多因素作用對某變量的影響情況，多變量的共同作用與因變量的相關(guān)性討論應(yīng)要用多元回歸解決，我們還是用果重、果長、果徑、種鱗數(shù)為自變量，種粒數(shù)為因變量進行分析，在分析過程中再刪除t檢驗值不顯著的自變量再進行回歸分析，以免影響回歸的顯著性，分析如下：

對表6的回歸參數(shù)表進行分析發(fā)現(xiàn)問題，截距、種鱗數(shù)的回歸分析的t檢驗值（多元回歸t值檢驗更具能反映該自變量的顯著程度）達不到顯著水平，對種鱗數(shù)做數(shù)據(jù)刪除再進行回歸分析，得出表7的結(jié)果，截距和果重也達不到t檢驗值的顯著水平，對果重數(shù)據(jù)進行刪除再進行回歸分析，果長與果徑相互作用對種粒數(shù)有極顯著的線性回歸關(guān)系，根據(jù)表8得出方程為＝－44.374－8.808xi（果長值）＋53.619xi（果徑值），并且該方程有線性回歸意義。

表6 回歸參數(shù)表

表8 回歸參數(shù)表

4 結(jié)論

4.1 通過多元回歸分析，發(fā)現(xiàn)果長和果徑直接影響這些標(biāo)志值相互作用的因素。一元回歸不能對總體標(biāo)志值相互對出粒數(shù)作用進行全面分析，還是以多元回歸分析更具有代表性。

4.2 對各標(biāo)志值的分析，得出馬尾松家系球果果長短與果徑大小，呈顯著線性相關(guān)，直接影響出種粒數(shù)。下一步要進行增加家系、營養(yǎng)狀況、溫度、積溫、采摘時間等方面研究。

［1］ 畢慶雨.數(shù)理統(tǒng)計［M］.北京：中國林業(yè)出版社，1994.