DOI QR코드

DOI QR Code

Computing Algorithm for Genetic Evaluations on Several Linear and Categorical Traits in A Multivariate Threshold Animal Model

범주형 자료를 포함한 다형질 임계개체모형에서 유전능력 추정 알고리즘

  • Published : 2004.04.30

Abstract

Algorithms for estimating breeding values on several categorical data by using latent variables with threshold conception were developed and showed. Thresholds on each categorical trait were estimated by Newton’s method via gradients and Hessian matrix. This algorithm was developed by way of expansion of bivariate analysis provided by Quaas(2001). Breeding values on latent variables of categorical traits and observations on linear traits were estimated by preconditioned conjugate gradient(PCG) method, which was known having a property of fast convergence. Example was shown by simulated data with two linear traits and a categorical trait with four categories(CE=calving ease) and a dichotomous trait(SB=Still Birth) in threshold animal mixed model(TAMM). Breeding value estimates in TAMM were compared to those in linear animal mixed model (LAMM). As results, correlation estimates of breeding values to parameters were 0.91${\sim}$0.92 on CE and 0.87${\sim}$0.89 on SB in TAMM and 0.72~0.84 on CE and 0.59~0.70 on SB in LAMM. As conclusion, PCG method for estimating breeding values on several categorical traits with linear traits were feasible in TAMM.

불연속 범주형 자료에 대한 잠재변수가 존재한다는 가정하에 임계값을 추정하고 잠재변수를 생성하며 생성된 잠재변수 및 기타 연속변량에 대한 관측치를 포함하는 다변량 임계개체모형을 설정하고 유전능력을 예측하기 위한 방법을 제시하였다. 각각의 범주형 조사 자료의 특성을 갖는 형질에 있어서 임계점의 추정은 추정 가능한 임계점에 대한 1차 미분값(gradient)과 2차 미분값(Hessian)을 이용한 Newton 방법을 이용하면 추정가능하며 지역모수인 육종가의 추정은 PCG 방법으로 구현 가능하다. 이러한 이론은 Quaas(2001)가 제시한 하나의 이산형 자료와 하나의 연속형 자료의 2변량 동시 분석방법을 확장하여 전개한 것이며 이때 잠재변수 및 임계점의 추정은 기타 형질의 잔차 회귀계수 및 상관을 고려해야 한다. 본 연구를 위한 모의실험은 2개의 연속변량으로 체중과 유량을 고려하였고 또 다른 2개의 불연속 변량인 분만난이도와 출생시 생존유무를 고려하여 4형질 동시 분석을 실시하였다. 임계모형에 의한 육종가 추정치의 정확도는 4개의 구간으로 분류되어 기록된 분만난이도의 경우에 91${\sim}$92%의 정확도를 보였고 이항분포인 분만시 생존유무에 대하여는 87~89%의 정확도를 보였다. 반면에 이들 범주형 자료를 선형으로 간주하고 분석한 선형 동물개체 혼합모형에서는 72${\sim}$84% 및 59${\sim}$70%으로 비교적 낮은 추정의 정확도를 보였다. 따라서 범주형 자료의 유전분석은 선형 혼합모형 보다 임계형 혼합모형이 크게 타당할 것으로 사료되었다.

Keywords

References

  1. Hendrson, C. R. 1973. Sire evaluation and genetic trends. In: Proceedings of the Animal Breeding and Genetics Symposium in Honour of J.L. Lush. American Society for Animal Science, Blackburgh, Champaign, Illinois.
  2. Janss, L. L. G. and Foulley, J. L. 1993. Bivariate analysis for one continuous and one dichotomous trait with unequal design matrices and an application to birth weight and calving difficulty. Livestock Prod. Sci. 33:183-195. https://doi.org/10.1016/0301-6226(93)90001-X
  3. Falconer, D. S. and Mackay, T. F. C. 1996. Introduction to Quantitative Genetics. Edinburgh, Longman. 4th Ed.
  4. Lee, D. H. 2002. Investigation of biases for variance components on multiple traits with varying number of categories in threshold models using Bayesian inferences. Asian-Aust, J. Anim, Sci. 15 (7):925-931. https://doi.org/10.5713/ajas.2002.925
  5. Lidauer, M., Stranden, I., Mantysaari, E. A., Poso, J. and Kettunen, A. 1999. Solving large test-day models by iteration on data and preconditioned conjugate gradient. J. Dairy Sci. 82:2788-2796.
  6. Meuwissen, T. H. E. and Luo, Z. 1992. Computing inbreeding coefficients in large populations. Genet. Sel. Evol. 24:305-313.
  7. Quaas, R. L 2001. Threshold model. Proceeding of a Symposium at the Organizational Meeting for a NCR Coordinating Committee pp. 56-58.
  8. Schaeffer, L. R. and Kennedy, B. W. 1986. Computing strategies for solving mixed model equations. J. Dairy Sci. 69:575-579. https://doi.org/10.3168/jds.S0022-0302(86)80441-6
  9. Shewchuk, J R. 1994. An introduction to the conjugate gradient method without the agonizing pain. Technical Report CMU-CS-94-125., Carnegie Mellon Univ. Pittsburgh, PA.
  10. Stranden I. and Lidauer, M. 1999. Solving large mixed linear models using preconditioned conjugate gradients iteration. J. Dairy Sci. 82:2779-2787.
  11. Tsuruta, S., Misztal, I. and Stranden, I. 2001. Use of preconditioned conjugate gradient algorithm as a generic solver for mixed-model equations in animal breeding applications. J. Animal Sci. 1166-1172.