Software-automatized individual lactation model fitting, peak and persistence and Bayesian criteria comparison for milk yield genetic studies in Murciano-Granadina goats

SPSS model syntax was defined and used to evaluate the individual performance of 49 linear and non-linear models to fit the lactation curve of 159 Murciano-Granadina goats selected for genotyping analyses. Lactation curve shape, peak and persistence were evaluated for each model using 3107 milk yield controls with an average of 3.78 ± 2.05 lactations per goat. Best fit (Adjusted R2) values (0.47) were reached by the five-parameter logarithmic model of Ali and Schaeffer. Three main possibilities were detected: non-fitting (did not converge), standard (Adjusted R2 over 75%) and atypical curves (Adjusted R2 below 75%). All the goats fitted for 38 models. The ability to fit different possible functional forms for each goat, which progressively increased with the number of parameters comprised in each model, translated into a higher sensitivity to explaining curve shape individual variability. However, for models for which all goats fitted, only moderate increases in explanatory and predictive potential (AIC, AICc or BIC) were found. The Ali and Schaeffer model reported the best fitting results to study the genetic variability behind goat milk yield and perhaps enhance the evaluation of curve parameters as trustable future selection criteria to face the future challenges offered by the goat dairy industry.