AA肉雞生產指標間的相關回歸分析

馬發順，秦星星，張凌源，荊學毅，元雪湞，梁秀麗

(1. 安陽工學院生物與食品工程學院，河南 安陽 455000；2. 河南省獸用生物制品研發與應用國際聯合實驗室，河南 安陽 455000；3. 河南省獸藥飼料監察所，河南 鄭州 450000)

AA肉雞是由美國艾拔益加種雞公司育成的白羽肉用四系配套雜交雞，具有生長發育快、出欄早、適應性強、飼養效益高、屠體品質好等特點[1]。目前有關AA肉雞生產性能的研究，主要針對的是不同飼養方式、不同光照制度和不同能量水平日糧等對生產性能和屠宰性能的影響[1-3]。楊志剛等[4]擬合了0～8周齡AA肉雞的生長曲線。馬發順等[5]建立了1～6周齡AA肉雞的生長模型。王剛等[6]運用主成分分析法對AA肉雞肉質性狀進行了綜合評價。相關分析是研究2個或2個以上變量間相關關系的統計分析方法，而典型相關分析是利用綜合變量之間的相關關系來反映2組變量之間整體相關性的多元統計分析方法，它可以識別并量化2組變量之間的聯系[7]。關于AA肉雞生產性能間的相關研究，龔爭等[8]進行了AA肉雞屠宰性狀與肌肉品質的簡單相關分析和典型相關分析。程郁昕等[9]運用回歸分析法對AA肉雞的胸肌、腿肌重等進行了活體估測。馬發順等[10]對AA肉雞的11個生產指標進行了簡單相關分析，并對其中的部分指標建立了回歸方程。關于AA肉雞多個生產指標間的典型相關分析目前未見報道，也沒有在典型相關分析基礎上進行回歸分析的報道。因此，本研究擬采用簡單相關、典型相關和回歸分析的方法對AA肉雞的9項生產指標進行統計分析，以揭示這些生產指標間的內在聯系，并用回歸方程定量描述這些指標之間的關系，從而進行生產水平的快速綜合評價，為飼養管理的科學化、精細化提供參考依據。

1 材料與方法

1.1 試驗動物與飼養管理

以河南省滑縣某養殖場6棟雞舍16.524萬只AA肉雞為試驗動物。該肉雞場采用厚墊料地面平養，0～3周齡人工喂料，3周齡后自動化喂料，自由飲水；按AA肉雞管理規程組織飼養管理和免疫接種。0～3周齡日糧：組成為玉米55.2%、豆粕32%、魚粉2%、菜粕4%、磷酸氫鈣1.5%、石粉1%、食鹽0.3%、油3%、添加劑1%，營養水平為代謝能12.93 MJ/kg，粗蛋白22.37%，鈣1.08%，有效磷0.49%。4～5周齡日糧：組成為玉米57.7%、豆粕27%、魚粉2%、菜粕4%、棉粕3%、磷酸氫鈣1.3%、石粉1.2%、食鹽0.3%、油2.5%、添加劑1%，營養水平為代謝能13.14 MJ/kg，粗蛋白20.11%，鈣0.96%，有效磷0.44%。

1.2 生產指標的測定與計算

按周齡記錄整理每棟雞舍的原始數據，獲得9項生產指標。每棟雞舍入欄數為27 540只。

各項生產指標：周齡(x1)為從出殼之日起按周計齡(周)，存欄數(x2)為1周內每日存欄數的平均值(只)，成活率(x3)即(存欄數/入欄數)×100%，周耗料量(x4)為本周累計耗料量與上周累計耗料量之差(kg)，累計只耗料量(x5)為每只AA肉雞耗料量的累計量(g)，體重(x6)為周末空腹重(g)，周增重(x7)為本周平均體重與上周平均體重之差(g)，周效率(x8)為每周飼料轉化率，即x2·x7/(1 000·x4)×100%，累計效率(x9)為從出殼至報告期累計飼料轉化率，即x6/x5×100%。

1.3 數據處理與統計分析

2 結果與分析

2.1 AA肉雞生產指標的描述統計

AA肉雞1～5周齡各項生產指標描述統計結果見表1。從表1可見，隨著周齡增長，存欄數(x2)和成活率(x3)遞減，而變異程度遞增；隨著周耗料量(x4)的增加，累計只耗料量(x5)也呈增長趨勢，二者變異程度在4周齡之前呈增大趨勢，5周齡時均有所下降；體重(x6)呈上升趨勢，周增重(x7)在4周齡前呈上升趨勢，5周齡時有所下降，二者的變異程度均存在波動性，并且波動走勢一致；周效率(x8)變化趨勢與周增重(x7)一致，變異程度在4周齡前呈上升趨勢，4周齡后下降；累計效率(x9)在4周齡前呈上升趨勢，4周齡后下降，略有波動，其變異程度變化與指標變化趨勢基本一致。

表1 AA肉雞生產指標描述統計結果

2.2 AA肉雞生產指標間的簡單相關

9項生產指標間的簡單相關系數見表2。

從表2可知，36個簡單相關系數中有30個達到極顯著水平(P<0.01)，有6個不顯著(P<0.05)；在極顯著水平的相關關系中有16個為正相關，有14個為負相關。根據生產實際情況，把9項生產指標劃分為因類指標(x1、x2、x3、x4、x5)和果類指標(x6、x7、x8、x9)；從因類指標內部來看，x1與x4、x5之間，x2與x3之間，x4與x5之間，均存在正相關關系(P<0.01)，其余指標之間均為負相關關系(P<0.01)；從果類指標內部來看，x6與x7之間，x8與x9之間，均存在正相關關系(P<0.01)；從因類指標與果類指標之間來看，x6、x7與x1、x4、x5之間均有正相關關系(P<0.01)，x6、x7與x2、x3之間均有負相關關系(P<0.01)，x8、x9與x2、x3之間均有正相關關系(P<0.01)，x8、x9與x4、x5之間均有負相關關系(P<0.01)。

2.3 AA肉雞生產指標間的典型相關

AA肉雞9項生產指標按照因類指標(x1、x2、x3、x4、x5)和果類指標(x6、x7、x8、x9)進行典型相關分析，獲得的典型相關系數及顯著性，詳見表3。

表3 AA肉雞9項生產指標間的典型相關系數及顯著性

從表3可知，4對典型變量中前3對典型變量的相關系數均達到極顯著水平(P<0.01)，第4對僅達到顯著水平(P<0.05)，因此選取前3對典型變量較為合適。AA肉雞的9項生產指標的變異能被3對達到極顯著水平的典型變量所解釋的比例見表4。

表4 AA肉雞各項生產指標的變異能被3對典型變量解釋的比例

從表4可見，5個因類指標的變異能被典型變量U1、U2、U3所解釋的比例分別為0.690 6、0.147 1、0.070 4，能被相對的典型變量V1、V2、V3所解釋的比例分別為0.685 1、0.143 1、0.031 5；4個果類指標的變異能被典型變量V1、V2、V3所解釋的比例分別為0.446 3、0.481 2、0.053 5，能被相對的典型變量U1、U2、U3所解釋的比例分別為0.442 8、0.467 9、0.023 9。可見，9個生產指標的變異主要是由第1對和第2對典型變量所解釋。3對達到極顯著水平的典型變量的構成為：

U1=0.806 1x1+9.179 9x2-9.064 0x3-0.006 4x4+0.329 4x5，

U2=1.515 7x1+0.352 3x2-0.440 6x3-1.581 9x4-0.240 7x5，

U3=-0.916 4x1+65.652 6x2-62.178 7x3-0.882 3x4+4.685 8x5，

V1=0.940 7x6+0.057 1x7-0.040 9x8-0.029 2x9，

V2=-0.020 3x6+0.208 8x7+0.782 2x8+0.172 7x9，

V3=1.943 7x6-2.138 6x7+0.374 6x8+0.319 5x9。

由典型變量表達式可知，U1主要概括了x2、x3，U2主要概括了x1、x4，U3主要概括了x2、x3；V1對x6概括性強，V2對x8概括性強，V3對x6、x7概括性強。所以，U1、U2、U3主要概括了x1、x2、x3、x4，V1、V2、V3主要概括了x6、x8。因類指標和果類指標間強的相關關系主要是x1、x2、x3、x4和x6、x8這些指標間的密切關系所造成。

2.4 AA肉雞生產指標間的回歸分析

AA肉雞的生產指標中與經濟效益密切相關的指標有體重x6、周效率x8和累計效率x9，因此將這3個果類指標作為因變量。x6主要被V1概括，其相對變量U1主要概括了x2、x3，而x6與x1、x4間存在正相關關系(P<0.01)，所以可以x6為因變量，以x1、x3、x4為自變量建立回歸方程。x8被V2概括，其相對變量U2主要概括了x1、x4，而x8與x5間存在負相關關系(P<0.01)，雖然x8與x1間相關不顯著(P>0.05)，但x1、x5間相關極顯著(P<0.01)，所以可以x8為因變量，以x1、x4、x5為自變量建立回歸方程。x9與x2、x8間存在強正相關關系(P<0.01)，與x5間存在負相關關系(P<0.01)，所以可以x9為因變量，以x2、x5、x8為自變量建立回歸方程。所建立的回歸方程為：

x6=1 695.898+455.938x1-20.976x3+0.003x4，

x8=57.012+101.455x1-0.016x4-0.008x5，

x9=-2 855.399+0.107x2+0.019x5+0.572x8。

3個回歸方程的擬合度(R2)分別為0.986、0.968、0.897；對3個回歸方程進行顯著性檢驗，均達到極顯著水平(P<0.01)。

3 討論

從描述統計結果來看，AA肉雞9項生產指標的變異程度存在較大差別，周耗料量(x4)、累計只耗料量(x5)、周效率(x8)和累計效率(x9)變異程度較大，存欄數(x2)、成活率(x3)、體重(x6)、周增重(x7)變異程度較小。可能是雞場規范化的飼養管理使x2、x3較穩定，肉雞的生長發育較為整齊，所以x6、x72項指標也較穩定；可能各棟雞舍并不是以周齡為飼養周期嚴格控制取料量，從而造成各周齡用料量數據誤差較大，這是x4、x5變異程度大的主要原因，進而影響到x8、x9的變異程度。

從簡單相關分析來看，9項生產指標間存在錯綜復雜的相關關系，在達到極顯著的相關系數中有正相關，也有負相關，這與馬發順等[10]的研究結果一致。x6與x1、x4、x5存在正相關關系，與x2、x3存在負相關關系；x8與x2、x3存在正相關關系，與x4、x5存在負相關關系；x9與x2、x3、x8存在正相關關系，與x4、x5存在負相關關系；這些相關關系絕大多數在回歸分析中以回歸系數的表現進一步得到證實，只有x5與x9之間的負相關關系沒有被回歸分析所確認。在x9為因變量的回歸方程中，x5的回歸系數沒有達到顯著水平，所以回歸分析結果并不構成否定x5與x9之間存在負相關關系的理由。x8、x9與x1、x6、x7間相關系數均未達到顯著水平，這與馬發順等[10-11]的研究結果相同或相似，說明無論周效率還是累計效率均與周齡無關，也不受AA肉雞體重及周增重的影響。

典型相關分析顯示，果類指標與因類指標之間的相關性主要由x1、x2、x3、x4和x6、x8之間的密切關系所造成，這與馬發順等[11]主成分分析結果基本一致。結合生產實際綜合分析認為，以x6為因變量選取x1、x3、x4為自變量，以x8為因變量選取x1、x4、x5為自變量，以x9為因變量選取x2、x5、x8為自變量，建立回歸方程是合適的。在建立回歸方程時，為了得到更好的回歸效果，首先采用二次多項式逐步回歸法，結果不收斂；后來采用線性逐步回歸，得到了較為理想的回歸效果。在3個回歸方程中，x1、x3權數較大，而x4、x5權數較小，這與馬發順等[11]建立生產性能評定模型時對各項指標權重的認識有所不同，可能是因為研究的樣本不同和建立回歸方程的方法不同所造成。盡管所建立的3個回歸方程能夠較好地反映AA肉雞生產指標間的定量關系，但是否可以用其他方法建立其他形式的回歸方程，以更好地描述這些生產指標間的定量關系，有待以后進一步研究解決。

4 小結

AA肉雞9項生產指標間存在錯綜復雜的相關關系，在36個簡單相關系數中有30個達到極顯著水平，有6個不顯著；在30個達到極顯著水平的相關系數中有16個為正相關，14個為負相關。周效率和累計效率均不受周齡、體重、周增重的影響。5個因類指標與4個果類指標間3個典型相關系數為0.996 0、 0.986 1、0.668 9；2類指標間的相關關系主要是x1、x2、x3、x4與x6、x8間的密切關系所造成。生產指標間的3個回歸方程為：x6=1 695.898+455.938x1-20.976x3+0.003x4(R2=0.986)；x8=57.012+101.455x1-0.016x4-0.008x5(R2=0.968)；x9=-2 855.399+0.107x2+0.019x5+0.572x8(R2=0.897)。此研究結果可以充分反映AA肉雞9項生產指標之間的定量關系。