This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Adding separate color a list of genes in plotMA

Hello all and happy NEW YEAR. I am trying to create either plotMA, DESeq2 or VolcanoPlot such that in addition to the original plot (Usually in black and red) I can supply a list of genes (present in the plot) from the DESeq2 result to appear as separate colors (e.g green) in the plot. I haven't seen a good answer to address this, help. See my original DESeq2 command bellow.

library(DESeq2)
coldata <- read.csv("Conditions2.txt", sep = '\t')
cts <- as.matrix(read.csv("counts", sep = '\t', row.names = 1, header = TRUE))
dds <- DESeqDataSetFromMatrix(countData = cts, colData = coldata, design = ~ Condition)
dds <- DESeq(dds)
res <- results(dds) 
write.table(res, file = "NEGvsZC.txt", sep = '\t')
pdf("NEGvsZC.pdf")
dev.off()
deseq2 plotma volcanoplot r rna-seq
• 1,462 views • link
written by jaqx008 jaqx008

If external libraries are an option, I'd recommend trying Kevin Blighe's EnhancedVolcano package. That makes highly customizable volcano plots.

• 0 views • link
written by Ram Ram

Thanks for your recommendation. I am not sure what you mean by external libraries.

• 0 views • link
written by jaqx008 jaqx008

I was referring to anything that's not a part of basic/default R packages (graphics, stats, etc are basic default R packages, ggplot2, ComplexHeatmap, EnhancedVolcano, dplyr etc are external libraries)

• 0 views • link
written by Ram Ram

0 answers

No answers yet.

Log in to answer this question.

More posts like this