Entering edit mode
3.7 years ago
CY ▴ 710
DESeq2 uses GLM with negative binomial distibution to model the regressor for variance shrinkage. It is obvious that the expression variance is different for gene with different expression level and a ordinary linear model is not suited. I am wonder if it is possible to model this regressor using LOESS instead of GLM?
You can ask the maintainer for an experts opinion at BioC (https://support.bioconductor.org/) but I would avoid the term
obviousunless you are 100% sure what you are talking about and have experts knowledge on why things are modeled/used as they are.