Chciałbym wygładzić szereg czasowy, aby uniknąć fałszywego fluktuacji / błędu. Innymi słowy, chcę wykonać bardzo lokalne lokalne wygładzanie.

Natknąłem się na rollmean i rollmedian w pakiecie zoo, ale wpadłem na problem, ponieważ mój wektor miał w nim NA. Potem przeczytałem gdzieś, że te funkcje zoo używają biegu i w tym tkwi problem.

== przykłady ==

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)

Pierwsza linia zwraca 2, ale zwróciłaby NA jeślina.rm = TRUE nie został uwzględniony. Powraca druga liniaError 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). Nie ma możliwości dodania argumentu na.rm do linii.

Jak mogę zostać zabity do obsługi NA? Przy okazji, rollmean zwraca wektor, który jest poprawny do NA, a następnie zwraca NA dla każdej późniejszej wartości.