readme.aireml
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
readme.aireml [2014/09/05 16:10] – shogo | readme.aireml [2024/03/25 18:22] (current) – external edit 127.0.0.1 | ||
---|---|---|---|
Line 2: | Line 2: | ||
===== Summary ===== | ===== Summary ===== | ||
- | A modification of REMLF90 for estimating variances with the Average-Information algorithm. Initially written by Shogo Tsuruta in 03/ | + | A modification of REMLF90 for estimating variances with the Average-Information algorithm. Initially written by Shogo Tsuruta in 03/99-07/99. AIREMLF90 uses a second derivative REML algorithm 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. The final results will be saved in " |
\\ | \\ | ||
- | 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 models, AI-REML may fail to converge when the covariance matrix is close to non-positive definite. Adjust sensitivity of the program by setting the appropriate tolerance or setting good starting values. The final results will be saved in " | ||
\\ | \\ | ||
See PREGSF90 with genotypes (SNP) for options. | See PREGSF90 with genotypes (SNP) for options. | ||
Line 10: | Line 9: | ||
===== Options ===== | ===== Options ===== | ||
< | < | ||
- | OPTION conv_crit 1d-12 | + | OPTION conv_crit 1d-10 |
</ | </ | ||
- | Convergence criterion (default 1d-10). | + | Convergence criterion (default 1d-12). |
< | < | ||
- | OPTION maxrounds | + | OPTION maxrounds |
</ | </ | ||
- | Maximum rounds (default 5000). When the number < 2, the program calculates BLUP without iterating REML. | + | Maximum rounds (default 5000). When n = 0, the program calculates BLUP without iterating REML and provides some statistics (-2logL, AIC, SE for (co)variances, |
< | < | ||
- | OPTION EM-REML | + | OPTION EM-REML |
</ | </ | ||
- | Run EM-REML (REMLF90) for first 10 rounds to get initial variances within the parameter space (default 0). | + | Run EM-REML (REMLF90) for first n rounds to get initial |
< | < | ||
OPTION use_yams | OPTION use_yams | ||
Line 45: | Line 44: | ||
OPTION missing -999 | OPTION missing -999 | ||
</ | </ | ||
- | Specify the missing value (default 0). | + | Specify the missing value (default 0) in integer. |
< | < | ||
OPTION constant_var 5 1 2 ... | OPTION constant_var 5 1 2 ... | ||
Line 58: | Line 57: | ||
OPTION hetres_pos 10 11 | OPTION hetres_pos 10 11 | ||
</ | </ | ||
- | Specify the position | + | Specify the column positions |
< | < | ||
OPTION hetres_pol 4.0 0.1 0.1 | OPTION hetres_pol 4.0 0.1 0.1 | ||
Line 73: | Line 72: | ||
OPTION hetres_pos 10 11 12 13 | OPTION hetres_pos 10 11 12 13 | ||
</ | </ | ||
- | Specify the position | + | Specify the column positions |
"10 10" or "10 11" could be linear for first and second traits.\\ | "10 10" or "10 11" could be linear for first and second traits.\\ | ||
"11 11" or "12 13" could be quadratic. | "11 11" or "12 13" could be quadratic. | ||
Line 89: | Line 88: | ||
< | < | ||
- | As an alternative of SE, calculate SD for function of (co)variances by repeated sampling of parameters estimates from their asymptotic multivariate normal distribution, | + | As an alternative of SE, calculate SD for function of (co)variances by repeated sampling of parameters estimates from their asymptotic multivariate normal distribution, |
\\ | \\ | ||
''< | ''< | ||
Line 97: | Line 96: | ||
A formula to calculate a function of (co)variances to estimate SD. All terms of the function should be written with no spaces.\\ | A formula to calculate a function of (co)variances to estimate SD. All terms of the function should be written with no spaces.\\ | ||
\\ | \\ | ||
- | Each term of the function corresponds to (co)variance elements and could include any random effects (G, PE, ...) and residual (R) (co)variances.\\ | + | Each term of the function corresponds to (co)variance elements and could include any random effects (G) and residual (R) (co)variances.\\ |
\\ | \\ | ||
Notation is with reference to the effect number and the trait number ('' | Notation is with reference to the effect number and the trait number ('' | ||
Line 125: | Line 124: | ||
< | < | ||
Set the number of samples to calculate SE for function of (co)variances.\\ | Set the number of samples to calculate SE for function of (co)variances.\\ | ||
- | default value 5000 | + | default value 10000 |
< | < | ||
Indicate to store in file samples of (co)variances function for postprocessing (histogram, etc.) | Indicate to store in file samples of (co)variances function for postprocessing (histogram, etc.) | ||
+ | |||
+ | ===== Tricks ===== | ||
+ | When the covariance matrix is close to non-positive definite, the AIREMLF90 may not converge. | ||
+ | There are two options you might want to try: | ||
+ | |||
+ | 1. change the tolerance value (xx) in the option: | ||
+ | |||
+ | OPTION tol xx | ||
+ | |||
+ | to a very strict value (e.g., 1d-20) or a lenient value (1d-06). | ||
+ | |||
+ | 2. use an option to use EM-REML inside AI-REML: | ||
+ | |||
+ | OPTION EM-REML xx | ||
+ | |||
+ | where xx is the number of iterations for EM-REML you expect to get a good starting value for AI-REML. After running xx rounds with EM-REML, the AIREMLF90 program will automatically switch from EM-REML to AI-REML using the last estimate from EM-REML as a starting value for AI-REML. | ||
+ |
readme.aireml.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1