You'll want something like:

cnts = counts(dds, normalized=T)
d = data.frame(cnts=cnts[some_gene_of_interest, ], group=c(rep("male", 4), rep("female", 4)))
library(ggplot2)
ggplot(d, aes(x=group, y=cnts)) + geom_boxplot(aes(color=group))

As an aside, there's a plotCounts() function that comes with DESeq2 that's also handy.

Hi there, I have something similar to that but mine is just a count table in a .txt file that has six columns as follows:

Gene Id WT_1 WT_2 WT_3 MUT_1 MUT_2

GENE1 321 5454 6565 564 535
GENE2 ..... ...... ...... ..... .....

What I would like to do is just boxplot the data like this (see pic).

Where F1 F2 F3... have to be the name of my genes.

How can i do that?? Thanks in advance