Principal Component Analysis of samples based on the 500 most variant features on VST- or rlog-counts (if use of DESeq2)
PCAPlot(counts.trans, group, n = min(500, nrow(counts.trans)),
col = c("lightblue", "orange", "MediumVioletRed", "SpringGreen"),
outfile = TRUE)
a matrix a transformed counts (VST- or rlog-counts)
factor vector of the condition from which each sample belongs
number of features to keep among the most variant
colors to use (one per biological condition)
TRUE to export the figure in a png file
A file named PCA.png in the figures directory with a pairwise plot of the three first principal components