RESTRICTED M A X I M U M LIKELIHOOD TO E S T I M A T E GENETIC P A R A M E T E R S - IN PRACTICE
|
|
- Ira Owen
- 5 years ago
- Views:
Transcription
1 RESTRICTED M A X I M U M LIKELIHOOD TO E S T I M A T E GENETIC P A R A M E T E R S - IN PRACTICE K. M e y e r Institute of Animal Genetics, Edinburgh University, W e s t M a i n s Road, Edinburgh EH9 3JN, Scotland, U.K. S U M M A R Y An efficient computing strategy for estimating variance and covariance components between and within random effects by Restricted Maximum Likelihood is described. It pertains to data with one random factor and equal design matrices for all traits or traits recorded on different sets of animals. INTRODUCTION Until recently, the estimation of variance (and covariance c o m p o n e n t s f r o m animal breeding d a t a has relied a l m o s t exclusively on Henderson's (1953) m e t h o d s for unbalanced m i xed models. M e t h o d III especially has found w i d e s p r e a d use, greatly aided by t h e availability of a 'general purpose' c o m p u t e r p r o g r a m tailored t o w a r d s the estimation of genetic p a r a m e t e r s (Harvey, 1960 and 1977). Interest has g r o w n since in m a x i m u m likelihood (ML) and related procedures w h i c h yield estimators w i t h m u c h m o r e desirable properties (Harville, 1977). The so-called Restricted M a x i m u m Likelihood (REML), d e veloped by P a t t e r s o n and T h o m p s o n (1971), has b e c o m e accepted as t h e prefered m e t h o d to e s t i m a t e (co)variance c o m p o n e n t s in animal breeding. T h o m p s o n (1982) discussed R E M L for estimating genetic p a r a m e t e r s in different cases. In spite of an extensive theoretical treatment, R E M L has b e e n put to comparatively little u s e in practice. Algorithms are complicated a n d c o m putational r e q u i r e m e n t s are extensive, especially for multivariate analyses or models w i t h m o r e than one r a n d o m factor. Although R E M L or analogous procedures (MINQUE, MIVQUE) are n o w available in s o m e general statistical c o m p u t e r packages (e.g. GENSTAT, Robinson et a/., 1982, or SAS), t h ese are often n o t suitable for the analysis of large d a t a sets f r o m livestock i m p r o v e m e n t schemes. This paper describes a c o m putational a p p r o a c h to R E M L estimation of (co)variance c o m p o n e n t s b e t w e e n and within genetic groups, considering a simple m o d e l but possibly large a m o u n t s of data, as appropriate, for instance, to the analysis of dairy cattle data. M O D E L O F ANALYSIS The computing strategy described allows for a linear m o del w i t h : One "major" fixe d e ffe c t. This c a n be simply an overall m e a n or, on the o t her extreme, a s y s t e m a t i c environmental factor w i t h a very large n u m b e r of levels, as, for instance, h e r d - y e a r - s e a s o n effects for dairy cattle data. All levels of this effect are absorbed and e s t i m a t e s are not obtained. Several additional ("minor") fixe d e ffe c ts (optional). These factors are a s s u m e d to h a v e relatively f e w levels and estimates of their effects are determined. For dairy records this could be, for example, sire country of origin or m o n t h of calving within season. Several covariables (optional). A linear or higher order regression on each covariable can be fitted. E s t i m a t e s of regression coefficients are obtained. For t h e dairy e x a m p l e this could be age at calving fitted as a linear and quadratic covariable. 454
2 One random e ffe c t. This is c o m m o n l y s o m e genetic grouping, e.g. sires or families. Optionally, s o m e of the r a n d o m effect levels can be t r e a t e d as fixed, i.e. records pertaining to t h ese levels contribute to t h e estimation of the within but n o t b e t w e e n g r o u p (co)variance c o m ponents. This option is of particular relevance in t h e analysis of dairy cattle d a t a w h e r e young, unselected sires are used c o n t e m p o r a n e o u s l y to proven, selected sires. T o avoid bias due to selection only the variance b e t w e e n y o u n g sires should be utilised to e s t i m a t e the additive genetic variance. H o w e v e r, for such d a t a the m o d e l of analysis often includes h e r d - y e a r - s e a s o n s as fixed effects. Considering y o u n g sire records only then results in v e r y small subclass sizes while including records for proven sires, t r e a t e d as fixed, increases subclass sizes a n d improves the c o n n e c t e d n e s s in the d a t a and thus enhances the accuracy of estimation (Meyer, 1983; V a n Vleck, 1985). For each trait, b o t h r a n d o m effects and residual errors are a s s u m e d to be identically distributed w i t h m e a n zero. All covariances a m o n g errors a n d b e t w e e n errors a n d r a n d o m effects are a s s u m e d to be zero. A specific covariance structure a m o n g the r a n d o m effect levels due to relationships b e t w e e n animals can optionally be t a k e n into a c c o u n t by including the n u m e r a t o r relationship matrix. C O M P U T I N G STRATEGY Largely the s a m e c o mputational steps are required for univariate analyses a n d multivariate analyses w h e r e the w i t h i n g r oup ("error") covariances b e t w e e n traits are zero. The latter e n c o m p a s s e s t w o special cases -. Firstly, different traits are m e a s u r e d on distinct sets of animals. Secondly, the design matrices for all traits are equal. T h e n a transformation to canonical scale is feasible w h i c h leaves all traits uncorrelated, b o t h genetically and phenotypically (Thompson, 197 6). Figure 1 s u m m a r i z e s the sequence of steps of analysis. Absorption I T h e first step is to set u p the m i x e d m o d e l equations (MME) (Henderson, 1973) for the r a n d o m effect a n d a n y additional fixed effects a n d covariables, absorbing the major fixed effect. Ordering d a t a according to levels of this effect and r a n d o m levels w i t h i n subclasses, this can be d o n e very efficiently for one level at a time. A suitable p r o g r a m m i n g strategy is outlined by Schaeffer (1975). This yields the coefficient matrix, a right hand side (RHS) for each trait a n d the s u m s of squares and crossproducts (SS/CP) within subclasses. Absorption II If additional fixed effects or covariables are fitted, these are a b s o r b e d simultaneously into the r a n d o m effect levels, after all major fixed effect subclasses h a v e been processed. This requires the direct inverse of a matrix of order equal to the n u m b e r of additional fixed effects levels and regression coefficients. The inverse t o g e t h e r w i t h o t h e r parts of the M M E required to backsolve for the a b s o r b e d effects are saved. Residual SS/ CP are adjusted for the variation removed. Absorption III If a n y levels of the r a n d o m effect are to be t r e a t e d as fixed these are a b s o r b e d now. Again an inverse of order equal to t h e n u m b e r of levels to be a b s o r b e d is required, information to obtain backsolutions is saved an d the residual SS/CP are adjusted. 455
3 Relationship M a t r i x In m i x e d m o d e l methodology, the relationships b e t w e e n animals are usually t a k e n into a c c o u n t by adding the inverse of the n u m e r a t o r relationship matrix, w e i g h t e d according to a function of the within and b e t w e e n (co)variances, to the coefficient matrix for t h e r a n d o m effects. This is feasible for a large n u m b e r of r a n d o m levels as the inverse can be set up directly f r o m a list of pedigree information (Henderson, 1976). For an iterative estimation procedure, however, use of a n equivalent model, as suggested by Quaas (1984), is m o r e efficient. L e t u d e n o t e the vector of r a n d o m effects for a trait, Z the design matrix relating it to the d a t a vector ^ a n d A the relationship matrix. The variance matri^c of u is2 t h e n A o 0 a n d its contribution to the variance of y is ZAZ'c>b w i t h ob the variance b e t w e e n r a n d o m effects. D e c o m p o s i n g A into LL' jjjives the latter as ZLL'Z'o 2. Hence defining a n e w design m a t r i x Z =ZL yields an equivalent 8 m o d e l w h i c h incorporates the relationships b e t w e e n animals but has V(u)=Io 2. L is a l o w e r triangular matrix w h i c h also has s o m e interpretation in describing the gene flow b e t w e e n generations. In computational t e r m s this requires to pre- and postmultiply the coefficient m a t r i x for r a n d o m effects w i t h L' a n d L, respectively, a n d to premultiply each RHS by L'. L c a n be set up f r o m a list of pedigree information one column at a t i m e as described by Quaas (197 6). The necessary multiplications can be carried out likewise, considering only the n o n - z e r o elements of each c o l u m n of L. Due to the triangular f o r m of L, m o s t of the n e w M M E can be stored in the original arrays (Meyer, 1985a). Tridiagonalisation T h e major c o m p u t a t i o n a l r e q u i r e m e n t in estimating (co)variance c o m p o n e n t s by R E M L is the inversion a matrix of order equal to a multiple to the n u m b e r of r a n d o m effect levels for each r o und of iteration. E v e n for a multiple of one (uni- or multivariate 'canonical' analysis) resources required can be considerable. Fortunately, in this case, a t r a n s f o r m a t i o n can eliminate the n e e d for a direct inverse. Q u a a s and S m ith (1981, cit. Taylor et a/., 1985) s u g g e s t e d to apply a series of Householder transformations to t h e M M E for r a n d o m effects w h i c h leaves the resulting coefficient matrix tridiagonal. R a n d o m effects solutions can t h e n be obtained in linear time and, moreover, the trace of the inverse of t h e coefficient matrix can be a c c u m u l a t e d indirectly, utilising t e r m s w h i c h arise in obtaining t h e solutions. A detailed outline of the procedure can be found in S m i t h a n d Graser (1985) and Taylor et a/. (1985). Transforming the M M E, coefficient matrix and all RHS, requires just s o m e w h a t m o r e time t h a n to invert a matrix of the s a m e size once. For N r a n d o m levels, N-2 Householder matrices are required, the i-th m a t r i x characterised by a vector of length N-i. T h ese vectors are s a ved to allow a backtransformation. Estimation In the univariate analysis o n e RHS is considered at a time. Variances b e t w e e n a n d within r a n d o m effects are e s t i m a t e d using an E M - a l g o r i t h m (Dempster et a/., 1977) on the tridiagonal scale t o g e t h e r w i t h a reparameterisation to speed up convergence, as suggested by T h o m p s o n a n d M e y e r (1985), iterating until the c h a n g e in e s timates b e t w e e n rounds reaches a specified minimum. At convergence, r a n d o m effect solutions are t r a n s f o r m e d to the original scale. If applicable, backsolutions for the fixed levels of the r a n d o m effect and subsequently a n y additional fixed effects a n d regression coefficients are obtained. 456
4 W i t h a M e t h o d of Scoring algorithm, large sample s t a n d a r d errors of the e s t i m a t e s w o u l d arise as a b y - p r o d u c t f r o m t h e inverse of the information matrix. To determine the information m a t r i x requires the trace of the square of the inverse of the coefficient matrix. W i t h the EM-algorithm, a good approximation of this quantity, a n d hence the sampling errors, can be obtained by numerical differentiation (Thompson, pers. comm.), requiring computational effort equivalent to less t h a n half a r o und of iteration. Using a t r a n s f o r m a t i o n to canonical scale, the multivariate analysis of q traits w i t h equal design matrices reduces to q corresponding univariate analyses (Meyer, 1985b). For an E M - a l g o r i t h m e s t i m a t e s can t h e n again be obtained on the tridiagonal scale w i t h o u t inverting the coefficient matrix; see Q u a a s and S m i t h (1981, cit. Taylor e t at., 1985) for a detailed description. A reparameterisation to increase the speed of c o n v e r g e n c e analogous to the univariate case w a s disregarded here to exploit the property of the E M - a l g o r i t h m t h a t it forces estimates to be w i t h i n the p a r a m e t e r space (Harville, 1977). A t convergence, backsolutions on the canonical scale are obtained as a b o v e for one trait at a time and t r a n s f o r m e d to the original scale subsequently. Again numerical differentiation (Thompson, pers. comm.) is utilised to d e t e r m i n e additional traces required in the information matrix. This requires t i m e equivalent to about q /2 multivariate (or q2/2 univariate) rounds of iteration. A p p r o x i m a t e sampling variances of the b e t w e e n and within c o m p o n e n t s o n the canonical scale are d e t e r m i n e d using formulae for the e l e m e n t s of (twice) the information m a t r i x given by M e y e r (1985b). T r a n s forming these then gives a half stored s y m m e t r i c m a t r i x of order q(q+1), describing the (co)variance structure of e s timates on the original scale. A p p r o x i m a t e sampling errors of genetic parameters, derived f r o m the e s timated (co)variance c o mponents, are t h e n calculated using formulae obtained by statistical differentiation. If different traits are m e a s u r e d on differen t sets of animals, d a t a h a v e to be split according to traits a n d t h e first 4 steps (if applicable) be carried o u t for each seperately. The resulting M M E are t h e n combined, yielding a coefficient matrix of order n u m b e r of traits * n u m b e r of r a n d o m levels to be inverted for each r o u n d of iteration. Variances and covariances b e t w e e n and variances within r a n d o m effects are e s t i m a t e d using a M e t h o d of Scoring t y p e algorithm (Schaeffer et a/., 1978). At c o n v e r g e n c e backsolutions for a b s o r b e d effects are obtained as a b ove for one trait at a time. C O M M E N T S The c o m p u t i n g strategy presented has b e e n u s e d to analyse dairy cattle d a t a w i t h up to 22 traits and m o r e t h a n 600 sires. T h o u g h requiring s o m e care in manipulating input a n d o u t p u t files f r o m the various steps, it has proven to be e a s y - t o - u s e a n d efficient. All information specifying the model of analysis and r u n - t i m e options (convergence criteria, trait(s) to be analysed a n d optionally starting values) are specified interactively in free format. T h e univariate analysis a p p e a r e d to be robust against starting values quite different f r o m the eventual estimates. It has b e e n found g o o d practice to use univariate variance estimates, together w i t h g u e s s e d or zero covariances, as starting values for a multivariate analysis, gradually increasing the n u m b e r of traits considered simultaneously. P r o g r a m s are w r i t t e n in F O R T R A N 77 and self-contained; listings can be m a d e available on request. 457
5 ACKNOW LEDGEMENTS This study w a s supported by the Agricultural and Food Research Council. REFERENCES D E MPSTER, A.P., LAIRD,N.M. and RUBIN, D.B M a x i m u m likelihood f r o m incomplete d a t a w i t h t h e E M algorithm. J. Roy. Stat. Soc. B 39,1-22. HARVEY, W.R L e a s t squares analysis of data w i t h unequal subclass numbers. USDA, ARS HARVEY, W.R Users guide for LSML76. M i x e d m o d e l least squares and m a x i m u m likelihood c o m p u t e r program. USDA, ARS. HARVILLE, D.A M a x i m u m likelihood approaches to variance c o m p o n e n t estimation and to related problems. J. Amer. Stat. Ass. 72, HENDERSON, C.R Estimation of variance and covariance components. Biometrics 9, HENDERSON, C.R Sire evaluation and genetic trends. Proc. Anim. Breed. Genet. Symp. in Honour o f Dr.J.L. Lush, Amer. Soc. Anim. Sci., pp HENDERSON, C.R A simple m e t h o d for c o m p u t i n g t h e inverse of a n u m e r a t o r relationship matrix used in prediction of breeding values. Biometrics 32, MEYER, K M a x i m u m likelihood procedures for estimating genetic p a r a m e t e r s for later lactations in dairy cattle. J. Dairy Sci. 66, MEYER, K. 1985a. Restricted M a x i m u m likelihood estimation of variance c o m p o n e n t s accounting for relationships b e t w e e n animals. Unpubl. man. MEYER, K. 1985b. M a x i m u m likelihood estimation of variance c o m p o n e n t s for a multivariate m i x e d m o d e l w i t h equal design matrices. Biometrics 41, PATTERSON, H.D. and T H O M P S O N, R Recovery of inter-block information w h e n block sizes are unequal. Biometrika 58, QUAAS, R.L C o m p u t i n g the diagonal elements and inverse of a large n u m e r a t o r relationship matrix. Biometrics 32, QUAAS, R.L Linear prediction. In BLUP School Handbook, AGBU, Univ. N e w England, Armidale, pp ROBINSON, D.L., T H O M P S O N, R. and DIGBY, P.G.N R E M L - a p r o g r a m for the analysis of n o n - o r t h o g o n a l d a t a by Restricted M a x i m u m Likelihood. In Proc. COMPSTAT Symp. : Short communications, A u g u s t SCHAEFFER, L.R Dairy sire evaluation fo r milk and fa t production. Dept. Anim. Poultry Sci., Uni. Guelph. SCHAEFFER, L.R., WILTON, J.W. and THOMPSON, R Simultaneous estimation of variance a n d covariance c o m p o n e n t s f r o m multitrait m i x e d m o del equations. Biom etrics 34, SMITH, S.P. and GRASER, H.U Estimating variance c o m p o n e n t s in a class of models by Restricted M a x i m u m Likelihood. J. Dairy Sci. in press. TAYLOR, J.F., BEAN, B., M A R S H A L L, C.E. and SULLIVAN, J.J Genetic and environmental c o m p o n e n t s of s e m e n production traits of artificial insemination Holstein bulls. J. Dairy Sci. 68, THOMPSON, R Estimation of quantitative genetic parameters. In Proc. Int. Conf. Quantit. Genet., Ames, Iowa,pp THOMPSON,R M e t h o d s of estimation of genetic parameters. In Proc. 2nd World Congr. Genet. Appl. Livest. Prod., Madrid, V, THOMPSON, R. and MEYER, K Estimation of variance c o m p o n e n t s : W h a t is missing in the E M algorithm? S u b m i t t e d to J. S tatist. Comput. V A N VLECK, L.D Including records of daughters of selected bulls in estimation of sire c o m p o n e n t s of variance. J. Dairy Sci. 68,
6 Figure 1 : Schematic outline of computational steps for the analysis of d a t a w i t h one r a n d o m factor and equal design matrices for all traits or traits m e a s u r e d on different animals. DATA I ABSORPTION I SET UP MME ABSORBING THE MAJOR FIXED EFFECT ABSORPTION II MME ABSORBING ALL FIXED EFFECTS ABSORPTION III MME FOR RANDOM LEVELS ONLY UNIVARIATE ANALYSIS A * f MULTIVARIATE ANALYSIS FOR EQUAL DESIGN MATRICES 459
Alternative implementations of Monte Carlo EM algorithms for likelihood inferences
Genet. Sel. Evol. 33 001) 443 45 443 INRA, EDP Sciences, 001 Alternative implementations of Monte Carlo EM algorithms for likelihood inferences Louis Alberto GARCÍA-CORTÉS a, Daniel SORENSEN b, Note a
More informationSummary INTRODUCTION. Running head : AI-REML FOR EQUAL DESIGN MATRICES. K. Meyer
Running head : AI-REML FOR EQUAL DESIGN MATRICES An average information Restricted Maximum Likelihood algorithm for estimating reduced rank genetic covariance matrices or covariance functions for animal
More informationEstimating genetic parameters using
Original article Estimating genetic parameters using an animal model with imaginary effects R Thompson 1 K Meyer 2 1AFRC Institute of Animal Physiology and Genetics, Edinburgh Research Station, Roslin,
More informationUse of sparse matrix absorption in animal breeding
Original article Use of sparse matrix absorption in animal breeding B. Tier S.P. Smith University of New England, Anirreal Genetics and Breeding Unit, Ar!nidale, NSW 2351, Australia (received 1 March 1988;
More informationVARIANCE COMPONENT ESTIMATION & BEST LINEAR UNBIASED PREDICTION (BLUP)
VARIANCE COMPONENT ESTIMATION & BEST LINEAR UNBIASED PREDICTION (BLUP) V.K. Bhatia I.A.S.R.I., Library Avenue, New Delhi- 11 0012 vkbhatia@iasri.res.in Introduction Variance components are commonly used
More informationApproximation of Sampling Variances and Confidence Intervals for Maximum Likelihood Estimates of Variance Components. Abstract
Running head : Approximate REML sampling variances Approximation of Sampling Variances and Confidence Intervals for Maximum Likelihood Estimates of Variance Components K. Meyer 1 and W.G. Hill Institute
More informationMIXED MODELS THE GENERAL MIXED MODEL
MIXED MODELS This chapter introduces best linear unbiased prediction (BLUP), a general method for predicting random effects, while Chapter 27 is concerned with the estimation of variances by restricted
More informationModification of negative eigenvalues to create positive definite matrices and approximation of standard errors of correlation estimates
Modification of negative eigenvalues to create positive definite matrices and approximation of standard errors of correlation estimates L. R. Schaeffer Centre for Genetic Improvement of Livestock Department
More informationMixed-Model Estimation of genetic variances. Bruce Walsh lecture notes Uppsala EQG 2012 course version 28 Jan 2012
Mixed-Model Estimation of genetic variances Bruce Walsh lecture notes Uppsala EQG 01 course version 8 Jan 01 Estimation of Var(A) and Breeding Values in General Pedigrees The above designs (ANOVA, P-O
More informationLecture 9 Multi-Trait Models, Binary and Count Traits
Lecture 9 Multi-Trait Models, Binary and Count Traits Guilherme J. M. Rosa University of Wisconsin-Madison Mixed Models in Quantitative Genetics SISG, Seattle 18 0 September 018 OUTLINE Multiple-trait
More informationEstimating Variances and Covariances in a Non-stationary Multivariate Time Series Using the K-matrix
Estimating Variances and Covariances in a Non-stationary Multivariate ime Series Using the K-matrix Stephen P Smith, January 019 Abstract. A second order time series model is described, and generalized
More informationGenetic Parameter Estimation for Milk Yield over Multiple Parities and Various Lengths of Lactation in Danish Jerseys by Random Regression Models
J. Dairy Sci. 85:1596 1606 American Dairy Science Association, 2002. Genetic Parameter Estimation for Milk Yield over Multiple Parities and Various Lengths of Lactation in Danish Jerseys by Random Regression
More informationChapter 12 REML and ML Estimation
Chapter 12 REML and ML Estimation C. R. Henderson 1984 - Guelph 1 Iterative MIVQUE The restricted maximum likelihood estimator (REML) of Patterson and Thompson (1971) can be obtained by iterating on MIVQUE,
More informationContrasting Models for Lactation Curve Analysis
J. Dairy Sci. 85:968 975 American Dairy Science Association, 2002. Contrasting Models for Lactation Curve Analysis F. Jaffrezic,*, I. M. S. White,* R. Thompson, and P. M. Visscher* *Institute of Cell,
More informationEvaluation of Autoregressive Covariance Structures for Test-Day Records of Holstein Cows: Estimates of Parameters
J. Dairy Sci. 88:2632 2642 American Dairy Science Association, 2005. Evaluation of Autoregressive Covariance Structures for Test-Day Records of Holstein Cows: Estimates of Parameters R. M. Sawalha, 1 J.
More informationChapter 11 MIVQUE of Variances and Covariances
Chapter 11 MIVQUE of Variances and Covariances C R Henderson 1984 - Guelph The methods described in Chapter 10 for estimation of variances are quadratic, translation invariant, and unbiased For the balanced
More informationMaternal Genetic Models
Maternal Genetic Models In mammalian species of livestock such as beef cattle sheep or swine the female provides an environment for its offspring to survive and grow in terms of protection and nourishment
More informationComparison of computing properties of derivative and derivative-free algorithms in variance component estimation by REML.
April, 000 Comparison of computing properties of derivative and derivative-free algorithms in variance component estimation by REML By Ignacy Misztal University of Illinois, Urbana, Illinois 0, USA 0 0
More informationGenetic parameters for various random regression models to describe total sperm cells per ejaculate over the reproductive lifetime of boars
Published December 8, 2014 Genetic parameters for various random regression models to describe total sperm cells per ejaculate over the reproductive lifetime of boars S. H. Oh,* M. T. See,* 1 T. E. Long,
More informationModels with multiple random effects: Repeated Measures and Maternal effects
Models with multiple random effects: Repeated Measures and Maternal effects 1 Often there are several vectors of random effects Repeatability models Multiple measures Common family effects Cleaning up
More informationRECENT DEVELOPMENTS IN VARIANCE COMPONENT ESTIMATION
Libraries Conference on Applied Statistics in Agriculture 1989-1st Annual Conference Proceedings RECENT DEVELOPMENTS IN VARIANCE COMPONENT ESTIMATION R. R. Hocking Follow this and additional works at:
More informationPREDICTION OF BREEDING VALUES FOR UNMEASURED TRAITS FROM MEASURED TRAITS
Libraries Annual Conference on Applied Statistics in Agriculture 1994-6th Annual Conference Proceedings PREDICTION OF BREEDING VALUES FOR UNMEASURED TRAITS FROM MEASURED TRAITS Kristin L. Barkhouse L.
More informationMultiple Trait Evaluation of Bulls for Calving Ease
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Papers and Publications in Animal Science Animal Science Department February 1984 Multiple Trait Evaluation of Bulls
More informationLinear Models for the Prediction of Animal Breeding Values
Linear Models for the Prediction of Animal Breeding Values R.A. Mrode, PhD Animal Data Centre Fox Talbot House Greenways Business Park Bellinger Close Chippenham Wilts, UK CAB INTERNATIONAL Preface ix
More informationMixed-Models. version 30 October 2011
Mixed-Models version 30 October 2011 Mixed models Mixed models estimate a vector! of fixed effects and one (or more) vectors u of random effects Both fixed and random effects models always include a vector
More informationUnivariate and multivariate parameter estimates for milk production
. Roslin Original article Univariate and multivariate parameter estimates for milk production traits using an animal model. II. Efficiency of selection when using simplified covariance structures PM Visscher
More informationLecture 5: BLUP (Best Linear Unbiased Predictors) of genetic values. Bruce Walsh lecture notes Tucson Winter Institute 9-11 Jan 2013
Lecture 5: BLUP (Best Linear Unbiased Predictors) of genetic values Bruce Walsh lecture notes Tucson Winter Institute 9-11 Jan 013 1 Estimation of Var(A) and Breeding Values in General Pedigrees The classic
More informationMultiple random effects. Often there are several vectors of random effects. Covariance structure
Models with multiple random effects: Repeated Measures and Maternal effects Bruce Walsh lecture notes SISG -Mixed Model Course version 8 June 01 Multiple random effects y = X! + Za + Wu + e y is a n x
More informationShould genetic groups be fitted in BLUP evaluation? Practical answer for the French AI beef sire evaluation
Genet. Sel. Evol. 36 (2004) 325 345 325 c INRA, EDP Sciences, 2004 DOI: 10.1051/gse:2004004 Original article Should genetic groups be fitted in BLUP evaluation? Practical answer for the French AI beef
More informationImpact of Using Reduced Rank Random Regression Test-Day Model on Genetic Evaluation
Impact of Using Reduced Rank Random Regression Test-Day on Genetic Evaluation H. Leclerc 1, I. Nagy 2 and V. Ducrocq 2 1 Institut de l Elevage, Département Génétique, Bât 211, 78 352 Jouy-en-Josas, France
More informationGenetic parameters for female fertility in Nordic dairy cattle
Genetic parameters for female fertility in Nordic dairy cattle K.Muuttoranta 1, A-M. Tyrisevä 1, E.A. Mäntysaari 1, J.Pösö 2, G.P. Aamand 3, J-Å. Eriksson 4, U.S. Nielsen 5, and M. Lidauer 1 1 Natural
More informationAnimal Models. Sheep are scanned at maturity by ultrasound(us) to determine the amount of fat surrounding the muscle. A model (equation) might be
Animal Models 1 Introduction An animal model is one in which there are one or more observations per animal, and all factors affecting those observations are described including an animal additive genetic
More informationLongitudinal random effects models for genetic analysis of binary data with application to mastitis in dairy cattle
Genet. Sel. Evol. 35 (2003) 457 468 457 INRA, EDP Sciences, 2003 DOI: 10.1051/gse:2003034 Original article Longitudinal random effects models for genetic analysis of binary data with application to mastitis
More informationPrediction of Future Milk Yield with Random Regression Model Using Test-day Records in Holstein Cows
9 ` Asian-Aust. J. Anim. Sci. Vol. 19, No. 7 : 9-921 July 26 www.ajas.info Prediction of Future Milk Yield with Random Regression Model Using Test-day Records in Holstein Cows Byoungho Park and Deukhwan
More informationGenetic Heterogeneity of Environmental Variance - estimation of variance components using Double Hierarchical Generalized Linear Models
Genetic Heterogeneity of Environmental Variance - estimation of variance components using Double Hierarchical Generalized Linear Models L. Rönnegård,a,b, M. Felleki a,b, W.F. Fikse b and E. Strandberg
More informationAnimal Model. 2. The association of alleles from the two parents is assumed to be at random.
Animal Model 1 Introduction In animal genetics, measurements are taken on individual animals, and thus, the model of analysis should include the animal additive genetic effect. The remaining items in the
More informationEstimation of covariance components between
Original article Estimation of covariance components between one continuous and one binary trait H. Simianer L.R. Schaeffer 2 Justus Liebig University, Department of Animal Breeding and Genetics, Bismarckstr.
More informationLikelihood Methods. 1 Likelihood Functions. The multivariate normal distribution likelihood function is
Likelihood Methods 1 Likelihood Functions The multivariate normal distribution likelihood function is The log of the likelihood, say L 1 is Ly = π.5n V.5 exp.5y Xb V 1 y Xb. L 1 = 0.5[N lnπ + ln V +y Xb
More informationReduced Animal Models
Reduced Animal Models 1 Introduction In situations where many offspring can be generated from one mating as in fish poultry or swine or where only a few animals are retained for breeding the genetic evaluation
More informationGenetic Parameters for Stillbirth in the Netherlands
Genetic Parameters for Stillbirth in the Netherlands Arnold Harbers, Linda Segeren and Gerben de Jong CR Delta, P.O. Box 454, 68 AL Arnhem, The Netherlands Harbers.A@CR-Delta.nl 1. Introduction Stillbirth
More informationpopulation when only records from later
Original article Estimation of heritability in the base population when only records from later generations are available L Gomez-Raya LR Schaeffer EB Burnside University of Guelph, Centre for Genetic
More informationBest unbiased linear Prediction: Sire and Animal models
Best unbiased linear Prediction: Sire and Animal models Raphael Mrode Training in quantitative genetics and genomics 3 th May to th June 26 ILRI, Nairobi Partner Logo Partner Logo BLUP The MME of provided
More information3. Properties of the relationship matrix
3. Properties of the relationship matrix 3.1 Partitioning of the relationship matrix The additive relationship matrix, A, can be written as the product of a lower triangular matrix, T, a diagonal matrix,
More informationSolving Large Test-Day Models by Iteration on Data and Preconditioned Conjugate Gradient
Solving Large Test-Day Models by Iteration on Data and Preconditioned Conjugate Gradient M. LIDAUER, I. STRANDÉN, E. A. MÄNTYSAARI, J. PÖSÖ, and A. KETTUNEN Animal Production Research, Agricultural Research
More informationContemporary Groups for Genetic Evaluations
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Papers and Publications in Animal Science Animal Science Department January 1987 Contemporary Groups for Genetic
More informationSingle and multitrait estimates of breeding values for survival using sire and animal models
Animal Science 00, 75: 15-4 1357-798/0/11300015$0 00 00 British Society of Animal Science Single and multitrait estimates of breeding values for survival using sire and animal models T. H. E. Meuwissen
More informationBest linear unbiased prediction when error vector is correlated with other random vectors in the model
Best linear unbiased prediction when error vector is correlated with other random vectors in the model L.R. Schaeffer, C.R. Henderson To cite this version: L.R. Schaeffer, C.R. Henderson. Best linear unbiased
More informationRestricted maximum likelihood estimation for animal models using derivatives of the likelihood
Original article Restricted maximum likelihood estimation for animal models using derivatives of the likelihood K Meyer, SP Smith Animal Genetics and Breeding Unit, University of New England, Armidale,
More informationRANDOM REGRESSION IN ANIMAL BREEDING
RANDOM REGRESSION IN ANIMAL BREEDING Course Notes Jaboticabal, SP Brazil November 2001 Julius van der Werf University of New England Armidale, Australia 1 Introduction...2 2 Exploring correlation patterns
More informationASPECTS OF SELECTION FOR PERFORMANCE IN SEVERAL ENVIRONMENTS WITH HETEROGENEOUS VARIANCES
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Papers and Publications in Animal Science Animal Science Department 2-3-1987 ASPECTS OF SELECTION FOR PERFORMANCE
More informationChapter 5 Prediction of Random Variables
Chapter 5 Prediction of Random Variables C R Henderson 1984 - Guelph We have discussed estimation of β, regarded as fixed Now we shall consider a rather different problem, prediction of random variables,
More informationQuantitative characters - exercises
Quantitative characters - exercises 1. a) Calculate the genetic covariance between half sibs, expressed in the ij notation (Cockerham's notation), when up to loci are considered. b) Calculate the genetic
More informationLecture 7 Correlated Characters
Lecture 7 Correlated Characters Bruce Walsh. Sept 2007. Summer Institute on Statistical Genetics, Liège Genetic and Environmental Correlations Many characters are positively or negatively correlated at
More informationRepeated Records Animal Model
Repeated Records Animal Model 1 Introduction Animals are observed more than once for some traits, such as Fleece weight of sheep in different years. Calf records of a beef cow over time. Test day records
More informationPrediction of breeding values with additive animal models for crosses from 2 populations
Original article Prediction of breeding values with additive animal models for crosses from 2 populations RJC Cantet RL Fernando 2 1 Universidad de Buenos Aires, Departamento de Zootecnia, Facultad de
More informationA simple method of computing restricted best linear unbiased prediction of breeding values
Original article A simple method of computing restricted best linear unbiased prediction of breeding values Masahiro Satoh Department of Animal Breeding and Genetics, National Institute of Animal Industry,
More informationFull conjugate analysis of normal multiple traits with missing records using a generalized inverted Wishart distribution
Genet. Sel. Evol. 36 (2004) 49 64 49 c INRA, EDP Sciences, 2004 DOI: 10.1051/gse:2003050 Original article Full conjugate analysis of normal multiple traits with missing records using a generalized inverted
More informationof Nebraska - Lincoln
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Papers and Publications in Animal Science Animal Science Department January 1993 Sequential Transformation for Multiple
More informationEstimates of genetic parameters for total milk yield over multiple ages in Brazilian Murrah buffaloes using different models
Estimates of genetic parameters for total milk yield over multiple ages in Brazilian Murrah buffaloes using different models R.C. Sesana 1, F. Baldi 1, R.R.A. Borquis 1, A.B. Bignardi 1, N.A. Hurtado-Lugo
More informationChapter 19. Analysis of longitudinal data -Random Regression Analysis
Chapter 19 Analysis of longitudinal data -Random Regression Analysis Julius van der Werf 1 Introduction In univariate analysis the basic assumption is that a single measurement arises from a single unit
More informationHeterogeneity of variances by herd production level and its effect on dairy cow and sire evaluation
Retrospective Theses and Dissertations 1989 Heterogeneity of variances by herd production level and its effect on dairy cow and sire evaluation Keith George Boldman Iowa State University Follow this and
More informationProcedure 2 of Section 2 of ICAR Guidelines Computing of Accumulated Lactation Yield. Computing Lactation Yield
of ICAR Guidelines Computing of Accumulated Lactation Yield Table of Contents 1 The Test Interval Method (TIM) (Sargent, 1968)... 4 2 Interpolation using Standard Lactation Curves (ISLC) (Wilmink, 1987)...
More informationLecture 11: Multiple trait models for QTL analysis
Lecture 11: Multiple trait models for QTL analysis Julius van der Werf Multiple trait mapping of QTL...99 Increased power of QTL detection...99 Testing for linked QTL vs pleiotropic QTL...100 Multiple
More informationThe equivalence of the Maximum Likelihood and a modified Least Squares for a case of Generalized Linear Model
Applied and Computational Mathematics 2014; 3(5): 268-272 Published online November 10, 2014 (http://www.sciencepublishinggroup.com/j/acm) doi: 10.11648/j.acm.20140305.22 ISSN: 2328-5605 (Print); ISSN:
More informationA reduced animal model with elimination of quantitative trait loci equations for marker-assisted selection
Original article A reduced animal model with elimination of quantitative trait loci equations for marker-assisted selection S Saito H Iwaisaki 1 Graduate School of Science and Technology; 2 Department
More informationMultiple-Trait Across-Country Evaluations Using Singular (Co)Variance Matrix and Random Regression Model
Multiple-rait Across-Country Evaluations Using Singular (Co)Variance Matrix and Random Regression Model Esa A. Mäntysaari M Agrifood Research Finland, Animal Production, SF 31600 Jokioinen 1. Introduction
More informationMIVQUE and Maximum Likelihood Estimation for Multivariate Linear Models with Incomplete Observations
Sankhyā : The Indian Journal of Statistics 2006, Volume 68, Part 3, pp. 409-435 c 2006, Indian Statistical Institute MIVQUE and Maximum Likelihood Estimation for Multivariate Linear Models with Incomplete
More informationSimulation Study on Heterogeneous Variance Adjustment for Observations with Different Measurement Error Variance
Simulation Study on Heterogeneous Variance Adjustment for Observations with Different Measurement Error Variance Pitkänen, T. 1, Mäntysaari, E. A. 1, Nielsen, U. S., Aamand, G. P 3., Madsen 4, P. and Lidauer,
More informationCovariance functions and random regression models for cow weight in beef cattle
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Papers and Publications in Animal Science Animal Science Department January 2004 Covariance functions and random
More informationREDUCED ANIMAL MODEL FOR MARKER ASSISTED SELECTION USING BEST LINEAR UNBIASED PREDICTION. R.J.C.Cantet 1 and C.Smith
REDUCED ANIMAL MODEL FOR MARKER ASSISTED SELECTION USING BEST LINEAR UNBIASED PREDICTION R.J.C.Cantet 1 and C.Smith Centre for Genetic Improvement of Livestock, Department of Animal and Poultry Science,
More informationREML Variance-Component Estimation
REML Variance-Component Estimation In the numerous forms of analysis of variance (ANOVA) discussed in previous chapters, variance components were estimated by equating observed mean squares to expressions
More informationEstimation of genetic parameters of egg production traits of laying hens by restricted maximum likelihood applied
Original article Estimation of genetic parameters of egg production traits of laying hens by restricted maximum likelihood applied to a multiple-trait reduced animal model B Besbes 4 V Ducrocq JL Foulley
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Genotype by Environment Interaction and Genetic Correlations Among Parities for Somatic Cell Count and Milk Yield Citation for published version: BANOS, G & SHOOK, GE 1990,
More informationBLUP without (inverse) relationship matrix
Proceedings of the World Congress on Genetics Applied to Livestock Production, 11, 5 BLUP without (inverse relationship matrix E. Groeneveld (1 and A. Neumaier ( (1 Institute of Farm Animal Genetics, Friedrich-Loeffler-Institut,
More informationHeritability, Reliability of Genetic Evaluations and Response to Selection in Proportional Hazard Models
J. Dairy Sci. 85:1563 1577 American Dairy Science Association, 00. Heritability, Reliability of Genetic Evaluations and Response to Selection in Proportional Hazard Models M. H. Yazdi,* P. M. Visscher,*
More informationIMPLEMENTATION ISSUES IN BAYESIAN ANALYSIS IN ANIMAL BREEDING. C.S. Wang. Central Research, Pfizer Inc., Groton, CT 06385, USA
IMPLEMENTATION ISSUES IN BAYESIAN ANALYSIS IN ANIMAL BREEDING C.S. Wang Central Research, Pfizer Inc., Groton, CT 06385, USA SUMMARY Contributions of Bayesian methods to animal breeding are reviewed briefly.
More informationBayes factor for testing between different structures of random genetic groups: A case study using weaning weight in Bruna dels Pirineus beef cattle
Genet. Sel. Evol. 39 (007) 39 53 39 c INRA, EDP Sciences, 006 DOI: 10.1051/gse:006030 Original article Bayes factor for testing between different structures of random genetic groups: A case study using
More informationEM-REML estimation of covariance parameters in Gaussian mixed models for longitudinal data analysis
Genet. Sel. Evol. 32 (2000) 129 141 129 c INRA, EDP Sciences EM-REML estimation of covariance parameters in Gaussian mixed models for longitudinal data analysis Original article Jean-Louis FOULLEY a, Florence
More information2.2 Selection on a Single & Multiple Traits. Stevan J. Arnold Department of Integrative Biology Oregon State University
2.2 Selection on a Single & Multiple Traits Stevan J. Arnold Department of Integrative Biology Oregon State University Thesis Selection changes trait distributions. The contrast between distributions before
More informationDistinctive aspects of non-parametric fitting
5. Introduction to nonparametric curve fitting: Loess, kernel regression, reproducing kernel methods, neural networks Distinctive aspects of non-parametric fitting Objectives: investigate patterns free
More informationAn indirect approach to the extensive calculation of relationship coefficients
Genet. Sel. Evol. 34 (2002) 409 421 409 INRA, EDP Sciences, 2002 DOI: 10.1051/gse:2002015 Original article An indirect approach to the extensive calculation of relationship coefficients Jean-Jacques COLLEAU
More informationECM approaches to heteroskedastic mixed models with constant variance ratios
Original article ECM approaches to heteroskedastic mixed models with constant variance ratios JL Foulley Station de génétique quantitative et appliquée, Institut national de la recherche agronomique, 78352
More information5. Best Linear Unbiased Prediction
5. Best Linear Unbiased Prediction Julius van der Werf Lecture 1: Best linear unbiased prediction Learning objectives On completion of Lecture 1 you should be able to: Understand the principle of mixed
More informationLecture 32: Infinite-dimensional/Functionvalued. Functions and Random Regressions. Bruce Walsh lecture notes Synbreed course version 11 July 2013
Lecture 32: Infinite-dimensional/Functionvalued Traits: Covariance Functions and Random Regressions Bruce Walsh lecture notes Synbreed course version 11 July 2013 1 Longitudinal traits Many classic quantitative
More informationGENERALIZED LINEAR MIXED MODELS: AN APPLICATION
Libraries Conference on Applied Statistics in Agriculture 1994-6th Annual Conference Proceedings GENERALIZED LINEAR MIXED MODELS: AN APPLICATION Stephen D. Kachman Walter W. Stroup Follow this and additional
More informationLikelihood calculations to evaluate experimental designs to estimate genetic variances
Running head : LIKELIHOOD EXPERIMENTAL DESIGN Likelihood calculations to evaluate experimental designs to estimate genetic variances Karin Meyer Animal Genetics and Breeding Unit 1, University of New England,
More informationThe use of independent culling levels and selection index procedures in selecting future sires for artificial insemination
Retrospective Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 1971 The use of independent culling levels and selection index procedures in selecting future sires for
More informationSelection on selected records
Selection on selected records B. GOFFINET I.N.R.A., Laboratoire de Biometrie, Centre de Recherches de Toulouse, chemin de Borde-Rouge, F 31320 Castanet- Tolosan Summary. The problem of selecting individuals
More informationG-BLUP without inverting the genomic relationship matrix
G-BLUP without inverting the genomic relationship matrix Per Madsen 1 and Jørgen Ødegård 2 1 Center for Quantitative Genetics and Genomics Department of Molecular Biology and Genetics, Aarhus University
More informationGenetic grouping for direct and maternal effects with differential assignment of groups
Original article Genetic grouping for direct and maternal effects with differential assignment of groups RJC Cantet RL Fernando 2 D Gianola 3 I Misztal 1Facultad de Agronomia, Universidad de Buenos Aires,
More informationEstimation in Generalized Linear Models with Heterogeneous Random Effects. Woncheol Jang Johan Lim. May 19, 2004
Estimation in Generalized Linear Models with Heterogeneous Random Effects Woncheol Jang Johan Lim May 19, 2004 Abstract The penalized quasi-likelihood (PQL) approach is the most common estimation procedure
More informationMULTIVARIATE ESTIMATION OF GENETIC PARAMETERS QUO VADIS? Karin Meyer
Proc. Assoc. Advmt. Anim. Breed. Genet. 19:71 78 MULTIVARIATE ESTIMATION OF GENETIC PARAMETERS QUO VADIS? Karin Meyer Animal Genetics and Breeding Unit *, University of New England, Armidale, NSW 2351
More informationThis book is dedicated to Professor Dr G. K. Constantinescu, founder of the modern animal husbandry science in Romania, originator of the National
This book is dedicated to Professor Dr G. K. Constantinescu, founder of the modern animal husbandry science in Romania, originator of the National Animal Husbandry Institute (1926) and initiator of the
More informationGenetic variation of traits measured in several environments. II. Inference
Original article Genetic variation of traits measured in several environments. II. Inference on between-environment homogeneity of intra-class correlations C Robert JL Foulley V Ducrocq Institut national
More informationBias and sampling covariances of estimates of variance components due to maternal effects
Original article Bias and sampling covariances of estimates of variance components due to maternal effects K Meyer Edinburgh University, Institute for Cell, Animal and Population Biology, West Mains Road,
More informationApplication of Householder's transformations and the QL algorithm to REML estimation of variance components
Application of Householder's transformations and the QL algorithm to REML estimation of variance components K.V. KonstantinClv and G.J. Erasmus Department of Animal Science, University of the Orange Free
More informationA link function approach to heterogeneous variance components
Original article A link function approach to heterogeneous variance components Jean-Louis Foulley Richard L. Quaas b Thaon d Arnoldi a Station de génétique quantitative et appliquée, Institut national
More informationROBUST ESTIMATION AND OUTLIER DETECTION FOR REPEATED MEASURES EXPERIMENTS. R.M. HUGGINS
77 ROBUST ESTIMATION AND OUTLIER DETECTION FOR REPEATED MEASURES EXPERIMENTS. 1. INTRODUCTION. R.M. HUGGINS Standard methods of testing for repeated measures experiments and other mixed linear models are
More informationLecture 5 Basic Designs for Estimation of Genetic Parameters
Lecture 5 Basic Designs for Estimation of Genetic Parameters Bruce Walsh. jbwalsh@u.arizona.edu. University of Arizona. Notes from a short course taught June 006 at University of Aarhus The notes for this
More information