Aplicación de modelos lineales mixtos para estudiar la producción de leche en vacas de la raza Lucerna

Application of linear mixed models to study milk yield per day in cows of the Lucerna breed

A Colombia se le considera como un país muy diverso y con múltiples opciones productivas a nivel ganadero. La raza colombiana Lucerna se caracteriza por presentar buenos desempeños productivos y reproductivos en el trópico bajo de Colombia, debido a su buena rusticidad. Existen pocos estudios que permitan identificar las ventajas productivas que se pueden obtener, al utilizar este tipo de animales. El estudio de la producción de leche por día es una herramienta importante para evaluar la eficiencia de un hato lechero; esta variable es afectada por factores genéticos y no genéticos. Entre tanto, la evaluación de Condición Corporal (CC) por ejemplo, ha sido propuesta como herramienta en el manejo nutricional, reproductivo y de salud en hatos bovinos, por ser un buen indicador del estado nutricional de la vaca y por ende, del nivel de desempeño que pueden tener las hembras para producir leche durante la lactancia. En la presente investigación se utilizaron los registros productivos de 197 bovinos colombianos sintéticos de la laza Lucerna y sus cruces, procedentes de una empresa agropecuaria ubicada en el Valle del Cauca, Colombia. La variable producción de leche por día fue modelada en función del componente racial, la edad de la vaca al momento del parto y la condición corporal. Para la construcción de los modelos fue utilizada una estrategia incremental, iniciando con un modelo lineal clásico con efectos fijos, agregando posteriormente efectos aleatorios y terminando con la modelación de la varianza y las correlaciones entre las observaciones para cada animal. Para identificar el modelo con el mejor ajuste, se utilizaron la prueba de razón de verosimilitud y el criterio de información de Akaike. El modelo que mejor se ajustó para explicar la producción de leche por día fue el que incluyó la variable condición corporal, la edad de la vaca al momento del parto y un término cuadrático de la edad. El modelo estimado acompañó bien la evolución de la producción de leche por día a lo largo de las lactancias, adicionalmente, el coeficiente de correlación de Pearson entre el valor real y el valor estimado fue de 0.792. Los supuestos del modelo se cumplieron satisfactoriamente. Las mayores producciones de leche fueron observadas para vacas con condición corporal alta y edades al parto entre los 100 y 150 meses

Colombia is considered a very diverse country with multiple productive options at the livestock level. The Colombian race Lucerna is characterized by good productive and reproductive performance in the low tropic of Colombia, due to its good hardiness. There are few studies that allow to identify the productive advantages that can be obtained, when using this type of animals. The study of dairy milk yield is an important tool to evaluate the efficiency of a milk production system; this variable is affected by genetic and non-genetic factors. On the other hand, the evaluation of Body Condition, has been proposed as a tool in the nutritional, reproductive and health management of bovine herds, as it is a good indicator of the nutritional status of the cow and therefore of the performance level that females can have to produce milk during lactation. The productive records of 197 synthetic Colombian cattle of the Lucerna breed and their crosses were used, coming from a farm located in Valle del Cauca, Colombia. The variable dairy milk yield per parturition was modeled according to the racial component, age of the cow at the time of delivery, body condition and the interaction between these last two variables. For the construction of the models used in the description of the evaluated characteristic, an incremental strategy was used, starting with a classical linear model with fixed effects, adding random effects later and ending with the correlation estimate between the correlations for each animal. To identify the model with the best fit, the likelihood ratio test and the Akaike Information Criterion were used. The final model to the explain the dairy milk yield variation was the one that included the body condition, the age of the cow at the time of delivery, the interaction between these two variables and a random intercept. The estimated model accompanied well the evolution of milk production per day throughout the eight lactations. The Pearson correlation coefficient between the real value and the estimated value was 0.792. The assumptions of the model were satisfactorily fulfilled. Animals with low body condition had lower daily milk production with increasing age at the time of delivery. The highest milk yields were observed for cows with high body condition and calving ages between 100 and 150 months