недостающее значение при расчете бегущих медиан?

Я хотел бы сгладить временной ряд, чтобы избежать ложного дрожания / ошибки. Другими словами, я хочу сделать очень локальное надежное сглаживание.

Я столкнулся с 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 для каждого значения после этого.

Ответы на вопрос(2)

Ваш ответ на вопрос