如何用省略ggplot2的箭头连接多个点

Question

我使用vegan包对一些植被数据的时间序列进行了排序。由于排序图通常杂乱无章，有许多数据点，因此我提取了前两个排序轴的特征值，并取了每组的平均值。现在，每个站点只有一个点（总共11个站点）。为了仍然显示某些变化，我添加了具有标准偏差和95％置信区间的椭圆：

我要做的最后一件事是用箭头连接同一组（A，B或C）的点，以指示随时间变化的方向。所有运动都是从右到左。

我最初想在ordiarrow中使用vegan函数，但这仅在class为decorana时有效。我的课程是factor。

使用ggplot2似乎不是有效的选项，因为ordiellipse函数（创建椭圆）在这里不起作用。

用于绘制数据的代码：

install.packages("vegan")
library(vegan)
plot(Ord_KIKKER, type = "n", main = "Kikkervalleien",
     xlab = "DCA1 Eigenvalue = 0.62", ylab = "DCA2 Eigenvalue = 0.39")


points(ORD_KIKKER, cex = 2, pch = 19,
       col = c("black", "black", "black", "red","red", "green", "green", "green", "blue", "blue", "blue"))

自从我在此处发布了精简的数据集以来，生成的图看起来有点不同。

我的数据（Ord_KIKKER）：

structure(list(DCA1 = c(2.676616032, 0.361181861, -1.363464067, 
3.176862449, -0.087190269, 2.059548542, 0.167440366, -0.459090096, 
1.571536367, 0.309623788, -0.25787459), DCA2 = c(0.276788721, 
0.422077659, 0.181723453, 0.221610649, 0.940063655, -0.116083905, 
-0.539375059, -0.545053063, -0.06120542, -0.367148924, -1.679257818
), Unique = structure(c(1L, 5L, 8L, 2L, 9L, 3L, 6L, 10L, 4L, 
7L, 11L), .Label = c("2001A", "2001B", "2001C", "2001D", "2008A", 
"2008C", "2008D", "2018A", "2018B", "2018C", "2018D"), class = "factor"), 
    BLOCK = structure(c(1L, 1L, 1L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 
    4L), .Label = c("A", "B", "C", "D"), class = "factor")), .Names = c("DCA1", 
"DCA2", "Unique", "BLOCK"), class = "data.frame", row.names = c("2001A", 
"2008A", "2018A", "2001B", "2018B", "2001C", "2008C", "2018C", 
"2001D", "2008D", "2018D"))

Answer 1

vegan::ordiarrows()将有效，如果您仅为其赋予具有分数的变量：

ordiarrows(Ord_KIKKER[,1:2], Ord_KIKKER$BLOCK) # one way

但是，您还应该记住在初始图中使用asp=1来使轴的纵横比相等。

我无法进行全面测试，因为无法使用发布的数据来重现图形：如果向plot(Ord_KIKKER, ...)发送数据框，则不会得到普通图，而是所有变量相对于每个变量的面板图其他（pairs()图），也为type = "n"参数给出了错误。似乎您改用了一些非标准的图形工具，而且我不确定vegan::ordiarrows()的标准 R 图形是否可以与这些图形结合使用。