Differences

This shows you the differences between two versions of the page.

--- readme.reml [2012/05/28 14:11] – created shogo
+++ readme.reml [2012/05/29 17:06] – shogo
@@ Line 1: / Line 1: @@
-REMLF90 - Additions to PDF documentation
+====== REMLF90 ======
-IM 02/17/2003
+Variance component estimation program with EM-REML algorithm.
-Two new random types were added. Type:
+===== Summary =====
-user_file
-accepts a matrix from a file specified under keyword FILE in the following
-format:
-row col value
-row col value
-.....
-Values can be lower- or upper-stored but not both.
-Second type:
+===== Options =====
-user_file_inv
+<file>
-works as befoer except that the given matrix is inverted before being
+OPTION conv_crit1d-12
-used in mixed model equations.
+</file>
+Convergence criterion (default 1d-10).
+<file>
+OPTION maxrounds500
+</file>
+Maximum rounds (default 5000).
+<file>
+OPTION constant_var 5 1 2
+</file>
+: effect number\\
+: first trait number\\
+: and the second trait number\\
+implying the covariance between traits 1 and 2 for effect 5.
+<file>
+OPTION SNP_file snp
+</file>
+Specify the SNP file name to use genotype data.
+===== 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. If airemlf90 does not converge but remlf90 converges with the same dataset, rerun remlf90 with a small starting value. The estimate could be artifact.