Differences

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

--- readme.aireml [2014/09/05 16:16] – shogo
+++ readme.aireml [2014/11/25 11:18] – shogo
@@ Line 3: / Line 3: @@
 ===== Summary =====
 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 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 "AIREMLF90.log".
 \\
 See PREGSF90 with genotypes (SNP) for options.
@@ Line 128: / Line 126: @@
 <file>OPTION out_se_covar_function</file>
 Indicate to store in file samples of (co)variances function for postprocessing (histogram, etc.)
+===== 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 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 "AIREMLF90.log".