how to draw the volcano plot for DEGs in a different way
1
0
Entering edit mode
12 weeks ago

Dear guys,

May I know if there is some sample code to use to draw the plot like below?

I have DEGs results from limma, does it need to re run use DESeq2 to use plotMA to get this?

Thank you very much!

enter image description here

volcano-plot DEG • 386 views
ADD COMMENT
2
Entering edit mode
12 weeks ago
CTLong ▴ 90

From you Limma DEG results, you could use glimmaMA to generate a MA plot as demonstrated in this vignette (https://bioconductor.org/packages/devel/bioc/vignettes/Glimma/inst/doc/limma_edger.html).

Alternatively, given the DEG matrix, you could also useggplot2 with geom_point to draw the MA plot. Set the aesthetics so that the Y-axis corresponds log2FC and X axis as average log2 expression. Color the points based on if it is upregulated or downregulated (according to your threshold). The code will look something like this:

ggplot(LimmaResult, aes(x=averageExp, y=log2FC)) +
  geom_point(aes(color = Significance))
ADD COMMENT
0
Entering edit mode

Thank you very much! This is super helpful! Appreciate!

ADD REPLY

Login before adding your answer.

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