adehabitatHS; Предупреждение в presel() fxn
Я использую adehabitatHS, чтобы сделать анализ среды обитания третьего порядка на данных ошейника GPS. Все мои данные прочитаны нормально и, похоже, правильно отформатированы. Однако когда я бегу
pks<-prepksel(slope, hr, cp)
Warning message:
In rbind(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, :
number of columns of result is not a multiple of vector length (arg 2)`
Я побежал след и получил следующее
traceback()
17: getClassDef(Class, where)
16: getClass(cl)
15: checkSlotAssignment(object, name, value)
14: `slot<-`(`*tmp*`, what, value = slot(from, what))
13: asMethod(object)
12: as(from, "SpatialPoints")
11: initialize(value, ...)
10: initialize(value, ...)
9: new("SpatialPointsDataFrame", as(from, "SpatialPoints"), data =
from@data,
coords.nrs = from@coords.nrs)
8: asMethod(object)
7: as(x, "SpatialPointsDataFrame")
6: dv[!is.na(slot(dv, "data")[, 1]), ]
5: dv[!is.na(slot(dv, "data")[, 1]), ]
4: FUN(X[[i]], ...)
3: lapply(1:ncol(hr), function(i) {
dv <- hr[, i]
dv <- dv[!is.na(slot(dv, "data")[, 1]), ]
dat[over(dv, geometry(sa)), ]
})
2: do.call("rbind", lapply(1:ncol(hr), function(i) {
dv <- hr[, i]
dv <- dv[!is.na(slot(dv, "data")[, 1]), ]
dat[over(dv, geometry(sa)), ]
}))
1: prepksel(slope, hr, cp)`
Между предупреждением и трассировкой я обнаружил, что проблема была в моем файле SpatialPixelsDataFrame, поэтому я вывел его резюме, однако не смог найти ничего, что могло бы вызвать ошибку в функции rbind. Помимо отличных фактических значений, данные выглядели так же, как примеры данных из puech.
Любая помощь будет оценена. Я могу опубликовать резюме из hr или любого другого фрагмента кода или результатов, этот пост уже становился длинным. Благодарю.