清洁矢量与发生在与像下面的向量提交系列

Question

元件，

vec01 <- c(1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 1, 2, 1, 
      2, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 2, 2, 
      1, 2, 2, 1, 1, 2, 3, 4, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 1, 1, 1, 
      2, 1, 2, 3, 4, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 
      1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1) 

问题1：怎样去除异常下面突出显示：

vec01 <- c(1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 1, 2, 1, 
      2, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, *2*, *2*, 
      1, 2, *2*, 1, 1, 2, 3, 4, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 1, 1, 1, 
      2, 1, 2, 3, 4, *2*, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 
      1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1) 

异常的定义：该元件需要是部分系列1,2，...以上加粗标记的部分

问题2：如何识别删除异常后的系列组，每个序列属于一个组i。即，像

result <- structure(list(vec = c(1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 3L, 
           4L, 5L, 6L, 7L, 8L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
           1L, 1L, 1L, 2L, 1L, 2L), 
         group = c(1L, 1L, 2L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 8L, 
            8L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 10L, 10L, 11L, 
            11L, 12L, 13L, 14L, 15L, 15L, 16L, 16L)), 
         .Names = c("vec", "group"), 
         row.names = c(NA, 30L), class = "data.frame") 

Answer 1

有趣的问题取代is.unsorted，这里是另一种解决方案。它看起来值正在增加，并建立相应的理想（无异常）序列vec02。那么这只是比较vec01和vec02的问题。

is.incr <- c(FALSE, diff(vec01) == 1) 
vec02 <- rep(1, length(vec01)) + sequence(rle(is.incr)$lengths) * is.incr 
vec  <- vec01[vec01 == vec02] 
result <- data.frame(vec = vec, group = cumsum(vec == 1)) 

Answer 2

一个输出清洁序列（问题1）：

m <- vec01[1]==1 
for (i in seq(2,length(vec01))) 
    m[i] <- vec01[i]==1 || vec01[i]==vec01[i-1]+1 && m[i-1] 
vec01 <- vec01[m] 

现在使你的期望结构（感谢@statquant为cumsum()想法）：

data.frame(vec=vec01, group=cumsum(c(1,diff(vec01)!=1))) 

Answer 3

这里是问题2（问题1，如果你最后删除全部TRUE）

library(data.table) #load data.table because syntax is nice (matter of pers taste) 
DT = data.table(vec01) 
DT[,counter:=ifelse(vec01==1,1,0)] #identify each sequence starting with one 
DT[,counter:=cumsum(counter)] #trick to give a diff ID to each seq so we can use by 
DT[,flag:=is.unsorted(vec01),by=counter] #check sorting for each sequence 

编辑：由f(vec01)与f = function(x){!(x==Reduce(max,x,accumulate=T))}