How to add Log2FoldChange as an argument for plotMA in DESeq2
2
0
Entering edit mode
5.6 years ago

Hello,

I would like to add an argument to the plotMA function in DESeq2.  Specifically, I would like to highlight points on the MA plot that are not only alpha<0.1, but also log2FoldChange <-1 or >1.  Is there a way to achieve this?  I tried the following "ifelse" statement but it did not work.

plotMA(res,alpha=0.1,col=ifesle(log2FoldChange>0.75 | log2FoldChange<-0.75, "red", "grey"), cex=2, ylim=c(-2,2))

thank you!

DESeq2 • 2.6k views
ADD COMMENT
0
Entering edit mode
5.6 years ago
vivekbhr ▴ 630

You can subset the DESeq result object by abs(logFC) > 1, then plot by plotMA. That's the fastest..

ADD COMMENT
0
Entering edit mode
5.6 years ago

Thanks for the reply, Vivek.  I'd like to retain the non-significant points in my plot as grey.  By subsetting I only plot those that meet my criteria. 

 

ADD COMMENT
0
Entering edit mode

You can manually make the MA plot with all the genes, then color the genes of interest.

ADD REPLY

Login before adding your answer.

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