Я использовал пакет spatialEco
для определения тренда, значения p и значения тау (raster.kendal). Он давал p-значение до 0,2, поэтому я хочу сделать prewhitnning.
, в то время как я использую функцию kendall для prewhitnning. Это дает ошибку как для формата ASCII, так и для формата GTIFF. (Error1). когда я выбираю файл вместо стека, он выдает другую ошибку (error2)
ascidata <- list.files(path = "F:/ShareBoth/../ascii/", full.names = TRUE)
length(ascidata)
asciistack <- stack(ascidata)
prewit <- kendall(asciistack, p.value = TRUE, prewhiten = TRUE)
tifdata <- list.files(path = "F:/ShareBoth/../tif/", full.names = TRUE)
length(tifdata)
tifstack <- stack(tifdata)
prewit <- kendall(tifstack, p.value = TRUE, prewhiten = TRUE)
, что дает error1: Error in y[!is.na(y)] : object of type 'S4' is not subsettable
tifdata <- list.files(path = "F:/ShareBoth/../tif/", full.names = TRUE)
length(tifdata)
tifstack <- stack(tifdata)
prewit <- kendall(tifstack$layer1, p.value = TRUE, prewhiten = TRUE)
выдает вторую ошибку error2: Error in model.frame.default(formula = y ~ 1, na.action = function (object, :object is not a matrix