張強龍，余忠祥★，閻明毅，吳 森

（1.青海大學畜牧獸醫(yī)科學院，青海 西寧 810016；2.青海省畜牧獸醫(yī)科學院，青海 西寧 810016）

“歐拉羊”是在海拔3000m以上的青藏高原惡劣氣候條件下培育形成的地方良種，對高寒牧區(qū)嚴酷的自然環(huán)境有著很強的適應能力。自2002年在青海省河南縣歐拉羊中心產區(qū)實施本品種選育以來，生產性能得到大幅度提高，已形成較為穩(wěn)定的群體數量，并于2018年被審定列為國家畜禽遺傳資源。在畜牧業(yè)生產中，家畜體尺數據直接反映家畜的體格大小、體軀結構、生長發(fā)育狀況以及各部位之間相對發(fā)育關系，也間接反映家畜個體或群體的生產性能、生理機能、抗病能力和對外界環(huán)境的適應能力等，在育種工作中有著重要的意義。但作為青藏高原寶貴的動物遺傳資源歐拉羊，常測的體尺主要有體高、體長、胸寬、胸深、尻寬和胸圍等，部分體尺之間有很強的相關性，而在各體尺之間的主次關系、主成分及每個主成分的代表體尺方面的分析尚無相關報道。毛學榮、蒲小劍等通過對歐拉羊生長發(fā)育規(guī)律的研究表明，在初生～6月齡是歐拉羊生長發(fā)育最快的階段，以后生長強度呈遞減趨勢，歐拉羊性公、母羊前期生長發(fā)育個體間差異不顯著。鑒于此，本文應用多元統(tǒng)計分析方法，對384只6月齡羔羊進行體尺和體重的聚類和主成分分析，以期選出歐拉羊代表性的測量指標，為歐拉羊今后的早期選種選育提供參考。

1 材料與方法

1.1 供試羊來源

2020年6月，在青海省黃南州河南蒙古族自治縣歐拉羊繁育基地，隨機抽取健康無病、發(fā)育良好的384只6月齡歐拉羊羔羊。

1.2 測定指標

測定的指標包括：體重（X1）、體高（X2）、體長（X3）、胸寬（X4）、胸深（X5）、尻寬（X6）、胸圍（X7）。主要測定工具是電子秤、測杖、軟尺等。

1.3 統(tǒng)計分析

對測量數據應用Microsoft Excel 2007進行整理，并用SPSS 19.0軟件進行聚類與主成分分析，分別得出相關系數矩陣、聚類分析圖、提取主成分的特征值、貢獻率、累計貢獻率及主成分載荷矩陣。

2 結果與分析

2.1 體尺與體重的相關系數矩陣

由表1可知，6月齡羔羊6個所測體尺指標均與體重呈極顯著正相關（P<0.01），其中，胸圍與體重的相關系數最大達到0.586，體高次之為0.505，其他依次為胸寬、體長、尻寬和胸深；其余各體尺指標間都達到了極顯著水平（P<0.01）。因此，可根據各測定指標間的相關性，把相似程度較高的歸為一類。

表1 羔羊體重與各體尺指標的相關系數矩陣

2.2 聚類分析結果

利用系統(tǒng)聚類法，對各所測指標根據相似性進行了歸類，聚類分析的樹狀圖見圖1。由圖1可知，胸寬、尻寬和胸深聚為一類，后又與體重聚為一大類；體長先于體高聚為一類，后又與胸圍聚為一大類；為了進一步定量地確定各指標的貢獻及其相互關系，對這些指標進行了主成分分析。

圖1 羔羊體尺、體重的聚類分析樹狀圖

2.3 主成分分析結果

累積貢獻率代表各復合性狀相對所有復合性狀對遺傳方差貢獻的百分率，特征向量表示對復合性狀貢獻的大小，其絕對值反映了各性狀對該主成分作用的大小和性質。本文根據累計貢獻率大于70%的要求，選取前3個主成分，3個主成分的累計貢獻率為78.176%，具有一定的代表性，其主成分載荷矩陣見表2。經SPSS統(tǒng)計分析，各成分體尺指標的特征向量系數見表3。

表2 各體尺指標的主成分載荷矩陣

表3 各成分體尺指標的特征向量系數

由表2可知，第一主成分特征根為3.592，貢獻率為51.309%，對第一主成分貢獻較大的依次為體重（X1）、胸圍(X4)、體高(X2)、尻寬（X6）、體長(X3)、胸深（X5）、胸寬（X4），可稱為重量因子，說明該6個體尺指標對體重的影響均較大。第二主成分特征根為1.216，貢獻率為17.378%，其主要正相關因素是胸寬（X4）、尻寬（X6），可稱為寬度因子，負相關因素是體高（X2）、胸深（X5），其絕對值較大，說明胸寬和尻寬越寬的羊，體高和胸深較小。第三主成分特征根為0.664，貢獻率為9.489%，其主要正相關因素是體長（X3），可稱為長度因子，負相關因素胸圍（X7）的絕對值較大，說明體長越長的羊，胸圍較小。

由表3可得3個主成分的函數表達式分別為：

3 討論與結論

聚類與主成分分析在其他羊品種上相對報道較多，而在歐拉羊上尚無相關報道，通過對6月齡歐拉羊進行體尺和體重的聚類和主成分分析，旨在為歐拉羊今后的選育、鑒定和評價提供了選育的方向和重點。本文基于6月齡歐拉羊7個主要性狀進行了聚類分析，結果可大致劃分為特征類（體重、胸圍）、寬度類和長度類，這也與主成分分析的結果基本一致。利用主成分分析，提取的3個主成分主要為重量因子、寬度因子和長度因子，3個主成分的累計貢獻率為78.176%，說明在以后歐拉羊的體尺分析中還應考慮其它體尺性狀。

利用聚類與主成分分析對6月齡歐拉羊的體重、體尺指標進行研究是一個有效的探索和嘗試，應用此方法，不僅可以對歐拉羊的性狀指標進行分類，選出代表性的指標，而且對后期開展選種選育及生長發(fā)育體尺指標的測定提供依據。