(Viene de la página anterior)
Página anterior | ![]() Volver al principio del trabajo | Página siguiente ![]() |
Ambas rectas pasan por el
origen, y los intervalos de confianza calculados para sus pendientes no
incluyen el punto 0,25 que sería el valor ideal al pensar que el CMTM es
simplemente el valor promedio de los CMT individuales de los cuatro cuartos. De
este modo, las estimaciones del CMTM resultaron ser estadísticamente inferiores
al promedio del CMT, particularmente en vaquillas. También el grado de ajuste
de la curva fue menor en este tipo de animal.
Las cifras de la segunda
columna de ambos cuadros son fácilmente interpretables: cuando la suma de CMT es
cero, el CMTM es cero. En vacas (Cuadro 3), cuando la suma de los CMT es 4, es
decir un cuarto en promedio con nota 1, la estimación del CMTM da 0,9. Cuando
es 8, da un valor de 1,8, un tanto inferior al promedio simple de 2.
b) Modelo logit
La segunda alternativa de
estimación se realizó mediante el modelo logit. En este caso la función
estimada general tomó la forma:
pj = 1 / (1 +
EXP(b SCMT - ck)), con:
b = 0,5794; c1 =
2,6341; c2 = 4,2041 y c3 = 6,0633; R2 = 0,82
para vaquillas, y
b = 0, 6242; c1 =
2,3686; c2 = 4,1423; c3 = 6,3403 y c4 =
11,8273; R2 = 0,76 para vacas
A partir de estas curvas
logísticas se calcularon las probabilidades de cada uno de los CMTM, usando los
parámetros del modelo:
Prob (CMTM = 0) = 1 / (1 +
EXP(b SCMT - c1))
Prob (CMTM = x) = 1 / (1 +
EXP(b SCMT - cx+1)) - 1 / (1 + EXP(b SCMT - cx))
Para valores de x entre 1 y
2; y
Prob ( CMTM = 3) = 1 / (1 +
EXP(b SCMT - c3))
Estos valores multiplicados
por 100 son los que se presentan en las columnas
La capacidad predictiva de
este modelo parece ser mejor y da más información que el modelo lineal
anterior. Así, en vaquillas (Cuadro 2), cuando la suma de los CMT de los
cuartos es 0, el modelo logit predice un CMTM de 0 con probabilidad 93,3 % y un
CMTM 1 con probabilidad 5,2 %. De este modo, la predicción desde SCMT = 0 puede
resumirse afirmando que existe un 98,5 % de probabilidad de tener un CMTM de 0
ó 1. Los valores así calculados conforman las dos últimas columnas de ambos
cuadros.
Observando estas columnas
puede advertirse que, agrupando de dos en dos los valores de CMTM contiguos, la
probabilidad resultó siempre bastante alta, es decir existe una alta certeza de
obtener dichos resultados. Al mismo tiempo esta información confirma la
consistencia del modelo, ya que los valores contiguos se desplazan de izquierda
a derecha a medida que
Algunos valores
particulares resultan interesantes de destacar, especialmente para las sumas
extremas de CMT. En vaquillas para una suma de CMT de 16, el valor más alto
posible de la variable, la probabilidad de tener un CMTM de 2 ó 3 es de 99,4%.
En vacas, para una suma de cero CMT la probabilidad de tener un CMTM entre 0 y
1 es de 98,4%. Análogamente, en vacas para una suma de 16 se tiene una
probabilidad de 97,5% de lograr un CMT entre 3 y 4.
Si se analizan los valores
intermedios de las sumas de los CMT también las probabilidades son altas, (
Estimación de los rangos
celulares a partir de SCMT
a) Modelo lineal
Las curvas de regresión
lineal que pasan por el origen para los RC a partir de la suma de los CMT de
los cuartos resultaron ser:
RC = 0,74 SCMT; R2
= 0,26 para vaquillas, y
RC = 0,74 SCMT; R2
= 0,47 para vacas.
En este caso, aunque los
dos coeficientes de regresión resultaron significativos al 5%, el grado de ajuste
es bajo, especialmente en vaquillas, de modo que las estimaciones no parecen
ser muy precisas bajo este modelo lineal. Cabe destacar, sin embargo, el hecho
que las pendientes resultaron numéricamente iguales para ambos tipos de animal.
Los valores estimados de estas curvas se presentan en la segunda columna de los
Cuadros 4 y
Cuadro 4. Estimación de los rangos de recuentos celulares (RC) en función de la
suma de los cuartos individuales (SCMT), usando una función lineal y el modelo
logit para vaquillas.
Table 4. Estimation of somatic cell recounts (RC) in dairy heifers as a
function of the addition of individual quarter CMT scores (SCMT), using a
linear function and a logit model.
|
SCMT |
RC1 |
Probabilidades estimadas para valores de RC usando el modelo logit RC |
Valores RC más probables |
Probabilidad |
|||||||||
|
|
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
|
|
0 |
0,0 |
18,6 |
16,5 |
20,4 |
19,1 |
12,1 |
7,2 |
3,6 |
1,6 |
0,6 |
0,3 |
|
74,6 |
|
1 |
0,7 |
14,1 |
13,8 |
19,3 |
20,6 |
14,6 |
9,4 |
4,9 |
2,2 |
0,8 |
0,4 |
|
68,3 |
|
2 |
1,5 |
10,5 |
11,2 |
17,3 |
21,0 |
16,9 |
11,9 |
6,5 |
3,0 |
1,1 |
0,6 |
|
67,1 |
|
3 |
2,2 |
7,7 |
8,8 |
14,8 |
20,4 |
18,7 |
14,6 |
8,6 |
4,1 |
1,5 |
0,8 |
|
68,5 |
|
4 |
3,0 |
5,6 |
6,8 |
12,2 |
18,8 |
19,7 |
17,2 |
11,0 |
5,5 |
2,0 |
1,2 |
|
66,7 |
|
5 |
3,7 |
4,1 |
5,1 |
9,7 |
16,5 |
19,5 |
19,5 |
13,8 |
7,3 |
2,8 |
1,6 |
|
69,3 |
|
6 |
4,4 |
3,0 |
3,8 |
7,5 |
13,9 |
18,4 |
20,9 |
16,8 |
9,6 |
3,8 |
2,2 |
|
69,5 |
|
7 |
5,2 |
2,1 |
2,8 |
5,7 |
11,2 |
16,5 |
21,4 |
19,5 |
12,4 |
5,2 |
3,1 |
|
75,0 |
|
8 |
5,9 |
1,5 |
2,0 |
4,3 |
8,8 |
14,2 |
20,7 |
21,8 |
15,5 |
7,0 |
4,3 |
|
79,2 |
|
9 |
6,7 |
1,1 |
1,5 |
3,2 |
6,8 |
11,6 |
19,0 |
23,0 |
18,7 |
9,2 |
5,9 |
|
81,5 |
|
10 |
7,4 |
0,8 |
1,1 |
2,3 |
5,1 |
9,3 |
16,6 |
23,2 |
21,7 |
11,8 |
8,1 |
|
73,3 |
|
11 |
8,1 |
0,6 |
0,8 |
1,7 |
3,8 |
7,2 |
14,0 |
22,1 |
24,2 |
14,8 |
11,0 |
|
75,1 |
|
12 |
8,9 |
0,4 |
0,6 |
1,2 |
2,8 |
5,5 |
11,3 |
20,0 |
25,5 |
18,0 |
14,7 |
|
78,2 |
|
13 |
9,6 |
0,3 |
0,4 |
0,9 |
2,0 |
4,1 |
8,9 |
17,3 |
25,7 |
21,0 |
19,4 |
|
83,4 |
|
14 |
10,4 |
0,2 |
0,3 |
0,6 |
1,5 |
3,0 |
6,8 |
14,4 |
24,5 |
23,5 |
25,2 |
|
87,6 |
|
15 |
11,1 |
0,1 |
0,2 |
0,5 |
1,1 |
2,2 |
5,1 |
11,5 |
22,2 |
25,0 |
32,0 |
|
90,7 |
|
16 |
11,8 |
0,1 |
0,1 |
0,3 |
0,8 |
1,6 |
3,8 |
9,0 |
19,2 |
25,3 |
39,7 |
|
93,2 |
Página anterior | ![]() Volver al principio del trabajo | Página siguiente ![]() |
Trabajos relacionados
Ver mas trabajos de Agricultura y Ganaderia |
|
Nota al lector: es posible que esta página no contenga todos los componentes del trabajo original (pies de página, avanzadas formulas matemáticas, esquemas o tablas complejas, etc.). Recuerde que para ver el trabajo en su versión original completa, puede descargarlo desde el menú superior.