我有一个数字序列,我想检查并拒绝不以类似水平增加的数字。
data <- c(1, 2, 5, 6, 6.25, 6.49, 6.75, 7.01, 8, 9, 11)
例如,对于此处的数据,我想增加 0.25+/-0.1 的子集数字,并拒绝不遵循此规则的数字。在这种情况下,子集序列将是 (6, 6.25, 6.49, 6.75, 7.01)。
这是获取要保留的索引的丑陋方法。如果您对其他截止值感兴趣,请更改0.35
myfun <- function(D) {
index <- unique(c(which(abs(diff(D)) < 0.35), which(abs(diff(D)) < 0.35)+1))
return(sort(index))
}
调用函数以获取所需的答案
data[myfun(data)]
# [1] 6.00 6.25 6.49 6.75 7.01
另一个测试
test <- c(1, 1.25, 2, 2.25, 3, 4.5, 4.75, 5, 5.45, 5.65)
test[myfun(test)]
# [1] 1.00 1.25 2.00 2.25 4.50 4.75 5.00 5.45