Analysis of culling probability in dairy buffalo using survival models

In order to contribute to the genetic breeding programs of buffaloes, this study aimed to determine the influence of environmental effects on the stayability (ST) of dairy female Murrah buffalo in the herd. Data from 1016 buffaloes were used. ST was defined as the ability of the female to remain in the herd for 1, 2, 3, 4, 5 or 6 years after the first calving. Environmental effects were studied by survival analysis, adjusted to the fixed effects of farm, year and season of birth, class of first-lactation milk yield and age at first calving. The data were analyzed using the LIFEREG procedure of the SAS program that fits parametric models to failure time data (culling or ST = 0), and estimates parameters by maximum likelihood estimation. Breeding farm, year of birth and first-lactation milk yield significantly influenced (P < 0.0001) the ST to the specific ages (1 to 6 years after the first calving). Buffaloes that were older at first calving presented higher probabilities of being culled 1 year after the first calving, without any effect on culling at older ages. Buffaloes with a higher milk yield at first calving presented a lower culling probability and remained for a longer period of time in the herd. The effects of breeding farm, year of birth and first-lactation milk yield should be included in models used for the analysis of ST in buffaloes.