Identifying DEGs using lfcshrink()
Entering edit mode
Zainab • 0
Last seen 2 hours ago
United States

Hi, I have a question about using lfcshrink() vs results() to identify DEGs.

I'm currently identifying upregulated DEGs as having padj<0.01 and log2FC>1. Using unshrunken results, I am getting around 800 DEGs which is higher then I expected. A lot of these DEGs come from genes with low counts, so if I apply lfcshrink() and use these shrunken log2FC to identify DEGs, I end up with around 400 DEGs.

It seems like lfcshrink() is used for ranking DEGs, can it also be used to identify DEGs? I'm currently filtering low counts by CPM. Is there anything else I can use to clean up my data?

DESeq2 lfcshrink • 292 views
Entering edit mode
Last seen 5 hours ago
United States

Yes, you can use lfcShrink with svalue=TRUE to identify DEG. See the apeglm vignette for some discussion about s-value compared to q-value (padj).


Login before adding your answer.

Traffic: 443 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6