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
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?

##Exporting results to CSV

write.csv(as.data.frame(resOrdered),
file="CTRL_results.csv")

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?

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!