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 или любого другого фрагмента кода или результатов, этот пост уже становился длинным. Благодарю.

0 ответов

Другие вопросы по тегам