基于薈萃分析研究肉禽非植酸磷需要量的最佳擬合模型

張嘉杰，信愛國，劉亞萍，王潤蓮

摘 要：本研究基于薈萃分析，采用3種常用模型對非植酸磷水平與肉雞、鴨、鵝脛骨灰分含量的劑量-效應試驗結果進行擬合比較，確定選擇最佳擬合模型，為肉禽非植酸磷需要量的評定提供依據。共收集32組肉禽非植酸磷水平與脛骨灰分含量的劑量-效應試驗數據，分別用二次曲線、線性折線及二次折線模型進行擬合，通過對不同模型的平均擬合指數及最佳擬合模型出現頻度來確定最優擬合方程。結果表明，三種模型對32組數據的擬合度范圍分別為：二次曲線模型65.7%～100.0%，二次折線模型4.6%～100.0%，折線模型4.6%～99.8%，擬合度均值分別為92.1%、89.5%、88.6%，最佳擬合出現頻度分別為17/32、18/32、15/32。綜上所述，以脛骨灰分含量為效應值來評定肉禽非植酸磷需要量時，建議首選二次曲線模型，次選二次折線模型。

關鍵詞：非植酸磷水平；肉禽；脛骨灰分；擬合模型

中圖分類號：S831.5文獻標識碼：B文章編號：1673-1085（2023）05-0001-07

磷是家禽生長必需的礦物元素之一，不僅是骨骼重要的結構物質，而且在細胞膜平衡、物質代謝、機體內環境的穩定中起重要作用。家禽缺乏磷會表現出精神低迷、食欲下降、站立不穩、跛行、行走困難、甚至癱瘓不起等臨床癥狀[1]，因此，滿足家禽的磷需要非常重要。由于植物性飼料中有70%以上的磷是以植酸磷的形式存在，很難被利用[2，3]，生產上家禽的磷需要量是指非植酸磷（Non-phytate phosphorus，NPP）的需要量[4]。許多動物營養工作者為確定家禽的非植酸磷需要量進行了大量研究，多采用劑量-效應試驗，通過不同的效應指標和擬合模型得出不同品種家禽在不同生長階段的非植酸磷需要量。如Yan等[5]以趾骨灰分含量作為效應指標，通過斷線模型得出1～21日齡肉仔雞非植酸磷適宜水平為0.39%；王晉晉[6]以日增重為效應指標，采用響應面回歸模型得出4～6周齡白羽肉雞非植酸磷適宜水平為0.56%；劉松柏[7]綜合多種指標，采用斷線模型得出1～21日齡肉仔雞非植酸磷適宜水平為0.39%；代述均[8]以脛骨灰分含量為效應指標，采用斷線模型得出14日齡和21日齡肉鴨的非植酸磷適宜水平分別為0.40%及0.47%。不同的方法適合于不同的數據分布規律[9]。本研究采用薈萃分析的方法，收集并整理肉雞、肉鴨和肉鵝三種禽類有關非植酸磷需要量研究的文獻，選用其中非植酸磷水平對脛骨灰分含量有顯著影響的數據，采用目前常見的3種模型進行擬合分析比較，選擇最佳擬合模型[10-12]，為今后評定禽類非植酸磷的需要量提供參考。

1 材料及方法

1.1 非植酸磷水平及效應指標的確定

本研究所選擇的試驗研究需符合三個條件：1）試驗動物為肉雞、肉鴨和肉鵝三種家禽之一；2）采用單因子試驗設計，以非植酸磷水平為試驗因子；3）非植酸磷水平梯度大于或等于3個。很多試驗表明，脛骨灰分含量可作為敏感指標來評價禽類磷需要量[13]。因此，本研究選取脛骨灰分含量作為效應指標。

1.2 試驗數據收集與整理

通過中國知網和Science Direct兩大學術平臺檢索近25年發表的關于非植酸磷需要量的期刊或學位論文，關鍵詞采用“禽類”“肉雞”“肉鴨”“肉鵝”“非植酸磷”及“非植酸磷需要量”。根據查閱結果，共收集到28篇相關文獻，選用了符合要求的32組非植酸磷水平與脛骨灰分含量存在顯著劑量-效應關系的數據，根據組、文獻、品種、性別、日齡、非植酸磷水平和脛骨灰分含量整理成數據表，見表1。

1.3 擬合曲線

選用國內外文獻經常使用的3種數學模型，包括二次曲線模型、線性折線模型、二次折線模型進行比較[9]。

1.4 數據統計分析

選用SPSS 17.0統計軟件對32組非植酸磷水平及其對應的脛骨灰分含量數據進行擬合分析，通過對不同擬合模型的擬合指數R2及最佳擬合方程出現的次數或幾率確定最佳擬合曲線。

2 結果與分析

2.1 不同模型的擬合度比較

表2表示3種模型對32組數據的擬合度。由表2可知，二次曲線模型擬合度范圍為65.7%～100.0%，折線模型擬合度范圍為4.6%～99.8%，二次折線模型擬合度范圍為4.6%～100.0%。

2.2 不同模型的擬合度平均值及最佳擬合次數或幾率

表3為不同模型的擬合度（R2）均值（%）及最佳擬合次數或幾率。由表3可知，以脛骨灰分含量為效應指標，3種不同模型的平均擬合度從高到低依次為：二次曲線模型（92.1%）＞二次折線模型（89.5%）＞折線模型（88.6%）；其最佳擬合模型出現的數據組，二次曲線、折線模型及二次折線模型分別為17組、15組、18組，最佳擬合模型出現頻度從高到低依次為：二次折線模型（18/32）>二次曲線模型（17/32）>折線模型（15/32）。

2.3 二次曲線模型擬合肉用雞、鴨、鵝非植酸磷需要量

根據上述擬合度最高的二次曲線模型對32組數據進行模擬計算，得出肉用雞、鴨、鵝非植酸磷需要量，見表4。其中肉雞21日齡內的研究較多，其擬合非植酸磷需要量范圍為0.39%～0.55%。肉鴨、鵝的數據較少，非植酸磷需要量如下：北京鴨3～6周齡非植酸磷需要量為0.36%，櫻桃谷肉鴨1～35日齡非植酸磷需要量為0.47%～0.59%；東北肉鵝1～28日齡非植酸磷需要量為0.52%，陽江鵝1～70日齡非植酸磷需要量為0.41%。

3 討論

許多動物營養工作者已通過試驗得出了不同品種禽類在不同生長階段的非植酸磷需要量。然而，現代家禽品種繁多，其生長性能、骨骼結構特點以及肉品質各異，所以家禽磷需要量的評定是一項不斷持續和更新的工作，其研究方法也需要不斷探討和改進。

劑量-效應試驗是評定動物營養需要量最常用的方法之一，通過試驗可以得到成對的劑量-效應數據，采用不同數學模型擬合分析，找出效應反應值隨著劑量水平變化的規律，得出擬合方程，從而確定營養需要量。目前擬合分析最常用的是二次曲線形式的拋物線模型（y = b0 + b1x + b2x2），通常適用于效應指標對營養素攝入不足或過量的反應是對稱的情況。折線模型[y = m （x < r）; y = m+ k（r-x） （x > r）]是營養需要量評定應用最廣泛的一種模型，要求在一定劑量水平范圍內，效應反應值隨營養素水平變化呈線性或曲線反應（即線性折線或二次折線反應），超過一定閾值后趨于穩定，因此可根據該閾值判斷最佳劑量水平。由于不同的模型適合于不同分布特點的數據，所以在試驗中模型的選擇對于能否得出準確的需要量評定結果至關重要。

近年來有一些研究采用了薈萃分析的方法，如李獻華[14]研究殼寡糖對肉雞生長性能的影響，通過整合相關文獻數據，對殼寡糖提高肉雞生長性能的效果進行定量分析，得出飼糧中添加適宜劑量殼寡糖可顯著提高平均日增重和降低料肉比，為生產中合理使用功能性寡糖提供參考依據。Wang等[15]匯總53個賴氨酸劑量-效應的肉雞試驗數據，通過9個模型的擬合分析及比較，得出肉雞賴氨酸營養需要量評定S線模型、二次曲線和酶動力模型效果較好，其擬合度分別為85.51%、87.45%、85.58%；最佳擬合出現的頻度分別為17/53、14/53、13/53。陳麗欣等[16]采用7種常用模型對肉雞、鴨飼糧鈣水平與脛骨灰分含量的劑量效應數據進行擬合比較，得出二次曲線、線性折線和二次折線模型的擬合度和最佳擬合出現的頻度高于其余模型，其擬合度分別為85.5%、71.5%、67.3%；最佳擬合出現的頻度分別為23/45、14/45、11/45，該三種模型的擬合效果優于其他模型。本研究收集28篇于1999～2022年期間發表的禽類非植酸磷需要量的研究論文，從中選用32組非植酸磷水平與脛骨灰分含量呈顯著性變化的劑量-效應數據，采用目前常用的3種擬合模型進行分析及對比。結果表明，3種模型的擬合度變異都比較大，說明不同的數據適合不同的模型，其平均擬合度依次為：二次曲線模型（92.1%）＞二次折線模型（89.5%）＞線性折線模型（88.6%）；各模型最佳擬合出現的頻度依次為：二次折線模型（18/32）＞二次曲線模型（17/32）＞線性折線模型（15/32）。從整體來看，二次曲線模型擬合效果最佳，二次折線模型次之，線性折線模型最差。因此，在家禽非植酸磷營養需要量評定的試驗中，建議采用二次曲線模型，其次是二次折線模型。

本試驗用擬合度最高的二次曲線擬合32組數據中肉雞、鴨、鵝非植酸磷需要量，分別得出肉用雞、鴨、鵝在不同飼養階段的非植酸磷需要量，可為相同或相近階段的肉禽飼糧配制和科學飼養提供重要的借鑒和參考。

4 結論

（1）以脛骨灰分含量為效應指標，研究選用3種擬合模型進行擬合分析比較，其擬合度范圍分別為：二次曲線模型65.7%～100.0%，折線模型4.6%～99.8%，二次折線模型4.6%～100.0%；擬合均值從高到低依次為：二次曲線模型（92.1%）＞二次折線模型（89.5%）＞線性折線模型（88.6%）；其最佳擬合出現頻度從高到低依次為：二次折線模型（18/32），二次曲線模型（17/32），線性折線模型（15/32）。

（2）根據脛骨灰分含量作為效應值來評定家禽非植酸磷需要量時，建議首選二次曲線模型，次選二次折線模型。

參考文獻：

[1] 羅玲萍，崔恒敏.雞鴨鈣磷缺乏癥的比較病理學[J]. 四川農業大學學報，1995，13（01）：96-100.

[2] 王成章，王恬. 飼料學[M].第二版.北京：中國農業出版社，2011：260.

[3] 陳雅湘，方熱軍. 磷的吸收代謝及影響因素研究進展[J].湖南飼料，2019（02）：30-32+48.

[4] 潘文.1～21日齡麒麟雞非植酸磷需要量及其植酸酶利用研究[D]. 湛江：廣東海洋大學，2016.

[5] Yan F， Keen C， Zhang K， et al. Comparison of methods to evaluate bone mineralization[J]. The Journal of Applied Poultry Research， 2005， 14（3）： 492-498.

[6] 王晉晉.0～6周齡白羽肉雞鈣磷需要量的研究[D].鄭州：河南工業大學，2010.

[7] 劉松柏.肉仔雞對飼料標準磷利用率、可利用磷需要量及小腸磷吸收機制研究[D].北京：中國農業科學院，2012.

[8] 代述均.飼糧非植酸磷水平對肉鴨生產性能、盲腸微生物及脛骨質量的影響[D].成都：四川農業大學，2019.

[9] Pesti G M， Vedenov D， Cason J A， et al. A comparison of methods to estimate nutritional requirements from experimental data[J]. British Poultry Science， 2009， 50（01）： 16-32.

[10] 何雪萍，信愛國，王潤蓮，等.肉鴨和蛋鴨蛋白質需要量評定擬合模型比較研究[J].中國家禽，2022，44（10）：50-55.

[11] 范海丹，信愛國，王潤蓮，等.肉鴨及生長期蛋鴨代謝能需要量評定擬合模型比較研究[J].廣東飼料，2022，31（08）：25-28.

[12] 陳麗欣，王潤蓮，譚鈞元，等.肉雞賴氨酸需要量最佳擬合模型薈萃分析研究[J].廣東飼料，2021，30（10）：47-50.

[13] 陳祖鴻，汪水平，彭祥偉，等.肉鴨鈣和磷營養研究進展[J].中國家禽，2014，36（08）：43-47.

[14] 李獻華.殼寡糖對肉仔雞生長性能影響的Meta分析[J].家禽科學，2022（03）：5-11.

[15] Wang R， Alhotan R A， Billard L， et al. A comparison of equations to depict and predict nutritional kinetics in fast-growing broiler chickens[J]. The Journal of Applied Poultry Research， 2019， 28（01）： 124-135.

[16] 陳麗欣，信愛國，王潤蓮，等.肉用雞、鴨鈣營養需要量評定擬合模型比較研究[J].家禽科學，2022（04）：5-11.

Studies on the Best Fitting Model for Non-phytate Phosphorus Requirement of Meat-Type Poultry Based on Meta-analysis

ZHANG Jiajie1， XIN Aiguo2， LIU Yaping1， WANG Runlian1

（1. Department of Animal Science， Guangdong Ocean University， Zhanjiang? 524088， China;

2. Yunnan Academy of Animal Science and Veterinary Sciences， Kunming? 650224， China）

Abstract：? The purpose of this study was to investigate the best fitting model used to fit the dose-response data between different non-phytate phosphorus （NPP） levels and tibia ash of poultry by 3 commonly used models based on meta-analysis， in an effort to find the optimum model for evaluating the scientific NPP level for chickens， ducks and goose in the future. A total of 32 sets of dose-response data with the significant effect of different NPP levels on tibia ash of chickens， ducks and goose were fitted with 3 models including quadratic regression model， broken line model， and broken quadratic model. The fitting performance （R2） and the incidence of models with the best fit were used to determine the best fitting model. The R2 ranges of 3 models were 65.7%～100.0%， 4.6%～100.0% and 4.6%～99.8%. The mean of R2 were 92.1%， 89.5% and 88.6%， respectively. The incidence of models with the best fit was 17/32 for quadratic regression model， 18/32 for broken quadratic model， and 15/32 for broken line model， respectively. In conclusion， to use the tibia ash concentration as a response value to determine NPP requirement for chickens， ducks and goose， it is recommended that the quadratic regression model is the first choice， followed by broken quadratic model.

Keywords：? Non-phytate phosphorus levels; Poultry; Tibia ash; Fitting model