Dopad operace levostranné dislokace slezu metodou omentopexe na užitkovost a reprodukci dojnic holfštýnského typu Jan Šterc1 , Jaroslav Marek2 1
e-mail:
1
Křídlovecká 16 603 00 Brno
[email protected], 2
[email protected]
Klíčová slova: levostranná dislokace slezu, dojivost, laktační cyklus, nelineární regrese, ANOVA, kontingenční tabulky
1
Úvod
Poznámky k zamyšlení a pro korektury: 1. V literatuře [4] jsou krávy děleny do skupin O1, O2, O3. Při korekturách by asi bylo vhodné označení skupin krav A, B, C, D vhodně přeznačit. 2. Pro porovnání dojivosti léčených a neléčených krav jsem jako zdroj pro průměrnou hodnotu v jednotlivých laktacích u neléčených krav použil výpočet ze souboru uzitkovost vyrazene09. Výpočet byl proveden postupem popsaným v článku [3]. Nemůžete získat z databáze přímo hodnoty celkového nádoje každé krávy ve všech laktací u krav vyřazených v r. 2009, resp. v jiných letech? 3. Poznatky shrnuté v závěru lze pomocí statistických testů ještě ověřit. Testování je poměrně pracné a má jistá úskalí. Pokud se rozhodneme pro jejich aplikaci, musíme společně precizně zformulovat zadání pro vlastní testování.
Léčené krávy byly klasifikovány do 4 skupin. skupina vlastnost A neúspěšně léčená kráva B dokončila laktaci a zabřezla C dokončila laktaci a nezabřezla D nedokončila laktaci 1
2 2.1
Vliv na zařazení léčené krávy do dalšího chovu Analýza přežívání neléčených dojnic v chovu
Z dat (uzitkovost vyrazene09.txt) obsahujících údaje o 278 dojnicích vyřazených v průběhu roku 2009 v kravínech firmy . . . ., zemědělská a.s. jsme získali strukturu dojnic v jednotlivých chovech. Viz histogram na obrázku 1.
Obr. 1: Histogram s počtem dojnic vyřazených v jednotlivých laktačních cyklech. V následující tabulce určíme pravděpodobnosti přežití, resp. vyřazení dojnice, v 1. až 8. laktačním cyklu. Konkrétně P1 je pravděpodobnost, že dojnice zařazená do chovu bude z chovu vyřazena v uvedeném laktačním cyklu. P2 je pravděpodobnost, že dojnice v daném laktačním cyklu nedojde do dalšího laktačního cyklu. laktace počet krav pravděpodobnost P1 pravděpodobnost P2 1 74 0.2662 0.2662 2 84 0.3022 0.4118 3 58 0.2086 0.4833 4 38 0.1367 0.6129 5 15 0.0540 0.6250 6 3 0.0108 0.3333 7 5 0.0180 0.8333 8 1 0.0036 1.0000 Tab. 1: pravděpodobnosti vyřazení kravy z chovu v jednotlivých laktacích. Na základě těchto údajů můžeme zkonstruovat aproximaci distribuční funkce a funkce přežití. 2
Obr. 2: distribuční funkce, tj. pravděpodobnost vyřazování.
Obr. 3: funkce přežití.
3
2.2
Analýza přežívání léčených dojnic v chovu
41 léčených krav bylo ve třídách a laktačních cyklech dle následující tabulky: P laktace/skupina A B C D 1 3 6 0 2 11 2 2 5 3 3 13 3 1 7 3 0 11 4 0 0 1 2 3 5 0 0 1 1 2 P 6 18 8 8 41 Tab. 2: klasifikace léčených krav Z tab. 2 určíme pravděpodobnosti vyřazování operovaných krav z chovu. laktace pravděpodobnost P2L 1 1 - 6/11 = 0.4545 2 1 - 5/13 = 0.6154 3 1 - 7/11 = 0.3636 4 1 - 0/3 = 1.0000 5 1 - 0/2 = 1.0000 Tab. 3: pravděpodobnosti vyřazení operované kravy z chovu v jednotlivých laktacích. Výsledky z tab. 1 a 2 porovnáme v části 4 Závěr.
3 3.1
Vliv operace na dojivost vyléčené krávy Analýza dojivosti neléčených dojnic v chovu
Z dat (soubor uzitkovost vyrazene09.txt) obsahujících údaje o 278 dojnicích vyřazených v průběhu roku 2009 v kravínech firmy . . . ., zemědělská a.s. jsme získali odhad průměrné dojivosti krav v různých laktacích, viz následující tabulka. Pro odhad dojivosti byl použit následující funkční předpis y(x) = θ1 +
2θ2 θ3 , pro x ∈ h1, 325i , π 4(x − θ4 )2 + θ32
(1)
kde x je pořadové číslo označující den laktace. Viz metoda popsaná v článku [3].
4
Obr. 4: závislost dojivosti na dnu laktace. laktace prům. dojivost 1 10502.66 2 11436.35 3 11547.94 4 10321.60 5 11113.31 6 10505.08 7 10774.86 Tab. 4: průměrná dojivost neléčených krav Průměrnou dojivost získanou z celého chovu porovnáme s dojivostí léčených krav.
3.2
Analýza dojivosti léčených dojnic
Průměrná dojivost ve skupinách B a C byla následující laktace prům. dojivost sk. B prům. dojivost sk. C prům. dojivost sk. B+C 1 10288.17 NaN 10288.17 2 11828.40 10629.33 11378.75 3 11622.43 9787.33 11071.90 4 NaN NaN 11213.50 5 NaN NaN NaN Tab. 5: průměrná dojivost léčených krav ve skupinách B, C a B+C Symbol NaN je použit v těch skupinách, kdy nešlo vzhledem k malému počtu dat (žádná nebo jedna kráva) ve skupině a laktaci určit průměrnou hodnotu.
5
4
Závěr – hypotézy
Porovnáním pravděpodobností vyřazení neléčených a léčených krav v jednotlivých laktacích uvedených v tab. 1 (poslední sloupec) a v tab. 3 zjišťujeme, že vyřazování léčených krav je vyšší. Pravděpodobnost je vyšší cca o 20% v první a druhé laktaci a cca o 15% ve třetí laktaci. Porovnáním průměrné dojivosti v různých laktačních cyklech uvedených v tab. 4 a 5 zjišťujeme, že dojivost neléčených a léčených krav je shodná.
5
Literatura [1 ] Anděl, J.: Základy matematické statistiky, Univerzita Karlova v Praze, Praha 2002 [2 ] Kubáček, L., Kubáčková, L.: Statistika a metrologie, UP Olomouc, 2000 [3 ] Zelinková, G., Marek, J.: Dependence of milk yield, fat: protein ratio and somatic cell counts. Proceedings of XI. Middle-European Buiatrics Congress, 17. 6.–19. 6. 2010, Brno, Czech republic [4 ] Haloun, T.: Vliv vybraných chirurgických zákroků u dojnic na jejich následnou užitkovost a reprodukční ukazatele, zpráva, 2010
6
