Pre-weaning performance evaluation of a multibreed Aberdeen Angus × Nellore population using different genetic models

Avaliação do desempenho na pré desmama de uma população bovina multirracial Aberdeen Angus × Nelore utilizando diferentes modelos genéticos

BREEDING, GENETIC AND REPRODUCTION

Pre-weaning performance evaluation of a multibreed Aberdeen Angus × Nellore population using different genetic models

Avaliação do desempenho na pré desmama de uma população bovina multirracial Aberdeen Angus × Nelore utilizando diferentes modelos genéticos

Jader Silva Lopes^I,III; Paulo Roberto Nogara Rorato^II; Tomás Weber^III; Ronyere Olegário de Araújo^III; Mariana de Almeida Dornelles^IV; Juliana Grigolleto Comin^V

^IDepartamento de Zootecnia - Universidade Federal de Santa Maria - Santa Maria, RS, Brasil, CEP: 97105-970

^IIDepartamento de Zootecnia da UFSM

^IIICurso de Pós-Graduação em Zootecnia da UFSM. Bolsista CAPES

^IVCurso de Pós-Graduação em Zootecnia da UFSM. Bolsista CNPq

^VCurso de Graduação em Zootecnia da UFSM. Bolsista PIBIC/CNPq

ABSTRACT

This work aimed at estimating the genetic effects that affect the pre-weaning performance of animals from multibreed crosses. In order to do so, it was used information of the weight at weaning of 79,521 animals, sired by 1,020 bulls and 61,898 cows from Aberdeen Angus and Nellore breeds and from many genetic groups resulted from their crosses. Five genetic models were tested: model 1, containing the fixed breed genetic effects (additive direct and maternal effects, heterozygote direct and maternal effects, epystatic direct and maternal effects, joint additive direct and maternal effects); model 2, equal to model 1, excluding direct and maternal joint additive effects; model 3, equal to model 1, excluding direct and maternal epystatic effects; model 4, equal to model 1, excluding direct and maternal epystactic effect and direct and maternal joint additive effects; and model 5, equal to model 1, excluding direct and maternal heterozygotic effects, direct and maternal epystatic effects and direct and maternal joint additive effects. The models were analyzed by the following methods: least square means method, ridge regression method, and the restricted maximum likelihood method. The dominant additive models usually used for genetic evaluations do not give a good description of the pre-weaning performance variations, making it necessary to add the heterozygote and epystatic effects; the joint additive effects do not significantly improve the adjustment of the analysis model and the heterozygote effects are efficient in representing a quadratic breed additive effect, in addition to insert an unnecessary bias assigned to multicollinearity related to the joint additive effects.

Key Words: additive effect, complementarity among breeds, epistasis, heterozygous, multicollinearity, weaning weight

RESUMO

Com o objetivo de estimar os efeitos genéticos que influenciam o desempenho pré-desmama de animais produtos de cruzamentos multirraciais, foram utilizadas informações do peso a desmama de 79.521 bezerros, filhos de 1.020 touros e 61.898 vacas das raças Aberdeen Angus e Nelore e de diversos grupos genéticos resultantes do cruzamento entre elas. Foram testados cinco modelos genéticos: o modelo 1, contendo os efeitos genéticos de raça fixos (aditivos diretos e maternos, heterozigóticos diretos e maternos, epistáticos diretos e maternos e aditivos conjuntos diretos e maternos); o modelo 2, igual ao modelo 1, porém sem os efeitos aditivos conjuntos diretos e maternos; o modelo 3, igual ao modelo 1, porém sem os efeitos epistáticos diretos e maternos; o modelo 4, igual ao modelo 1, sem os efeitos epistáticos diretos e maternos e os aditivos conjuntos diretos e maternos; e o modelo 5, igual ao modelo 1, sem os efeitos heterozigóticos diretos e maternos, espistáticos diretos e maternos e aditivos conjuntos diretos e maternos. Os modelos foram analisados pelos métodos dos quadrados mínimos, de regressão de cumeeira e de máxima verossimilhança restrita. Os modelos aditivo-dominantes, usualmente utilizados em avaliações genéticas, não descrevem adequadamente as variações no desempenho pré-desmama, devendo ser adicionados os efeitos heterozigóticos e epistáticos; os efeitos aditivos conjuntos não melhoram significativamente os ajustes dos modelos de análise e os efeitos heterozigóticos são eficientes para representar um efeito quadrático aditivo de raça, além de inserirem um viés desnecessário atribuído a multicolineariedade relacionada aos efeitos aditivos conjuntos.

Palavras-chave: complementaridade entre raças, efeito aditivo, epistasia, heterozigose, multicolineariedade, peso à desmama

Introduction

In Brazil, as in other countries, the use of multibreed populations for commercial purposes has been growing, given the possibility to increase the efficiency of production determined by heterosis and complementarity between breeds used in crossing, enhanced by the selection of superior animals.

The response to selection is proportional to the accuracy of prediction of breeding values (Falconer & Mackay, 1996). On the other hand, the accurate estimation of genetic value of individuals subjected to a breeding program depends, in large part, on the effects considered in the statistical model used for evaluation of the animals.

According to Cardoso & Tempelman (2004), a model that includes additive and non-additive genetic fixed effects, in addition to additive random genetic effects, may be suitable for genetic evaluation of multiracial populations. However, Kinghorn (1993), Fries (1996), Sharma et al. (2000) and Pimentel et al. (2006) emphasized the importance of effects such as epistasis, heterozygosity and complementarity between breeds, in addition to additives.

The most used method to derive prediction equations is the least-squares method. However, when there are strong linear relationships between independent variables (multicollinearity), the simple estimation of individual regression coefficients by least squares tends to be unstable, often with high standard error, which can lead to erroneous inferences (Bergmann & Hohenboken, 1995).

The ridge regression method (Hoerl & Kennard, 1970) is one of the alternative methods of estimation that provide a more informative analysis of the data in the presence of multicollinearity.

The purpose of this study was to evaluate the importance of including non-additive genetic effects in the analysis model of records of pre-weaning performance of a multiracial population of cattle, derived from crosses between the Aberdeen Angus and Nelore breeds.

Material and Methods

The original database contained performance information of 121,241 calves sired by 1,359 bulls and 84,465 cows of Aberdeen Angus and Nellore breeds by crossing between them different genetic group, collected from 1986 to 2002 from 75 farms located in the states of Goiás, Minas Gerais, Mato Grosso do Sul, Paraná, Rio Grande do Sul, São Paulo and Tocantins.

The data were structured using the SAS application (2001) to allow the file to be used for estimation of genetic and environmental effects. In the data edition, contemporary group variables were created, considering animals born on the same farm (1-75) in the same year (1986-2002) and in the same season (1-4), belonging to the same gender and to the same management group until weaning (1-388), totalizing 2,439 contemporary groups; age at weaning in classes (IDC), corresponding to the weaning age in months (5-11); and age of the cow at calving in classes (IVC), indicating the age of the cow in years (3-15). It was eliminated data from product of cows whose type of reproduction did not result from artificial insemination or controlled breeding from bulls with less than 10 offspring, contemporary groups with less than five components, weaning weight (PD) outside the range of ± 2.5 standard deviations compared to the average population, genetic groups with fewer than 100 animals or observed in only one farm, farms with only one single genetic group, and cows at less than three years of age and at more than fifteen years of age when delivery occurred. After editing, the working file was composed of the weaning weights of 79,521 calves sired by 1,020 bulls and 61,898 cows, totalizing 135,051 animals in the relationship matrix.

The coefficients for direct additive (AD) and maternal additive of breed (AM) effects were defined by the contribution of genes of the Nellore breed in the genetic composition of each individual. The coefficients for joint additive direct (ACD) and maternal (ACM) effects between the Aberdeen Angus and Nellore breeds were calculated as ACD=AD*(1-AD) and ACM=AM*(1-AM), respectively, theorizing a quadratic additive effect.

The direct and maternal heterozygosity were calculated as: h_ij = α α + α α, in which α and α denote the proportion of the 'i' gene race in the mother and father of the animal, respectively.

The effects of direct and maternal epistasis were calculated as the average heterozygosity present in the gametes that generate each individual or as the average heterozygosity in the parents of an individual (Fries et al., 2000).

For heterozygous and maternal epistatic effects, genetic information regarding the maternal grandparents of the individuals was necessary. When this information was not available, it was considered that the mothers were produced by inter se mating.

Five models were tested using three different methods of analysis: least squares method (MQM) and ridge regression with λ = 0.05 (RC) using SAS statistical package (2001), and Restricted Maximum Likelihood Method (REML) using the Multiple Trait Derivative Free Restricted Maximum Likelihood (MTDFREML) application described by Boldman et al. (1995).

For the least squares and ridge regression methods, the statistical model (model 1) can be described as: Y_ijkl = µ+ GC_i + IDC_j + IVC_k + β₁AD + β₂AM + β₃HD + β₄HM + β₅ED + β₆EM + β₇ACD + β₈ACM + ε_ijkl, in which Y_ijkl = observation of weight to weaning of the l^th animal; µ = overall mean of the trait; GC_i = effect of the i^th contemporary group; IDC_j = effect of j^th age class at weaning; IVC_k = effect of k^th age class of the cow at calving; β₁, β₂, β₃, β₄, β₅, β₆, β₇ and β₈ = regression coefficients of direct and maternal additive genetic, direct and maternal heterozygous, direct and maternal epistatic and joint additive direct and maternal effects, respectively, and ε_ijkl = random error associated with each observation, NID (0, δ²).

The ridge regression method (Hoerl, 1962) consists on the addition of a scalar (λ) to the elements of the main diagonal of the coefficient matrix, in order to break the linear dependence found between its columns. The estimator becomes: = (X'X+λI)^-1X'y

Originally, Hoerl & Kennard (1970) proposed a graphical method, namely ridge trace, to choose the λ value. By this method, the parameter estimates are plotted against various values of λ, starting from zero. The value of λ is then chosen, from which it is possible to see stabilization in the values of the estimates. However, high values of λ lead to the decrease of the regression coefficients and of standard errors of these estimates, which could affect the estimation of the model effects, especially those which are not involved in cases of multicollinearity.

The procedure implemented using the SAS statistical package, as described by Freund & Littell (2000), was the following:

_SAS= [R_xx (I+D_λ)]^-1R_Xy, in which R_XX= correlation matrix of the explanatory variables; R_Xy = correlation vector of the explanatory variables with the response variable; and D_λ = diagonal matrix with λ values on the diagonal.

The parameters and breeding values were obtained by REML with convergence criterion of 1x10^-9, under the model described in the matrix form: y = Xβ + Wg + Z₁a + Z₂m + e, in which y = vector of observations of the weights to weaning; X = incidence matrix associated with fixed environmental effects (GC, IDC and IDV), β = vector of solutions for fixed environmental effects; W = incidence matrix for fixed genetic effects (direct and maternal additive, direct and maternal heterozygous, direct and maternal epistatic and direct and maternal joint additive); g = vector of solutions for fixed genetic effects; Z₁ = incidence matrix associated with random direct genetic effect of each animal; a = vector of solutions to the random direct additive genetic effects; Z₂ = incidence matrix associated with random maternal additive genetic effect of each animal; m = vector of solutions for random maternal additive genetic effects; and e = vector of random waste.

Briefly, the models were the following: model 1 = as described above; model 2 = model 1 without the direct and maternal joint additive effects; model 3 = model 1 without the direct and maternal epistatic effects; model 4 = model 1 without the direct and maternal epistatic, direct and maternal joint additive effects; and model 5 = model 1 without the direct and maternal heterozygous, direct and maternal epistatic and joint additive direct and maternal effects.

The efficiency of the methodologies and the comparison among models were performed by analysis of the coefficients of determination (R²), the F test, the C(p) statistic of Mallows, the variance inflation factors and the magnitude of correlations (Pearson correlation) and ranking (Spearman correlation), in addition to the likelihood ratio test.

To test whether the inclusion of parameters for fixed genetic effects effectively improved the efficiency of the model in terms of adjustment, the F test was undertaken, providing evidence against the simplest model (x-1) if the F value is greater than the tabulated values for a distribution as F(gl_x-1 - gl_x, gl_x, α).

The F value was calculated using the following formula, as described by Weisberg (1980). , in which SQR_x and SQR_x-1 are the sums of squared residuals of models x and x-1, respectively; and gl_x and gl_x-1 are the degrees of freedom of model residues x and x-1, respectively, in which x represents models 2, 3, 4 or 5.

As described by Freund & Littell (2000), the C(p) statistic of Mallows (1973), a measure of error variance plus the bias introduced by the exclusion of a variable of type, is calculated as it follows: C (p) = [SQR (p) / MSE] - (N - 2p) + 1, in which MSE is the mean square error of the full model; SQR (p), the sum of squares of a model containing a subset of p variables; and N, the number of observations. When it is observed that C (p) > (p +1) for a model containing explanatory p variables, there is evidence of bias by exclusion of an important variable in the model.

The identification of the problem of multicollinearity in regression analysis is usually done by examining the correlations between pairs of explanatory variables. However, in some cases, associations between three or more variables cannot be detected by examining correlations in pairs. A more efficient way of diagnosing multicollinearity is the examination of variance inflation factors (FIV), defined by: FIV (xi) = , in which is the squared multiple correlation coefficient that results from the regression of variable x_i against all other explanatory variables included in the model.

The likelihood ratio test was used to verify if the models differ statistically by the difference between the values of -2log L. This test is based on a chi-square distribution with g degrees of freedom and error probability of 5%, in which g is the difference in numbers of estimated parameters in the compared models (Dobson, 1990).

Results and Discussion

In this experiment, the animals were in average at 217.63 days of age, with an average weight of 177.69 kg, an their dams calved them at an average weight of 5.59 years. About 80% of the animals were weaned from six to eight months of age, and approximately 25% of mothers were four years of age at birth, indicating that this was a fairly young flock.

The working file was first examined by the five models submitted to the method of least squares. It can be checked the significance of all covariate fixed genetic effects on weaning weight (P<0.01), except for the direct additive effect on the model 1 and of joint additive direct effect in model 3 (Table 1). The models explained about 60% of the variation in weaning weight, whereby model 5, which contained only additive fixed genetic effects, showed the lowest setting, with a loss of about seven percentage points in the coefficient of determination (R²) confirmed by the high estimate for the C(p) statistic.

However, the wide variation in the regressors is a signal for possible collinearity among the independent variables, which was confirmed by analysis of correlations between these variables, whose correlation coefficients ranged from -0.12 (between direct additive effects and maternal joint additive) to 0.98 (between epistatic and maternal heterozygous) with an average of 0.45.

Due to multicollinearity, it was employed the ridge regression method, indicated by Hoerl & Kennard (1970) as one of alternative methods of estimation that provide a more informative analysis of the data in the presence of multicollinearity (Table 2).

From analysis of the ridge trace cart (Hoerl & Kennard, 1970) we chose the value of λ as 0.05, a minimal value that provided the breakdown of multicollinearity, reducing the variance inflation factors, and not dramatically, the regressors and their standard errors.

Analyzing the comparative values of F, it is observed a loss of adjustment when a couple of effects is removed. Only the comparison value between models 1 and 2 was not significant, which indicates that withdrawal of joint additive genetic effects.

The variance inflation factor is due to increased variance of the regressors, as determined by the correlation between independent variables. When the variance inflation factor is greater than 10, it indicates the occurrence of multicollinearity (Freund & Wilson, 1998). According to Dias et al. (2003), this procedure is performed because, in the presence of multicollinearity, statistical tests may fail to detect significant differences among the factors. Thus, it may be noted that the effects with greater involvement in cases of multicollinearity are the heterozygous and joint additives, and that the application of λ = 0.05 by the ridge regression method was enough to substantially reduce the values of variance inflation factors (Table 3).

In an attempt to verify the possible effects of multicollinearity and distinctions imposed by the use of different models in multiracial genetic evaluation, we estimated genetic and phenotypic covariances used later for prediction of breeding values at the studied population.

Direct heritabilities and their standard errors estimates ranged from 0.29 ± 0.020 to 0.35 ± 0.022 (Table 4), showing the possibility of genetic gain in this population through the use of mass selection for weaning weight. Lower values were reported by Everling et al. (2001) for products from the cross between the Angus and Nellore breeds (0.23), by Ferraz Filho et al. (2002) for the Tabapuã breed (0.23) and by Weber (2008) for the Aberdeen Angus breed (0.24). However, heritabilities higher than those found in this study were described by Kaps et al. (2000), as 0.53, and by Souza et al. (2007) as 0.56 for animals of the Aberdeen Angus breed and by Lopes et al. (2008) for the Nellore breed in the three Southern States of Brazil (from 0.41 to 0.47).

Maternal heritability values indicated less variability in the maternal additive random effect rather than in the direct one, with estimates ranging from 0.15 ± 0.015 to 0.26 ± 0.019. These figures are evidence of the importance of maternal effect on the development of the calves at pre-weaning phase.

On the other hand, the values of genetic correlation between direct and maternal additive genetic effects (r_am) show antagonism between these two groups of genes. Negative values for r_am are not uncommon in the available literature; authors like Marques et al. (2000), Everling et al. (2001), Corrêa et al. (2006) and Souza et al. (2007) reported values ranging from -0.83 to -0.17.

According to Crump et al. (1997), problems with maternal effects and their correlations may be related to the structure of the analyzed data, depending on the number of generations covered by the available kinship structure. In this study, only 12% of female born products took her offspring with accompanied performance, which may add maternal bias in the estimates. This is because the genes responsible for this effect are expressed only in females, and therefore the assessment of their importance may depend on the existence, in the data set, of a considerable number of female relatives with progeny evaluated.

The values of the likelihood ratio test indicate the existence of significant differences between the tested models (models 4 and 5; models 3 and 4; models 2 and 3; and models 1 and 2). However, comparisons between models 1 and 2 and between models 3 and 4 showed minor differences with lower values in the likelihood ratio test. The difference between each pair of models are the joint additive effects, and these minor differences show the importance of this small effect of variation in weaning weight in this population.

It is observed a reduction in the residual variance in the application of model 3, however, the high values of variance inflation factor for the heterozygous and joint additive effects present in the model, combined with the large variation in estimates of fixed genetic effects (Table 5) among the different models used, confirm the existence of multicollinearity, apart from the problems of estimation by REML in the presence of this type of interaction among estimators.

The correlations of magnitude (Pearson's) and of rank (Spearman) of breeding values, predicted by different models, ranged from 0.911 to 0.989 and from 0.819 to 0.978, respectively. The average accuracy of breeding values estimated by model 5 was 26% lower than those estimated by model 1 (0.55 and 0.74, respectively). The major similarities between the genetic values were observed between models 1 and 2, confirming the results obtained by estimates of C(p) statistic, F test and likelihood ratio test.

The results of this study show that the inclusion of heterozygous effects in the analysis model (model 4) promotes an increase in the coefficient of determination in relation to the additive model (model 5). Very small changes occurred when joint additive effects were added, reinforced by the negative fact that these effects are strongly involved in cases of multicollinearity.

The inclusion of epistatic effects in model 2 promoted substantial gain in the adjustment. These effects affect negatively the performance of pre-weaning calves in the studied population, in accordance with Trematore et al. (1998), Fries et al. (2000) and Pimentel et al. (2006). This model was very similar to the more complete model (model 1), with the advantage of not containing multicollinearities that may add bias to the estimates.

Deleterious effects were observed for epistasis, according to Arthur et al. (1999), Fries et al. (2000), Roso et al. (2005) and Pimentel et al. (2006), whereby the first authors, studying a Hereford-Brahman population, were those who reported the most significant effect, where they noticed a difference of -61.5 ± 19 kg in weight at weaning adjusted to 240 days of age, in favor of animals with minimal epistasis. Fries et al. (2000) explain that, when breeds cross, genes are forced to interact and to cooperate with other genes with which they are not accustomed to. A crossed animal must then be out of harmony, and it is expected that epistasis, if important, has a negative effect.

The genetic additive direct effects negative suggest that yield potential increase by increasing the contribution of genes of the Aberdeen Angus breed; but rather, the positive value for maternal additive effects indicate higher performance for offspring of cows with a higher proportion of Nelore genes. This behavior may be due to the fact that most herds of this population have been created in the Southeast and Midwest regions of Brazil, where there are known major environmental constraints and where the Nellore breed gets better performances.

McMorris & Wilton (1986), studying biologically different genetic groups , noticed that heavier cows and/or major milk producers increase their nutritional needs during the dry and lactation periods. Similarly, Ferrel & Jenkins (1985), studying the interaction between the morphological type of the cows and the environment, reported that cattle require energy for maintenance, growth, pregnancy and lactation, and that the requirements for each of these conditions vary with the type of cattle. Exclusively under pasture without supplemental feeding, the size of the cow may become a very important factor in the efficiency of production systems, since 70-75% of total energy requirements are for maintenance functions. These maintenance requirements vary more often than the requirements for other functions and they seem to be more associated with the genetic potential for production measures (growth rate and milk production), indicating that animals with high genetic potential for production, as is the Aberdeen Angus breed, may have fewer advantages or present more disadvantages in restrictive environment.

Conclusions

The dominant additive models, commonly used in genetic evaluations, do not correctly describe the changes in pre-weaning performance of the cattle in a multiracial population, and the heterozygous and epistatic effects should be included in the analysis models. The inclusion of joint additive effects does not substantially improve the adjustment promoted by the analysis models, besides inserting a bias attributed to multicollinearity. The inclusion of the heterozygous effects in the models is efficient in representing the additive quadratic effect of the breed. The study of the relationship between the independent variables must always precede the choice of a genetic analysis model so that will be estimated effects and breeding values of low accuracy, generating expectations of genetic progress that will not be confirmed.

Acknowledgments

To the companies GenSys Consultores Associados Ltda and Natura Genética Sul Americana for the cession of database. To the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), for the financial support.

Received July 2, 2009 and accepted March 9, 2010.

Corresponding author: [e-mail]

• ARTHUR, P.F.; HEARNSHAW, H.; STEPHENSON, P.D. Direct and maternal additive and heterosis effects from crossing Bos indicus and Bos taurus cattle: cow and calf performance in two environments. Livestock Production Science, v.57, p.231-241, 1999.
• BERGMANN, J.A.G.; HOHENBOKEN, W.D. Alternatives to least squares in multiple linear regression to predict production traits. Journal of Animal Breeding and Genetics, v.112, p.1-16, 1995.
• BOLDMAN, K.G.; KRIESE, L.A.; VAN VLECK, L.D. et al. A manual for use of MTDFREML. A set of program to obtain estimates of variances and covariances (DRAFT) Lincoln: Department of Agriculture/Agricultural Research Service, 1995. 120p.
• CARDOSO, F.F.; TEMPELMAN, R.J. Hierarchical Bayes multiple-breed inference with an application to genetic evaluation of a Nelore-Hereford population. Journal of Animal Science, v.82, p.1589-1601, 2004.
• CORRÊA, M.B.B.; DIONELLO, N.J.L.; CARDOSO, F.F. Estimativa de parâmetros genéticos, componentes de (co)variância e tendências genéticas e fenotípicas para características produtivas pré-desmama em bovinos Devon no Rio Grande do Sul. Revista Brasileira de Zootecnia, v.35, n.3, p.997-1004, 2006 (supl.
• CRUMP, R.E.; HALEY, C.S.; THOMPSON, R. et al. Individual animal model estimates of genetic parameters for performance test traits of male and female Landrace pigs tested in a commercial nucleus herd. Journal of Animal Science, v.65, p.275-283, 1997.
• DOBSON, A.J. An introduction to generalized linear models Melbourne: Chapman and Hall, 1990. 174p.
• DIAS, L.T.; EL FARO, L.; ALBUQUERQUE, L.G. Estimativas de herdabilidade para perímetro escrotal de animais da raça Nelore. Revista Brasileira de Zootecnia, v.32, n.6, p.1878-1882, 2003 (suppl. 2).
• EVERLING, D.M.; FERREIRA, G.B.B.; RORATO, P.R.N. et al. Estimativas de herdabilidade e correlação genética para características de crescimento na fase de pré-desmama e medidas de perímetro escrotal ao sobreano em bovinos Angus-Nelore. Revista Brasileira de Zootecnia, v.30, n.6, p.2002-2008, 2001 (suppl.
• FALCONER, D.S.; MACKAY, T.F.C. Introduction to quantitative genetics Essex: Longman, 1996. 464p.
• FERRAZ FILHO, P.B.; RAMOS, A.A.; SILVA, L.O.C. et al. Tendência genética dos efeitos direto e materno sobre os pesos à desmama e pós-desmama de bovinos da raça Tabapuã no Brasil. Revista Brasileira de Zootecnia, v.31, n.2, p.635-640, 2002.
• FERREL, C.L.; JENKINS, T.G. Cow type and the nutritional environment: nutritional aspects. Journal of Animal Science, v.61, n.3, p.725-741, 1985.
• FREUND, R.J.; LITTELL, R.C. SAS System for regression 3.ed. Cary: SAS Institute Inc., 2000. 235p.
• FREUND, R.J.; WILSON, W.J. Regression analysis: Statistical modeling of a response variable. San Diego: Academic Press, 1998. 444p.
• FRIES, L.A. Cruzamentos em gado de corte. In.: SIMPÓSIO SOBRE PECUÁRIA DE CORTE: PRODUÇÃO DE NOVILHOS DE CORTE, 4., 1996, Piracicaba. Anais... (lecture) Piracicaba: FEALQ/ESALQ, 1996. (CD-ROM).
• FRIES, L.A.; JOHNSTON, D.J.; HEARNSHAW, H. et al. Evidence of epistatic effects on weaning weight in crossbred beef cattle. Asian-Australasian Journal of Animal Sciences, v.13, p.242, 2000 (suppl. B).
• HOERL, A.E. Application of Ridge analysis to regression problems. Chemical Engineering Progress, v.58, p.54-59, 1962.
• HOERL, A.E.; KENNARD, R.W. Ridge regression: Biased estimation for onorthogonal problems. Technometrics, v.12, n.1, p.55-67, 1970.
• KAPS, M.; HERRING, W.O.; LAMBERSON, W.R. Genetic and environmental parameters for traits derived from the Brody Growth curve and their relationships with weaning weight in Angus cattle. Journal of Animal Science, v.78, p.1436-1442, 2000.
• KINGHORN, B.P. Theory of breed utilization. In: Design of livestock breeding programs: short course in animal Breeding. Armidale: AGBU, 1993. p.187-204.
• LOPES, J.S.; RORATO, P.R.N.; WEBER, T. et al. Efeito da interação genótipo x ambiente sobre o peso ao nascimento, aos 205 e aos 550 dias de idade de bovinos da raça Nelore na região Sul do Brasil. Revista Brasileira de Zootecnia, v.37, n.1, p.54-60, 2008.
• MALLOWS, C.P. Some comments on C(p). Technometrics, v.15, p.661-675, 1973.
• MARQUES, L.F.A.; PEREIRA, J.C.C.; OLIVEIRA, H.N. et al. Análise de características de crescimento da raça Simental. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.52, n.5, p.527-533, 2000.
• McMORRIS, M.R.; WILTON, J.W. Breeding systems, cow weight and milk yield effects on various biological variables in beef production. Journal of Animal Science, v.63, p.1361-1372, 1986.
• PIMENTEL, E.C.G.; QUEIROZ, S.A.; CARVALHEIRO, R. et al. Estimativas de efeitos genéticos em bezerros cruzados por diferentes modelos e métodos de estimação.Revista Brasileira de Zootecnia, v.35, n.3, p.1020-1027, 2006 (suppl.
• ROSO, V.M.; SCHENKEL, F.S.; MILLER, S.P. et al. Additive, dominance, and epistatic loss effects on preweaning weight gain of crossbred beef cattle from different Bos taurus breeds. Journal of Animal Science, v.83, p.1780-1787, 2005.
• SHARMA, B.S.; PRABHAKARAN, V.T.; PIRCHNER, F. Gene action and heterosis in lifetime traits of Friesian x Sahiwal crosses. Journal of Animal Breeding and Genetics, v.117, p.319-330, 2000.
• SOUZA, J.C.; CAMPOS, L.T.; FREITAS, J.A. et al. Parâmetros genéticos dos pesos ao nascer e aos 205 dias de idade em animaisda raça Angus no Brasil. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 44., 2007, Jaboticabal. Anais... Jaboticabal: Sociedade Brasileira de Zootecnia, 2007. (CD-ROM).
• STATISTICAL ANALYSIS SYSTEM - SAS. User's guide Version 8. Cary, 2001. 956p.
• TREMATORE, R.L.; ALENCAR, M.M.; BARBOSA, P.B. et al. Estimativas de efeitos aditivos e heteróticos para características de crescimento pré-desmama em bovinos Charolês-Nelore. Revista Brasileira de Zootecnia, v.27, n.1, p.87-94, 1998.
• WEBER, T. Parâmetros genéticos e tendências genéticas e fenotípicas para características produtivas para uma população da raça Aberdeen Angus. 2008. 67f. Dissertação (Mestrado em Zootecnia) - Universidade Federal de Santa Maria, Santa Maria.
• WEISBERG, S. Applied linear regression New York: Wiley Series in probability and mathematical statistics, 1980. 283p.

Publication Dates

• Publication in this collection
  25 Nov 2010
• Date of issue
  Nov 2010

History

• Received
  02 July 2009
• Accepted
  09 Mar 2010