2

我目前正在一个包含 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 = 'horizo​​ntal', legend.position = 'top') print(g)

因为我不允许发布图片,所以我已经与情节建立了 gyazo 链接: https ://gyazo.com/b8706a3063243d5d8e19bb3f8b2eb520

被黑眼圈包裹的是形状空间中的化石。与其他化石相比,如何修改我的绘图代码以增加化石的大小?我试过g <- g + scale_size_manual了,但这不会起作用。非常感谢任何帮助。

图片

4

0 回答 0