RELAX does not reproduce its own estimates at default settings, and its output does not show it
Selection-intensity estimates from RELAX are reported as a point value of K with a likelihood-ratio P. We report that, at default settings and on data of ordinary size, the program does not reproduce its own fits. Of 27 enzyme entries refitted under two optimiser configurations, none reproduced its log-likelihood to within 0.01 units; the median change was 103 units, the largest over 3,400, and…
We haven't written up this one. bioRxiv has the full story — the link below goes straight to it.
