我目前正在一个包含 2000 个观测值和 21 个变量的大型数据集上运行 PCA(主成分分析),其中 20 个变量是我对化石所做的测量,最后一个是关于物种的。2000 个观察结果中有 7 个是我在绘制 PCA 时想要增加尺寸的化石。我目前正在运行的内容:
log.ir <- log(Size_adjusted_measurements[, 1:20])
ir.species <- Size_adjusted_measurements[, 21] ir.pca <- prcomp(log.ir, center = TRUE, scale. = TRUE)
现在绘制 PCA:
g <- ggbiplot(ir.pca, obs.scale = 1, var.scale = 1, groups = ir.species, labels = FALSE)
g <- g + scale_colour_discrete(name = "Modern") g <- g + theme(legend.direction = 'horizontal', legend.position = 'top') print(g)
因为我不允许发布图片,所以我已经与情节建立了 gyazo 链接: https ://gyazo.com/b8706a3063243d5d8e19bb3f8b2eb520
被黑眼圈包裹的是形状空间中的化石。与其他化石相比,如何修改我的绘图代码以增加化石的大小?我试过g <- g + scale_size_manual
了,但这不会起作用。非常感谢任何帮助。