An Approximate Procedure for Determining Solutions and Prediction Errors for Random Effects in Mixed Model Equations with a Large Numerator Relationship Matrix

1985 ◽  
Vol 56 (8) ◽  
pp. 685-689
Author(s):  
Kenji TOGASHI ◽  
Hisato TAKEDA
2021 ◽  
Vol 12 ◽  
Author(s):  
Mohammad Ali Nilforooshan ◽  
Dorian Garrick

Reduced models are equivalent models to the full model that enable reduction in the computational demand for solving the problem, here, mixed model equations for estimating breeding values of selection candidates. Since phenotyped animals provide data to the model, the aim of this study was to reduce animal models to those equations corresponding to phenotyped animals. Non-phenotyped ancestral animals have normally been included in analyses as they facilitate formation of the inverse numerator relationship matrix. However, a reduced model can exclude those animals and obtain identical solutions for the breeding values of the animals of interest. Solutions corresponding to non-phenotyped animals can be back-solved from the solutions of phenotyped animals and specific blocks of the inverted relationship matrix. This idea was extended to other forms of animal model and the results from each reduced model (and back-solving) were identical to the results from the corresponding full model. Previous studies have been mainly focused on reduced animal models that absorb equations corresponding to non-parents and solve equations only for parents of phenotyped animals. These two types of reduced animal model can be combined to formulate only equations corresponding to phenotyped parents of phenotyped progeny.


1985 ◽  
Vol 36 (3) ◽  
pp. 527 ◽  
Author(s):  
H-U Graser ◽  
K Hammond

A multiple-trait mixed model is defined for regular use in the Australian beef industry for the estimation of breeding values for continuous traits of sires used non-randomly across a number of herds and/or years. Maternal grandsires, the numerator relationship matrix, appropriate fixed effects, and the capacity to partition direct and maternal effects are incorporated in this parent model. The model was fitted to the National Beef Recording Scheme's data bank for three growth traits of the Australian Simental breed, viz 200-, 365- and 550-day weights. Estimates are obtained for the effects of sex, dam age, grade of dam, age of calf and breed of base dam. The range in estimated breeding value is reported for each trait, with 200-day weight being partitioned into 'calves' and 'daughters' calves', for the Simmental sires commonly used in Australia. Estimates of the fixed effects were large, and dam age, grade of dam and breed of base dam had an important influence on growth to 365 days of age. The faster growth of higher percentage Simmental calves to 200 days continued to 550 days. Estimates of genetic variance for the traits were lower than reported for overseas populations of Simmental cattle, and the genetic covariance between direct and maternal effects for 200-day weight was slightly positive.


2004 ◽  
Vol 3 (1) ◽  
pp. 1-21 ◽  
Author(s):  
R. Mark Thallman ◽  
Kathryn J Hanford ◽  
Stephen D Kachman ◽  
L. Dale Van Vleck

Gametic models for fitting breeding values at QTL as random effects in outbred populations have become popular because they require few assumptions about the number and distribution of QTL alleles segregating. The covariance matrix of the gametic effects has an inverse that is sparse and can be constructed rapidly by a simple algorithm, provided that all individuals have marker data, but not otherwise. An equivalent model, in which the joint distribution of QTL breeding values and marker genotypes is considered, was shown to generate a covariance matrix with a sparse inverse that can be constructed rapidly with a simple algorithm. This result makes more feasible including QTL as random effects in analyses of large pedigrees for QTL detection and marker assisted selection. Such analyses often use algorithms that rely upon sparseness of the mixed model equations and require the inverse of the covariance matrix, but not the covariance matrix itself. With the proposed model, each individual has two random effects for each possible unordered marker genotype for that individual. Therefore, individuals with marker data have two random effects, just as with the gametic model. To keep the notation and the derivation simple, the method is derived under the assumptions of a single linked marker and that the pedigree does not contain loops. The algorithm could be applied, as an approximate method, to pedigrees that contain loops.


Genetics ◽  
1997 ◽  
Vol 145 (4) ◽  
pp. 1243-1249 ◽  
Author(s):  
Piter Bijma ◽  
Johan A M Van Arendonk ◽  
Henk Bovenhuis

Under gynogenetic reproduction, offspring receive genes only from their dams and completely homozygous offspring are produced within one generation. When gynogenetic reproduction is applied to fully inbred individuals, homozygous clone lines are produced. A mixed model method was developed for breeding value and variance component estimation in gynogenetic families, which requires the inverse of the numerator relationship matrix. A general method for creating the inverse for a population with unusual relationships between animals is presented, which reduces to simple rules as is illustrated for gynogenetic populations. The presence of clones in gynogenetic populations causes singularity of the numerator relationship matrix. However, clones can be regarded as repeated observations of the same genotype, which can be accommodated by modifying the incidence matrix, and by considering only unique genotypes in the estimation procedure. Optimum gynogenetic sib family sizes for estimating heritabilities and estimates of their accuracy were derived and compared to those for conventional full-sib designs. This was done by means of a deterministic derivation and by stochastic simulation using Gibbs sampling. Optimum family sizes were smallest for gynogenetic families. Only for low heritabilities, there was a small advantage in accuracy under the gynogenetic design.


Forests ◽  
2020 ◽  
Vol 11 (11) ◽  
pp. 1169
Author(s):  
Gary R. Hodge ◽  
Juan Jose Acosta

Research Highlights: An algorithm is presented that allows for the analysis of full-sib genetic datasets using generalized mixed-model software programs. The algorithm produces variance component estimates, genetic parameter estimates, and Best Linear Unbiased Prediction (BLUP) solutions for genetic values that are, for all practical purposes, identical to those produced by dedicated genetic software packages. Background and Objectives: The objective of this manuscript is to demonstrate an approach with a simulated full-sib dataset representing a typical forest tree breeding population (40 parents, 80 full-sib crosses, 4 tests, and 6000 trees) using two widely available mixed-model packages. Materials and Methods: The algorithm involves artificially doubling the dataset, so that each observation is in the dataset twice, once with the original female and male parent identification, and once with the female and male parent identities switched. Five linear models were examined: two models using a dedicated genetic software program (ASREML) with the capacity to specify A or other pedigree-related functions, and three models with the doubled dataset and a parent (or sire) linear model (ASREML, SAS Proc Mixed, and R lme4). Results: The variance components, genetic parameters, and BLUPs of the parental breeding values, progeny breeding values, and full-sib family-specific combining abilities were compared. Genetic parameter estimates were essentially the same across all the analyses (e.g., the heritability ranged from h2 = 0.220 to 0.223, and the proportion of dominance variance ranged from d2 = 0.057 to 0.058). The correlations between the BLUPs from the baseline analysis (ASREML with an individual tree model) and the doubled-dataset/parent models using SAS Proc Mixed or R lme4 were never lower than R = 0.99997. Conclusions: The algorithm can be useful for analysts who need to analyze full-sib genetic datasets and who are familiar with general-purpose statistical packages, but less familiar with or lacking access to other software.


2020 ◽  
pp. 1-37
Author(s):  
Tal Yarkoni

Abstract Most theories and hypotheses in psychology are verbal in nature, yet their evaluation overwhelmingly relies on inferential statistical procedures. The validity of the move from qualitative to quantitative analysis depends on the verbal and statistical expressions of a hypothesis being closely aligned—that is, that the two must refer to roughly the same set of hypothetical observations. Here I argue that many applications of statistical inference in psychology fail to meet this basic condition. Focusing on the most widely used class of model in psychology—the linear mixed model—I explore the consequences of failing to statistically operationalize verbal hypotheses in a way that respects researchers' actual generalization intentions. I demonstrate that whereas the "random effect" formalism is used pervasively in psychology to model inter-subject variability, few researchers accord the same treatment to other variables they clearly intend to generalize over (e.g., stimuli, tasks, or research sites). The under-specification of random effects imposes far stronger constraints on the generalizability of results than most researchers appreciate. Ignoring these constraints can dramatically inflate false positive rates, and often leads researchers to draw sweeping verbal generalizations that lack a meaningful connection to the statistical quantities they are putatively based on. I argue that failure to take the alignment between verbal and statistical expressions seriously lies at the heart of many of psychology's ongoing problems (e.g., the replication crisis), and conclude with a discussion of several potential avenues for improvement.


Sign in / Sign up

Export Citation Format

Share Document