Заполнить NA во временном ряду только ограниченное число

Есть ли способ, которым мы можем заполнить NA s в zoo или xts объект с ограниченным номером NA s вперед. Другими словами, как заполнить NA s до 3 последовательных NA s, а затем сохранить NA s от 4-го значения до действительного числа.Заполнить NA во временном ряду только ограниченное число

Нечто подобное.

library(zoo) 
x <- zoo(1:20, Sys.Date() + 1:20) 
x[c(2:4, 6:10, 13:18)] <- NA 
x 

2014-09-20 2014-09-21 2014-09-22 2014-09-23 2014-09-24 2014-09-25 2014-09-26 
     1   NA   NA   NA   5   NA   NA 
2014-09-27 2014-09-28 2014-09-29 2014-09-30 2014-10-01 2014-10-02 2014-10-03 
     NA   NA   NA   11   12   NA   NA 
2014-10-04 2014-10-05 2014-10-06 2014-10-07 2014-10-08 2014-10-09 
     NA   NA   NA   NA   19   20

Желаемая выход, будет что-то с переменной п = 3

2014-09-20 2014-09-21 2014-09-22 2014-09-23 2014-09-24 2014-09-25 2014-09-26 
     1   1   1  1   5   5  5 
2014-09-27 2014-09-28 2014-09-29 2014-09-30 2014-10-01 2014-10-02 2014-10-03 
     5   NA   NA   11   12   12  12 
2014-10-04 2014-10-05 2014-10-06 2014-10-07 2014-10-08 2014-10-09 
     12   NA   NA   NA   19   20

Я попробовал много комбинации с na.locf(x, maxgap = 3) и т.д., без особого успеха. Я могу создать цикл, чтобы получить желаемый результат, мне было интересно, существует ли векторизованный способ достижения этого.

fillInTheBlanks <- function(v, n=3) { 
    result <- v 
    counter0 <- 1 
    for(i in 2:length(v)) { 
    value <- v[i] 
    if (is.na(value)) { 
     if (counter0 > n) { 
     result[i] <- v[i] 
     } else { 
     result[i] <- result[i-1] 
     counter0 <- counter0 + 1 
     } } 
    else { 
     result[i] <- v[i] 
     counter0 <- 1 
    } 
    } 
    return(result) 
}

Благодаря

источник

2014-09-19 Kate K

Добавление некоторых сценариях использования, когда мы имеем qtrly данные и мы знаем, что данные будут хорошими в течение следующих 3 месяцев и могут составлять максимум еще 3 месяца, но все, что выходит за пределы приемлемого предела, должно сделать данные действительно NA и не должно заполнять его до бесконечного типа сценариев. –

Вот еще один способ:

l <- cumsum(! is.na(x)) 
c(NA, x[! is.na(x)])[replace(l, ave(l, l, FUN=seq_along) > 4, 0) + 1] 
# [1] 1 1 1 1 5 5 5 5 NA NA 11 12 12 12 12 NA NA NA 19 20

редактировать: мой предыдущий ответ требуется, чтобы x нет дубликатов. В текущем ответе нет.

тесты

x <- rep(x, length.out=1e4) 

plourde <- function(x) { 
    l <- cumsum(! is.na(x)) 
    c(NA, x[! is.na(x)])[replace(l, ave(l, l, FUN=seq_along) > 4, 0) + 1] 
} 

agstudy <- function(x) { 
    unlist(sapply(split(coredata(x),cumsum(!is.na(x))), 
      function(sx){ 
      if(length(sx)>3) 
       sx[2:4] <- rep(sx[1],3) 
      else sx <- rep(sx[1],length(sx)) 
      sx 
      })) 
} 

microbenchmark(plourde(x), agstudy(x)) 
# Unit: milliseconds 
#  expr min  lq median  uq max neval 
# plourde(x) 5.30 5.591 6.409 6.774 57.13 100 
# agstudy(x) 16.04 16.249 16.454 17.516 20.64 100

источник

2014-09-19 20:54:37

Отлично, работает быстро. Спасибо. –

Действительно быстрее, чем слияние! – user3032689

без использования na.locf, но идея состоит в том, чтобы разделить ваши XTS группы лиц, не являющихся отсутствующих значений, то для каждой группы заменить только первое 3 значения (после отсутствия misssing один) с первое значение. Это цикл, но поскольку он применяется только к группе, он должен быть быстрее, чем простой цикл по всем значениям.

zz <- 
unlist(sapply(split(coredata(x),cumsum(!is.na(x))), 
     function(sx){ 
     if(length(sx)>3) 
      sx[2:4] <- rep(sx[1],3) 
     else sx <- rep(sx[1],length(sx)) 
     sx 
     })) 
## create the zoo object since , the latter algorithm is applied only to the values 
zoo(zz,index(x)) 

2014-09-20 2014-09-21 2014-09-22 2014-09-23 2014-09-24 2014-09-25 2014-09-26 2014-09-27 2014-09-28 2014-09-29 2014-09-30 2014-10-01 2014-10-02 
     1   1   1   1   5   5   5   5   NA   NA   11   12   12 
2014-10-03 2014-10-04 2014-10-05 2014-10-06 2014-10-07 2014-10-08 2014-10-09 
     12   12   NA   NA   NA   19   20

источник

2014-09-19 19:19:36 agstudy

Работает как шарм. Спасибо. Прежде чем принять ответ, вы увидите, есть ли какой-либо векторный метод его выполнения. –

И еще одна мысль, что, если я не пропустил что-то, кажется, действительно:

na_locf_until = function(x, n = 3) 
{ 
    wnn = which(!is.na(x)) 
    inds = sort(c(wnn, (wnn + n+1)[which((wnn + n+1) < c(wnn[-1], length(x)))])) 
    c(rep(NA, wnn[1] - 1), 
    as.vector(x)[rep(inds, c(diff(inds), length(x) - inds[length(inds)] + 1))]) 
} 
na_locf_until(x) 
#[1] 1 1 1 1 5 5 5 5 NA NA 11 12 12 12 12 NA NA NA 19 20

источник

2014-09-19 22:12:58

Заполнить NA во временном ряду только ограниченное число

ответ

Смежные вопросы