我修改了虹膜数据,以提供我想要做的事情的一个例子。我在虹膜数据中添加了一个额外的列,如下面的链接所示。这个额外的专栏还有一些基于器官的小组。
然后我做PCA并绘制它。我希望在图中有两组按颜色和形状区分。器官通过四种形状进行区分,但除了图中所示的蓝色外,物种不按颜色区分。
library(ggbiplot)
data(iris)
#### add an extra column with organ information
iris$organ<-c(rep("leaf",50),rep("root",50),rep("shrub",25),rep("petal",25))
# log transform
log.ir <- log(iris[, 1:4])
#grouping by species
ir.species <- iris[, 5]
#grouping by organs
ir.organ <- iris[,6]
##PCA
ir.pca <- prcomp(log.ir, center = TRUE,scale. = TRUE)
###PCA plots
g <- ggbiplot(ir.pca, obs.scale = 1, var.scale = 1,groups = ir.species)
g+geom_point(aes(shape=factor(ir.organ)))