gravatar for tianshenbio

2 hours ago by

In the following analysis using deseq2

dds_out <- DESeq(dds)
res <- results(dds_out)
summary(res)
out of 19287 with nonzero total read count
adjusted p-value < 0.1
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

link

modified 1 hour ago

written
2 hours ago
by

tianshenbio40



Source link