李金輝 陳佳欣 嚴(yán) 慧 劉 潔 王里彥 田星哲 楊彩虹 劉月琴 紀(jì)守坤* 張英杰*

(1.河北農(nóng)業(yè)大學(xué)動物科技學(xué)院，保定 071000；2.石家莊市農(nóng)林科學(xué)研究院，石家莊 050041)

營養(yǎng)物質(zhì)表觀消化率是評價飼糧營養(yǎng)價值以及動物消化能力的綜合指標(biāo)[1]，是影響肉羊生長性能的重要因素。多種因素如動物品種[2]、飼糧組成[3]、環(huán)境[4]等均會影響肉羊?qū)︼暭Z營養(yǎng)物質(zhì)的消化，因此肉羊在不同飼養(yǎng)條件和生理狀況下營養(yǎng)物質(zhì)表觀消化率差別較大，對肉羊的營養(yǎng)物質(zhì)表觀消化率進(jìn)行準(zhǔn)確估測，是提供與肉羊?qū)嶋H營養(yǎng)需要相匹配的飼糧營養(yǎng)水平，實現(xiàn)肉羊精準(zhǔn)飼養(yǎng)的重要保障。

肉羊飼糧營養(yǎng)物質(zhì)表觀消化率常用體內(nèi)法或體外法進(jìn)行測定[5]，但2種方法操作繁瑣、時效性差，難以對肉羊整體的飼糧營養(yǎng)物質(zhì)表觀消化率進(jìn)行快速、系統(tǒng)地評價，因此前人嘗試通過構(gòu)建預(yù)測模型來解決這一問題。Morrison[6]認(rèn)為營養(yǎng)物質(zhì)表觀消化率主要受飼糧酸性洗滌纖維(ADF)含量的影響；Adesogan等[7]認(rèn)為飼糧粗蛋白質(zhì)(CP)含量能較好地預(yù)測營養(yǎng)物質(zhì)表觀消化率；Nousiainen等[8]發(fā)現(xiàn)飼糧不可消化中性洗滌纖維含量可預(yù)測肉羊的營養(yǎng)物質(zhì)表觀消化率；劉潔等[9]發(fā)現(xiàn)飼糧某營養(yǎng)物質(zhì)表觀消化率與該營養(yǎng)物質(zhì)在飼糧中的含量呈正相關(guān)，與中性洗滌纖維(NDF)含量呈負(fù)相關(guān)；趙江波[10]研究發(fā)現(xiàn)，總能表觀消化率(GED)、粗蛋白質(zhì)表觀消化率(CPD)與飼糧NDF、ADF含量顯著相關(guān)，中性洗滌纖維表觀消化率(NDFD)、酸性洗滌纖維表觀消化率(ADFD)則與飼糧總能(GE)、CP含量顯著相關(guān)。以上研究基于各自研究結(jié)果構(gòu)建了肉羊?qū)︼暭Z營養(yǎng)物質(zhì)表觀消化率的預(yù)測方程，但這些基于試驗測定數(shù)據(jù)構(gòu)建的預(yù)測模型雖準(zhǔn)確性高但均存在普適性差的問題，產(chǎn)生這一現(xiàn)象的原因是試驗研究測定較準(zhǔn)確但普遍存在納入影響因素不足的問題。NRC(2007)[11]和CNCPS(2004)[12]針對該問題采用納入多項前人研究結(jié)果建立數(shù)據(jù)庫，進(jìn)而構(gòu)建預(yù)測模型以提高模型的普適性，本課題組利用該方法也成功建立了舍飼綿羊干物質(zhì)采食量(DMI)的預(yù)測模型并解析了其影響因素[13]，為本研究提供了方法學(xué)的借鑒。

本研究通過搜集國內(nèi)外肉羊飼糧營養(yǎng)物質(zhì)表觀消化率的相關(guān)文獻(xiàn)數(shù)據(jù)，構(gòu)建肉羊飼糧營養(yǎng)物質(zhì)表觀消化率及其影響因素的數(shù)據(jù)庫，使用多元回歸法建立預(yù)測方程，采用方差分解法解析其影響因素，以期建立普適性較強(qiáng)的營養(yǎng)物質(zhì)表觀消化率預(yù)測模型，為我國肉羊飼糧營養(yǎng)物質(zhì)表觀消化率的預(yù)測提供基礎(chǔ)數(shù)據(jù)。

1 材料與方法

1.1 文獻(xiàn)檢索策略

由3名研究人員分別在Web of Science、PubMed、中國知網(wǎng)等文獻(xiàn)數(shù)據(jù)平臺檢索2006年1月1日至2021年1月1日期間發(fā)表的關(guān)于肉羊飼糧營養(yǎng)物質(zhì)表觀消化率的相關(guān)中、英文文獻(xiàn)，中文文獻(xiàn)檢索以“羊”和“表觀消化率”為關(guān)鍵詞，英文文獻(xiàn)檢索分別以“sheep”&“apparent digestibility”和“goat”&“apparent digestibility”為關(guān)鍵詞。

1.2 納入標(biāo)準(zhǔn)

文獻(xiàn)納入標(biāo)準(zhǔn)：1)試驗研究對象為肉羊的文獻(xiàn)；2)飼糧營養(yǎng)成分組成包括GE、代謝能(ME)以及干物質(zhì)(DM)、CP、NDF、ADF、鈣(Ca)和磷(P)含量等基本信息的文獻(xiàn)；3)肉羊生長性能包括體重(BW)、相對體成熟度(RM)、平均日增重(ADG)和DMI等信息的文獻(xiàn)。

1.3 排除標(biāo)準(zhǔn)

文獻(xiàn)排除標(biāo)準(zhǔn)：1)綜述類文獻(xiàn)；2)重復(fù)檢索到的文獻(xiàn)；3)肉羊處于非正常狀態(tài)的文獻(xiàn)，包括患病、使用藥物和限飼等；4)沒有明確肉羊品種的文獻(xiàn)。

1.4 文獻(xiàn)數(shù)據(jù)的篩選、提取和數(shù)據(jù)庫的構(gòu)建

根據(jù)嚴(yán)格的納入和排除標(biāo)準(zhǔn)對檢索到的文獻(xiàn)通讀復(fù)篩后，提取以下信息：1)作者、發(fā)表年份、肉羊的品種、羊只數(shù)、BW、RM、ADG和DMI；2)飼糧ME、DM、CP、NDF、ADF、Ca和P含量等信息。納入文獻(xiàn)無成年體重時通過查閱相關(guān)文獻(xiàn)補齊數(shù)據(jù)，計算肉羊的RM(RM=BW/成年BW)；使用公式計算肉羊的代謝體重(MBW，MBW=BW0.75)。構(gòu)建包含參數(shù)MBW、RM、ADG、DMI，飼糧ME、DM、CP、NDF和ADF含量以及營養(yǎng)物質(zhì)表觀消化率[干物質(zhì)表觀消化率(DMD)、GED、CPD、NDFD和ADFD]在內(nèi)的數(shù)據(jù)庫。

1.5 肉羊飼糧營養(yǎng)物質(zhì)表觀消化率預(yù)測模型的構(gòu)建

使用Excel 2016對構(gòu)建的數(shù)據(jù)庫進(jìn)行整理，進(jìn)一步剔除參數(shù)缺失的數(shù)據(jù)，使用R-4.0.0軟件，通過cor程序分析各指標(biāo)之間的相關(guān)性，P<0.05表示存在顯著相關(guān)性。使用lm函數(shù)以MBW、RM、ADG、DMI以及飼糧DM、ME、CP、NDF和ADF含量為自變量，DMD、GED、CPD、NDFD和ADFD為因變量構(gòu)建營養(yǎng)物質(zhì)表觀消化率預(yù)測模型：

Y=a+b1X1+b2X2+……bnXn。

預(yù)測模型中：Y為因變量；X為自變量(參數(shù))；a為常數(shù)項；b為對應(yīng)參數(shù)的系數(shù)。

使用vif函數(shù)計算模型中各參數(shù)的方差膨脹因子(VIF)，當(dāng)VIF>10認(rèn)為各參數(shù)間存在共線性，兩參數(shù)間存在共線性時，依據(jù)該參數(shù)與因變量之間的相關(guān)性大小去除相關(guān)性較低的參數(shù)。

1.6 肉羊飼糧營養(yǎng)物質(zhì)表觀消化率預(yù)測模型準(zhǔn)確性的評估

使用平均偏差(mean bias,MB)和預(yù)測均方根誤差(root of mean square prediction error,RMSPE)評估營養(yǎng)物質(zhì)表觀消化率預(yù)測模型的預(yù)測偏差。

式中：MB表示平均偏差；n表示預(yù)測值和觀測值的數(shù)量；Pi表示預(yù)測值；Oi表示觀測值；i表示預(yù)測值和觀測值的次序。

式中：RMSPE表示預(yù)測均方根誤差；n表示預(yù)測值和觀測值的數(shù)量；Pi表示預(yù)測值；Oi表示觀測值；i表示預(yù)測值和觀測值的次序。

1.7 方差分解

首先使用scale函數(shù)對納入預(yù)測模型的數(shù)據(jù)進(jìn)行中心化和標(biāo)準(zhǔn)化處理，消除不同參數(shù)量綱的影響，隨后采用多元回歸模型進(jìn)行模型構(gòu)建。依據(jù)方差分解結(jié)果獲得營養(yǎng)物質(zhì)表觀消化率預(yù)測模型的方差解釋率、各個參數(shù)對于預(yù)測結(jié)果的貢獻(xiàn)度以及動物因素(MBW、RM、ADG和DMI)和飼糧因素(ME、DM、CP、NDF和ADF含量)對營養(yǎng)物質(zhì)表觀消化率(DMD、GED、CPD、NDFD和ADFD)的貢獻(xiàn)度。

2 結(jié)果與分析

2.1 數(shù)據(jù)庫的基本特征

本研究共納入文獻(xiàn)99篇，數(shù)據(jù)4 179條，數(shù)據(jù)庫基本特征見表1。所構(gòu)建的數(shù)據(jù)庫中肉羊MBW、RM、ADG和DMI分別為(12.51±2.89) kg、(0.42±0.18)%、(179.33±71.54) g/d和(1 144.31±386.82) g/d；飼糧ME、DM、CP、NDF和ADF含量分別為(9.91±1.86)MJ/kg、(86.44±10.00)%、(13.89±2.94)%、(41.25±11.56)%和(25.28±7.74)%；DMD、GED、CPD、NDFD和ADFD分別為(66.13±8.63)%、(63.81±10.23)%、(64.18±9.31)%、(52.57±12.13)%和(46.11±12.14)%。

表1 數(shù)據(jù)庫的基本特征

2.2 數(shù)據(jù)庫中參數(shù)間相關(guān)性

由圖1可知，所構(gòu)建的數(shù)據(jù)庫中各參數(shù)均服從正態(tài)分布；DMD與ADG、飼糧CP含量顯著正相關(guān)(P<0.05)，與飼糧ME、NDF含量顯著負(fù)相關(guān)(P<0.05)；GED與飼糧CP含量顯著正相關(guān)(P<0.05)，與飼糧ADF含量顯著負(fù)相關(guān)(P<0.05)；CPD與MBW、RM、ADG和飼糧CP含量顯著正相關(guān)(P<0.05)；NDFD與ADG、飼糧CP和NDF含量顯著正相關(guān)(P<0.05)，與RM顯著負(fù)相關(guān)(P<0.05)；ADFD與ADG、飼糧CP和ADF含量顯著正相關(guān)(P<0.05)，與RM顯著負(fù)相關(guān)(P<0.05)。

MBW：代謝體重 metabolic weight；RM：相對體成熟度 relative body maturity；ADG：平均日增重 average daily gain；DMI：干物質(zhì)采食量 dry matter intake；DM：干物質(zhì) dry matter；ME：代謝能 metabolizable energy；CP：粗蛋白質(zhì) crude protein；NDF：中性洗滌纖維 neutral detergent fiber；ADF：酸性洗滌纖維 acid detergent fiber；DMD：干物質(zhì)表觀消化率 dry matter apparent digestibility；GED：總能表觀消化率gross energy apparent digestibility；CPD：粗蛋白質(zhì)表觀消化率crude protein apparent digestibility；NDFD：中性洗滌纖維表觀消化率neutral detergent fiber apparent digestibility；ADFD：酸性洗滌纖維表觀消化率acid detergent fiber apparent digestibility。下圖同 the same as below。

2.3 肉羊飼糧營養(yǎng)物質(zhì)表觀消化率預(yù)測模型

由表2可知，DMD預(yù)測模型、CPD預(yù)測模型、NDFD預(yù)測模型和ADFD預(yù)測模型參數(shù)間均不存在共線性(VIF<10)，模型可分別解釋65.0%、67.0%、51.0%和47.0%的變異(P<0.05)，MB分別為0.332、-0.469、-0.073和0.098個百分點，RMSPE分別為7.1%、8.3%、17.5%和20.2%。構(gòu)建GED預(yù)測模型時，參數(shù)間存在共線性(VIF>10)，依據(jù)參數(shù)與因變量之間的相關(guān)性去除相關(guān)性較低的參數(shù)MBW、ADG、飼糧CP和ADF含量后，剩余參數(shù)間不存在共線性(VIF<10)，該模型可解釋95.0%的GED變異(P<0.05)，MB為6.346個百分點，RMSPE為23.1%。

2.4 肉羊飼糧營養(yǎng)物質(zhì)表觀消化率影響因素的方差分解

由圖2可知，本研究構(gòu)建的DMD預(yù)測模型可解釋肉羊飼糧65.0%的DMD變異，在可解釋變異中，動物因素可解釋66.6%的變異，飼糧因素可解釋33.4%的變異，其中MBW和DMI是影響DMD的最主要因素，分別解釋了35.7%和15.0%的變異；GED預(yù)測模型可解釋肉羊飼糧95.0%的GED變異，在可解釋變異中，動物因素可解釋57.7%的變異，飼糧因素可解釋42.3%的變異，其中RM和飼糧DM含量是影響GED的最主要因素，分別解釋了43.0%和17.5%的變異；CPD預(yù)測模型可解釋肉羊飼糧67.0%的CPD變異，在可解釋變異中，動物因素可解釋65.1%的變異，飼糧因素可解釋34.9%的變異，其中DMI、ADG和MBW是影響CPD的最主要因素，分別解釋了28.7%、18.6%和15.6%的變異；NDFD預(yù)測模型可解釋肉羊飼糧51.0%的NDFD變異，在可解釋變異中，動物因素可解釋43.1%的變異，飼糧因素可解釋56.9%的變異，其中飼糧CP含量、RM、飼糧DM含量和MBW是影響NDFD的最主要因素，分別解釋了21.0%、17.1%、16.8%和10.7%的變異；ADFD預(yù)測模型可解釋肉羊飼糧47.0%的ADFD變異，在可解釋變異中，動物因素可解釋41.9%的變異，飼糧因素可解釋58.1%的變異，其中飼糧CP、DM含量及MBW和RM是影響ADFD的最主要因素，分別解釋了16.8%、14.7%、14.6%和11.4%的變異。

表2 肉羊營養(yǎng)物質(zhì)表觀消化率預(yù)測模型

圖2 肉羊營養(yǎng)物質(zhì)表觀消化率影響因素的方差分解

3 討 論

3.1 肉羊飼糧營養(yǎng)物質(zhì)表觀消化率及其影響因素的數(shù)據(jù)庫

營養(yǎng)物質(zhì)表觀消化率受動物、飼糧、環(huán)境等多種因素的影響[14]，因此二次整合文獻(xiàn)數(shù)據(jù)構(gòu)建數(shù)據(jù)庫，通過系統(tǒng)、全面的數(shù)據(jù)再分析，能得出預(yù)測精度較高且更具普適性的營養(yǎng)物質(zhì)表觀消化率預(yù)測模型，NRC(2007)[11]和CNCPS(2004)[12]均采用這種方法來構(gòu)建數(shù)據(jù)庫進(jìn)行營養(yǎng)參數(shù)預(yù)測。本研究圍繞影響肉羊飼糧營養(yǎng)物質(zhì)表觀消化率的主要因素(動物因素和飼糧因素)構(gòu)建數(shù)據(jù)庫，建立了包含43個肉羊品種、4 179條數(shù)據(jù)的數(shù)據(jù)庫，其中包含國內(nèi)肉羊品種或雜交肉羊品種數(shù)據(jù)3 721條，國外肉羊品種數(shù)據(jù)458條，包括影響DMD、GED、CPD、NDFD和ADFD的潛在因素9個。本研究構(gòu)建的數(shù)據(jù)庫數(shù)據(jù)量多，涉及的影響因素范圍全面，并且大量納入我國的肉羊品種，因此建立的預(yù)測模型更適合預(yù)測我國肉羊的飼糧營養(yǎng)物質(zhì)表觀消化率。

3.2 肉羊飼糧營養(yǎng)物質(zhì)表觀消化率的主要影響因素

DMD可反映動物對飼糧營養(yǎng)物質(zhì)的消化程度。本研究構(gòu)建的數(shù)據(jù)庫中DMD最低為20.7%，最高可達(dá)89.8%，可見不同研究中肉羊飼糧DMD存在較大差別。通過數(shù)據(jù)的匯總分析發(fā)現(xiàn)，肉羊飼糧DMD平均值為66.1%，與大多數(shù)研究結(jié)果[15-16]相近，表明該值可大致反映肉羊整體對飼糧DM的利用情況。本研究通過數(shù)據(jù)庫構(gòu)建的DMD預(yù)測模型R2為0.65，這與劉潔等[9]基于動物試驗構(gòu)建的模型(R2=0.79)相比偏低，是由于動物試驗過程中對試驗條件進(jìn)行了良好的控制，用于構(gòu)建預(yù)測模型的數(shù)據(jù)測定較準(zhǔn)確，因而相較于通過數(shù)據(jù)庫構(gòu)建的預(yù)測模型而言R2更高，但基于數(shù)據(jù)庫構(gòu)建的預(yù)測模型納入的影響因素更多，涉及的范圍更加全面，因此本研究構(gòu)建的預(yù)測模型也更具普適性。本研究進(jìn)一步通過將數(shù)據(jù)庫中納入預(yù)測模型的影響因素方差分解表明，肉羊飼糧DMD主要受MBW和DMI的影響，其中肉羊MBW與DMD顯著正相關(guān)，這可能是隨MBW的增加，肉羊胃腸道發(fā)育更加完善，腸道運動能力、腸道激素釋放和瘤胃發(fā)酵功能增強(qiáng)[17]，提高了肉羊?qū)︼暭ZDM的消化率；肉羊DMI與DMD顯著負(fù)相關(guān)，是由于采食量增加導(dǎo)致消化道流通速率增加，飼糧難以得到充分消化，同時采食量增加使底物增加，消化酶與飼糧接觸不充分消化不徹底，使得DMD降低。

GED指飼糧中可消化利用的能量占飼糧所含能量總量的百分比，可用于評價飼糧的營養(yǎng)價值和肉羊的能量利用效率。前人研究發(fā)現(xiàn)肉羊在不同生長階段和生理狀況下，對不同飼糧的能量利用效率存在較大差別[11]，因此本研究通過匯總文獻(xiàn)數(shù)據(jù)建立的數(shù)據(jù)庫中GED覆蓋范圍較廣(39.8%～86.0%)。數(shù)據(jù)匯總的結(jié)果顯示，肉羊飼糧GED平均值為63.8%，這與趙敏孟等[18]報道的杜泊母羊(GED=62.6%)、張海容[20]報道的藏綿羊(GED=64.4%)和武渝瀚等[20]報道的云南半細(xì)毛羊(GED=63.2%)的數(shù)值相近，表明該結(jié)果可反映肉羊整體對飼糧能量的利用效率。本研究通過數(shù)據(jù)庫構(gòu)建的預(yù)測模型可解釋95.0%的GED變異，這與劉潔等[9]通過動物試驗構(gòu)建的預(yù)測模型(R2=0.96)相近，表明與通過動物試驗構(gòu)建的預(yù)測模型相比，基于數(shù)據(jù)庫建立的預(yù)測模型對肉羊GED同樣具有較高的解釋度，但本研究通過數(shù)據(jù)庫構(gòu)建的模型MB為6.346個百分點，表明存在高估GED的可能，這是由于目前關(guān)于GED預(yù)測模型的研究較少，符合構(gòu)建預(yù)測模型的數(shù)據(jù)量少導(dǎo)致其自由度低所致。前人通過動物試驗構(gòu)建的肉羊飼糧GED預(yù)測模型將飼糧作為主要影響因素：趙江波[10]構(gòu)建的肉用綿羊精料GED的預(yù)測模型認(rèn)為，GED與飼糧CP、NDF和ADF含量顯著負(fù)相關(guān)，與GE顯著正相關(guān)；劉潔等[9]構(gòu)建的肉用綿羊飼糧GED和有效能的預(yù)測模型認(rèn)為，GED與飼糧GE顯著正相關(guān)，與飼糧NDF含量顯著負(fù)相關(guān)。本研究通過數(shù)據(jù)庫將動物因素和飼糧因素相結(jié)合，通過影響因素的方差分解分析表明，飼糧DM含量和肉羊RM是影響GED最主要的因素，其中飼糧DM含量與GED顯著負(fù)相關(guān)，肉羊RM與GED顯著正相關(guān)。

CPD可反映肉羊?qū)︼暭ZCP的利用效率。本研究構(gòu)建的數(shù)據(jù)庫中CPD最低為30.9%，最高為90.0%，這比CNCPS(2004)[12]數(shù)據(jù)庫(CPD：52.4%～83.5%)覆蓋范圍更廣，這是由于CNCPS(2004)[12]主要納入的是生長期和成年期的綿羊，而本研究將山羊和綿羊不同生長階段的CPD均納入數(shù)據(jù)庫。數(shù)據(jù)匯總結(jié)果顯示，肉羊飼糧CPD平均值為64.2%，該結(jié)果低于NRC(2007)[11]中綿羊?qū)︼暭ZCP的利用效率(67.0%)，其原因也與納入數(shù)據(jù)庫的肉羊品種和生長階段有關(guān)。本研究通過數(shù)據(jù)庫構(gòu)建的CPD預(yù)測模型可解釋67.0%的變異，明顯高于CNCPS(2004)[12]基于糞便CP含量構(gòu)建的預(yù)測模型(R2=0.34)。前人研究表明飼糧CP含量、飼喂量及肉羊品種、性別等均會影響飼糧中CP的消化[21]，但主要影響因素尚不明確。本研究通過將數(shù)據(jù)庫中納入預(yù)測模型的影響因素方差分解表明，動物因素(DMI、MBW和ADG)是影響肉羊飼糧CPD的主要因素，其中DMI與CPD顯著負(fù)相關(guān)，這是由于肉羊DMI增加，刺激消化道運動能力增強(qiáng)，食糜排空速度加快，縮短了飼糧在消化道內(nèi)的滯留時間，使其與蛋白酶作用不充分[22]，因而CPD降低。肉羊MBW、ADG與CPD顯著正相關(guān)，目前關(guān)于體重和日增重對CPD影響的研究未見報道，但前人研究發(fā)現(xiàn)幼齡動物隨年齡增長，消化道器官發(fā)育逐漸完善，消化酶的生成量和活性逐漸達(dá)到成年水平，對飼糧粗脂肪、CP和粗纖維的消化率提高[23]，因此羔羊CPD、MBW和ADG均會同步增加。

飼糧纖維物質(zhì)是反芻動物能量的重要來源之一，NDFD和ADFD可用來評價瘤胃發(fā)酵能力和飼糧營養(yǎng)價值。本研究通過搜集文獻(xiàn)構(gòu)建的數(shù)據(jù)庫中肉羊飼糧NDFD平均值(52.8 %)與CNCPS(2004)[12]數(shù)據(jù)庫(平均值：54.6%)相近，通過動物因素和飼糧因素建立的NDFD預(yù)測模型R2=0.51，與CNCPS(2004)[12]構(gòu)建的預(yù)測模型(R2=0.51)相同，表明本研究構(gòu)建的肉羊飼糧NDFD數(shù)據(jù)庫與CNCPS(2004)[12]數(shù)據(jù)庫相比，在納入大量國內(nèi)肉羊品種后具有相近的預(yù)測效果。目前關(guān)于肉羊飼糧ADFD預(yù)測模型的研究較少，本研究通過搜集匯總?cè)庋虿煌L階段飼糧ADFD的數(shù)據(jù)構(gòu)建了數(shù)據(jù)庫，結(jié)果顯示肉羊飼糧ADFD平均值為46.1%，數(shù)據(jù)的覆蓋范圍為9.4%～74.7%，可見不同研究中肉羊飼糧ADFD存在較大差別。本研究通過數(shù)據(jù)庫構(gòu)建的ADFD預(yù)測模型可解釋47.0%的變異，且模型MB(0.098個百分點)較小，可為我國肉羊飼糧ADFD的預(yù)測提供參考。本研究通過將數(shù)據(jù)庫中納入預(yù)測模型的影響因素解析發(fā)現(xiàn)：肉羊飼糧NDFD和ADFD的主要影響因素相同，分別為飼糧CP、DM含量及MBW、RM，其中飼糧CP含量與NDFD、ADFD顯著正相關(guān)，Broderick[24]研究發(fā)現(xiàn)，初產(chǎn)奶牛NDFD和ADFD隨飼糧CP含量升高而線性升高，與本研究結(jié)果一致，這可能是瘤胃中大多數(shù)CP降解菌都具有一定程度降解飼糧粗纖維的能力，飼糧CP含量升高可增強(qiáng)肉羊瘤胃發(fā)酵能力，進(jìn)而增強(qiáng)對纖維物質(zhì)的降解能力[25]，使得NDFD和ADFD升高。飼糧DM含量與NDFD、ADFD顯著負(fù)相關(guān)，這可能是飼糧DM含量降低導(dǎo)致肉羊DMI減少，瘤胃微生物轉(zhuǎn)而充分降解纖維物質(zhì)為其供能[26]，使得肉羊?qū)︼暭ZNDF和ADF的消化能力增強(qiáng)。本研究同時表明肉羊MBW與NDFD和ADFD顯著正相關(guān)，這可能是隨肉羊年齡的增長MBW增加，同時伴隨著瘤胃的快速發(fā)育，逐漸具備可發(fā)酵粗纖維的瘤胃內(nèi)環(huán)境和瘤胃菌群，提高了對飼糧中纖維物質(zhì)的降解能力，因而提高了NDFD和ADFD。前人研究結(jié)果也證實該推論：Nie等[27]研究發(fā)現(xiàn)牦牛隨年齡增長，食源從易消化的母乳轉(zhuǎn)變?yōu)殡y消化的草料，瘤胃纖維降解菌數(shù)量增長，并形成穩(wěn)定的微生物區(qū)系，對飼糧中粗纖維的降解能力增強(qiáng)。

3.3 本研究的局限性

肉羊品種、泌乳期、飼糧組成、營養(yǎng)水平、環(huán)境等因素均會影響飼糧營養(yǎng)物質(zhì)表觀消化率[14]，本研究通過數(shù)據(jù)庫構(gòu)建的GED預(yù)測模型可解釋度強(qiáng)，但DMD、CPD、NDFD和ADFD預(yù)測模型仍分別有35.0%、67.0%、49.0%和53.0%的變異未查明來源。研究表明環(huán)境[4]和飼養(yǎng)管理模式[28]也會影響營養(yǎng)物質(zhì)表觀消化率，這可能是本研究預(yù)測模型誤差的來源，但關(guān)于環(huán)境和飼養(yǎng)管理模式影響飼糧營養(yǎng)物質(zhì)表觀消化率的研究相對較少，難以滿足構(gòu)建數(shù)據(jù)庫的要求，隨今后納入數(shù)據(jù)庫該方面文獻(xiàn)數(shù)據(jù)和影響因素的增多，預(yù)測模型的準(zhǔn)確性將逐步提高。本研究通過數(shù)據(jù)庫構(gòu)建的預(yù)測模型大量納入了我國的肉羊品種，因此更適合預(yù)測我國肉羊的飼糧營養(yǎng)物質(zhì)表觀消化率，模型普適性強(qiáng)。

4 結(jié) 論

① 本研究通過數(shù)據(jù)庫構(gòu)建的肉羊飼糧DMD、GED、CPD、NDFD和ADFD預(yù)測模型可分別解釋65.0%、95.0%、67.0%、51.0%和47.0%的變異，MB分別為0.332、6.346、-0.469、-0.073和0.098個百分點，RMSPE分別為7.1%、23.1%、8.3%、17.5%和20.2%。

② 本研究通過數(shù)據(jù)庫構(gòu)建的預(yù)測模型明確了肉羊飼糧DMD最主要的影響因素是MBW和DMI；GED最主要的影響因素是RM和飼糧DM含量；CPD最主要的影響因素是DMI、ADG和MBW；NDFD和ADFD最主要的影響因素是飼糧DM、CP含量及MBW、RM。