本發(fā)明屬于犢牛粗飼料飼養(yǎng)，尤其涉及基于meta分析的犢牛粗飼料需要水平的評估方法。

背景技術(shù)：

1、近年來，關(guān)于犢牛是否應(yīng)該飼喂粗飼料的爭論受到了廣泛關(guān)注，盡管這方面的研究早在20世紀中期就已開始。美國的研究表明，在最近的20年間，斷奶前飼喂粗飼料的概率從86.7％下降到43.3％，這表明奶牛生產(chǎn)者的觀念發(fā)生了變化。支持者認為，飼喂粗飼料有利于犢牛的生長和瘤胃環(huán)境的改善。有些研究表明，在犢牛開食料之外添加粗飼料可提高采食量和飼喂效益、促進反芻、生長、增加瘤胃液ph值并使得瘤胃酸中毒得以緩解，對瘤胃壁的外觀，瘤胃的肌肉生長都有積極的改善。但在另一方面，反對飼喂粗飼料的人則認為，犢牛的瘤胃上皮發(fā)育不全，而且瘤胃容量有限，在牛乳飼喂階段無法有效利用粗飼料，應(yīng)飼喂精料來促進丁酸鹽和丙酸鹽的產(chǎn)生，這對斷奶前的瘤胃發(fā)育至關(guān)重要。然而，精料具有很高的適口性和易消化性，其中快速發(fā)酵碳水化合物會大量產(chǎn)生揮發(fā)性脂肪酸(vfas)，并可能通過影響反芻和唾液以及瘤胃中營養(yǎng)物質(zhì)的消化率而加重犢牛的瘤胃酸中毒。犢牛需要充足的能量和蛋白質(zhì)來維持生長，并在斷奶前達到推薦體重。發(fā)育成熟的腸道可提高飼料消化率，從而獲得更多的營養(yǎng)物質(zhì)來支持新陳代謝和犢牛生長。精料可隨時為瘤胃的生長和發(fā)育提供營養(yǎng)。然而，飼喂粗飼料會影響乳牛的dm、有機物(om)和總養(yǎng)分消化率。與不飼喂粗飼料的犢牛相比，飼喂粗飼料的犢牛dm、om和cp的總表觀消化率較低。此外，飼喂苜蓿干草(ah)的犢牛的cp消化率高于不飼喂粗飼料的開食性犢牛，但與飼喂小麥straw和甜菜漿的犢牛相似。然而，有些研究者指出，當(dāng)犢牛在飼喂精飼料的同時飼喂苜蓿干草或只飼喂精飼料時，在dm、om或cp消化率方面沒有觀察到明顯的差異。許多因素都與混雜結(jié)果有關(guān)，包括粗飼料來源、粗飼料水平、開食的物理形式、粗飼料飼喂方法和牛奶飼喂策略。最近，一些針對幼犢的研究發(fā)現(xiàn)，飼喂全混合日糧(tmr；定義為混合犢牛開食料和粗飼料)比自由選擇日糧(在不同桶中飼喂開食料和粗飼料)更理想。在改善瘤胃ph值和咀嚼活動方面，在犢牛開食料中添加來自ah的ndf比添加甜菜漿更有效，但不會同時影響開食料的攝入量或adg。盡管已經(jīng)使用元分析工具探討了生長中的犢牛飼喂粗飼料的問題，但關(guān)于犢牛日糧中是否應(yīng)包括粗飼料的問題仍然層出不窮(nrc，2001年)。盡管該領(lǐng)域的研究成果不斷增加，但支持或反對飼喂粗飼料的定量證據(jù)仍然很少。此外，理想的粗飼料水平尚不清楚，研究中使用了不同類型的粗飼料。試驗中使用的粗飼料包括：溴草干草、燕麥干草和紫花苜蓿。此外，不同粗飼料類型之間的比較在文獻中也很有限，盡管飼喂粗飼料對犢牛表現(xiàn)的影響可能存在差異。因此，本發(fā)明的目標是：(1)評估犢牛飼喂粗飼料對生長性能的影響；(2)確定粗飼料來源在粗飼料飼喂反應(yīng)異質(zhì)性中的作用；(3)根據(jù)報道粗飼料水平的文章中的數(shù)據(jù)，確定提高生長性能的最佳粗飼料飼喂水平。

技術(shù)實現(xiàn)思路

1、本發(fā)明的目的在于克服現(xiàn)有技術(shù)的缺點，提供了基于meta分析的犢牛粗飼料需要水平的評估方法，通過meta分析結(jié)合混合模型構(gòu)建犢牛生長性能和粗飼料水平的回歸曲線，尋找粗飼料最佳添加水平；首先采用meta分析并使用cochrane's?q檢驗估計效應(yīng)大小統(tǒng)計異質(zhì)性；隨后亞組分析以評估不同粗飼料來源和粗飼料供應(yīng)方法對結(jié)果異質(zhì)性的影響；最后采用混合模型預(yù)測粗飼料水平與生長性能之間的關(guān)系；利用此方法分析文獻后，結(jié)果表明，補充粗飼料可提高犢牛的生長性能；本發(fā)明是基于meta分析并探究粗飼料對犢牛生長性能的影響以及犢牛飼喂最佳粗飼料水平，從而精準的得到飼養(yǎng)犢牛時的飼養(yǎng)方法；可為犢牛精準飼養(yǎng)提供思路，以幫助犢牛飼養(yǎng)策略，操作簡便易行。

2、為實現(xiàn)以上技術(shù)效果，采用如下技術(shù)方案：

3、基于meta分析的犢牛粗飼料需要水平的評估方法，包括如下步驟：

4、步驟s1：文章檢索和納入

5、在文獻數(shù)據(jù)庫中搜索一定期限內(nèi)發(fā)表的研究論文；主要檢索詞包括但不限于：calves、forage、performance；

6、按照pico原則閱讀論文文章，按照納入標準，決定哪些研究可以納入；綜述、調(diào)查文獻以及未包含隨機對照試驗的文獻被排除；

7、所述納入標準為：所述研究論文的研究對象必須為犢牛；所述研究論文中試驗必須包括一個不飼喂任何粗飼料的對照組和一個或多個至少飼喂一種粗飼料的干預(yù)組；

8、步驟s2：數(shù)據(jù)提取

9、必要數(shù)據(jù)信息提?。好科芯空撐膶φ战M的犢牛平均日增重adg值、每篇研究論文干預(yù)組的犢牛平均日增重adg值、樣本數(shù)量n以及標準差sd，并統(tǒng)一平均日增重adg值的單位；構(gòu)建meta-analysis數(shù)據(jù)庫；

10、若某些研究論文中只是公開了干預(yù)組的標準誤差se，則需要使用標準誤差se換算成標準差sd；換算公式為：

11、步驟s3：數(shù)據(jù)統(tǒng)計

12、記錄步驟s2中的所有數(shù)據(jù)，然后導(dǎo)出到stata/mp?16.0中進行統(tǒng)計分析和繪圖；森林圖可簡單直觀了解粗飼料對犢牛某一指標的影響；做森林圖時，以均值差md和95％置信區(qū)間描述變量的效應(yīng)大?。挥捎诓煌芯空撐拇嬖谠囼灧椒?、動物數(shù)量、試驗設(shè)計差異，選擇隨機效應(yīng)模型；同時，i2統(tǒng)計量和卡方q檢驗用于描述結(jié)果異質(zhì)性；異質(zhì)性的大小使用i2值進行分類，小于50％被認為是低異質(zhì)性，大于50％被認為是高異質(zhì)性，表明不同研究論文結(jié)果差異較大；

13、如果粗飼料有15項或更多的試驗對比，且森林圖結(jié)果表示研究指標的總體效應(yīng)顯著且具有顯著異質(zhì)性，則進行粗飼料亞組分析，以探究該粗飼料亞組對指標的影響程度；

14、步驟s4：回歸分析

15、基于meta-analysis數(shù)據(jù)庫，把各個干預(yù)組的adg值作為因變量、粗飼料水平作為自變量、研究作為隨機效應(yīng)，生成數(shù)據(jù)庫；利用混合模型建立犢牛adg與粗飼料水平的回歸關(guān)系圖，以得到最佳粗飼料添加水平，使用r語言的lme4包里面的混合模型進行擬合；

16、所述粗飼料亞組為oh水平和straw水平；

17、步驟如下：1、選擇合適的函數(shù)回歸模型；2、使用混合模型消除研究間隨機效應(yīng)，得到模型預(yù)測值；3、利用預(yù)測值建立回歸模型并擬合曲線；

18、選取的模型包括線性混合模型和非線性混合模型：

19、線性混合模型函數(shù)表達式如下：yij＝b0+b1xij+si+bixij+eij；

20、非線性混合模型函數(shù)表達式如下：yij＝b0+b1xij+b2x2ij+si+bixij+eij；

21、其中，yij為因變量(第i項研究1，...m，第j個觀測值1，...ni)，b0為所有研究的總截距，b1為自變量的系數(shù)，b2為自變量的二次系數(shù)，xij為自變量的值，si為第i項研究的隨機效應(yīng)截距，bi為與第i項研究相關(guān)的隨機斜率，eij為殘差；

22、r2和估計殘差方差的平方根srrv用于衡量回歸模型擬合的優(yōu)良性，r2越大，srrv越小，表明模型擬合度越好。

23、進一步的，所述步驟s1中的數(shù)據(jù)庫包括：google?scholar、sciencedirect和x-mol數(shù)據(jù)庫。

24、進一步的，所述步驟s1中的一定期限為20年以上，所述納入文獻數(shù)據(jù)至少為15篇，干預(yù)組和對照組至少40組。

25、進一步的，所述步驟s1中的pico原則包括研究對象、干預(yù)措施、比較因素、干預(yù)措施的診療結(jié)果。

26、進一步的，所述步驟s1中的閱讀論文文章包括標題、摘要、研究方法和試驗結(jié)果。

27、進一步的，所述步驟s1中的納入標準還必須提供一些有關(guān)生長性能的基本數(shù)據(jù)，所述基本數(shù)據(jù)包括試驗動物的數(shù)量、平均值、標準誤差或標準差。

28、進一步的，所述步驟s3中使用microsoft?excel記錄步驟s2中的所有數(shù)據(jù)。

29、進一步的，所述步驟s4中：

30、粗飼料水平和adg關(guān)系式為：y＝0.663+4.416×10-3x-1.72×10-4x2，且達到最佳adg時，粗飼料水平為13％；

31、oh水平與adg的關(guān)系式為：y＝0.657+0.042×10-3x-2.826×10-3x2，達到adg峰值時所需要的水平為7％；

32、straw水平與adg的關(guān)系式為：y＝0.681+0.018x-9.419×10-4x2，達到adg峰值時所需要的水平為9％。

33、本發(fā)明的有益效果為：

34、本發(fā)明公開了基于meta分析的犢牛粗飼料需要水平的評估方法，通過meta分析結(jié)合混合模型構(gòu)建犢牛生長性能和粗飼料水平的回歸曲線，尋找粗飼料最佳添加水平；首先采用meta分析并使用cochrane's?q檢驗估計效應(yīng)大小統(tǒng)計異質(zhì)性；隨后亞組分析以評估不同粗飼料來源對結(jié)果異質(zhì)性的影響；最后采用混合模型預(yù)測粗飼料水平與生長性能之間的關(guān)系；利用此方法分析文獻后，結(jié)果表明，補充粗飼料可提高犢牛的生長性能；本發(fā)明是基于meta分析并探究粗飼料對犢牛生長性能的影響以及犢牛飼喂最佳粗飼料水平，從而精準的得到飼養(yǎng)犢牛時的飼養(yǎng)方法；可為犢牛精準飼養(yǎng)提供思路，以幫助犢牛飼養(yǎng)策略，操作簡便易行。