У меня есть данные о скоплениях гнездящихся птиц из карьеров и известняковых пастбищ, и я хочу отобразить их на диаграмме распределения НМДС вместе с параметрами окружающей среды в качестве векторов. Как мне убрать теги видов? Я использую веганский пакет для посвящения.
Я уже уменьшил количество отображаемых видов с помощью функции ordiselect
из пакета goeveg. Теги все еще перекрываются, но я не хочу еще больше сокращать виды.
С функцией orditkplot
я не могу получить результат, который выглядел бы так хорошо, как я хотел бы с точки зрения эстетики. Мне также не удалось отобразить только подмножество видов при использовании этой функции.
ord.nmds <- metaMDS(species, trace=FALSE)
sel <- ordiselect(species, ord.nmds, ablim = 0.5, method = "axes", freq = TRUE)
ef <- envfit(ord.nmds ~ cgr + L_wdl + str + Hab_Div, env_var, permu=999)
grp <- factor(rep(c("GRASS", "QUARRY"), each = 14))
cols <- c("gray70","gray30")
plot(ord.nmds, type="n", xlim = c(-1.1, 1.1), ylim = c(-0.9,0.9))
points(ord.nmds, display = "sites", pch = 15, col = rep(cols, each = 14), cex=1)
text(ord.nmds, display = "species", col="black", cex=0.55, select = sel)
plot(ef, p.max = 0.05, col ="black", cex = 0.8)
legend(x = 0.6, y = 1.05, legend = tools::toTitleCase(levels(grp)), fill = cols, bty = "n", cex = 0.8)
Пока что мой график выглядит так: