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 [2014/11/25 11:18] – shogo | readme.aireml [2014/11/25 11:29] – shogo | ||
---|---|---|---|
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/ |
+ | \\ | ||
\\ | \\ | ||
See PREGSF90 with genotypes (SNP) for options. | See PREGSF90 with genotypes (SNP) for options. | ||
Line 128: | Line 129: | ||
===== Tricks ===== | ===== Tricks ===== | ||
- | 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 | + | When the covariance matrix |
+ | 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 | ||
+ | |||
+ | OPTION EM-REML xx | ||
+ | |||
+ | where xx is the number of iterations for EM-REML you expect to get a good starting | ||
readme.aireml.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1