x=c(NA, 2, -3, NA, -5, 5, -7, -8, -9, -10, NA, -2, 2, -14, -15, NA, -17, 2, NA, -20)
我想选择前 10 个不是NA
. 所以我这样做:
head(x[!is.na(x)], 10)
现在我想知道有多少 NA 被忽略了:
sum(is.na(head(x[is.na(x)], 10)))
# [1] 5
这是不正确的:应该是 3
我们需要在实际向量的子集上而不是仅在 NA 数据上的 logicla 向量
sum(is.na(head(x, 10)))
#[1] 3
head(x[is.na(x)], 10)
#[1] NA NA NA NA NA
仅给出 NA,因为我们正在对 NA 元素进行子集化
根据评论中的逻辑,一种方法是
i1 <- !is.na(x)
i1[i1] <- seq_along(i1[i1])
sum(is.na(head(x, which(i1 == 10))))
这应该这样做:
sum(is.na(head(x, which(cumsum(!is.na(x)) == 10))))
给出的cumsum
逻辑向量的!is.na
给我们第十个非NA
数字的索引,这是我们想要的head
函数的输入(因为我们想要所有的值,直到我们得到第十个非NA
值)。