сводится к классической задаче «заполнение отсутствующих данных предыдущим значением».
вот тут много разных вариантов приведено. есть даже сравнение их производительности.
понятно, что на первом месте ffi-функция, написанная на компилируемом языке. но можно взять и следующий по производительности вариант — в ответах функция названа na.lomf()
.
т.е., основываясь на ваших данных, имеется вот такой вектор, полученный из исходного:
> a
[1] NA NA NA 1 1 1 NA NA NA NA 0 0
где значения меньше 1.8 заменены на единицы, а значения больше нуля заменены на нули:
> z <- c(-1,-1.4,-1.1,-1.8,-2.2,-2.5,-1.7,-1,-0.5,-0.1,0.2,0.4)
> a <- ifelse(z <= -1.8, 1, ifelse(z > 0, 0, NA))
> a
[1] NA NA NA 1 1 1 NA NA NA NA 0 0
применяя к нему ту самую функцию:
# last observation moved forward
# replaces all NA values with last non-NA values
na.lomf <- function(x) {
na.lomf.0 <- function(x) {
non.na.idx <- which(!is.na(x))
if (is.na(x[1L])) {
non.na.idx <- c(1L, non.na.idx)
}
rep.int(x[non.na.idx], diff(c(non.na.idx, length(x) + 1L)))
}
dim.len <- length(dim(x))
if (dim.len == 0L) {
na.lomf.0(x)
} else {
apply(x, dim.len, na.lomf.0)
}
}
получаем (почти) то, что вам требуется:
> na.lomf(a)
[1] NA NA NA 1 1 1 1 1 1 1 0 0
теперь оставшиеся NA
можно заменить на требующиеся вам нули каким-нибудь способом.
хотя можно обойтись и без этой замены, требующей доп. вычислений, если заменить первый элемент a
на ноль перед вызовом na.lomf()
:
> a[1] <- 0
или, более универсально, заменить на ноль только если там NA
:
> a[1] <- ifelse(is.na(a[1]), 0, a[1])
теперь получим уже ровно то, что требуется:
> a
[1] 0 NA NA 1 1 1 NA NA NA NA 0 0
> na.lomf(a)
[1] 0 0 0 1 1 1 1 1 1 1 0 0