一種通過牛體表溫度計(jì)算牛體溫的方法
【專利摘要】本發(fā)明提供了一種通過牛體表溫度計(jì)算牛體溫的方法，包括使用如模型1所示的混合線性模型計(jì)算牛直腸溫度，模型1：Y＝a+X0b0+Xsbs+Xtbt+Xibi+e；其中，Y為直腸溫度，X0為體表溫度，Xs為季節(jié)設(shè)計(jì)矩陣，bs為Xs的非隨機(jī)參數(shù)向量；Xt為時(shí)間設(shè)計(jì)矩陣，b0為X0的回歸系數(shù)，bt為Xt的非隨機(jī)參數(shù)向量；Xi為個(gè)體設(shè)計(jì)矩陣，bi為Xi的隨機(jī)參數(shù)向量；a代表截距，e為隨機(jī)向量誤差。本發(fā)明在將體表溫度對(duì)直腸溫度回歸分析時(shí)，以季節(jié)因素與時(shí)間點(diǎn)因素為固定效應(yīng)，以個(gè)體因素為隨機(jī)效應(yīng)。運(yùn)用該模型校正后的體表溫度與實(shí)際直腸溫度比較，各校正溫度與直腸溫度差值絕對(duì)數(shù)的平均數(shù)達(dá)到0.104℃，90％數(shù)據(jù)差值的絕對(duì)值不超過0.2℃，取得了理想的效果。
【專利說明】
一種通過牛體表溫度計(jì)算牛體溫的方法
[00011本申請(qǐng)要求2015年7月2日提交的中國專利申請(qǐng)201510383183.4的優(yōu)先權(quán)，將其通 過引用并入本文。
技術(shù)領(lǐng)域
[0002] 本發(fā)明涉及畜牧養(yǎng)殖領(lǐng)域，具體的涉及一種通過牛體表溫度計(jì)算牛體溫的方法。
【背景技術(shù)】
[0003] 體溫是機(jī)體重要的生理指標(biāo)，正常情況下，牛體溫會(huì)在相對(duì)恒定的范圍內(nèi)，隨生理 活動(dòng)變化而成一定規(guī)律性變化(Cooper-Prado et al. 2011)，比如發(fā)情、排卵、妊娠及分娩 等。然而，在牛體溫檢測(cè)的實(shí)踐中，大部分牛場(chǎng)目前依然采用水銀體溫計(jì)進(jìn)行直腸測(cè)溫，該 方法測(cè)溫效率低下，易引起奶牛疾病交叉?zhèn)鞑?，且不能?zhǔn)確實(shí)時(shí)獲取每頭牛的體溫?cái)?shù)據(jù)，需 要的人力也較多，耗時(shí)費(fèi)力，遠(yuǎn)遠(yuǎn)不能滿足現(xiàn)代規(guī)?；翀?chǎng)的管理要求(朱翔宇，2013)。隨 著生物傳感技術(shù)的發(fā)展，部分發(fā)達(dá)國家已將無線傳感技術(shù)應(yīng)用于奶牛體表測(cè)溫，并取得了 一定的研究成果。Bligh等1974年首次提出將無線遙測(cè)技術(shù)運(yùn)用于動(dòng)物體溫測(cè)定的構(gòu)想 (Bligh et al.l974);Kyle等與Suthar等（Kyle et al.l998;Suthar et al.2011)在奶牛 陰道植入含溫度采集、無線電接收設(shè)備的棒狀無線遙測(cè)系統(tǒng)，對(duì)奶牛陰道溫度進(jìn)行實(shí)時(shí)檢 測(cè)；Raul等與Miranda等（Raul et al.2006;Miranda et al.2009)將與會(huì)陰肌肉組織無排 斥反應(yīng)的無線溫度傳感器，植入奶牛會(huì)陰部肌肉，檢測(cè)會(huì)陰部溫度，進(jìn)行發(fā)情鑒定;AlZahal 等與Edouard等(AlZahal et al.2009,2010;Edouard et al.2011)分別通過食道將無線電 瘤胃丸投入瘤胃，檢測(cè)奶牛瘤胃溫度變化，結(jié)果顯示，無線電瘤胃丸能夠有效檢測(cè)奶牛體內(nèi) 溫度變化。
[0004] 上述方法雖不同程度地達(dá)到了測(cè)溫目的，但無論瘤胃丸還是皮下埋置都會(huì)對(duì)奶牛 造成一定傷害，均不符合動(dòng)物福利要求。近年來，紅外測(cè)溫不斷被用于豬體表溫度測(cè)量，然 而效果都不很理想，無法指導(dǎo)生產(chǎn)應(yīng)用?，F(xiàn)代醫(yī)學(xué)的發(fā)展，推動(dòng)了對(duì)人體體表溫度變化規(guī)律 的深入研究，并取得了一定階段性成果，然而，牛體表被毛較多，無毛裸露的體表部位較少， 并且檢測(cè)很不方便，體表溫度研究報(bào)道較少，僅見到關(guān)于熱應(yīng)激對(duì)奶牛體表溫度的影響(溫 雅俐，2011)和乳區(qū)表面溫度變化與乳房炎關(guān)系(Metzner et al.2014)等少量研究，這些結(jié) 果只揭示了部分環(huán)境因素及體表個(gè)別部位溫度與部分病理現(xiàn)象的關(guān)系，沒有深入揭示體表 溫度與體溫(肛溫)間的關(guān)系及其變化規(guī)律，對(duì)牛體表溫度變化規(guī)律缺乏足夠深入系統(tǒng)的研 究。
[0005] 因此需要提供一種通過牛體表溫度計(jì)算牛體溫(肛溫）的方法，為方便的測(cè)量牛體 溫提供技術(shù)手段。

【發(fā)明內(nèi)容】

[0006] 本發(fā)明的目的在于克服現(xiàn)有技術(shù)中存在的缺陷，提供一種能夠準(zhǔn)確簡便的通過牛 體表溫度計(jì)算牛體溫的方法。
[0007] 本發(fā)明提供了一種通過牛體表溫度計(jì)算牛體溫的方法，包括使用如模型1所示的 混合線性模型計(jì)算牛直腸溫度，
[0008] 模型1:模型1: Y = a+Xobo+Xsbs+Xtbt+Xibi+e ;模型1中，Y為直腸溫度，Xo為牛腿跖部 溫度，Xs為季節(jié)設(shè)計(jì)矩陣，bs為Xs的非隨機(jī)參數(shù)向量;Xt為時(shí)間設(shè)計(jì)矩陣，bo為Χο的回歸系數(shù)， bt為Xt的非隨機(jī)參數(shù)向量;Xi為個(gè)體設(shè)計(jì)矩陣，bi為Xi的隨機(jī)參數(shù)向量;a代表截距，e為隨機(jī) 向量誤差。
[0009] 可選的，所述牛腿跖部為牛后腿跖部。
[0010] 可選的，所述牛腿跖部為牛腿跖骨下部近籽骨上方。
[0011] 可選的，所述方法包括將體溫檢測(cè)裝置固定于牛腿跖部，實(shí)時(shí)檢測(cè)牛只溫度數(shù)據(jù)。
[0012] 可選的，所述方法包括利用體溫檢測(cè)裝置每小時(shí)檢測(cè)2-10次。
[0013] 本發(fā)明還提供了所述計(jì)算牛體溫的方法在牛機(jī)體生理情況檢測(cè)中的應(yīng)用。
[0014] 可選的，所述應(yīng)用包括檢測(cè)牛腿跖部溫度通過混合線性模型計(jì)算得到牛直腸溫 度，利用所述牛直腸溫度作為牛機(jī)體生理情況的指標(biāo)。
[0015] 可選的，牛機(jī)體生理情況包括發(fā)情、排卵、妊娠、分娩。
[0016] 具體的，可以將利用本發(fā)明所提供的方法計(jì)算得到的牛直腸溫度應(yīng)用于了解被測(cè) 動(dòng)物的生理情況，從而為指導(dǎo)畜牧生產(chǎn)提供依據(jù)。值得注意的是，其直接目的不是為了獲得 動(dòng)物疾病的診斷結(jié)果或健康狀況。
[0017] 本發(fā)明以檢測(cè)時(shí)間和季節(jié)因素作為固定效應(yīng)，個(gè)體因素作為隨機(jī)效應(yīng)對(duì)體表溫度 和直腸溫度數(shù)據(jù)進(jìn)行混合線性模型分析，得到了最佳回歸效果的混合線性模型:模型1:Y = a+Xobo+Xsbs+Xtbt+Xibi+e，運(yùn)用該模型校正后的體表溫度與實(shí)際直腸溫度比較，各校正溫度 與直腸溫度差值絕對(duì)數(shù)的平均數(shù)為0.104°C，90%數(shù)據(jù)差值的絕對(duì)值不超過0.2°C，說明該 擬合方程與本測(cè)溫系統(tǒng)結(jié)合，可以取得較理想的測(cè)溫結(jié)果。
【附圖說明】
[0018] 圖1為體表測(cè)溫裝置安裝示意圖；
[0019]圖2為體表測(cè)溫設(shè)備與水銀溫度計(jì)測(cè)溫一致性分析的Bland-Altman圖；
[0020] 圖3為不同季節(jié)牛一天內(nèi)體表溫度和直腸溫度變化曲線圖。
【具體實(shí)施方式】
[0021] 下面將通過【具體實(shí)施方式】對(duì)本發(fā)明進(jìn)行詳細(xì)說明。需要理解的是以下實(shí)施例的給 出僅是為了起到說明的目的，并不是用于對(duì)本發(fā)明的范圍進(jìn)行限制。本領(lǐng)域的技術(shù)人員在 不背離本發(fā)明的宗旨和精神的情況下，可以對(duì)本發(fā)明進(jìn)行各種修改和替換。
[0022] 實(shí)施例1
[0023]本實(shí)施例用于說明利用本發(fā)明所提供的方法計(jì)算牛體溫。
[0024] 1材料方法
[0025] 1.1試驗(yàn)時(shí)間與地點(diǎn)
[0026] 分別于2014年冬（2月2日至10日；winter，Win·)、夏（6月2日至10日；summer， Sum.)、秋（10月2日至10日；autumn，Aut.)三個(gè)季節(jié)，在河北石家莊某牛場(chǎng)進(jìn)行體表溫度與 直腸溫度測(cè)定，因?yàn)榇杭九c秋季氣候條件相似，所以，僅選擇秋季進(jìn)行溫度測(cè)定分析，以代 表春秋季測(cè)定效應(yīng)。
[0027] 1.2試驗(yàn)動(dòng)物及其飼養(yǎng)管理
[0028]每個(gè)季節(jié)分別隨機(jī)選取7頭健康、產(chǎn)后60至90天二胎正?？煞蔽鏖T塔爾牛作為實(shí) 驗(yàn)牛;試驗(yàn)牛采用散欄式飼養(yǎng)、全自動(dòng)TMR日糧飼喂，每天分別在2:00、6:00、12:00、18:00、 22:00投放飼料、保證自由采食和清潔飲水。
[0029] 1.3測(cè)溫設(shè)備及其測(cè)試與校正
[0030] 市售水銀溫度計(jì)1支，獸用水銀體溫計(jì)2支，體表溫度自動(dòng)測(cè)定系統(tǒng)，包括體表溫度 檢測(cè)裝置7套(該裝置采用接觸式測(cè)溫）、數(shù)據(jù)采集裝置1套(該裝置用于采集檢測(cè)裝置中存 儲(chǔ)的體溫?cái)?shù)據(jù)）。如圖1所示，該系統(tǒng)采用接觸式測(cè)溫方式，選擇牛體表被毛稀疏的后肢蹄腕 部為體表測(cè)溫部位，將蹄腕部肌肉、血管較豐厚部位（即牛腿跖部，牛腿跖骨下部近籽骨上 方)作為溫度傳感器固定位置，并根據(jù)該部位形狀設(shè)計(jì)了具有一定弧度的固定外殼，通過外 殼與腿部緊密貼合使傳感器能夠與肌體緊密接觸。同時(shí)，又在外殼中與肌體貼合一側(cè)的中 心設(shè)計(jì)了一個(gè)凹巢，溫度傳感器固定于凹巢內(nèi)，進(jìn)一步降低了溫度、風(fēng)力等外界環(huán)境因素對(duì) 測(cè)溫效果的影響，檢測(cè)獲得的溫度。
[0031]為使實(shí)驗(yàn)結(jié)果準(zhǔn)確科學(xué)，每次試驗(yàn)前將7個(gè)體表溫度檢測(cè)裝置置于恒溫箱中，進(jìn)行 測(cè)試實(shí)驗(yàn)。首先設(shè)定恒溫箱溫度從30°C升至40°C (每0.5h升溫1°C)，測(cè)溫設(shè)備每0.5h自動(dòng)測(cè) 溫一次。試驗(yàn)期間，使用膠帶將水銀溫度計(jì)固定在恒溫箱玻璃門內(nèi)側(cè)，在自動(dòng)測(cè)溫時(shí)間點(diǎn)， 同時(shí)讀取并記錄水銀溫度計(jì)顯示的恒溫箱內(nèi)溫度。該測(cè)試實(shí)驗(yàn)重復(fù)三次，證明體表溫度自 動(dòng)測(cè)定系統(tǒng)準(zhǔn)確可靠后，進(jìn)行后續(xù)體表溫度自動(dòng)測(cè)定實(shí)驗(yàn)。
[0032] 1.4體表溫度和體內(nèi)溫度的采集
[0033] 預(yù)試驗(yàn)：正式試驗(yàn)開始前兩天，將體表溫度（surface temperature，ST)檢測(cè)裝置 緊密固定于實(shí)驗(yàn)牛后腿蹄腕部（即牛后腿跖部），并將測(cè)溫間隔設(shè)定為半小時(shí)，每4個(gè)小時(shí)用 數(shù)據(jù)采集裝置收集一次溫度數(shù)據(jù);根據(jù)溫度采集結(jié)果進(jìn)行設(shè)備調(diào)試，同時(shí)進(jìn)行試驗(yàn)牛適應(yīng) 性訓(xùn)練。待設(shè)備測(cè)溫正常穩(wěn)定、試驗(yàn)牛適應(yīng)后，進(jìn)行正式試驗(yàn)。
[0034] 正式試驗(yàn)：自動(dòng)測(cè)溫裝置繼續(xù)沿用預(yù)試期的半小時(shí)測(cè)溫設(shè)置，連續(xù)測(cè)溫3天；同時(shí)， 采用水銀體溫計(jì)，每2小時(shí)測(cè)定一次試驗(yàn)牛直腸溫度(rectal temperature，RT)。測(cè)定直腸 溫度時(shí)，先將試驗(yàn)牛保定，待其休息10分鐘后將溫度計(jì)插入直腸，5分鐘后取出溫度計(jì)，讀取 數(shù)據(jù)，每頭牛每天共采集12個(gè)時(shí)間點(diǎn)的直腸溫度數(shù)據(jù)，也連續(xù)測(cè)溫3天。
[0035] 1.5數(shù)據(jù)處理與分析
[0036]設(shè)備實(shí)驗(yàn)室測(cè)試數(shù)據(jù)，使用SAS9.2進(jìn)行差異顯著性分析與線性回歸分析，并使用 Bland一Altman法進(jìn)行一致性評(píng)價(jià)。
[0037]每頭牛每天有48條體表溫度數(shù)據(jù)與12條直腸溫度數(shù)據(jù)。因直腸溫度兩個(gè)小時(shí)內(nèi)變 化較小，故取相鄰兩個(gè)時(shí)間點(diǎn)的直腸溫度平均值作為中間時(shí)間點(diǎn)的直腸溫度，并舍棄所有 非整點(diǎn)體表溫度數(shù)據(jù)，使每頭牛每天各有23個(gè)相同時(shí)間點(diǎn)的體表溫度與直腸溫度數(shù)據(jù)，三 個(gè)季節(jié)一共獲得1449條相同時(shí)間點(diǎn)的體表溫度與直腸溫度數(shù)據(jù)。采用SAS 9.2進(jìn)行體溫?cái)?shù) 據(jù)差異顯著性分析與協(xié)方差分析，同時(shí)采用TableCUrVe2D/3D進(jìn)行線性回歸分析與非線性 回歸分析，并進(jìn)行體表溫度校正分析。
[0038] 2結(jié)果與分析
[0039] 2.1體表溫度檢測(cè)裝置的測(cè)試分析
[0040]實(shí)驗(yàn)室測(cè)試結(jié)果表明，恒溫箱溫度在30~40°C范圍內(nèi)變化時(shí)，7個(gè)檢測(cè)裝置的測(cè)定 結(jié)果發(fā)生了同步變化，且差異不顯著，自動(dòng)測(cè)溫結(jié)果與恒溫箱內(nèi)實(shí)際溫度呈強(qiáng)相關(guān)（R2: 0.9987)。測(cè)溫結(jié)果一致性金標(biāo)準(zhǔn)分析得到的Bland-Altman圖顯示，7個(gè)測(cè)溫設(shè)備間的差異 僅為0 ±0.07°C (95%置信區(qū)間為-0.22~0.22°C，圖2)，可以忽略不計(jì)。說明該裝置的準(zhǔn)確 性和穩(wěn)定性較好，可用于后續(xù)體表溫度研究。
[0041 ] 2.2體表溫度與直腸溫度的季節(jié)及測(cè)定日差異分析
[0042]跖部體表溫度與直腸溫度統(tǒng)計(jì)分析結(jié)果(表1、2)顯示，直腸溫度變化幅度較小，各 季節(jié)間、各測(cè)定日間均無顯著差異(P>〇.05)，說明生理?xiàng)l件下，直腸溫度基本恒定。各測(cè)定 日間體表溫度(表2)差異也不顯著(p>0.05)，說明自動(dòng)測(cè)溫裝置性能較為穩(wěn)定。然而，不同 于直腸溫度，跖部體表溫度測(cè)定結(jié)果季節(jié)差異顯著(P<〇.05)，其中，夏季溫度最高，冬季溫 度最低，各季節(jié)的測(cè)定結(jié)果均顯著低于直腸溫度測(cè)定結(jié)果(P < 〇. 〇 5)。
[0043] 表1體表溫度與直腸溫度統(tǒng)計(jì)分析結(jié)果
[0044] 單位:°C
[0045]
[0046] 注:A、B、C表示同行數(shù)據(jù)的差異顯著性，字母相同者表示差異不顯著(p>0.05)，字 母不同者表示差異顯著(P<〇.05);
[0047] 表2相同季節(jié)三天體表溫度與直腸溫度情況
[0048]
[0050]注:A、B、C表示同列數(shù)據(jù)的差異顯著性，字母相同者表示差異不顯著(p>0.05)，字 母不同者表示差異顯著(P<〇.05);
[0051 ] 2.3妬部體表溫度和直腸溫度的晝夜差異分析
[0052] 表3體表溫度與直腸溫度晝夜統(tǒng)計(jì)分析
[0053]
[0054] 注:A、B表示同列數(shù)據(jù)的差異顯著性，字母相同者表示差異不顯著(p>0.05)，字母 不同者表示差異顯著(P<〇.05);
[0055] 實(shí)驗(yàn)牛不同季節(jié)各測(cè)定時(shí)間點(diǎn)溫度曲線（圖3)與統(tǒng)計(jì)分析結(jié)果(表3)表明，體表溫 度受時(shí)間因素影響顯著并呈現(xiàn)明顯晝夜變化規(guī)律，5:00至7:00為一天中的最低值，隨后逐 漸升高，到13:00到15:00左右達(dá)到一天的最高值，隨后逐漸降低，直至最低溫度。直腸溫度 變化幅度較小，但變化規(guī)律類似。
[0056] 2.4體表溫度與直腸溫度回歸分析
[0057]假設(shè)體表溫度與直腸溫度間存在線性關(guān)系，將體表溫度對(duì)直腸溫度進(jìn)行回歸分 析?？紤]到單因素方差分析提示不同測(cè)量時(shí)間對(duì)體表測(cè)溫結(jié)果有顯著影響。因此，以季節(jié)因 素與測(cè)量時(shí)間為固定效應(yīng)，個(gè)體為隨機(jī)效應(yīng)，對(duì)體表溫度和直腸溫度數(shù)據(jù)進(jìn)行混合線性模 型分析，建立的混合線性模型如下：
[0058] Y = a+Xob()+Xsbs+Xtbt+Xibi+e其中，Y為直腸溫度，Χο為牛腿跖部溫度（即體表溫度）， Xs為季節(jié)設(shè)計(jì)矩陣，bs為Xs的非隨機(jī)參數(shù)向量;Xt為時(shí)間設(shè)計(jì)矩陣，bo為Χο的回歸系數(shù)，bt為Xt 的非隨機(jī)參數(shù)向量;Xi為個(gè)體設(shè)計(jì)矩陣，bi為Xi的隨機(jī)參數(shù)向量;a代表截距，e為隨機(jī)向量誤 差。采用模型對(duì)體表溫度和直腸溫度進(jìn)行混合線性分析，各參數(shù)的估計(jì)值見表4,模型解釋 的變異為0.562。
[0059] 表4同時(shí)考慮季節(jié)因素測(cè)定時(shí)間因素為固定效應(yīng)的混合線性模型結(jié)果
[0060]
[00611 2.5直腸溫度與校正后體表溫度對(duì)比
[0062]運(yùn)用混合線性模型得出的回歸方程對(duì)體表測(cè)溫?cái)?shù)據(jù)進(jìn)行矯正，各校正溫度與直腸 溫度差值絕對(duì)數(shù)的平均數(shù)為0.104°c，90%數(shù)據(jù)差值的絕對(duì)值不超過0.2°C，說明該擬合方 程與本測(cè)溫系統(tǒng)結(jié)合，可以取得較理想的測(cè)溫結(jié)果。
[0063] 2.6交叉驗(yàn)證
[0064]為驗(yàn)證上述線性模型的預(yù)測(cè)效果，將全部數(shù)據(jù)隨機(jī)平均分為10等份，隨機(jī)選取9份 建立混合線性模型，剩下1份使用建立的模型進(jìn)行預(yù)測(cè)，依次重復(fù)10次，得到每個(gè)體表溫度 對(duì)應(yīng)的預(yù)測(cè)直腸溫度。最后對(duì)預(yù)測(cè)的直腸溫度與真實(shí)直腸測(cè)溫?cái)?shù)據(jù)做相關(guān)分析以評(píng)價(jià)模型 的預(yù)測(cè)效能?？紤]時(shí)間點(diǎn)因素模型的相關(guān)系數(shù)為0.644,明顯高于不考慮時(shí)間點(diǎn)因素模型的 相關(guān)系數(shù)ο. 569。經(jīng)SAS9.3差異顯著性分析后發(fā)現(xiàn)，校正后的體表溫度與真實(shí)直腸溫度間差 異不顯著(P>〇.05)，說明該方法校正后得到了理想的效果，同時(shí)也說明牛體表溫度與直腸 溫度之間存在良好的線性關(guān)系?？梢酝ㄟ^將后續(xù)測(cè)量得到的牛腿跖部溫度帶入混合線性模 型計(jì)算獲得牛直腸溫度，完成牛體溫的測(cè)量。即可以通過體表溫度直接獲得直腸溫度數(shù)據(jù)。
[0065] 對(duì)比例1
[0066] 本對(duì)比例僅考慮季節(jié)為固定效應(yīng)。
[0067] 以季節(jié)因素為固定效應(yīng)，個(gè)體為隨機(jī)效應(yīng)，對(duì)體表溫度和直腸溫度數(shù)據(jù)進(jìn)行混合 線性模型分析。建立的混合線性模型如下：
[0068] Υ = &+Χ(Λ()+Χ1+ΧΛ^Θ其中，Y為直腸溫度，Χο牛腿跖部溫度（即體表溫度），X S為季 節(jié)設(shè)計(jì)矩陣，bs為Xs的非隨機(jī)參數(shù)向量，即固定效應(yīng);Xi為個(gè)體設(shè)計(jì)矩陣，bo為Χο的回歸系數(shù)， 匕為&的隨機(jī)參數(shù)向量，即隨機(jī)效應(yīng);a代表截距，e為隨機(jī)向量誤差。
[0069] 采用模型對(duì)體表溫度和直腸溫度進(jìn)行混合線性分析，得到模型的參數(shù)估計(jì)見表5， 模型解釋的變異為0.399。
[0070]表5僅考慮季節(jié)效應(yīng)因素為固定效應(yīng)的混合線性模型結(jié)果 「00711
[0072]上述方法獲得的擬合度不如實(shí)施例1。
[0073]雖然，上文中已經(jīng)用一般性說明及具體實(shí)施方案對(duì)本發(fā)明作了詳盡的描述，但在 本發(fā)明基礎(chǔ)上，可以對(duì)之作一些修改或改進(jìn)，這對(duì)本領(lǐng)域技術(shù)人員而言是顯而易見的。因 此，在不偏離本發(fā)明精神的基礎(chǔ)上所做的這些修改或改進(jìn)，均屬于本發(fā)明要求保護(hù)的范圍。
【主權(quán)項(xiàng)】
1. 一種通過牛體表溫度計(jì)算牛體溫的方法，其特征在于，包括使用如模型1所示的線性 方程計(jì)算牛直腸溫度， 1? 型 1 ： Y = a+Xobo+Xsbs+Xtbt+Xibi+e ； 模型1中，γ為直腸溫度，xQ為牛腿跖部溫度，xs為季節(jié)設(shè)計(jì)矩陣，bs為Xs的非隨機(jī)參數(shù)向 量;Xt為時(shí)間設(shè)計(jì)矩陣，bo為Χο的回歸系數(shù)，bt為Xt的非隨機(jī)參數(shù)向量;Xi為個(gè)體設(shè)計(jì)矩陣，bi 為Xi的隨機(jī)參數(shù)向量;a代表截距，e為隨機(jī)向量誤差。2. 根據(jù)權(quán)利要求1所述的方法，其特征在于，所述牛腿跖部為牛后腿跖部。3. 根據(jù)權(quán)利要求1或2所述的方法，其特征在于，所述牛腿跖部為牛腿跖骨下部近籽骨 上方。4. 根據(jù)權(quán)利要求3所述的測(cè)量方法，其特征在于，所述方法包括將體溫檢測(cè)裝置固定于 牛腿跖部，實(shí)時(shí)檢測(cè)牛只溫度數(shù)據(jù)。5. 根據(jù)權(quán)利要求4所述的測(cè)量方法，其特征在于，所述方法包括利用體溫檢測(cè)裝置每小 時(shí)檢測(cè)2-10次。6. 權(quán)利要求1-5中任意一項(xiàng)所述的方法在牛機(jī)體生理情況檢測(cè)中的應(yīng)用。7. 根據(jù)權(quán)利要求6所述的應(yīng)用，其特征在于，所述應(yīng)用包括檢測(cè)牛腿跖部溫度通過混合 線性模型計(jì)算得到牛直腸溫度，利用所述牛直腸溫度作為牛機(jī)體生理情況的指標(biāo)。8. 根據(jù)權(quán)利要求6或7所述的應(yīng)用，其特征在于，牛機(jī)體生理情況包括發(fā)情、排卵、妊娠、 分娩。
【文檔編號(hào)】A61B5/01GK106037683SQ201610515375
【公開日】2016年10月26日
【申請(qǐng)日】2016年7月1日
【發(fā)明人】王棟, 蔡勇, 寇紅祥, 陳曉麗
【申請(qǐng)人】中國農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所