森林のさまざまな層 (樹冠と下層) における昆虫の豊富さに関する 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))
結果はとてもいいです:
しかし、いくつかのバイプロットを結合してより大きなグラフにするため、ファミリーのラベル (通常は「種のラベル」と呼ばれます) が小さすぎます。cex=1.5inを使用して変更しようとしましbiplot()たが、関数のデフォルトによって上書きされるようです。を設定して空のバイプロット フレームを作成すると、関数type=c("none", "none")を使用してより大きな種のラベルを追加できます。text()
text(family.pca, display = "species", cex = 1.0, col = "black")
しかし、その後、矢印をプロットに追加する方法がわかりません（そこに矢印を入れたいと思います...）。
この場合の解決策を知っている人はいますか？回答をお待ちしております。