Eu gostaria de suavizar uma série de tempo para evitar jitter / erro espúrio. Em outras palavras, eu quero fazer algumas suavizações robustas muito locais.

Eu me deparei com rollmean e rollmedian no pacote de zoológico mas colidiu com um problema porque meu vetor teve um NA nisto. Eu li então em algum lugar que aquelas funções do jardim zoológico usam runmed e nisso está o problema.

== exemplos ==

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)

A primeira linha retorna 2, mas teria retornado NA sena.rm = TRUE não foi incluído. A segunda linha retornaError 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). Não há como adicionar um argumento na.rm à linha.

Como posso ser treinado para lidar com o NA? By the way, rollmean retorna um vetor que está correto até o NA e, em seguida, retorna NA para cada valor posterior.