R и Zoo: скользящая средняя с оборванными данными?

Question

Я хочу рассчитать скользящие средние, чтобы заполнить записи NA известными позициями 3, 5 и 1. Как я могу это сделать с зоопарком пакета в R?

Входной

> library(zoo) 
> hh <- c(NA, NA, NA, 3, NA, 5, NA, 1, NA, NA, NA, NA) 

Сбой

неудачно с rollmean

> rollmean(hh,na.omit=TRUE,k=1) 
    [1] NA NA NA NA NA NA NA NA NA NA NA NA 
    > rollmean(hh,4, na.omit=TRUE,k=1) 
    [1] NA NA NA NA NA NA NA NA NA NA NA NA 
    > rollmean(hh,4, na.rm=TRUE,k=1) 
    [1] NA NA NA NA NA NA NA NA NA NA NA NA 

неудачно с Rollapply

> rollapply(hh, 4, function(x) mean(x)) 
    [1] NA NA NA NA NA NA NA NA NA 

Предназначенный выход что-то вроде

> COMMAND(hh, movingAverageNumber, function(x) mean(x)) 
    [1] 3 3.3 3.4 3 4 5 3 1 2 1.5 1.2 1 0.8 

Связанные

1. Using rollmean when there are missing values (NA)

Answer 1

В зависимости от того, что вы ищете, как насчет:

rollapply(hh, 5, mean, na.rm = TRUE) 
    [1] 3 4 4 3 3 3 1 1 

или

rollapply(hh, 4, mean, na.rm = TRUE) 
    [1] 3 3 4 4 3 3 1 1 NaN 

Answer 2

Предполагая, что вы хотите, прокатные mean для четырех значений (вы передаете k=1 к rollmean), вам нужно сказать mean уронить NA

library(zoo) 
rollapply(hh, 4, function(x) mean(x, na.rm=TRUE)) 

[1] 3 3 4 4 3 3 1 1 NaN 

Хотя это возвращает нечто иное, чем ваш желаемый результат. Вам нужно будет объяснить это дальше, если вы хотите что-то другое.