楊 帆，邱昌功

（陜西省畜牧技術(shù)推廣總站，陜西 西安710016）

奶牛產(chǎn)奶量和體細(xì)胞之間存在著負(fù)相關(guān)關(guān)系，即體細(xì)胞上升奶產(chǎn)量下降，體細(xì)胞下降奶產(chǎn)量上升。為了摸清體細(xì)胞和產(chǎn)奶量之間的定量關(guān)系，為標(biāo)準(zhǔn)化飼養(yǎng)尋找依據(jù)。我們利用陜西省DHI中心提供的日產(chǎn)量和體細(xì)胞，月產(chǎn)量和體細(xì)胞，年產(chǎn)量（即305d產(chǎn)奶量）和體細(xì)胞的相關(guān)數(shù)據(jù)研究和統(tǒng)計分析奶產(chǎn)量和體細(xì)胞之間的相關(guān)關(guān)系。

從三方面進(jìn)行分析。即測定日產(chǎn)奶量和體細(xì)胞的關(guān)系，月產(chǎn)奶量和體細(xì)胞的關(guān)系，305d產(chǎn)奶量和體細(xì)胞的關(guān)系。

1 測定日產(chǎn)奶量和體細(xì)胞的關(guān)系

1.1 隨機(jī)選取杜陽奶牛場2014年9月測定日產(chǎn)奶量和體細(xì)胞數(shù)進(jìn)行相關(guān)分析。

共測定了297頭牛的產(chǎn)奶量和相應(yīng)體細(xì)胞數(shù)（平均產(chǎn)奶量22.29kg，平均體細(xì)胞49.28萬／mL），計算相關(guān)系數(shù)，對相關(guān)系數(shù)進(jìn)行顯著性檢驗（t檢驗），經(jīng)計算，r＝－0.127，由于n＝297，所以自由度df＝n－2＝295，又因為相同自由度時，r0.05＝0.114，r0.01＝0.150，r0.05＜r＜r0.01，所以r呈顯著負(fù)相關(guān)，推翻零相關(guān)假設(shè)，即r＝－0.127＊。

計算相關(guān)系數(shù)顯著性。經(jīng)計算，得t＝2.181＊，因為t0.05＝1.96，t0.01＝2.58，t0.05＜t＜t0.01，所以差異顯著，和相關(guān)系數(shù)顯著性相吻合。

1.2 隨機(jī)選取綠源牛場2014年7月份測定日奶產(chǎn)量和體細(xì)胞的數(shù)據(jù)進(jìn)行計算。共測定了253頭牛的產(chǎn)奶量和相應(yīng)的體細(xì)胞數(shù)，平均產(chǎn)奶量21.44kg，平均體細(xì)胞68.32萬／mL。計算相關(guān)系數(shù)，并進(jìn)行t檢驗。

經(jīng)計算，r＝－0.25，因為df＝253－2＝251，r0.01＝0.164，0.25＞0.164，所以負(fù)相關(guān)極顯著，即r＝－0.25＊＊。

計算t值，t＝4.099，因為4.099＞t0.01＝2.576，所以差異極顯著，即t＝4.099＊＊，同相關(guān)系數(shù)顯著性相吻合。

1.3 選取常興牛場2013年3月測定日產(chǎn)奶量和體細(xì)胞數(shù)據(jù)進(jìn)行相關(guān)系數(shù)的計算。共測定了161頭牛，平均產(chǎn)奶量29.15 kg，平均體細(xì)胞12.42萬／mL。經(jīng)計算相關(guān)系數(shù)r＝0.026，因為r0.05＝0.134＞r＝0.026，所以正相關(guān)不顯著。計算t值，t＝0.327，因為t＝0.327＜t0.05＝1.960，所以差異不顯著，結(jié)果相符。

2 月產(chǎn)奶量和體細(xì)胞的關(guān)系

2.1 隨機(jī)選取農(nóng)得利牛場2013年的月均產(chǎn)奶量數(shù)據(jù)和月均體細(xì)胞數(shù)據(jù)進(jìn)行相關(guān)系數(shù)計算。全年共測定241頭牛，按月統(tǒng)計日均產(chǎn)奶量和平均體細(xì)胞數(shù)，共統(tǒng)計了12個月的完整數(shù)據(jù)。數(shù)據(jù)如表1。

計算相關(guān)系數(shù)，r＝－0.879＊＊，因為df＝n－2＝10，r0.01＝0.708，0.879＞0.708，所以負(fù)相關(guān)極顯著。計算t值，得t＝5.83＊＊，因為t0.01＝3.169，t＞t0.01，差異極顯著。結(jié)果和相關(guān)系數(shù)顯著性完全吻合。

2.2 隨機(jī)選取省奶牛中心2013年月均產(chǎn)奶量數(shù)據(jù)和月均體細(xì)胞數(shù)據(jù)進(jìn)行相關(guān)性分析。數(shù)據(jù)如表2。

表1 農(nóng)得利2013年月均日產(chǎn)奶量和體細(xì)胞數(shù)統(tǒng)計表

表2 奶牛中心2013年月均日產(chǎn)奶量和體細(xì)胞數(shù)統(tǒng)計表

經(jīng)計算，r＝－0.399，df＝10，r0.05＝0.576，因為0.399＜0.576，所以r呈負(fù)相關(guān)且相關(guān)不顯著。計算t值，t＝1.376，因為t＜t0.05＝2.228，所以差異不顯著，結(jié)果和r顯著性相吻合。

2.3 2012年全省參測牛場月均產(chǎn)奶量和體細(xì)胞數(shù)據(jù)間的相關(guān)分析。數(shù)據(jù)如表3。

經(jīng)計算，r＝－0.214，df＝10，r0.05＝0.576，因為r＝0.214＜r0.05＝0.576，所以負(fù)相關(guān)不顯著。將r值帶入公式②求t值，t＝0.675，因為t＜t0.05＝2.228，所以差異不顯著，結(jié)果和r顯著性相吻合。

3 305天產(chǎn)奶量（即通常說的年產(chǎn)量）和體細(xì)胞數(shù)的相關(guān)分析。

3.1 咸陽西安草灘第二奶牛場2014年305d產(chǎn)奶量和體細(xì)胞的相關(guān)分析。數(shù)據(jù)如表4。

表3 陜西參測牛場2012年月均產(chǎn)奶量和體細(xì)胞數(shù)統(tǒng)計表

表4 西安草灘奶牛場2014年305天產(chǎn)奶量和體細(xì)胞統(tǒng)計表

經(jīng)計算，r＝－0.365，因為df＝9時，r0.05＝0.602，r＜r0.05，所以負(fù)相關(guān)不顯著。計算t值，t＝1.175，因為t0.05＝2.262＞t所以差異不顯著，結(jié)果和相關(guān)系數(shù)的顯著性相吻合。

3.2 農(nóng)得利奶牛場2013年305d產(chǎn)奶量和體細(xì)胞的相關(guān)分析。數(shù)據(jù)見表5。

經(jīng)計算，r＝－0.999＊＊，因為df＝10時，r0.01＝0.708，r＞r0.01，所以相關(guān)極顯著。計算t值，t＝70.66＊＊，因為t0.01＝3.169＜t所以差異極顯著，結(jié)果和相關(guān)系數(shù)顯著性相吻合。

3.3 2012年全省參測牛場305d產(chǎn)奶量和體細(xì)胞的相關(guān)分析。數(shù)據(jù)見表6。

表5 農(nóng)得利2013年305天產(chǎn)奶量和體細(xì)胞統(tǒng)計表

經(jīng)計算，r＝－0.08，因為df＝10時，r0.05＝0.576，r＜r0.05，所以負(fù)相關(guān)不顯著。計算t值，t＝0.253，因為t0.05＝2.228＞t所以差異不顯著，結(jié)果和相關(guān)系數(shù)顯著性相吻合。

4 產(chǎn)奶量和體細(xì)胞相關(guān)系數(shù)匯總和分析

4.1 匯總

測定日均產(chǎn)奶量、月均產(chǎn)奶量、305d均產(chǎn)奶量和體細(xì)胞數(shù)的相關(guān)系數(shù)及顯著性見表7。

4.2 分析

當(dāng)平均體細(xì)胞高于某個臨界值時（15～20萬／ml），體細(xì)胞和奶產(chǎn)量呈負(fù)相關(guān)；體細(xì)胞低于此臨界值時則相關(guān)不明顯，甚至出現(xiàn)不顯著正相關(guān)（如常興牛場），這個閾值應(yīng)該是健康奶牛體細(xì)胞的生理指標(biāo)值。超越正常生理指標(biāo)值時就會出現(xiàn)負(fù)相關(guān)，超出越多負(fù)相關(guān)越顯著（見表7）。

表7 奶產(chǎn)量和體細(xì)胞數(shù)相關(guān)系數(shù)匯總表＊

5 結(jié)論

由相關(guān)系數(shù)的計算結(jié)果可見，在體細(xì)胞高于正常生理指標(biāo)時，產(chǎn)奶量和體細(xì)胞呈負(fù)相關(guān)關(guān)系，當(dāng)體細(xì)胞值超越到一定范圍時，則二者之間的關(guān)系呈強(qiáng)負(fù)相關(guān)或極強(qiáng)負(fù)相關(guān)（見表7），這個范圍在49.28萬／ml左右，其產(chǎn)量與體細(xì)胞呈強(qiáng)負(fù)相關(guān)，超過68.32萬／ml左右時，呈極強(qiáng)負(fù)相關(guān)。這一結(jié)論與測定牛群的胎次有很大關(guān)系，當(dāng)一胎牛占的比例大時，范圍會降低，甚至40多萬時就表現(xiàn)出極強(qiáng)負(fù)相關(guān)，當(dāng)二、三胎牛較多時，范圍會提高，甚至40多萬才表現(xiàn)強(qiáng)負(fù)相關(guān)。

［1］ 山田淳三.統(tǒng)計方法在畜牧上的應(yīng)用［M］.上海：上?？茖W(xué)技術(shù)出版社，1965.

［2］ 周亞平，劉琴.乳體細(xì)胞與產(chǎn)奶量、乳成分的關(guān)系研究［J］.中國奶牛，2011（4）：40－42.

［3］ 張慧林，等.牛乳中體細(xì)胞與產(chǎn)奶量和乳成分的相關(guān)分析［J］.西北農(nóng)業(yè)學(xué)報，2010，19（4）：1－4.