ordiellipse中的多边形面积是NaN - 为什么?

时间:2017-04-16 17:28:33

标签: r nan vegan mds

我正在尝试将椭圆添加到我在R上用Vegan包创建的NMDS图上,但是虽然代码没有出现错误,但没有多边形被绘制到我的图形上。使用summary()函数后,我发现多边形的面积是NaN,因此没有绘制多边形的原因。我不确定为什么我没有区域 - 这与我的数据有关吗?

我的数据可在此处找到:https://docs.google.com/spreadsheets/d/1uxWbKAvhdVqnorIMXURvYLrDZuoqejJpUsc9N6wSDxA/edit?usp=sharing

在三种类型的栖息地中进行了三次横断面 - 内陆森林,森林边缘和不安的栖息地。每只蜻蜓和豆娘都被看到了。

我的R代码如下:

OdonateNMDSdata <- read.csv(file.choose(), header=TRUE)
Odonaterownames <- row.names(OdonateNMDSdata) <- c("Interior", "Edge", "Disturbed")
library(vegan)
OdonateNMDS <- metaMDS(OdonateNMDSdata, k=2)
ordiplot(OdonateNMDS,type="n")
orditorp(OdonateNMDS,display="species",col="red",air=0.01)
orditorp(OdonateNMDS,display="sites",cex=1.25,air=0.01)
Ellipse <- ordiellipse(OdonateNMDS, groups=Odonaterownames, kind = "ehull", draw="polygon", col="blue", cex=0.7, conf=0.95)
summary(Ellipse)

由于

1 个答案:

答案 0 :(得分:1)

你有三个点,你想绘制三个椭圆,每个点一个。每个椭圆需要不止一个点(即使是两个点,封闭的椭圆也是连接点的直线)。

然而,似乎在包围椭圆(kind = "ehull")时我们将NaN作为单点椭圆的区域,而对于其他kind s,我们将该区域指定为{ {1}}一点。我会改变它。