Caracterización de la curva de lactancia de bovinos Siboney con modelos no lineales mixtos

Abstract

En bovinos productores de leche la curva de lactancia (CL) se ha analizado con modelos no lineales (MNL), y recientemente se ha implementado la metodología de MNL mixtos. Los objetivos fueron: a) ajustar un MNL a la CL a 240 días; b) comparar MNL vs MNL mixtos; c) estimar producción al pico (PMAX) y total a 240 días (PTOTAL), así como días para alcanzar la producción máxima (DPMAX). Se analizaron 15,324 observaciones de producción diaria (kg), correspondientes a primera lactancia de 2,809 vacas Siboney, con fecha de parto de 2000 a 2012 en 28 hatos. Se evaluaron cinco MNL: Wood, Wiltmink, Cobby, Brody y Sikka. Se realizaron dos análisis con el procedimiento NLMIXED de SAS: ANA1, incluyó sólo el efecto aleatorio de los residuales; y ANA2, incluyó el efecto aleatorio correspondiente a un coeficiente de regresión, más los residuales. La selección del modelo con mejor ajuste se realizó con: error de predicción promedio; varianza del error de predicción; estadístico Durbin Watson; coeficiente de determinación (R2); criterios de información Akaike (AIC) y Bayesiano (BIC). Los MNL mixtos mostraron el mejor ajuste con base en: a) los R2 del ANA2 fueron superiores, en un intervalo de 9 a 12 %; b) por la reducción en la varianza residual en el ANA2; y, c) los criterios AIC y BIC, mostraron mejor ajuste en los MNL del ANA2. El modelo de mejor ajuste en todo el estudio fue WOD del ANA2, con estimaciones para PMAX, DPMAX y PTOTAL de 7.71 kg, 40 días y 1,653 kg, respectivamente.Nonlinear models (NLM) have been used to evaluate dairy cattle lactation curves (LC), and nonlinear mixed models are increasingly applied. The most appropriate model for estimating LC components in Siboney cattle was identified by applying a NLM for LC to 240 d; comparing a non-mixed NML vs a mixed NML; and estimating peak production (PMAX), total production at 240 d (PTOTAL) and days to maximum production (DPMAX). The database consisted of 15,324 records for daily milk production (kg) corresponding to first lactation in 2,809 Siboney cows (calving date in 2000-2012) in 28 herds in Cuba. Five NLM were evaluated: Wood, Wiltmink, Cobby, Brody and Sikka. Two SAS NLMIXED analyses were run: ANA1, including only the random effect of residuals; and ANA2, including a regression coefficient random effect plus residuals. Selection of the best fitting model was done based on average prediction error; prediction error variance; the Durbin Watson statistic; the determination coefficient (R2); the Akaike information criterion (AIC); and the Bayesian information criterion (BIC). The mixed NLM (ANA2) had the best fit: a) its R2 was higher (9 to 12 %); b) residual variance was notably reduced; and c) the AIC and BIC criteria exhibited a better fit. The model with the overall best fit was the ANA2 WOD, with estimates of 7.71 kg for PMAX, 40 d for DPMAX and 1,653 kg for PTOTAL.