Number of DEGs with adj pvalue < 0.05
1
0
Entering edit mode
GlycineMax • 0
@user-24655
Last seen 27 days ago
United States

Hello, after using DESeq2 to analyze my dataset, I called the following code to count the number of genes with adj pvalue<0.05. However, when I export the data, the number of genes with adj pvalue<0.05 in the spreadsheet is much lower than that (125). I exported right after the above code, so no more data filtering/transformation was applied.


> ##using p=0.05
> res05 <- results(dds, alpha=0.05)
> summary(res05)

out of 45412 with nonzero total read count
adjusted p-value < 0.05
LFC > 0 (up)       : 2626, 5.8%
LFC < 0 (down)     : 2197, 4.8%
outliers [1]       : 270, 0.59%
low counts [2]     : 7801, 17%
(mean count < 13)
[1] see 'cooksCutoff' argument of ?results
[2] see 'independentFiltering' argument of ?results

> 
> sum(res05$padj < 0.05, na.rm=TRUE)
[1] 4823

Anybody has any idea why that is happening?

Thank you!

DESeq2 • 195 views
ADD COMMENT
0
Entering edit mode

The relevant part is missing => the export code.

ADD REPLY
0
Entering edit mode

ATpoint, sorry for the late reply. The export code is:

####p-values and adjusted p-values
#We can order our results table by the smallest p value:

resOrdered <- res_CTRL[order(res_CTRL$pvalue),]

summary(res_CTRL)

#How many adjusted p-values were less than 0.1?

sum(res_CTRL$padj < 0.05, na.rm=TRUE)

##Exporting results to CSV

write.csv(as.data.frame(resOrdered), 
          file="CTRL_results.csv")
ADD REPLY
0
Entering edit mode
@mikelove
Last seen 6 hours ago
United States

I don't know how you are counting adj p-value < 0.05 outside of R, could it be there is a mistake on that end?

ADD COMMENT
0
Entering edit mode

I was able to figure out, it was something wrong in my code, since I was running several treatments in one code, some files were overwriting.

Thank you all for helping!

ADD REPLY

Login before adding your answer.

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