• Title/Summary/Keyword: REML method

Search Result 15, Processing Time 0.025 seconds

Comparison between REML and Bayesian via Gibbs Sampling Algorithm with a Mixed Animal Model to Estimate Genetic Parameters for Carcass Traits in Hanwoo(Korean Native Cattle) (한우의 도체형질 유전모수 추정을 위한 REML과 Bayesian via Gibbs Sampling 방법의 비교 연구)

  • Roh, S.H.;Kim, B.W.;Kim, H.S.;Min, H.S.;Yoon, H.B.;Lee, D.H.;Jeon, J.T.;Lee, J.G.
    • Journal of Animal Science and Technology
    • /
    • v.46 no.5
    • /
    • pp.719-728
    • /
    • 2004
  • The aims of this study were to estimate genetic parameters for carcass traits on Hanwoo(Korean Native Cattle) and to compare two different statistical algorithms for estimating genetic parameters. Data obtained from 1526 steers at Hanwoo Improvement Center and Hanwoo Improvement Complex Area from 1996 to 2001 were used for the analyses. The carcass traits considered in these studies were carcass weight, dressing percent, eye muscle area, backfat thickness, and marbling score. Estimated genetic parameters using EM-REML algorithm were compared to those by Bayesian inference via Gibbs Sampling to find out statistical properties. The estimated heritabilities of carcass traits by REML method were 0.28, 0.25, 0.35, 0.39 and 0.51, respectively and those by Gibbs Sampling method were 0.29, 0.25, 0.40, 0.42 and 0.54, respectively. This estimates were not significantly different, even though the estimated heritabilities by Gibbs Sampling method were higher than ones by REML method. Since the estimated statistics by REML method and Gibbs Sampling method were not significantly different in this study, it is inferred that both mothods could be efficiently applied for the analysis of carcass traits of cattle. However, further studies are demanded to define an optimal statistical method for handling large scale performance data.

Restricted maximum likelihood estimation of a censored random effects panel regression model

  • Lee, Minah;Lee, Seung-Chun
    • Communications for Statistical Applications and Methods
    • /
    • v.26 no.4
    • /
    • pp.371-383
    • /
    • 2019
  • Panel data sets have been developed in various areas, and many recent studies have analyzed panel, or longitudinal data sets. Maximum likelihood (ML) may be the most common statistical method for analyzing panel data models; however, the inference based on the ML estimate will have an inflated Type I error because the ML method tends to give a downwardly biased estimate of variance components when the sample size is small. The under estimation could be severe when data is incomplete. This paper proposes the restricted maximum likelihood (REML) method for a random effects panel data model with a censored dependent variable. Note that the likelihood function of the model is complex in that it includes a multidimensional integral. Many authors proposed to use integral approximation methods for the computation of likelihood function; however, it is well known that integral approximation methods are inadequate for high dimensional integrals in practice. This paper introduces to use the moments of truncated multivariate normal random vector for the calculation of multidimensional integral. In addition, a proper asymptotic standard error of REML estimate is given.

Methods and Techniques for Variance Component Estimation in Animal Breeding - Review -

  • Lee, C.
    • Asian-Australasian Journal of Animal Sciences
    • /
    • v.13 no.3
    • /
    • pp.413-422
    • /
    • 2000
  • In the class of models which include random effects, the variance component estimates are important to obtain accurate predictors and estimators. Variance component estimation is straightforward for balanced data but not for unbalanced data. Since orthogonality among factors is absent in unbalanced data, various methods for variance component estimation are available. REML estimation is the most widely used method in animal breeding because of its attractive statistical properties. Recently, Bayesian approach became feasible through Markov Chain Monte Carlo methods with increasingly powerful computers. Furthermore, advances in variance component estimation with complicated models such as generalized linear mixed models enabled animal breeders to analyze non-normal data.

Estimation of Variance Component and Environment Effects on Somatic Cell Scores by Parity in Dairy Cattle (젖소집단의 산차에 따른 체세포점수의 환경효과 및 분산성분 추정)

  • 조광현;나승환;서강석;김시동;박병호;이영창;박종대;손삼규;최재관
    • Journal of Animal Science and Technology
    • /
    • v.48 no.1
    • /
    • pp.39-48
    • /
    • 2006
  • This study utilized test day of somatic cell score data of dairy cattle from 2000 to 2004. The number of data used were 124,635 of first parity, 134,308 of second parity, 77,862 of third parity, 41,787 of forth parity and 37,412 of fifth parity. The data was analyzed by least square mean method using GLM to estimate the effects of calving year, age, lactation stage, parity and season on somatic cell score. Variance component estimation using test day model was determined by using expectation maximization algorithm- restricted maximum likelihood (EM-REML) analysis method. In each parity, somatic cell score was low for younger group and was relatively high in older groups. Likewise, for lactation stage, the score was low in early-lactation and high in late-lactation in first parity and second parity. Nevertheless, for the third, fourth and fifth parity, however, high somatic cell score was observed in mid-lactation. Generally, the score was high in the peak. Although in fourth and fifth parity, the score was low in late-lactation. Environmental effect of season, somatic cell score was generally low from September to November for all parities. The score was high between June and August when the milk production is usually low. The heritability in each parity were 0.05, 0.09, 0.10, 0.05 and 0.05 for parity 1, 2, 3, 4, 5, respectively. Genetic variance value was estimated to be high in second, third and fifth parity in early-lactation and to be low in first and forth parity.

Relationship of Somatic Cell Score and Udder Type Traits of Holstein Cattle (체세포점수와 홀스타인 유방형질간의 관계)

  • Choi, Tae Jeong;Seo, Kang Seok;Kim, Sidong;Park, Byung Ho;Choi, Je Kwan;Yoon, Ho Paek;Na, Seung Hwan;Son, Sam Kyu;Kwon, Oh Sub;Cho, Kwang Hyun
    • Journal of Animal Science and Technology
    • /
    • v.50 no.3
    • /
    • pp.285-292
    • /
    • 2008
  • Data were taken from the dairy herd improve- ment program from the year 2000, composed of 10,929 first lactation cows consisting of 290,144 test-day records and 37,723 udder type records. The objective of the study was to estimate genetic and phenotypic correlation between fore udder attachment, rear udder height, rear udder width, udder cleft, udder depth, and somatic cell score (SCS) and to calculate heritability of udder depth, front teat length and SCS in Holstein cattle in Korea. The variance component estima- tion using test day model was determined by a derivative-free algorithm-restricted maximum likeli- hood(DF-REML) analysis method. Generally phenotypic correlations were very low between udder traits and lactation SCS which varied from -0.03 to -0.06. Heritability of all type traits and SCS was smaller than 0.12. The results of this study would be applicable to SCS using linear genetic evaluation for future studies.

Estimation of Genetic Variance Components of Body Size Measurements in Hanwoo (Korean Cattle) Using a Multivariate Linear Model

  • Lee, Jung-Jae;Kim, Nae-Soo
    • Journal of Animal Science and Technology
    • /
    • v.52 no.3
    • /
    • pp.167-174
    • /
    • 2010
  • The objectives of this study were to quantify the combination values of the principal components and factors calculated using body measurements of Hanwoo (Korean Cattle) and estimate their heritabilities. The technique of multivariate analysis was used to reduce a large number of variables to a smaller number of new variables and characterize cattle according to body shape. The analyses were performed using 1,979 cattle at 12 months of age and 936 cattle at 24 months of age. The data for the analyses was obtained from progeny tests performed on Korean Cattle for 6 years from 2003 to 2008. The phenotypic correlations among these traits were estimated to range from 0.32 to 0.90 at 12 months of age and from 0.21 to 0.82 at 24 months of age. The first principal components (PC1s) indicated a weighed average of overall body measurements, accounting for 99.91% of the total variation for both periods of test. The two first PCs had positive coefficients for all body measurements. The major sources of PC, such as chest girth (CG), body length (BL), rump height (RH), and wither height (WH) were similar for both test periods. The heritabilities for PC1, the first factor score (FS1), and the second factor score (FS2) were estimated by multivariate REML method. The estimated heritabilities for PC1, FS1, and FS2 were 0.33, 0.38, and 0.40, respectively, at 12 months of age and 0.26, 0.76, and 0.58 at 24 months of age. Further studies are needed to determine whether the heritabilities of FS1 and FS2 at 24 months of age were overestimated.

A mixed model for repeated split-plot data (반복측정의 분할구 자료에 대한 혼합모형)

  • Choi, Jae-Sung
    • Journal of the Korean Data and Information Science Society
    • /
    • v.21 no.1
    • /
    • pp.1-9
    • /
    • 2010
  • This paper suggests a mixed-effects model for analyzing split-plot data when there is a repeated measures factor that affects on the response variable. Covariance structures are discussed among the observations because of the assumption of a repeated measures factor as one of explanatory variables. As a plausible covariance structure, compound symmetric covariance structure is assumed for analyzing data. The restricted maximum likelihood (REML)method is used for estimating fixed effects in the model.

Power analysis of testing fixed effects with two way classification (이원혼합모형에서 고정효과 유의성검정에 대한 검정력 분석)

  • 이장택
    • The Korean Journal of Applied Statistics
    • /
    • v.10 no.1
    • /
    • pp.177-187
    • /
    • 1997
  • This article considers the power performance of the tests in unbalanced two way mixed linear models with one fixed factor. The generalized least squares (GLS) F statistic testing no differences among the effects of the levels of the fixed factor is estimated using Henderson's method III, minimum norm quadratic unbiased estimator (MINQUE) with prior guess 1, maximum likelihood (ML) and resticted maximum likelihood (REML). We investigate the power performance of these test statistics. It can be shown, through simulation, that the GLS F statistics using four estimators produce similar type I error rates and power performance.

  • PDF

Power Comparison of EGLS Test Statistic for Fixed Effects with Arbitrary Distributions

  • Lee, Jang-Taek
    • Communications for Statistical Applications and Methods
    • /
    • v.10 no.1
    • /
    • pp.11-18
    • /
    • 2003
  • Quite often normality assumptions are not satisfied in practical applications. In this paper, an estimated generalized least squares(EGLS) analysis are considered in two way mixed linear models with arbitrary types of distributions for random effects. We investigate the power performance of EGLS analysis based on Henderson's method III, ML, REML and MINQUE(1). The power performances depend on the imbalance of design, on the actual values of ratio of variance components, and on the skewness and kurtosis parameters of the underlying distributions slightly. Results of our limited simulation study suggest that the EGLS F-statistics using four estimators and arbitrary distributions produce similar type I error rates and power performance.

Heritability Estimates under Single and Multi-Trait Animal Models in Murrah Buffaloes

  • Jain, A.;Sadana, D.K.
    • Asian-Australasian Journal of Animal Sciences
    • /
    • v.13 no.5
    • /
    • pp.575-579
    • /
    • 2000
  • First lactation records of 683 Murrah buffaloes maintained at NDRI, Karnal which were progeny of 84 sires used for comparing the heritability estimates of age at first calving, first lactation milk yield and first service period under single and multiple trait models using restricted maximum likelihood (REML) method of estimation under an individual animal model. The results indicated that the heritability estimates may vary under single and multiple trait models depending upon the magnitude of genetic and environmental correlation among the traits being considered. Therefore, a single or multiple trait model is recommended for estimation of variance components depending upon the goal of breeding programme. However, there may not be any advantage of considering a trait with zero or near zero heritability and having no or very low genetic correlation with other traits in the model. Lower heritability estimates of part lactation yield (120-day milk yield) implied that there may not be any advantage of considering this trait in place of actual 305-day milk yield, whereas, comparable heritability estimates of predicted 305-day milk yield suggested that it could be used for sire evaluation to reduce the cost of milk recording under field conditions.