Hi, Very simple question but I couldn't find any answer online. I'm currently analyzing some RNAseq data using DESeq 2 and produced PCA plot using the plotPCA funtion. I want to add a simple title to my plot but I can't figure it out;

The current code I'm using that works without a title is this

plotPCA(vsdata, intgroup=c("infection","CellType"))

I tried using arguments like main = "my title" as follows

plotPCA(vsdata, intgroup=c("infection","CellType"), main = "my title" )

or adding a line like title("infection") which works for plotMDS when I worked with limma-voom

plotPCA(vsdata, intgroup=c("infection","CellType"))
title("infection")

However neither of those actually produced a title for my PCA plot.

Thanks in advance

The function uses ggplot2, so you need + ggtitle('title'). Please take a ggplot2 tutorial for details.