Hi all, I'm using DESeq2 for the first time and ran into some problems during the LFC shrinkage step. I got this error when I tried to run the process function "lfcShrink()":

using 'normal' for LFC shrinkage, the Normal prior from Love et al (2014). Note that type='apeglm' and type='ashr' have shown to have less bias than type='normal'. See ?lfcShrink for more details on shrinkage type, and the DESeq2 vignette. Reference: https://doi.org/10.1093/bioinformatics/bty895 Error in lfcShrink(ddshi, contrast = c("surgery", "SCI", "Sham"), res = hires) : LFC shrinkage type='normal' not implemented for designs with interactions

But the thing is, when I plotted the dispersion of the genes, this is what I got in the output:

Honestly, it doesn't look all that dispersed at all, can I just simply skip this step?

You don't need LFC shrinkage if you are only interested in getting a set of genes with FDR bound.

Re: the error above, see more details in the vignette. You can use type="apeglm" or "ashr" instead.