Multi-Dimensional Scaling plot of samples based on the 500 most variant features (for edgeR analyses)
MDSPlot(dge, group, n = min(500, nrow(dge$counts)),
gene.selection = c("pairwise", "common"), col = c("lightblue", "orange",
"MediumVioletRed", "SpringGreen"), outfile = TRUE)
a DGEList
object
vector of the condition from which each sample belongs
number of features to keep among the most variant
"pairwise"
to choose the top features separately for each pairwise comparison between the samples or "common"
to select the same features for all comparisons. Only used when method="logFC"
colors to use (one per biological condition)
TRUE to export the figure in a png file
A file named MDS.png in the figures directory