недостающее значение при расчете бегущих медиан?
Я хотел бы сгладить временной ряд, чтобы избежать ложного дрожания / ошибки. Другими словами, я хочу сделать очень локальное надежное сглаживание.
Я столкнулся с rollmean и rollmedian в пакете зоопарка, но столкнулся с проблемой, потому что мой вектор содержал NA. Затем я где-то читал, что эти функции зоопарка используют runmed и в этом заключается проблема.
== Примеры ==
median(c(1,1,1,2,2,2,7,NA,1,2,3,10,10,10),na.rm = TRUE)
runmed(c(1,1,1,2,2,2,7,NA,1,2,3,10,10,10),k=3)
Первая строка возвращает 2, но вернула бы NA, еслиna.rm = TRUE
не был включен. Вторая строка возвращаетError in runmed(c(1, 1, 1, 2, 2, 2, 7, NA, 1, 2, 3, 10, 10, 10), k = 3) :
NA/NaN/Inf in foreign function call (arg 1)
, Нет способа добавить аргумент na.rm в строку.
How can I get runmed to handle the NA? Кстати, rollmean возвращает вектор, который является правильным вплоть до NA, а затем возвращает NA для каждого значения после этого.