This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to highlight DE genes in plotMD ?

How to highlight the DE genes with FDR =0.10 using plotMD as we usually see in volcanoplot? Example figures are given below ![volcanoplot][1] [1]: https://ibb.co/fEZejx ![maplot][1] [1]: https://ibb.co/mZ3dcH How to create volcanoplot in edgeR using coefficient of contrasts.

edger plotmd volcanoplot mirnaseq
• 5,680 views • link
updated by raulgcova97 • raulgcova97 • written by Björn • Björn

Could you give use a reproducible example of your input dataset please ?

• 0 views • link
written by Nicolas Rosewick Nicolas Rosewick

plotMD doesn't use ggplot graphing system as backend. It uses plotwithhighlights as engine and probably supports values function to highlight genes of interest.

• 0 views • link
written by cpad0112 cpad0112

2 answers

I don't understand how you used ggplot on plotMD. It doesn't seem to be compatible with it.

• 0 views • link
written by valizad2 • valizad2

Hi I managed to add DE tags to the plotMD using the code I published on Kaggle: Kaggle DE analysis

Here is an example of the code used:

fit <- glmQLFit(d, design)
qlftest<-glmQLFTest(fit, coef=i)
results <- as.data.frame(topTags(qlftest,n = Inf))
plotMD(qlftest,main="Mean-Difference Control vs T15")
text(results$logCPM,results$logFC,labels = rownames(results),col="black",cex=0.5,pos=3)

Here is the result:

plotMD

Also here is the original code which I took as reference (Author: Charleen D. Adams): Differential expression (ER-/ER+) for 125 breast cancers

• 0 views • link
written by raulgcova97 • raulgcova97

Log in to answer this question.

More posts like this