我正在尝试在 vhf 跟踪数据上计算 MCP 和 KDE,我的数据框有三列,Individual、Lat 和 Long。
我正在尝试按照 adahabitatHR 小插图创建一个 SpatialPointsDataFrame,其中包含多个人的甚高频数据(https://cran.r-project.org/web/packages/adehabitatHR/vignettes/adehabitatHR.pdf),但使用我自己的数据集。
我加载我的数据
vhf <- read.csv("vhfdatacrom.csv")`
str(vhf)
数据帧':186 obs。3 个变量: $ 个体:因子 w/ 12 个级别 "PM01","PM03",..: 2 1 2 2 1 1 1 1 1 4 ... $ Lat : 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(,"package") [1] "sp"
> clusthr(idsp)
tmp[1, ] 中的错误:下标越界
我已经尽我所能,但无法让它运行,任何帮助将不胜感激。
谢谢