我希望从第一个元素开始提取向量的递增子序列。例如,从这个向量:
a = c(2, 5, 4, 0, 1, 6, 8, 7)
...我想回复:
res = c(2, 5, 6, 8)
。
我以为我可以使用循环,但我想避免它。 sort
的另一次尝试:
a = c(2, 5, 4, 0, 1, 6, 8, 7)
ind = sort(a, index.return = TRUE)$ix
mat = (t(matrix(ind))[rep(1, length(ind)), ] - matrix(ind)[ , rep(1, length(ind))])
mat = ((mat*upper.tri(mat)) > 0) %*% rep(1, length(ind)) == (c(length(ind):1) - 1)
a[ind][mat]
基本上我对输入向量进行排序并检查索引是否验证条件“右侧没有索引较低”,这意味着事先没有更大的值。
但它似乎有点复杂,我想知道是否有更容易/更快的解决方案,或R中的预建功能。
由于
答案 0 :(得分:10)
一种可能性是找到向量的累积最大值,然后提取唯一元素:
unique(cummax(a))
# [1] 2 5 6 8
答案 1 :(得分:1)
另一个答案更好,但我做了这个迭代功能,也可以。它通过使所有连续的差异> 0
increasing <- function (input_vec) {
while(!all(diff(input_vec) > 0)){
input_vec <- input_vec[c(1,diff(input_vec))>0]
}
input_vec
}