The glm_degs function filters axial genes using qval and log-likelihood after calculating them(spateo-release/spateo/tools/glm.py line 133). After reviewing relevant materials, I found that the smaller the absolute value of log-likelihood, the better the fitting degree.
However, when the glm_degs function filters axial genes based on log-likelihood , it only keep the glm test results whose log-likelihood is less than the (default: -2000), which appears to be incorrect. Similarly, the sorting direction of log-likelihood also seems to be reversed(spateo-release/spateo/tools/glm.py line 128).
Could you please confirm this point?
Thank you !
