How to generate a dataframe of deseq2 result excluding 'low counts' and 'outliers' ?
1
0
Entering edit mode
11 months ago
tianshenbio ▴ 110

In the following analysis using deseq2

dds_out <- DESeq(dds)
res <- results(dds_out)
summary(res)
out of 19287 with nonzero total read count
LFC > 0 (up)       : 3502, 18%
LFC < 0 (down)     : 4187, 22%
outliers [1]       : 0, 0%
low counts [2]     : 2992, 16%
(mean count < 1)
[1] see 'cooksCutoff' argument of ?results
[2] see 'independentFiltering' argument of ?results


How can I save the res output without the 'outliers' and 'low counts'? Since I need to get the whole set of genes that are analyzed to generate the DE genes as the background gene set in my subsequent GO enrichment analysis using goseq2

deseq2 RNA-Seq • 584 views
0
Entering edit mode

I think the p-value or FDR or both for those should be set to NA. There are no outliers in your example.

0
Entering edit mode
11 months ago

[1] see 'cooksCutoff' argument of ?results

[2] see 'independentFiltering' argument of ?results

By the way, I would get out of the habit of running results() like that unless you are absolutely certain of the contrast that is being performed. In doing it this way, you are also avoiding the log fold-change shrinkage step.

Kevin