Published in

Cambridge University Press, Journal of Agricultural Science, 5(152), p. 829-842, 2013

DOI: 10.1017/s0021859613000798

Links

Tools

Export citation

Search in Google Scholar

Comparison of growth curves of two genotypes of dairy goats using nonlinear mixed models

This paper is available in a repository.
This paper is available in a repository.

Full text: Download

Green circle
Preprint: archiving allowed
Red circle
Postprint: archiving forbidden
Red circle
Published version: archiving forbidden
Data provided by SHERPA/RoMEO

Abstract

SUMMARYThe objective of the current study was to assess the use of nonlinear mixed model methodology to fit the growth curves (weightv.time) of two dairy goat genotypes (Alpine, +A and Saanen, +S). The nonlinear functions evaluated included Brody, Von Bertalanffy, Richards, Logistic and Gompertz. The growth curve adjustment was performed using two steps. First, random effectsu1,u2andu3were linked to the asymptotic body weight (β1), constant of integration (β2) and rate constant of growth (β3) parameters, respectively. In addition to a traditional fixed-effects model, four combinations of models were evaluated using random variables: all parameters associated with random effects (u1,u2andu3), onlyβ1andβ2(u1andu2), onlyβ1andβ3(u1andu3) and onlyβ1(u1). Second, the fit of the best adjusted model was refined by using the power variance and modelling the error structure. Residual variance ($σ _e^2 $) and the Akaike information criterion were used to evaluate the models. After the best fitting model was chosen, the genotype curve parameters were compared. The residual variance was reduced in all scenarios for which random effects were considered. The Richards (u1andu3) function had the best fit to the data. This model was reparameterized using two isotropic error structures for unequally spaced data, and the structure known in the literature as SP(MATERN) proved to be a better fit. The growth curve parameters differed between the two genotypes, with the exception of the constant that determines the proportion of the final size at which the inflection point occurs (β4). The nonlinear mixed model methodology is an efficient tool for evaluating growth curve features, and it is advisable to assign biologically significant parameters with random effects. Moreover, evaluating error structure modelling is recommended to account for possible correlated errors that may be present even when using random effects. Different Richard growth curve parameters should be used for the predominantly Alpine and Saanen genotypes because there are differences in their growth patterns.