AA肉雞生產(chǎn)指標(biāo)間的相關(guān)回歸分析

馬發(fā)順，秦星星，張凌源，荊學(xué)毅，元雪湞，梁秀麗

(1. 安陽工學(xué)院生物與食品工程學(xué)院，河南 安陽 455000；2. 河南省獸用生物制品研發(fā)與應(yīng)用國際聯(lián)合實(shí)驗(yàn)室，河南 安陽 455000；3. 河南省獸藥飼料監(jiān)察所，河南 鄭州 450000)

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

1 材料與方法

1.1 試驗(yàn)動(dòng)物與飼養(yǎng)管理

以河南省滑縣某養(yǎng)殖場6棟雞舍16.524萬只AA肉雞為試驗(yàn)動(dòng)物。該肉雞場采用厚墊料地面平養(yǎng)，0～3周齡人工喂料，3周齡后自動(dòng)化喂料，自由飲水；按AA肉雞管理規(guī)程組織飼養(yǎng)管理和免疫接種。0～3周齡日糧：組成為玉米55.2%、豆粕32%、魚粉2%、菜粕4%、磷酸氫鈣1.5%、石粉1%、食鹽0.3%、油3%、添加劑1%，營養(yǎng)水平為代謝能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%，營養(yǎng)水平為代謝能13.14 MJ/kg，粗蛋白20.11%，鈣0.96%，有效磷0.44%。

1.2 生產(chǎn)指標(biāo)的測定與計(jì)算

按周齡記錄整理每棟雞舍的原始數(shù)據(jù)，獲得9項(xiàng)生產(chǎn)指標(biāo)。每棟雞舍入欄數(shù)為27 540只。

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

1.3 數(shù)據(jù)處理與統(tǒng)計(jì)分析

2 結(jié)果與分析

2.1 AA肉雞生產(chǎn)指標(biāo)的描述統(tǒng)計(jì)

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

表1 AA肉雞生產(chǎn)指標(biāo)描述統(tǒng)計(jì)結(jié)果

2.2 AA肉雞生產(chǎn)指標(biāo)間的簡單相關(guān)

9項(xiàng)生產(chǎn)指標(biāo)間的簡單相關(guān)系數(shù)見表2。

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

2.3 AA肉雞生產(chǎn)指標(biāo)間的典型相關(guān)

AA肉雞9項(xiàng)生產(chǎn)指標(biāo)按照因類指標(biāo)(x1、x2、x3、x4、x5)和果類指標(biāo)(x6、x7、x8、x9)進(jìn)行典型相關(guān)分析，獲得的典型相關(guān)系數(shù)及顯著性，詳見表3。

表3 AA肉雞9項(xiàng)生產(chǎn)指標(biāo)間的典型相關(guān)系數(shù)及顯著性

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

表4 AA肉雞各項(xiàng)生產(chǎn)指標(biāo)的變異能被3對典型變量解釋的比例

從表4可見，5個(gè)因類指標(biāo)的變異能被典型變量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個(gè)果類指標(biāo)的變異能被典型變量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個(gè)生產(chǎn)指標(biāo)的變異主要是由第1對和第2對典型變量所解釋。3對達(dá)到極顯著水平的典型變量的構(gòu)成為：

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。

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

2.4 AA肉雞生產(chǎn)指標(biāo)間的回歸分析

AA肉雞的生產(chǎn)指標(biāo)中與經(jīng)濟(jì)效益密切相關(guān)的指標(biāo)有體重x6、周效率x8和累計(jì)效率x9，因此將這3個(gè)果類指標(biāo)作為因變量。x6主要被V1概括，其相對變量U1主要概括了x2、x3，而x6與x1、x4間存在正相關(guān)關(guān)系(P<0.01)，所以可以x6為因變量，以x1、x3、x4為自變量建立回歸方程。x8被V2概括，其相對變量U2主要概括了x1、x4，而x8與x5間存在負(fù)相關(guān)關(guān)系(P<0.01)，雖然x8與x1間相關(guān)不顯著(P>0.05)，但x1、x5間相關(guān)極顯著(P<0.01)，所以可以x8為因變量，以x1、x4、x5為自變量建立回歸方程。x9與x2、x8間存在強(qiáng)正相關(guān)關(guān)系(P<0.01)，與x5間存在負(fù)相關(guān)關(guān)系(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個(gè)回歸方程的擬合度(R2)分別為0.986、0.968、0.897；對3個(gè)回歸方程進(jìn)行顯著性檢驗(yàn)，均達(dá)到極顯著水平(P<0.01)。

3 討論

從描述統(tǒng)計(jì)結(jié)果來看，AA肉雞9項(xiàng)生產(chǎn)指標(biāo)的變異程度存在較大差別，周耗料量(x4)、累計(jì)只耗料量(x5)、周效率(x8)和累計(jì)效率(x9)變異程度較大，存欄數(shù)(x2)、成活率(x3)、體重(x6)、周增重(x7)變異程度較小?？赡苁请u場規(guī)范化的飼養(yǎng)管理使x2、x3較穩(wěn)定，肉雞的生長發(fā)育較為整齊，所以x6、x72項(xiàng)指標(biāo)也較穩(wěn)定；可能各棟雞舍并不是以周齡為飼養(yǎng)周期嚴(yán)格控制取料量，從而造成各周齡用料量數(shù)據(jù)誤差較大，這是x4、x5變異程度大的主要原因，進(jìn)而影響到x8、x9的變異程度。

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

典型相關(guān)分析顯示，果類指標(biāo)與因類指標(biāo)之間的相關(guān)性主要由x1、x2、x3、x4和x6、x8之間的密切關(guān)系所造成，這與馬發(fā)順等[11]主成分分析結(jié)果基本一致。結(jié)合生產(chǎn)實(shí)際綜合分析認(rèn)為，以x6為因變量選取x1、x3、x4為自變量，以x8為因變量選取x1、x4、x5為自變量，以x9為因變量選取x2、x5、x8為自變量，建立回歸方程是合適的。在建立回歸方程時(shí)，為了得到更好的回歸效果，首先采用二次多項(xiàng)式逐步回歸法，結(jié)果不收斂；后來采用線性逐步回歸，得到了較為理想的回歸效果。在3個(gè)回歸方程中，x1、x3權(quán)數(shù)較大，而x4、x5權(quán)數(shù)較小，這與馬發(fā)順等[11]建立生產(chǎn)性能評(píng)定模型時(shí)對各項(xiàng)指標(biāo)權(quán)重的認(rèn)識(shí)有所不同，可能是因?yàn)檠芯康臉颖静煌徒⒒貧w方程的方法不同所造成。盡管所建立的3個(gè)回歸方程能夠較好地反映AA肉雞生產(chǎn)指標(biāo)間的定量關(guān)系，但是否可以用其他方法建立其他形式的回歸方程，以更好地描述這些生產(chǎn)指標(biāo)間的定量關(guān)系，有待以后進(jìn)一步研究解決。

4 小結(jié)

AA肉雞9項(xiàng)生產(chǎn)指標(biāo)間存在錯(cuò)綜復(fù)雜的相關(guān)關(guān)系，在36個(gè)簡單相關(guān)系數(shù)中有30個(gè)達(dá)到極顯著水平，有6個(gè)不顯著；在30個(gè)達(dá)到極顯著水平的相關(guān)系數(shù)中有16個(gè)為正相關(guān)，14個(gè)為負(fù)相關(guān)。周效率和累計(jì)效率均不受周齡、體重、周增重的影響。5個(gè)因類指標(biāo)與4個(gè)果類指標(biāo)間3個(gè)典型相關(guān)系數(shù)為0.996 0、 0.986 1、0.668 9；2類指標(biāo)間的相關(guān)關(guān)系主要是x1、x2、x3、x4與x6、x8間的密切關(guān)系所造成。生產(chǎn)指標(biāo)間的3個(gè)回歸方程為：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)。此研究結(jié)果可以充分反映AA肉雞9項(xiàng)生產(chǎn)指標(biāo)之間的定量關(guān)系。