要注意PCA,princomp,prcomp的input,行是样本,列是特征(如基因)
但是princomp只能适用于特征数小于样本数的input,所以不适用于表达谱数据。
library(ggbiplot)
data(wine)
wine.pca <- prcomp(wine, scale. = TRUE)
ggbiplot(wine.pca, obs.scale = 1, var.scale = 1,
groups = wine.class, ellipse = TRUE, circle = TRUE) +
scale_color_discrete(name = '') +
theme(legend.direction = 'horizontal', legend.position = 'top')