DESeq2 with too many rejections
1
0
Entering edit mode
Raymond ▴ 20
@raymond-14020
Last seen 4.9 years ago

Hi,

  I had a groups of samples with 12 different conditions. I followed the DESeq2 tutorial in Bioconductor, and then plot the rejections:

 

{plot(metadata(res)$filterNumRej, 
     type="b", ylab="number of rejections",
     xlab="quantiles of filter")
lines(metadata(res)$lo.fit, col="red")
abline(v=metadata(res)$filterTheta)}

The figure shows that there are a lot of rejections, does this make sense?


https://imgur.com/a/lwtvGOF

Filter vs no filter


```{r}
resNoFilt <- results(dds, independentFiltering=FALSE)
addmargins(table(filtering=(res$padj < .1),
                 noFiltering=(resNoFilt$padj < .1)))
```
         noFiltering
filtering FALSE  TRUE   Sum
    FALSE 19068     0 19068
    TRUE    845 14266 15111
    Sum   19913 14266 34179

 

 

Thanks & regards,

Raymond

deseq2 • 562 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 15 hours ago
United States

Take a look at the PCA plot. Also see the section of the vignette on LFC threshold tests.

ADD COMMENT

Login before adding your answer.

Traffic: 482 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6