我正在尝试创建一系列关于森林不同层（树冠和林下）昆虫丰度的 PCA 双图。我正在使用 rda() 函数执行此操作。为了创建情节，我使用了 vegan 的biplot()函数并对其进行了一些修改：
Family=read.table("Family2.txt", header=T)
strata=read.table("Strata.txt", header=T)
family.pca=rda(Family)
with(strata, levels(Strata))
biplot(family.pca, type=c("text", "none"), col=c("black", "black"), xlab="",
ylab="")
title(xlab="PC1 (86.8%)", ylab="PC2 (9.7%)", mgp=c(2.2, 2.2, 0))
points=c(16, 1)
colour=c("black", "black")
with(strata, points(family.pca, display = "sites", col = colour, pch =
points))
with(strata, legend("topright", legend = levels(Strata), bty = "n", col =
colour, pch = points, pt.bg = points))
结果非常好：
但是因为我会将几个双图连接到一个更大的图表中，所以家庭的标签（通常称为“物种标签”）太小了。我尝试使用cex=1.5in更改它们biplot()，但它似乎被函数的默认值覆盖。当我通过设置创建一个空的双标图框时type=c("none", "none")，我可以使用该函数添加更大的物种标签text()，例如
text(family.pca, display = "species", cex = 1.0, col = "black")
但后来我不知道如何将箭头添加到情节中（我真的很想把箭头放在那里......）。
有人知道这种情况的解决方案吗？非常感谢答案。