readme.reml
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.reml [2013/07/17 16:19] – shogo | readme.reml [2014/11/25 11:16] – shogo | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== REMLF90 ====== | ====== REMLF90 ====== | ||
- | Variance component estimation program with EM-REML algorithm. | ||
- | See PREGSF90 with genotypes (SNP) for options. | ||
===== Summary ===== | ===== Summary ===== | ||
- | The final results will be saved in " | + | Variance component estimation program with EM-REML algorithm. |
+ | \\ | ||
+ | See PREGSF90 with genotypes (SNP) for options. | ||
===== Options ===== | ===== Options ===== | ||
Line 15: | Line 15: | ||
</ | </ | ||
Maximum rounds (default 5000). | Maximum rounds (default 5000). | ||
+ | < | ||
+ | OPTION sol se | ||
+ | </ | ||
+ | Store solutions and se. | ||
+ | < | ||
+ | OPTION residual | ||
+ | </ | ||
+ | y-hat and residuals will be included in " | ||
+ | < | ||
+ | OPTION missing -999 | ||
+ | </ | ||
+ | Specify missing observations (default 0). | ||
+ | < | ||
+ | OPTION use_yams | ||
+ | </ | ||
+ | Run the program with YAMS (modified FSPAK). The computing time can be dramatically improved. | ||
< | < | ||
OPTION constant_var 5 1 2 | OPTION constant_var 5 1 2 | ||
Line 28: | Line 44: | ||
===== Does remlf90 always converge? ===== | ===== Does remlf90 always converge? ===== | ||
- | Even when there is no variance, remlf90 will estimate a positive | + | When the expected variance |
+ | The remlf90 program will always estimate a positive variance (within the parameter space) even if there is no variance (using wrong data or format and/or a completely wrong model). | ||
+ | If airemlf90 does not converge but remlf90 converges with the same data set and the same model, rerun remlf90 with a small starting value to check the estimate because it could be artifact. |
readme.reml.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1