對(duì)瀘寧雞黑羽烏骨家系生長發(fā)育性狀的聚類分析

熊 斌 ，鄭 璐 ，陳兆賢 ，敖學(xué)成 ，熊應(yīng)龍 ，方元麗

（1.四川省冕寧縣畜牧局，四川 冕寧 615600；2.西昌學(xué)院，四川 西昌 615000；3.四川省涼山州畜牧局，四川 西昌 615000；4.四川省涼山州原生農(nóng)業(yè)綜合開發(fā)有限責(zé)任公司，四川 冕寧 615602）

1 材料與方法

1.1 分析材料 從瀘寧雞擴(kuò)繁群后代中選擇黑羽烏骨雞組成家系繁育群，按公母1∶10比例采用人工授精配種，在相同孵化飼養(yǎng)條件下，孵化出雛后在脫溫舍飼養(yǎng)35d，35 d后轉(zhuǎn)到育成舍飼養(yǎng)到90 d，90 d后采用林間放養(yǎng)，飼養(yǎng)到19周齡，然后隨機(jī)抽樣測(cè)定相關(guān)指標(biāo)作為分析材料。

1.2 測(cè)定內(nèi)容與分析方法 按農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)NY/T 823-2004家禽生產(chǎn)性能名詞術(shù)語和度量統(tǒng)計(jì)方法，測(cè)定體尺體重指標(biāo)：胸深（X1）、胸寬（X2）、龍骨長（X3）、脛長（X4）、盆骨寬（X5）、脛圍（X6）、體斜長（X7）、體重（X8）；利用 Excel和 SPSS軟件，分別進(jìn)行性狀描述分析、相關(guān)分析、主成分分析，獲取主成分值進(jìn)行聚類分析，分類后按樣本分類作類間性狀方差比較分析。

2 結(jié)果與分析

2.1 體重體尺性狀的統(tǒng)計(jì)分析 體重體尺數(shù)據(jù)列表1。從表1可看出各性狀平均值的參數(shù)指標(biāo)與中位數(shù)均較接近，表明性狀均值有代表性。各性狀峰度值在-0.02～0.8之間，偏度值在-0.22～0.51之間，各性狀最大值與最小值經(jīng)處理相差比率在36.10%～183.40%之間，說明各性狀分布離差較大，群體中存在性狀的明顯差異，表型值遺傳差異明顯，存在明顯的可選擇性。

表1 樣本體重體尺性狀的統(tǒng)計(jì)分析 cm、g

2.2 主成分分析

2.2.1 性狀相關(guān)分析。樣本相關(guān)系數(shù)列表2，從表2可看出，除胸寬與脛長的相關(guān)性較弱，龍骨長與脛長的相關(guān)性顯著（P＜0.05）外，其余性狀間的相關(guān)系數(shù)均達(dá)到極顯著水平（P＜0.01）。說明體重與體尺相關(guān)性狀間呈正相關(guān)關(guān)系，選擇體尺大的將會(huì)影響體重性狀和相關(guān)性狀的提高。

表2 樣本相關(guān)系數(shù)

2.2.2 主成分荷載與特征值貢獻(xiàn)率。從表3可知，按不低于85%累計(jì)特征值貢獻(xiàn)率算，前4個(gè)主成分的

表3 主成分特征值及貢獻(xiàn)率

特征值分別是4.810、0.907、0.682和0.557，方差貢獻(xiàn)率分別是60.121%、11.343%、8.529%和6.957%，累計(jì)貢獻(xiàn)率達(dá)到86.949%。主成分分析可以將多指標(biāo)問題轉(zhuǎn)化為較少的新的指標(biāo)，從而獲得許多重要信息，并包含樣本的主要信息。此例提取的4個(gè)主成分已包含了樣本主要信息，可取代樣本的綜合指標(biāo)。

2.2.3 依據(jù)主成分分析原理，可用提取的主成分荷載進(jìn)行主成分因子判斷，主成分荷載系數(shù)詳見表4。

表4 主成分荷載系數(shù)

第一主成分中8個(gè)性狀荷載系數(shù)在0.451～0.923之間，其中最大的荷載系數(shù)是體重（0.923），其次是體斜長（0.884）。當(dāng)體重大時(shí)相應(yīng)的體斜長等指標(biāo)會(huì)更大，體重會(huì)更重，因此第一主成分可定義為體斜長生長發(fā)育因子。

第二主成分中突出較大的荷載系數(shù)是脛長（0.825），表明當(dāng)脛長越長時(shí)，個(gè)體體重會(huì)更大，因此可定義為脛長因子。

第三主成分中較大的荷載系數(shù)是胸寬（0.537），其次是脛圍，當(dāng)胸寬較大時(shí)，胸圍相應(yīng)會(huì)更大，其個(gè)體體重也會(huì)更大，可定義為胸寬因子。

第四主成分中較大的荷載系數(shù)是并列的胸寬（0.362）和脛長（0.310），當(dāng)胸寬和脛長同時(shí)較大時(shí)，將對(duì)個(gè)體體重產(chǎn)生影響。鑒于第二主成分為脛長因子、第三主成分為胸寬因子，可將第四主成分定義為胸寬脛長互補(bǔ)因子。

綜上分析可知，4個(gè)主成分代表的主要因子——體斜長生長發(fā)育因子、脛長因子、胸寬因子和胸寬脛長互補(bǔ)因子，是影響黑羽瀘寧雞的主要性狀因素。

2.3 聚類分析 利用所獲得的主成分值進(jìn)行聚類分析（見聚類成員表5），可將樣本聚為4個(gè)大類，按分類性狀進(jìn)行方差分析（結(jié)果列表6）。從表6可看出，各類間性狀中，龍骨長、脛圍類間差異不顯著（P＞0.05），胸深、胸寬、脛長、骨盆寬、體斜長類間差異達(dá)到極顯著（P＜0.01），體重類間差異顯著（P＜0.05）。分別命名為第1、2、3因子的性狀——體長、脛長、胸寬類間差異極顯著，說明與影響體重性狀的因子是吻合的。

表5 聚類成員

通過分類性狀排序比較（見表7）可看出，四個(gè)類中1類所有性狀為最小，排序均居4；2類、3類與4類的體重相似（平均僅差1.33g和23.63g），在2類中體斜長排序?yàn)?，表明和提取主成分的因子命名相吻合，因此按2、4類群性狀特點(diǎn)進(jìn)行選擇，個(gè)體體重會(huì)更大，而與此相關(guān)的性狀——體斜長、脛長、胸寬等指標(biāo)將會(huì)更大，可見具有選擇意義。

表6 分類性狀比較分析

表7 分類性狀排序比較

3 結(jié)論

通過對(duì)瀘寧雞黑羽烏骨家系后代進(jìn)行主成分和聚類分析，揭示出了瀘寧黑羽烏骨雞的選育依據(jù)，可初步得出以下結(jié)論：

3.1 利用提取主成分最大荷載提供的命名因子，可知影響瀘寧雞黑羽烏骨家系個(gè)體體重的主要性狀是體斜長、脛長和胸寬。

3.2 主成分分析所提取的4個(gè)主成分，完全包含了樣本的主要信息，可取代樣本的綜合指標(biāo)；以2、4類群性狀的表現(xiàn)特點(diǎn)作為選育依據(jù)，可提高以群體體重為主的綜合性狀。

[1] 楊朝武，朱慶，蔣小松.基于主成分分析的優(yōu)質(zhì)雞肉質(zhì)評(píng)價(jià)研究[J].生物數(shù)學(xué)學(xué)報(bào)，2007，22（3）：547-553.

[2] 章文波，陳紅艷.實(shí)用數(shù)據(jù)統(tǒng)計(jì)分析及SPSS 12.0應(yīng)用[M].北京：人民郵電出版社，2006.

[3] 吳高嶺，朱德鳳.數(shù)值分類法在玉米品種資源上的應(yīng)用[J].玉米科學(xué)，2000，8（1）：23-27.

[4] 白厚義.試驗(yàn)方法及統(tǒng)計(jì)分析[M].北京：中國林業(yè)出版社，2005.

[5] 陳國順.運(yùn)用聚類分析和主成分分析篩選豬的血清指標(biāo)[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào)，2005（6）：723-727.