我正在尝试创建一系列可发表的PCA关于森林不同层次(树冠和林下)昆虫丰度的双图。我使用rda()函数完成此操作。为了创建图,我使用了素食的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))结果相当不错:

但是,因为我将把几个双图连接到一个更大的图中,所以科的标签(通常称为“物种标签”)太小了。我试图在biplot()中使用cex=1.5更改它们,但似乎被函数的默认设置覆盖了。当我通过设置type=c("none", "none")创建一个空的双线图框架时,我可以使用text()函数添加更大的物种标签,例如
text(family.pca, display = "species", cex = 1.0, col = "black")但是我不知道如何将箭头添加到图中(我真的希望在那里有箭头...)。
有没有人知道解决这个问题的办法?答案是非常感谢的。
https://stackoverflow.com/questions/47661984
复制相似问题