我正在尝试根据vhf跟踪数据计算MCP和KDE,我的数据框有三列,分别为:个人,纬度和经度。
我正在尝试遵循adahabitatHR小插图,尝试创建一个SpatialPointsDataFrame,其中包含多个人(https://cran.r-project.org/web/packages/adehabitatHR/vignettes/adehabitatHR.pdf)的vhf数据,但是使用我自己的数据集。
vhf <- read.csv("vhfdatacrom.csv")`
str(vhf)
data.frame':186磅。 3个变量: $个人:有12个等级的系数“ PM01”,“ PM03”,..:2 1 2 2 1 1 1 1 1 4 ... $纬度:num 54.2 54.2 54.2 54.2 54.2 ... $ Long:数字-7.44 -7.42 -7.45 -7.44 -7.42 ...
> coords.data.xy <- vhf[c("Long", "Lat")]
> id <- vhf[c("Individual")]
> idsp <- data.frame(id)
> coordinates(idsp) <- coords.data.xy
> class(idsp)
[1]“ SpatialPointsDataFrame” attr(,“包裹”) [1]“ sp”
> clusthr(idsp)
tmp [1,]中的错误:下标超出范围
我已经竭尽全力,但是无法运行它,任何帮助将不胜感激。
谢谢