DESeq2: thousands of NA for pvalues?
2
1
Entering edit mode
aishu.jp ▴ 30
@aishujp-15144
Last seen 4.2 years ago

I have a RNA seq data which I am trying to identify DEGs. Dataset contains three sets of Untreated - 2 replicates and Treated - 2 replicates. My DEseq2 results shows padj values NA and some p values and padj values as zero. A total of 3000 genes shows such result.

I tried using res <- results(dds, independentFiltering = F) and res <- results(dds, cooksCutoff=FALSE) but I am unable to understand what difference both of these filtering makes

Is it good to use independentFiltering = F or cooksCutoff=FALSE for 3000 genes which shows padj values NA and some p values and padj values as zero

Can anyone please help me how to solve this.

deseq2 rna-seq DEGs cooks cutoff • 25k views
ADD COMMENT
4
Entering edit mode
@mikelove
Last seen 10 hours ago
United States

The mechanism that sets adjusted p-values or p-values to NA is listed in the FAQ:

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#why-are-some-p-values-set-to-na

ADD COMMENT

Login before adding your answer.

Traffic: 573 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