Foram utilizados arquivos com 104.101, 55.063 e 60.782 registros de pesos corporais, do nascimento aos 630 de idades, de bovinos da raça Guzerá, da Associação Brasileira dos Criadores de Zebu (ABCZ), referentes ao período de 1975 a 2001, para estimar componentes de (co)variância e parâmetros genéticos, de acordo com três abordagens de análise, cada uma com um arquivo, utilizando diferentes modelos estatísticos e a metodologia REML. Nestas abordagens foram obtidos valores das estimativas de variância e parâmetros genéticos de modelos unicaracterísticos; modelos uni e bicaracterísticos e, de modelos de regressão aleatória (MRA), respectivamente. Na primeira, com quatro modelos unicaraterísticos, observou-se que o modelo 1, considerado completo, por incluir os efeitos genéticos direto e materno (GM) e os efeitos de ambiente permanente materno e residual, não diferiu significativamente (P<0,05), pelo teste razão de verossimilhança do modelo 2, que exclui o GM. As herdabilidades diretas (h2) estimadas nos modelos 1 e 2, foram muito semelhantes. Os valores de h2 cresceram da primeira idade até a segunda idade, mantiveram os valores até o desmame e depois cresceram. As estimativas de variância genética materna foram baixas, principalmente antes da desmama. Na segunda usando diferentes modelos uni e bicaracterísticos, o comportamento das estimativas foi semelhante à primeira abordagem. Os valores de h2 das análises uni e bicaracterísticas, para os pesos ajustados às idades padrão de 120(P120), 205(P205), 365(P365) e 550(P550) dias de idades foram 0,15; 0,10; 0,17; 0,14 e 0,14; 0,10; 0,16; 0,15, respectivamente. As correlações genéticas diretas (ra) semelhantes nos modelos 1 e 2 para P120/P205, P120/P365, P120P550, P205/P365, P205/P550 e P365/P550 foram 0,80; 0,54; 0,54; 0,74; 0,62 e 0,95 e, 0,80; 0,54; 0,53; 0,74; 0,62 e 0,95 respectivamente. Na comparação entre os modelos 1 e 2, das análises uni e bicaracterísticas, o modelo 2, que exclui GM, proporcionou ajuste similar ao modelo 1. Na análise de regressão aleatória foram testados dez modelos com diferentes estruturas de variâncias residuais e combinações de ordens (k) dos efeitos aleatórios. Pelos critérios de informação de Akaike e Bayesiano de Schwarz, os MRA mais adequados foram os modelos Reg666-r10 (RegkAkCkQ-r10) e Reg653-r10 (RegkAkCkQ-r10), respectivamente, para os efeitos genético aditivo direto (kA), de ambiente permanente de animal (kC) e de ambiente permanente materno (kQ), com as variâncias residuais divididas em dez (r10) classes. As estimativas de h2 de peso ao nascer (PN), 205, 365 e 550 dias de idades, foram de 0,13, 0,43, 0,46 e 0,48, para o modelo Reg666-r10; para o modelo mais parcimonioso (Reg653-10), foram de 0,13, 0,46, 0,54 e 0,56, respectivamente. As correlações fenotípicas, genéticas, de ambientes permanentes animal e materno foram todas positivas e similares entre os MRA. No modelo mais parcimonioso, Reg653-10, as ra do PN/205, PN/365, PN/550, 205/365, 205/550 e 365/550 foram, 0,40, 0,40, 0,42, 0,68, 0,74 e 0,81, respectivamente. Os modelos com homogeneidade de variância foram inadequados. A variância do resíduo dividido em dez classes distintas de variâncias foi mais adequada para modelar a variação residual. Na comparação das três abordagens de análise foi constatado que o MRA constitui alternativa importante para análise de dados de crescimento de bovinos.

Three different files containing 104,101, 55,063 and 60,782 body weight records from birth to 630 days of age, collected from 1975 to 2001, from Guzera cattle belonging to the Brazilian Association of Zebu Breeders (ABCZ) were used in this study. The data were analyzed according to three different approaches, corresponding to one of those files, in order to estimate the (co)variance components and the genetic parameters by different statistical models and the REML methodology: variance estimates and genetic parameters of single-trait models, single and two-trait models, and of random regression models (RRM). Regarding the first analysis, where four single-trait models were used, it was observed that model 1, the complete one which included the direct genetic and maternal (GM) effects, as well as the maternal permanent environmental effect and residual, did not statistically (P<0.05) differ by the likelihood ratio test from model 2, which excluded the GM. Direct heritabilities (h2) estimated by models 1 and 2 were quite similar, and increased from first to second age, remained until weaning, and then increased. Estimates of maternal genetic variance were low, especially before weaning. The second analysis, which used different single and two-trait models, resulted in similar behavior of estimates to the first analysis. The h2 estimated by single and two-trait analyses, for weight records adjusted to 120(P120), 205(P205), 365(P365) and 550(P550) days of age were 0.15, 0.10, 0.17, 0.14 and, 0.14, 0.10, 0.16, 0.15, respectively. The direct genetic correlations (ra) by models 1 and 2 for P120/P205, P120/P365, P120/P550, P205/P365, P205/P550 and P365/P550 were 0.80, 0.54, 0.54, 0.74, 0.62, 0.95 and, 0.80, 0.54, 0.53, 0.74, 0.62, 0.95, respectively. When comparing models 1 and 2, by the single and two-trait analyses, it was verified that model 2, which excluded the GM, was equivalent to the complete model. In the random regression analysis ten models presenting different residual variance structures and order combinations (k) of random effects were evaluated. According to the Akaike Information Criterion and the Schwarz Bayesian Information Criterion, the most suitable RRM were the models Reg666-r10 (RegkAkCkQ-r10) and Reg653-r10 (RegkAkCkQ-r10), respectively, for the direct additive effect (kA), the animal permanent environmental effect (kC) and the maternal permanent environmental effect (kQ), with the residue variances divided into ten distinct variance classes (r10). The h2 estimates for weights at birth (PN), 205, 365 and 550 days of age were 0.13, 0.43, 0.46, 0.48 and 0.13, 0.46, 0.54, 0.56 for the Reg666-r10 model and for the most parsimonious model (Reg653-10), respectively. The phenotypic, genetic, animal permanent environmental and maternal permanent environmental correlations were all positive and similar for all RRM. Considering the most parsimonious model, Reg653-10, the ra between PN/205, PN/365, PN/550, 205/365, 205/550 and 365/550 were 0.40, 0.40, 0.42, 0.68, 0.74 and 0.81, respectively. Models with variance homogeneity were inadequate. The residual variance divided into ten distinct variance classes was the most suitable to model the residual variation. In the comparison of the three approaches of analysis it was evidenced that the RRM constitutes important alternative for analysis of growth data of bovines.