how_to_run_pure_gblup
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revision | Next revisionBoth sides next revision | ||
how_to_run_pure_gblup [2019/06/07 00:40] – yutaka | how_to_run_pure_gblup [2019/06/07 01:34] – yutaka | ||
---|---|---|---|
Line 15: | Line 15: | ||
$$ | $$ | ||
- | where $\mathbf{y}$ is a vector of observations, | + | where $\mathbf{y}$ is a vector of observations, $\mathbf{X}$ is a design matrix relating the fixed effects to the observations (typically $\mathbf{1}$), $\mathbf{b}$ is a vector of fixed effects (typically a single $\mu$), $\mathbf{u}$ is a vector of additive genetic effects, $\mathbf{e}$ is a vector of residual effects. We assume $\mathrm{var}(\mathbf{y})=\mathrm{var}(\mathbf{u})+\mathrm{var}(\mathbf{e})=\mathbf{G}\sigma_{u}^{2}+\mathbf{R}\sigma_{e}^{2}$ where $\mathbf{G}$ is the genomic relationship matrix, $\mathbf{R}$ is a diagonal matrix (typically $\mathbf{I}$), |
$$ | $$ | ||
Line 38: | Line 38: | ||
\right] | \right] | ||
$$ | $$ | ||
- | Note that $\mathbf{R}^{-1}=\mathbf{I}/ | + | Note that $\mathbf{R}^{-1}=\mathbf{I}/ |
===== A way to build MME in blupf90 ===== | ===== A way to build MME in blupf90 ===== | ||
Line 44: | Line 44: | ||
The blupf90 program is designed to solve the animal model that has the inverse of the numerator relationship matrix ($\mathbf{A}^{-1}$). Also, this program is extended to perform the single-step GBLUP analysis including the inverse of a subset pedigree matrix for genotyped animals ($\mathbf{A}_{22}^{-1}$) as well as $\mathbf{G}^{-1}$. | The blupf90 program is designed to solve the animal model that has the inverse of the numerator relationship matrix ($\mathbf{A}^{-1}$). Also, this program is extended to perform the single-step GBLUP analysis including the inverse of a subset pedigree matrix for genotyped animals ($\mathbf{A}_{22}^{-1}$) as well as $\mathbf{G}^{-1}$. | ||
- | Regardless of the above GBLUP model, the program literally creates the following MME. | + | Regardless of whether you need the pedigree relationships, the program literally creates the following MME. |
$$ | $$ |
how_to_run_pure_gblup.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1