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()