Can anyone help me with adding labels to significant differentially regulated genes from DESeq plotMA in R? I generate MA plot with a command -
plotMA(res, col = ifelse(res$padj >=0.1, "black", "red"), cex=0.8)
but am not sure how to add labels
Many thanks, Pavel
Hi,
I am not familiar with the MA plot from DESeq but adding text based on a subset of values would probably work out for you. Here is an example for labeling only the significant ones with value < 0.05. If you can get access to the fold change values and the average intensity (the x axis) you can do:
xx
Agree with Ilian, the label can be added using "text(xxx)". I am sure Vasek can help you.
Hi Liuyang and Ilian, many thanks for your suggestions. I will try it. Best, Pavel
P.S. of course Vašek :)
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics
- Statistical Software
- R