data.frame(またはdata.table)で、NAを以前の最も近い非NA値で「転送」します。(の代わりにdata.frame
)ベクトルを使用する簡単な例は次のとおりです。
> y <- c(NA, 2, 2, NA, NA, 3, NA, 4, NA, NA)
私はそのようなものfill.NAs()
を構築することを可能にする関数を望みますyy
:
> yy
[1] NA NA NA 2 2 2 2 3 3 3 4 4
この操作を多くの(合計〜1 Tb)小さなサイズdata.frame
(〜30-50 Mb)で繰り返す必要があります。ここで、行はNAであり、そのエントリはすべてです。問題に取り組む良い方法は何ですか?
私が作った醜い解決策はこの関数を使っています:
last <- function (x){
x[length(x)]
}
fill.NAs <- function(isNA){
if (isNA[1] == 1) {
isNA[1:max({which(isNA==0)[1]-1},1)] <- 0 # first is NAs
# can't be forward filled
}
isNA.neg <- isNA.pos <- isNA.diff <- diff(isNA)
isNA.pos[isNA.diff < 0] <- 0
isNA.neg[isNA.diff > 0] <- 0
which.isNA.neg <- which(as.logical(isNA.neg))
if (length(which.isNA.neg)==0) return(NULL) # generates warnings later, but works
which.isNA.pos <- which(as.logical(isNA.pos))
which.isNA <- which(as.logical(isNA))
if (length(which.isNA.neg)==length(which.isNA.pos)){
replacement <- rep(which.isNA.pos[2:length(which.isNA.neg)],
which.isNA.neg[2:max(length(which.isNA.neg)-1,2)] -
which.isNA.pos[1:max(length(which.isNA.neg)-1,1)])
replacement <- c(replacement, rep(last(which.isNA.pos), last(which.isNA) - last(which.isNA.pos)))
} else {
replacement <- rep(which.isNA.pos[1:length(which.isNA.neg)], which.isNA.neg - which.isNA.pos[1:length(which.isNA.neg)])
replacement <- c(replacement, rep(last(which.isNA.pos), last(which.isNA) - last(which.isNA.pos)))
}
replacement
}
関数fill.NAs
は次のように使用されます。
y <- c(NA, 2, 2, NA, NA, 3, NA, 4, NA, NA)
isNA <- as.numeric(is.na(y))
replacement <- fill.NAs(isNA)
if (length(replacement)){
which.isNA <- which(as.logical(isNA))
to.replace <- which.isNA[which(isNA==0)[1]:length(which.isNA)]
y[to.replace] <- y[replacement]
}
出力
> y
[1] NA 2 2 2 2 3 3 3 4 4 4
動作しているようです... しかし、男、それは醜いです!助言がありますか?
fill
ていますR
tidyr::fill()
。
roll=TRUE
たと思いますdata.table
。