-4

我有一个数字序列,我想检查并拒绝不以相似级别增加的数字。

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)。

4

1 回答 1

-1

这是获取要保留的索引的一种丑陋方式。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
于 2017-10-05T19:44:41.073 回答