I am trying to draw box plots to show change in expression value between two conditions using deseq2.
d <- plotCounts(dds, gene="gene_id", intgroup="condition", returnData = TRUE) ggplot(d, aes(x=condition, y=log2(count))) + geom_boxplot()
Is this enough to show the differential expression of a gene between two conditions?
I would like to draw boxplots with logFC values on Y axis for a single gene in two conditions. Can please let me know how to do that?