Simulation Study on Covariance Component Estimation for Two Binary Traits in an Underlying Continuous Scale

E. A. Mntysaari ¹, R. L. Quaas ¹, and Y. T. Gröhn ¹

¹ Cornell University, Ithaca, NY 14853

The usefulness of the variance and covariance component estimation methods based on a threshold model was studied in a multiple-trait situation with two binary traits. Estimation equations that yield marginal maximum likelihood estimates of variance components on the underlying continuous variable scale and point estimates of location parameters with empirical Bayesian properties are described. Methods were tested on simulated data sets that were generated to exhibit three different incidences, 25, 15, and 5%. Results were compared with analyses of the same data sets with a REML method based on normal distribution and a linear model. Heritabilities and residual correlations calculated from discrete observations were transformed to underlying parameters.

In estimation of heritabilities, all methods performed equally well at all incidence levels and with no detectable bias. As suggested by threshold theory, the genetic correlation was accurately estimated directly from the observations without any need of correction for incidence. Marginal maximum likelihood estimates of genetic correlations were similar to linear model estimates; discrepancies from the true parameters were consistent with both methods. In estimation of residual correlations, the method with the linear model approach yielded satisfactory estimates only at the highest incidence level, 25%. For 5% incidence, the uncorrected estimate of residual correlation was 50% less than the true value, and after correction for incidence, the parameter was overestimated by 90%. The estimates of residual correlation from the threshold model were regarded fair, except at the lowest level of incidence, where the estimate was 27% higher than the true value. Results indicated that when an accurate estimate of residual correlation is needed, the marginal maximum likelihood estimates are superior to the estimates calculated with the linear model. Using correction for the incidence level for residual correlation did not work well except at the highest incidence level.

Key Words: binary traits • variance components • Bayesian methods

Submitted on March 2, 1990
Accepted on June 6, 1990

This article has been cited by other articles:

				H. Vehvilainen, A. Kause, C. Quinton, H. Koskinen, and T. Paananen Survival of the Currently Fittest: Genetics of Rainbow Trout Survival Across Time and Space Genetics, September 1, 2008; 180(1): 507 - 516. [Abstract] [Full Text] [PDF]

				E. Carlen, U. Emanuelson, and E. Strandberg Genetic evaluation of mastitis in dairy cattle using linear models, threshold models, and survival analysis: a simulation study. J Dairy Sci, October 1, 2006; 89(10): 4049 - 4057. [Abstract] [Full Text] [PDF]

				A. C. Sorensen, P. Madsen, M. K. Sorensen, and P. Berg Udder health shows inbreeding depression in danish holsteins. J Dairy Sci, October 1, 2006; 89(10): 4077 - 4082. [Abstract] [Full Text] [PDF]

				H. Mesa, T. J. Safranski, K. M. Cammack, R. L. Weaber, and W. R. Lamberson Genetic and phenotypic relationships of farrowing and weaning survival to birth and placental weights in pigs J Anim Sci, January 1, 2006; 84(1): 32 - 40. [Abstract] [Full Text] [PDF]

				M. Koivula, E. A. Mantysaari, E. Negussie, and T. Serenius Genetic and Phenotypic Relationships Among Milk Yield and Somatic Cell Count Before and After Clinical Mastitis J Dairy Sci, February 1, 2005; 88(2): 827 - 833. [Abstract] [Full Text] [PDF]

				J. Lassen, M. Hansen, M. K. Sorensen, G. P. Aamand, L. G. Christensen, and P. Madsen Genetic Relationship Between Body Condition Score, Dairy Character, Mastitis, and Diseases Other than Mastitis in First-Parity Danish Holstein Cows J Dairy Sci, November 1, 2003; 86(11): 3730 - 3735. [Abstract] [Full Text] [PDF]