Search

We used a class of stochastic differential equations (SDE) to model the evolution of cattle weight that, by an appropriate transformation of the weight, resulted in a variant of the Ornstein-Uhlenbeck model. In previous works, we have dealt with estimation, prediction, and optimization issues for this class of models. However, to incorporate individual characteristics of the animals, the average transformed size at maturity parameter α and/or the growth parameter β may vary randomly from animal to animal, which results in SDE mixed models. Obtaining a closed-form expression for the likelihood function to apply the maximum likelihood estimation method is a difficult, sometimes impossible, task. We compared the known Laplace approximation method with the delta method to approximate the integrals involved in the likelihood function. These approaches were adapted to allow the estimation of the parameters even when the requirement of most existing methods, namely having the same age vector of observations for all trajectories, fails, as it did in our real data example. Simulation studies were also performed to assess the performance of these approximation methods. The results show that the approximation methods under study are a very good alternative for the estimation of SDE mixed models. [ABSTRACT FROM AUTHOR]

In early work we have studied a class of stochastic differential equation (SDE) models, for which the Gompertz and the Bertalanffy-Richards stochastic models are particular cases, to describe individual growth in random environments, and applied it to model cattle weight evolution using real data. We have started to work on these type of models considering that the model parameters are fixed, i.e. the same for all the animals. Aiming to incorporate variability among individuals, we consider that the model parameters can be random variables, resulting in SDE mixed models. In additon, here we consider SDE mixed models, allowing the parameters to be random and propose to incorporate each animal's genetic characteristics considering the transformed animal's weight at maturity to be a function of its genetic values. The main objective is to extend the SDE mixed model to the more realistic case where the individual's genetic value becomes an important component in the estimated growth curve. For the estimation of the model parameters we have used maximum likelihood estimation theory. Estimates and asymptotic confidence intervals of the parameters are presented. A comparison with SDE non-mixed model and SDE mixed model without the inclusion of genetic characteristics is shown with the conclusion that the incorporation of some genetic characteristics in the model parameters improves the estimation of the animal's growth curve. [ABSTRACT FROM AUTHOR]

Stochastic differential equations (SDE) appropriately describe a variety of phenomena occurring in random environments, such as the growth dynamics of individual animals. Using appropriate weight transformations and a variant of the Ornstein–Uhlenbeck model, one obtains a general model for the evolution of cattle weight. The model parameters are α , the average transformed weight at maturity, β , a growth parameter, and σ , a measure of environmental fluctuations intensity. We briefly review our previous work on estimation and prediction issues for this model and some generalizations, considering fixed parameters. In order to incorporate individual characteristics of the animals, we now consider that the parameters α and β are Gaussian random variables varying from animal to animal, which results in SDE mixed models. We estimate parameters by maximum likelihood, but, since a closed-form expression for the likelihood function is usually not possible, we approximate it using our proposed delta approximation method. Using simulated data, we estimate the model parameters and compare them with existing methodologies, showing that the proposed method is a good alternative. It also overcomes the existing methodologies requirement of having all animals weighed at the same ages; thus, we apply it to real data, where such a requirement fails. [ABSTRACT FROM AUTHOR]