高程程，劉玲玲，何美升，王 軍，阿克瑪熱·阿樂，趙 茜，劉武軍

（新疆農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，烏魯木齊 830052）

新疆伊犁馬泌乳曲線影響因素的研究

高程程，劉玲玲，何美升，王 軍，阿克瑪熱·阿樂，趙 茜，劉武軍

（新疆農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，烏魯木齊 830052）

【目的】通過擬合伊犁母馬的最佳泌乳曲線來預(yù)測(cè)伊犁母馬產(chǎn)奶能力?！痉椒ā窟x取6～7歲、處于泌乳中期的伊犁母馬60匹為研究對(duì)象，利用4種模型（高斯模型、Nelder逆多項(xiàng)式模型、立方模型和六次多項(xiàng)式模型）和3種日間隔方式（間隔0、10、30d），對(duì)第3、4胎次的伊犁母馬的泌乳曲線進(jìn)行擬合，分析不同模型、日間隔和胎次3個(gè)因素對(duì)泌乳曲線擬合度及其參數(shù)值的影響。【結(jié)果】4種模型泌乳曲線擬合度變化范圍為0.077 7～0.827 0，不同模型極顯著影響參數(shù) b值、c值、擬合度 R2（P＜0.01），顯著影響參數(shù) a值（P＜0.05），日間隔、胎次水平對(duì)參數(shù) a值、b值、c值和擬合度均沒有顯著性影響（P＞0.05）。六次多項(xiàng)式模型擬合度最高，逆多項(xiàng)式回歸模型擬合度最低，立方模型擬合曲線更接近實(shí)際產(chǎn)奶?！窘Y(jié)論】該馬場(chǎng)中第3胎次伊犁馬間隔0d的產(chǎn)奶量的擬合效果較好，其中立方模型擬合曲線的泌乳潛力較高，達(dá)到泌乳高峰速度較慢，從泌乳高峰下降速度較慢，科學(xué)評(píng)定伊犁母馬的泌乳能力，確定立方模型為第3、4胎次伊犁馬最佳泌乳曲線模型，為伊犁馬的產(chǎn)奶性能的科學(xué)選育提供前期基礎(chǔ)。

伊犁馬；日產(chǎn)奶量；泌乳曲線；模型參數(shù)

0 引 言

【研究意義】伊犁馬，多產(chǎn)自新疆伊犁地區(qū)，屬于地方培育品種，具有良好的肉、乳兼用型體型［1］。馬奶作為馬的附屬產(chǎn)品，其包含多種微量元素和特殊功能因子，風(fēng)味獨(dú)特，主要的乳成分更接近人奶［2］。由于馬奶具有高營養(yǎng)價(jià)值特性，使得市場(chǎng)對(duì)馬奶的需求量逐漸增加，生產(chǎn)也隨之發(fā)展，發(fā)展乳用馬產(chǎn)業(yè)，具有很大的潛力。目前乳用馬品種少、母馬泌乳期短、產(chǎn)奶少，培育高產(chǎn)型乳用馬品種成為乳用馬產(chǎn)業(yè)化的重要出路。泌乳曲線模型是用來描述泌乳期內(nèi)產(chǎn)奶量隨時(shí)間變化的一個(gè)數(shù)學(xué)模型［3］，研究伊犁母馬的泌乳曲線影響因素，在伊犁母馬的泌乳規(guī)律性變化、產(chǎn)奶量預(yù)測(cè)、高產(chǎn)乳用馬的品系培育等方面具有重要意義?！厩叭搜芯窟M(jìn)展】馬的產(chǎn)奶量不僅受到遺傳因素的影響，更重要的是受到品種、年齡、胎次、體重、飼養(yǎng)水平、營養(yǎng)物質(zhì)、生理階段等其他因素的影響［4］。夸特馬從分娩后10～150d的產(chǎn)奶量從14.4下降到9.8 kg／d［5］。盧西塔諾馬產(chǎn)后第 31 d，達(dá)到泌乳高峰，為14.0 kg／d，之后開始逐漸下降，到180d時(shí)，產(chǎn)奶量下降到 7.57 kg／d［6］?！颈狙芯壳腥朦c(diǎn)】目前，國內(nèi)關(guān)于馬的泌乳規(guī)律變化的相關(guān)研究報(bào)道較少。以新疆伊犁州昭蘇馬場(chǎng)伊犁馬不同胎次的日產(chǎn)奶量為基礎(chǔ)，利用高斯模型、Nelder逆多項(xiàng)式模型（IQP模型）、立方模型（Cubic模型）及6次多項(xiàng)式模型等4種泌乳曲線模型對(duì)日產(chǎn)奶量進(jìn)行擬合曲線擬合，通過影響泌乳曲線的各參數(shù)及擬合度的方差分析，選擇確定最佳擬合泌乳曲線。【擬解決的關(guān)鍵問題】研究新疆伊犁州昭蘇馬場(chǎng)伊犁馬的日產(chǎn)奶量變化及泌乳曲線，為提高伊犁馬的泌乳性能提供基礎(chǔ)資料，同時(shí)為伊犁馬的品系培育和飼養(yǎng)水平等方面提供參考依據(jù)。

1 材料與方法

1.1 材 料

試驗(yàn)于2014年6～8月在新疆伊犁哈薩克自治州昭蘇縣馬場(chǎng)進(jìn)行，選擇第 3、4胎次的體格健康、精神狀況良好的伊犁母馬各30匹，以這60匹伊犁母馬的日產(chǎn)奶量為研究對(duì)象。

1.2 方 法

1.2.1 采集產(chǎn)奶量

根據(jù)馬的生理特點(diǎn)，每日分別在 11：00、13：00、15：00、17：00對(duì)伊犁母馬進(jìn)行人工擠奶，使用電子天平稱量并記錄每匹伊犁母馬在每個(gè)時(shí)間點(diǎn)的產(chǎn)奶量。

1.2.2 處理

為了能夠使制定出的校正系數(shù)更加精確而且能夠反應(yīng)實(shí)際，對(duì)收集的數(shù)據(jù)進(jìn)行如下處理：去除泌乳月數(shù)大于 5個(gè)月的記錄，去除泌乳月數(shù)不足2個(gè)月的記錄；去除日產(chǎn)奶量不足 1 kg的記錄。日產(chǎn)奶量的數(shù)據(jù)處理采用薩伊金公式。

其中W是一晝夜產(chǎn)奶量，G是12 h的全部產(chǎn)奶量，h是擠奶時(shí)間總和。

1.2.3 擬合曲線所用的日間隔

試驗(yàn)從伊犁母馬第二個(gè)泌乳月開始采集產(chǎn)奶量，根據(jù)不同日間隔分為5類：連續(xù)記錄即間隔 0d，間隔10d是每隔9 d測(cè)定次日的產(chǎn)奶量，間隔30d是每隔29 d測(cè)定次日的產(chǎn)奶量。

1.2.4 泌乳曲線的統(tǒng)計(jì)模型

泌乳曲線的擬合采用在奶牛上常使用模型中的四種，分別為高斯模型、Nelder逆多項(xiàng)式模型（IQP模型）、立方模型（Cubic模型）及六次多項(xiàng)式模型，列出模型的表達(dá)式。表1

表1 四種模型表達(dá)式Table 1 The expressions of four models

4個(gè)模型表達(dá)式中 X均表示泌乳天數(shù)，Y均為第 X d的產(chǎn)乳量；a、b、c、d、e、f、g均表示模型的參數(shù)，a是泌乳潛力，b是從泌乳高峰下降的速度，c是達(dá)到泌乳高峰的速度，exp是自然對(duì)數(shù)底數(shù)。分別利用以上4種模型對(duì)第3、4胎次伊犁馬進(jìn)行擬合泌乳曲線，使用MATLAB 2014軟件進(jìn)行分析，計(jì)算過程采用非線性（NLIN）回歸過程 DUD迭代算法。

1.2.5 模型初始值

模型Ⅰ參數(shù)的初始值由線性化的非線性模型獲得；模型Ⅱ的初始值是由MATLAB 2014軟件中的 NLIN過程自動(dòng)按所給范圍和步長搜索出誤差平方和SSE最小的參數(shù)初值組合；模型 Ⅲ、模型Ⅳ的初始值是通過選取 p觀察值代入模型，構(gòu)建一個(gè)包含 p個(gè)未知參數(shù)的m個(gè)方程組，從而求出方程解來當(dāng)作參數(shù)的迭代初始值。

1.2.6 擬合曲線過程及擬合效果評(píng)價(jià)

在擬合曲線時(shí)，由于胎次對(duì)產(chǎn)奶量可能會(huì)產(chǎn)生影響，研究單獨(dú)對(duì)第 3、4胎次母馬產(chǎn)奶數(shù)據(jù)進(jìn)行擬合。所有產(chǎn)奶量數(shù)據(jù)先用Excel進(jìn)行初步整理，最后用SPSS（17.0）非線性回歸子程序（Nonlinear Regression）進(jìn)行擬合。模型擬合效果用擬合度 R2進(jìn)行評(píng)價(jià)。

1.3 數(shù)據(jù)統(tǒng)計(jì)

采用Excel 2010軟件及 SPSS 17.0軟件進(jìn)行相關(guān)數(shù)據(jù)的處理，采用 MATLAB 2014b軟件進(jìn)行泌乳曲線的繪制及分析。

2 結(jié)果與分析

2.1 不同胎次泌乳期日產(chǎn)奶量統(tǒng)計(jì)

利用 Excel 2010和 SPSS 17.0軟件將伊犁母馬的7 200個(gè)日產(chǎn)奶量數(shù)據(jù)進(jìn)行記錄整理，計(jì)算得出不同胎次伊犁母馬日產(chǎn)奶量的平均值及標(biāo)準(zhǔn)差，結(jié)果是第3胎次伊犁母馬日產(chǎn)奶量的平均值是10.56 kg，最高是20.10 kg，最低為2.82 kg，第4胎次的平均值是 7.68 kg，最高是 15.60 kg，最低是1.29 kg。表2

2.2 最佳泌乳曲線模型的選擇

分別利用高斯函數(shù)模型、Nelder逆多項(xiàng)式模型、立方模型、六次多項(xiàng)式回歸模型對(duì)新疆伊犁馬第3和第 4胎次分別進(jìn)行了泌乳曲線的擬合。新疆伊犁母馬日產(chǎn)奶量的不同模型、日間隔及胎次的擬合結(jié)果表明，在間隔0d、間隔10d的測(cè)定條件下，第3、4胎次的模型比較均為模型Ⅳ＞模型Ⅲ＞模型Ⅰ＞模型Ⅱ，模型Ⅳ的擬合效果優(yōu)于其他三種模型，在間隔 30d的測(cè)定條件下，第3、4胎次的模型比較均為模型Ⅲ＞模型Ⅰ＞模型Ⅱ，模型Ⅲ的擬合效果較好，由表模型Ⅲ、模型Ⅳ的R2值達(dá)到0.8以上，剩余平方和SSE越大，則剩余標(biāo)準(zhǔn)差RMSE越?。荒Ｐ廷?、模型Ⅲ、模型Ⅳ的間隔30d的擬合度高于其他日間隔，間隔30d＞間隔10d＞間隔0d，模型Ⅰ的間隔0d的擬合度相對(duì)較高。表3

表2 泌乳期日產(chǎn)奶量統(tǒng)計(jì)（kg）Table 2 lactation milk yield statistics

表3 新疆伊犁馬不同模型的擬合結(jié)果Table 3 Comparison of different models fitting to the data of Xinjiang Yili horse

2.3 泌乳曲線擬合參數(shù)及擬合度方差分析的 F值及顯著性檢驗(yàn)

2.3.1 統(tǒng)計(jì)不同泌乳曲線的擬合參數(shù)和擬合度

對(duì)不同模型、日間隔及胎次的擬合參數(shù)和擬合度進(jìn)行統(tǒng)計(jì)，研究表明，模型Ⅲ、模型Ⅳ的擬合度較高，其中第 3胎次模型Ⅲ的間隔0、10、30d的擬合度分別是0.832 4、0.891 8、0.999 3，第 4胎次模型Ⅲ的間隔 0、10、30d的擬合度分別是0.803 9、0.861 1、0.913 6，第 3胎次模型Ⅳ的間隔0、10d的擬合度分別是0.875 1、0..896 3，第4胎次模型Ⅳ的間隔 0、10d的擬合度分別是0.839 3、0.924 6。表4

表4 不同模型、日間隔、胎次的泌乳曲線擬合參數(shù)及擬合度Table 4 The fitness and parameters for lactation curve in different models，interval days and parity

2.3.2 擬合曲線參數(shù)、擬合度方差分析的 F值及顯著性檢驗(yàn)

經(jīng)方差分析及 F值顯著性檢驗(yàn)，不同模型對(duì)模型參數(shù)a值達(dá)到極顯著性水平（P＜0.01），對(duì)模型參數(shù)b值、c值、泌乳曲線擬合度 R2達(dá)到顯著性水平（P＜0.05）。而日間隔、胎次水平對(duì)模型各參數(shù)值與泌乳曲線擬合度均沒有顯著性影響（P＞0.05）。表5

表5 不同模型、日間隔和胎次對(duì)泌乳曲線參數(shù)及擬合度方差分析的 F值及顯著性檢驗(yàn)Table 5 The F value and significance test of ANOVA of parameters and R2of lactation curve for different models，interval days and parity

2.4 最佳泌乳曲線繪制

2.4.1 第3胎次不同日間隔的最佳泌乳曲線

第3胎次泌乳伊犁母馬的120d的立方模型擬合曲線圖顯示，第60d左右，日產(chǎn)奶量達(dá)到最高峰，持續(xù)時(shí)間較長，第80d左右開始下降，相比其他模型擬合度較高，更接近伊犁母馬實(shí)際泌乳情況。圖1

圖1 最佳泌乳曲線擬合效果Fig．1 The diagram of optimal lactation curve fitting effect

第3胎次泌乳伊犁母馬的120d的立方模型擬合曲線圖顯示，該泌乳曲線的趨勢(shì)變化為緩慢上升到達(dá)頂峰再快速降低，在第7個(gè)間隔10d即第70d左右，產(chǎn)奶量達(dá)到最高峰，第90d左右開始下降。圖2

圖2 最佳泌乳曲線擬合效果Fig．2 The diagram of optimal lactation curve fitting effect

第3胎次泌乳伊犁母馬的120d的立方模型擬合曲線圖顯示，在第3個(gè)間隔30d與第4個(gè)間隔30d之間，產(chǎn)奶量處于最高峰，該階段屬于泌乳中期。圖3

圖3 最佳泌乳曲線擬合效果Fig．3 The diagram of optimal lactation curve fitting effect

2.4.2 第4胎次不同日間隔的最佳泌乳曲線

第4胎次泌乳伊犁母馬的120d的立方模型擬合曲線圖表明，在 30d左右時(shí)，開始進(jìn)入泌乳高峰期，在40d左右產(chǎn)奶量出現(xiàn)下降趨勢(shì)。圖4

圖4 最佳泌乳曲線擬合效果Fig．4 The diagram of optimal lactation curve fitting effect

第4胎次泌乳伊犁母馬的120d的6次多項(xiàng)式回歸模型擬合曲線圖表明，在第 6個(gè)間隔 10d即第60d左右，日產(chǎn)奶量達(dá)到最高峰，在第70d左右產(chǎn)奶量開始下降。圖5

圖5 最佳泌乳曲線擬合效果Fig．5 The diagram of optimal lactation curve fitting effect

第4胎次泌乳伊犁母馬的120d的立方模型擬合曲線圖顯示，在第 3個(gè)間隔 30日，產(chǎn)奶量處于最高峰，該階段屬于泌乳中期，隨后產(chǎn)奶量出現(xiàn)下降趨勢(shì)。圖6

3 討 論

3.1 模型對(duì)泌乳曲線的影響

目前泌乳曲線模型的擬合研究集中在荷斯坦牛上，主要是把非線性模型轉(zhuǎn)換成為線性以后，利用 SAS軟件中的非線性（NLIN）正割法（DUD）迭代求解。在 Yadav等［7］對(duì)奶牛的泌乳曲線研究中，發(fā)現(xiàn)逆多項(xiàng)式模型的擬合度比 Wood模型的擬合度要高，劉玲玲［8］利用 Wood、逆多項(xiàng)式及六次多項(xiàng)式回歸等3個(gè)模型擬合伊犁馬第4、5胎次泌乳曲線，結(jié)果表明，逆多項(xiàng)式模型為新疆伊犁馬最優(yōu)泌乳曲線模型。

研究所采用的4種泌乳曲線模型是目前國內(nèi)外常用的模型，不同泌乳曲線模型對(duì)擬合度影響極顯著（P＜0.01），不同模型擬合度變化范圍為0.077 7～0.827 0，六次多項(xiàng)式回歸模型擬合度最高，逆多項(xiàng)式回歸模型擬合度最低。不同模型對(duì)模型參數(shù) b值、c值、泌乳曲線擬合度 R2達(dá)到顯著性水平（P＜0.05），對(duì)模型參數(shù) a值達(dá)到極顯著性水平（P＜0.01）。與上述對(duì)不同品種牛泌乳曲線擬合模型的結(jié)果并不一致，說明模型對(duì)不同產(chǎn)奶家畜的泌乳曲線影響較大。由于馬屬動(dòng)物與奶牛在動(dòng)物學(xué)分類中的不同，以及馬的乳房結(jié)構(gòu)及生理特征與奶牛不同，形成不同的泌乳期，造成對(duì)選擇最優(yōu)泌乳曲線模型的結(jié)果不同。

3.2 胎次對(duì)泌乳曲線的影響

泌乳曲線的特點(diǎn)主要反映在其各參數(shù)值的差異，在泌乳曲線各參數(shù)中，a是泌乳潛力，b是從泌乳高峰下降的速度，c是達(dá)到泌乳高峰的速度。

Tekerli等［9］通過研究影響土耳其荷斯坦牛泌乳曲線參數(shù)的相關(guān)因素，結(jié)果表明，牛場(chǎng)、年份、產(chǎn)犢季節(jié)對(duì)參數(shù) c值有極顯著影響，胎次對(duì)產(chǎn)奶高峰日有顯著影響，牛場(chǎng)、年份、產(chǎn)犢季節(jié)、胎次對(duì)高峰產(chǎn)奶量均有極顯著影響。付雪峰等［10］用Wood、IQP、多項(xiàng)式回歸模型、AS和WIL 5種模型對(duì)2 843頭新疆褐牛泌乳曲線進(jìn)行擬合，擬合度變化趨勢(shì)基本上是：第二胎＞第三胎及三胎以上＞第一胎。研究通過3、4胎次的泌乳曲線方差分析得知，胎次對(duì)伊犁母馬泌乳曲線的參數(shù)值的影響不顯著，與前人結(jié)果并不一致。主要原因是馬與牛屬于不同動(dòng)物，牛的胎次水平對(duì)牛個(gè)體產(chǎn)奶性能影響較大，頭胎馬與經(jīng)產(chǎn)馬的泌乳能力與馬自身的生理因素有關(guān)，由于第 3和第 4胎次的馬身體發(fā)育穩(wěn)定，此年齡段的馬在不受外界環(huán)境變化的影響下，可仍舊保持相同水平的泌乳力。

3.3 日間隔對(duì)泌乳曲線的影響

采取日間隔是實(shí)際工作當(dāng)中節(jié)約時(shí)間的一種方法。王雅春等［11］利用 7種模型對(duì)新疆呼圖壁種牛場(chǎng)3個(gè)品種奶牛泌乳曲線擬合結(jié)果表明，當(dāng)間隔30d時(shí)，泌乳曲線擬合度比間隔10d有大幅度提高。研究中三種日間隔對(duì)模型各參數(shù)值與泌乳曲線擬合度均沒有顯著性影響（P＞0.05），與奶牛上的研究并不一致。主要原因是母馬的乳房結(jié)構(gòu)及生理特征與奶牛不同，造成馬與牛具有不同的泌乳期，馬的泌乳期是120到 180d，而牛的泌乳期是305 d，由于馬的泌乳期相對(duì)較短，采取不同日間隔，日間隔的跨度不影響整體泌乳變化。

泌乳曲線影響因素主要是樣本的數(shù)量和真實(shí)性，因此要較準(zhǔn)確的繪制伊犁馬的泌乳曲線，并預(yù)測(cè)其產(chǎn)奶水平，必須利用實(shí)際生產(chǎn)數(shù)據(jù)進(jìn)行擬合，并對(duì)影響泌乳曲線參數(shù)的相關(guān)因素進(jìn)行詳細(xì)分析，才能達(dá)到較好的效果。

3.4 泌乳曲線變化趨勢(shì)

通過立方模型、6次多項(xiàng)式回歸模型對(duì)第3、4胎次伊犁馬的不同間隔泌乳曲線的擬合趨勢(shì)能夠看出，兩個(gè)胎次泌乳曲線的波形趨勢(shì)是基本一致，都是先上升到最高點(diǎn)再下降至干奶期的過程。但是同第4胎比較，第3胎日產(chǎn)奶量上升較慢，最高日產(chǎn)奶量較高，高峰日持續(xù)時(shí)間長，泌乳后期產(chǎn)奶量下降較快；反映了第3胎伊犁馬泌乳持續(xù)力較強(qiáng)，而第4胎牛泌乳持續(xù)力較弱的特點(diǎn)。在鄧緋［12］的研究中，選取泌乳期、年齡和胎次相近、產(chǎn)駒時(shí)間接近的伊犁馬8匹，利用每10d記錄的日產(chǎn)奶量建立泌乳曲線，泌乳曲線為先上升后下降趨勢(shì)，下降趨勢(shì)較為緩慢。在姚新奎［3］研究的不同品種母馬泌乳曲線圖中，以泌乳天數(shù)為橫坐標(biāo)，日產(chǎn)奶量為縱坐標(biāo)繪制的泌乳曲線，同樣是先上升后下降的趨勢(shì)，在第80d時(shí)略有上升。論文利用立方模型對(duì)泌乳曲線進(jìn)行擬合，計(jì)算出日產(chǎn)奶量，繪制出泌乳曲線，結(jié)果與上述研究結(jié)果一致。

4 結(jié) 論

研究利用高斯模型、Nelder逆多項(xiàng)式模型（IQP模型）、立方模型（Cubic模型）及六次多項(xiàng)式模型對(duì)伊犁州昭蘇馬場(chǎng)的第 3、4胎次伊犁母馬的泌乳曲線進(jìn)行擬合，根據(jù) R2越大越好（更接近1）的原則，發(fā)現(xiàn)在伊犁州昭蘇馬場(chǎng)的第 3、4胎次伊犁母馬中，間隔0d立方模型的擬合效果最好，其泌乳潛力高，達(dá)到泌乳高峰速度慢，從泌乳高峰下降速度慢，更切合實(shí)際地評(píng)定伊犁母馬的泌乳能力；通過分析不同模型、胎次和日間隔對(duì)泌乳曲線的差異性，以及繪制擬合曲線與實(shí)際產(chǎn)奶的對(duì)比結(jié)果，發(fā)現(xiàn)影響伊犁母馬泌乳曲線的主要因素是模型的選擇。確定立方模型為新疆伊犁州昭蘇馬場(chǎng)伊犁馬的最佳泌乳曲線模型。

（

）

［1］葉再華.伊犁馬種質(zhì)特性研究［D］.烏魯木齊：新疆農(nóng)業(yè)大學(xué)碩士學(xué)位論文，2008. YE Zai－h(huán)ua.（2008）.Research of Germplasm characteristics of Yili horse［D］.Master Dissertation.Xinjiang Agricultural University，Urumqi.（in Chinese）

［2］Malacarne，M.，Martuzzi，F(xiàn).，Summer，A.，＆Mariani，P.（2002）.Protein and fat composition of mare＇s milk：some nutritional remarks with reference to human and cow＇s milk.International Dairy Journal，12（11）：869－877..

［3］王瑞軍，張志剛，吳宏軍，等.影響內(nèi)蒙古三河牛泌乳曲線參數(shù)的非遺傳因素［J］.中國畜牧獸醫(yī)，2010，37（6）：14－16. WANG Rui－jun，ZHANG Zhi－gang，WU Hong－jun，et al.（2010）.Study of Non－genetic Factors Affecting on Lactation Curve Traits of Inner Mongolia Sanhe Cattle［J］.China Animal Husbandry and Veterinary Medicine，37（6）：14－16.（in Chinese）

［4］姚新奎.伊犁馬、新吉馬及其雜交馬乳理化指標(biāo)、泌乳特性初步研究 ［D］.烏魯木齊：新疆農(nóng)業(yè)大學(xué)博士學(xué)位論文，2011. YAO Xin－kui.（2011）.Preliminary Study on Chemical－Physical Indicators of Mare｀s Milk and lactation Character of Yili Horse，New Kirgiz Horse and Hybrid Horse［D］.PhD Dissertation.Xinjiang Agricultural University，Urumqi.（in Chinese）

［5］Gibbs，P.G.，Potter，G.D.，Blake，R.W.，＆Mcmullan，W.C.（1982）.Milk production of quarter horse mares during 150days of lactation.Journal of Animal Science，54（3）：496－499.

［6］Santos，A.S.，＆Silvestre，A.M.（2008）.A study of lusitano mare lactation curve with wood＇s model.Journal of Dairy Science，91（2）：760－766.

［7］Yadav S B S，Sharma J S.（1985）.Functions for lactation curves in crossbred dairy cattle.Ind J Animal Science，55（1）：42.

［8］劉玲玲.新疆伊犁馬早期泌乳性狀及體型外貌與泌乳性能關(guān)系的研究［D］.烏魯木齊：新疆農(nóng)業(yè)大學(xué)碩士學(xué)位論文，2014. LIU Ling－ling.（2014）.Studies on Relationship of Early Lactation Traits and Body Image with Lactation Performance of YiLi Horse in Xinjiang［D］.Master Dissertation.Xinjiang Agricultural University，Urumqi.（in Chinese）

［9］Tekerli，M.，Akinci，Z.，Dogan，I.，＆Akcan，A.（2000）. Factors affecting the shape of lactation curves of holstein cows from the balikesir province of turkey.Journal of Dairy Science，83（6）：1，381－1，386.

［10］付雪峰，王亞春，郭俊青，等.新疆褐牛產(chǎn)奶量校正系數(shù)的制定［J］.畜牧獸醫(yī)學(xué)報(bào)，2010，41（5）：536－542. FU Xue－feng，WANG Ya－chun，GUO Jun－qing，et al.（2010）.Study of Adjustment Factors for Standardizing Milking Records in Xinjiang Brown Cattle［J］.Chinese Journal of Animal and Veterinary Sciences，41（5）：536－542.（in Chinese）

［11］王雅春，陳幼春，柏榮，等.奶牛泌乳曲線的擬合及其模型參數(shù)的遺傳分析［J］.畜牧獸醫(yī)學(xué)報(bào)，1999，30（5）：399－404. WANG Ya－chun，CHEN You－chun，BAI Rong，et al.（1999）.Genetic analysis of model＇s parameters and dairy lactation curve in cattle［J］.Chinese Journal of Animal and Veterinary Sciences，30（5）：399－404.（in Chinese）

［12］鄧緋.伊犁馬產(chǎn)奶量、乳成分測(cè)定與分析［J］.吉林農(nóng)業(yè)，2010，（7）：49－50. DENG Fei.（2010）.Analysis of milk yield and milk composition of Yili horse［J］.Jilin Agriculture，（7）：49－50.（in Chinese）.

Effects of Different Feed Models on Milk Production of Yili Horse

GAO Cheng－cheng，LIU Ling－ling，HE Mei－sheng，WANG Jun，Akemare·Ale，ZHAO Xi，LIU Wu－jun
（College of Animal Science，Xinjiang Agricultural University，Urumqi 830052，China）

【Objective】The research was designed to predict milking capacity by fitting the best lactation curves.【Method】60 mares of Yili horses were selected as samples from Zhaosu Stud Farm in Yili，Xinjiang. The lactation curve was fitted by 4 mathematical models（including the Gaussian model，Nelder inverse quadratic polynomial model，Cubic model and Six polynomial model）and 3 days interval（including 0day interval，10days interval，30days interval）on 3－4 parities＇Yili horses.The effects of models，day interval and parity on fitness and parameter of lactation curve were analyzed by ANOVA.【Result】The results showed that the fitness of 4 lactation curve models ranged from 0.077，7 to 0.827，0.The models showed significant effect on the fitness，the b value，the c value of lactation curve at P＜0.01，on the a value of lactation curve at P＜0.05.The day interval and parity showed no significant effect on the fitness，the a value，the b value，the c value of lactation curve.By the fitness of lactation curve，the sixth polynomial model was the largest，and the Nelder inverse quadratic polynomial model was the least，the cubic model was better near to the real milk production.【Conclusion】The study found the lactation curves of 0day interval on third parity Yili horse is betterr，and among them，the lactation potential of the cubic model fitting curve was higher，but the rising speed of the peak of lactation was lower，so was the dropping speed.The cubic model is more reliable to evaluate the milking capacity of Yili horse.So，the cubic model is the best model to curve lactation，which provided an initial basis for the later scientific breeding of Yili horses milk performance.

Yili horse；milk yield；lactation curve；model parameters

S821.4

A

1001－4330（2016）09－1743－08

10.6048／j.issn.1001－4330.2016.09.023

2016－04－23

國家科技支撐計(jì)劃項(xiàng)目課題（2012BAD44B01）；國家自然科學(xué)基金項(xiàng)目（31560620）

高程程（1990－），男，陜西人，碩士研究生，研究方向?yàn)閯?dòng)物遺傳育種與繁殖，（E－mail）515877292＠qq.com

（Cotresponding author）：劉武軍（1966－），女，新疆人，教授，博士，博士生導(dǎo)師，研究方向?yàn)閯?dòng)物遺傳育種與繁殖，（E－mail）lwj_ ws＠163.com

Fund project：National science and technology support project（2012BAD44B01）；The Project of National Natural Science Foundation of China（31560620）