Ich möchte eine Zeitreihe glätten, um unechte Jitter / Fehler zu vermeiden. Mit anderen Worten, ich möchte eine sehr lokale, robuste Glättung durchführen.

Ich bin im Zoo-Paket auf rollmean und rollmedian gestoßen, bin aber auf ein Problem gestoßen, weil mein Vektor eine NA enthielt. Ich habe dann irgendwo gelesen, dass diese Zoofunktionen runmed verwenden und darin liegt das Problem.

== Beispiele ==

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)

Die erste Zeile gibt 2 zurück, hätte aber NA zurückgegeben, wennna.rm = TRUE war nicht enthalten. Die zweite Zeile kehrt zurückError 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). Es gibt keine Möglichkeit, der Zeile ein na.rm-Argument hinzuzufügen.

Wie kann ich mit der NA fertig werden? Übrigens gibt rollmean einen Vektor zurück, der bis zur NA korrekt ist, und gibt danach für jeden Wert NA zurück.