readme.aireml
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
readme.aireml [2012/05/30 13:30] – shogo | readme.aireml [2014/06/04 13:07] – [Options] ignacio | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== AIREMLF90 ====== | ====== AIREMLF90 ====== | ||
- | \\ | ||
- | A modification of REMLF90 for estimating variances with the Average-Information algorithm. Initially written by Shogo Tsuruta in 03/ | ||
===== Summary ===== | ===== Summary ===== | ||
- | AIREMLF90 uses a second derivative REML algrithm with extra heuristics, as is described in Jensen et al. (1996-7). For most problems, it converges in far fewer rounds than EM REML as implemented in REMLF90. While typically REMLF90 takes 50-300 rounds to converge, AIREMLF90 converges in 5-15 rounds and to a higher accuracy. For selected | + | A modification of REMLF90 for estimating variances with the Average-Information algorithm. Initially written by Shogo Tsuruta in 03/ |
+ | \\ | ||
+ | AIREMLF90 uses a second derivative REML algrithm with extra heuristics, as is described in Jensen et al. (1996-7). For most models, it converges in far fewer rounds than EM-REML as implemented in REMLF90. While typically REMLF90 takes 50-300 rounds to converge, AIREMLF90 converges in 5-15 rounds and to a higher accuracy. For selected | ||
+ | \\ | ||
+ | See PREGSF90 with genotypes (SNP) for options. | ||
===== Options ===== | ===== Options ===== | ||
Line 19: | Line 21: | ||
</ | </ | ||
Run EM-REML (REMLF90) for first 10 rounds to get initial variances within the parameter space (default 0). | Run EM-REML (REMLF90) for first 10 rounds to get initial variances within the parameter space (default 0). | ||
+ | < | ||
+ | OPTION use_yams | ||
+ | </ | ||
+ | Run the program with YAMS (modified FSPAK). The computing time can be dramatically improved. | ||
< | < | ||
OPTION tol 1d-12 | OPTION tol 1d-12 | ||
</ | </ | ||
- | Tolerance (or precision) (default 1d-14) for positive definite matrix and g-inverse subroutines. Convergence may be much faster by changing this value. | + | Tolerance (or precision) (default 1d-14) for positive definite matrix and g-inverse subroutines.\\ |
+ | Convergence may be much faster by changing this value. | ||
< | < | ||
OPTION sol se | OPTION sol se | ||
Line 28: | Line 35: | ||
Store solutions and those s.e. | Store solutions and those s.e. | ||
< | < | ||
- | OPTION | + | OPTION |
</ | </ | ||
- | Set the missing value (default 0). | + | y-hat and residuals will be included in " |
+ | < | ||
+ | OPTION missing -999 | ||
+ | </ | ||
+ | Specify | ||
**Heterogeneous residual variances for a single trait** | **Heterogeneous residual variances for a single trait** | ||
Line 42: | Line 53: | ||
Initial values of coefficients for heterogeneous residual variances using //ln//(a0, a1, a2, ...) to make these values. | Initial values of coefficients for heterogeneous residual variances using //ln//(a0, a1, a2, ...) to make these values. | ||
- | **Heterogeneous residual variances for multiple traits** | + | **Heterogeneous residual variances for multiple traits**\\ |
Convergence will be very slow with multiple trait heterogeneous residual variances | Convergence will be very slow with multiple trait heterogeneous residual variances | ||
< | < | ||
Line 51: | Line 62: | ||
OPTION hetres_pos 10 11 12 13 | OPTION hetres_pos 10 11 12 13 | ||
</ | </ | ||
- | Specify the position of covariables (trait first). "10 10" or "10 11" could be linear for first and second traits, and "11 11" or "12 13" could be quadratic. | + | Specify the position of covariables (trait first). |
+ | "10 10" or "10 11" could be linear for first and second traits.\\ | ||
+ | "11 11" or "12 13" could be quadratic. | ||
< | < | ||
OPTION hetres_pol 4.0 4.0 0.1 0.1 0.01 0.01 | OPTION hetres_pol 4.0 4.0 0.1 0.1 0.01 0.01 | ||
</ | </ | ||
- | Initial values of coefficients for heterogeneous residual variances using //ln//(a0, a1, a2, ...) to make these values (trait first). "4.0 4.0" are intercept for first and second traits. "0.1 0.1" could be linear and "0.01 0.01" could be quadratic. To transform back to the original scale, use exp(a0+a1*X1+a2*X2). | + | Initial values of coefficients for heterogeneous residual variances using //ln//(a0, a1, a2, ...) to make these values (trait first).\\ |
+ | "4.0 4.0" are intercept for first and second traits.\\ | ||
+ | "0.1 0.1" could be linear and "0.01 0.01" could be quadratic.\\ | ||
+ | To transform back to the original scale, use exp(a0+a1*X1+a2*X2). | ||
+ | < | ||
+ | OPTION SNP_file snp | ||
+ | </ | ||
+ | Specify the SNP file name to use genotype data. | ||
+ | |||
+ | < | ||
+ | < | ||
+ | Set the number of samples to calculate SE for function of (co)variances.\\ | ||
+ | default value 5000 | ||
+ | < | ||
readme.aireml.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1