Method R and Restricted Maximum Likelihood (REML) were compared for estimating heritability (h2) and subsequent prediction of breeding values (a) with data subject to selection. A single-trait animal model was used to generate the data and to predict breeding values. The data originated from 10 sires and 100 dams and simulation progressed for 10 overlapping generations. In simulating the data, genetic evaluation used the underlying parameter values and sires and dams were chosen by truncation selection for greatest predicted breeding values. Four alternative pedigree structures were evaluated: complete pedigree information, 50% of phenotypes with sire identities missing, 50% of phenotypes with dam identities missing, and 50% of phenotypes with sire and dams identities missing. Under selection and with complete pedigree data, Method R was a slightly less consistent estimator of h2 than REML. Estimates of h2 by both methods were biased downward when there was selection and loss of pedigree information and were unbiased when no selection was practiced. The empirical mean square error (EMSE) of Method R was several times larger than the EMSE of REML. In a subsequent analysis, different combinations of generations selected and generations sampled were simulated in an effort to disentangle the effects of both factors on Method R estimates of h2. It was observed that Method R overestimated h2 when both the sampling that is intrinsic in the method and the selection occurred in generations 6 to 10. In a final experiment, BLUP(a) were predicted with h2 estimated by either Method R or REML. Subsequently, five more generations of selection were practiced, and the mean square error of prediction (MSEP) of BLUP(a) was calculated with estimated h2 by either method, or the true value of the parameter. The MSEP of empirical BLUP(a) using Method R was greater than the MSEP of empirical BLUP(a) using REML. The latter statistic was closer to prediction error variance of BLUP(a) than the MSEP of empirical BLUP(a) using Method R, indicating that empirical BLUP(a) calculated using REML produced accurate predictions of breeding values under selection. In conclusion, the variability of h2 estimates calculated with Method R was greater than the variability of h2 estimates calculated with REML, with or without selection. Also, the MSEP of EBLUP(a) calculated using estimates of h2 by Method R was larger than MSEP of EBLUP(a) calculated with REML estimates of h2.
-A Markov chain Monte Carlo (MCMC) algorithm to sample an exchangeable covariance matrix, such as the one of the error terms (R 0 ) in a multiple trait animal model with missing records under normal-inverted Wishart priors is presented. The algorithm (FCG) is based on a conjugate form of the inverted Wishart density that avoids sampling the missing error terms. Normal prior densities are assumed for the 'fixed' effects and breeding values, whereas the covariance matrices are assumed to follow inverted Wishart distributions. The inverted Wishart prior for the environmental covariance matrix is a product density of all patterns of missing data. The resulting MCMC scheme eliminates the correlation between the sampled missing residuals and the sampled R 0 , which in turn has the effect of decreasing the total amount of samples needed to reach convergence. The use of the FCG algorithm in a multiple trait data set with an extreme pattern of missing records produced a dramatic reduction in the size of the autocorrelations among samples for all lags from 1 to 50, and this increased the effective sample size from 2.5 to 7 times and reduced the number of samples needed to attain convergence, when compared with the 'data augmentation' algorithm.FCG algorithm / multiple traits / missing data / conjugate priors / normal-inverted Wishart
Contemporary groups (CG) are used in genetic evaluation to account for systematic environmental effects of management, nutritional level, or any other differentially expressed group effect; however, because the functional form of the distribution of those effects is unknown, CG serve as an approximation to a time-varying mean. Conversely, in semiparametric models, there is no need to assume any functional form for the time-varying effects. In this research, we present a semiparametric animal model (AMS) using the covariate day of birth (DOB) by means of penalized splines (P-splines), as an alternative to fitting CG. In the AMS, the functionality of the data on DOB is expressed by means of a Basic segmented polynomial line (B-spline) basis, and proper covariance matrices are used to reflect the serial correlation among the points of support (or knots) at different times. Three different covariance matrices that reflect either short- or long-range dependences among knots are discussed. Different models were fitted to birth weight data from Polled Hereford calves. Models compared were an animal model with CG, an animal model with CG and the covariate DOB nested within CG (CG + DOB), and P-splines with the first difference penalty matrix and three different AMS with 20, 40, 60, 80, or 120 knots. Models were compared using a modified Akaike information criterion (AICC), which was calculated as a byproduct of the estimation of variance components by REML using the expectation maximization algorithm. All three AMS had smaller (better) values of AICC than the regular model with CG, while producing almost the same ranking of predicted breeding values and similar average predicted error variance. In all AMS, the inference and all measures of comparison were similar when the number of knots was equal > or = 40. The model CG + DOB had analogous performance to the AMS, but at the expense of using more parameters. It is concluded that the use of penalized regression splines using a B-spline basis with proper covariance matrices is a competitive method to the fitting of CG into animal models for genetic evaluation, without having to assume any functional form for the covariate DOB.
Abstract. Growth and carcass data of Angus cattle were used to estimate heritabilities and genetic and environmental correlations between growth and carcass traits by means of a Bayesian data augmentation (DA) algorithm. Records were taken on 739 Angus steers from 31 sires, during 10 years of a designed progeny test. The cattle were entirely fed on grass during their lifelong. Growth traits evaluated were birth (BW), weaning (WW) and 18-month (W18) weights; and carcass traits were the weights of half the carcass (HCW), of hind "pistola" cut (HPW) and of three retail cuts (ECW). The model used for estimation was a multiple trait additive animal model. The prior densities used in the analyses were the multivariate normal for the fixed effects (with very large variances) and for the breeding values, and the inverted Wishart for the additive and environmental covariance matrices. The observed residual vector was augmented with sampled residuals for missing traits. The total number of samples drawn was 200,000. The heritabilities of growth traits increased with age at measure, and those of carcass traits were of sizeable magnitude. Whereas estimates of the genetic correlations were similar to those found in the literature for cattle fed on concentrates, environmental correlations were lower. Additive correlations between growth traits with either the HPW or ECW, were smaller than the correlations between growth characters and HCW.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.