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 [2012/05/29 14:30] – shogo | readme.reml [2012/05/29 14:47] – [Does remlf90 always converge?] ignacio | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== REMLF90 ====== | ====== REMLF90 ====== | ||
+ | Variance component estimation program with EM-REML algorithm. | ||
===== Summary ===== | ===== Summary ===== | ||
+ | |||
===== Options ===== | ===== Options ===== | ||
Line 15: | Line 17: | ||
OPTION constant_var 5 1 2 | OPTION constant_var 5 1 2 | ||
</ | </ | ||
- | | + | |
+ | 1: first trait number | ||
+ | 2: and the second trait number | ||
+ | | ||
< | < | ||
OPTION SNP_file snp | OPTION SNP_file snp | ||
Line 22: | Line 27: | ||
===== Does remlf90 always converge? ===== | ===== Does remlf90 always converge? ===== | ||
- | < | + | |
- | Even when there is no variance, remlf90 will estimate a positive variance (in the parameter space), which will be determined by the starting value. If you are not sure if there is any variance, use a small starting value like 0.1 or 0.01. If the estimate does not change, probably there is no variance, so check your parameter file, model, and/or data again. | + | Even when there is no variance, remlf90 will estimate a positive variance (in the parameter space), which will be determined by the starting value. |
- | </ | + | If you are not sure if there is any variance, use a small starting value like 0.1 or 0.01. |
+ | If the estimate does not change, probably there is no variance, so check your parameter file, model, and/or data again. | ||
+ | If airemlf90 does not converge but remlf90 converges with the same dataset, rerun remlf90 with a small starting value. | ||
+ | The estimate could be artifact. | ||
readme.reml.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1