如何从另一个向量中删除向量中的项目

Question

我有两个向量：

a = c(1,1,2,2,3,3,4,4) 
b = c(1)

我想从a 中删除b 的第一个 匹配项。因此，这里只删除了a 中的第一个1：

c = c(1,2,2,3,3,4,4)

a 中项目的顺序并不重要。

我试过这段代码：

a[a != b]
a[! a %in% b] 

两个结果都是：

[1] 2 2 3 3 4 4.

所有数字 1 都被删除。但是，我只想从a 中删除b 中的特定项目。

如果b = c(1, 1, 2)，那么我希望结果

[1] 2 3 3 4 4

a[-(1:3)]

上面的代码可能会导致[1] 2 3 3 4 4的结果。但是，我希望它可以更灵活。例如，当项目的顺序未知或随机时：

a = c(3,4,3,1,2,2,1,4)

我如何使用 R 来做到这一点？

Answer 1

vecsets 包可以执行标准集合操作，同时保留重复：

vecsets::vsetdiff( c(1,1,2,2,3,3,4,4), c(1) )
## [1] 1 2 2 3 3 4 4

vecsets::vsetdiff( c(1,1,2,2,3,3,4,4), c(1,1,2) )
## [1] 2 3 3 4 4

请注意，它将保留第一个参数的顺序。使用你的最后一个例子：

vecsets::vsetdiff( c(3,4,3,1,2,2,1,4), c(1,1,2) )
## [1] 3 4 3 2 4

Answer 2

从this answer 中获取灵感到我在评论中链接的一个问题，您可以使用data.table 包中的fsetdiff。
它将all 作为参数，这样可以避免只返回唯一值，就像setdiff 一样：

library(data.table)

# with your first example (b = c(1)):
unlist(fsetdiff(data.table(v1=a), data.table(v1=b), all = TRUE))
# v11 v12 v13 v14 v15 v16 v17 
#  1   2   2   3   3   4   4

# with second example (b = c(1, 1, 2)):
unlist(fsetdiff(data.table(v1=a), data.table(v1=b), all = TRUE))
# v11 v12 v13 v14 v15 
#  2   3   3   4   4

Answer 3

您可以使用which()

a = c(3, 4, 3, 1, 2, 2, 1, 4)
a
## [1] 3 4 3 1 2 2 1 4

b = 1

a[- which(a %in% b)[1]]
## [1] 3 4 3 2 2 1 4

案例b有两个元素：

b2 = c(1, 2)

sapply(seq_along(b1), function(x) a <<- a[- which(a == x)[1]])[[2]]
## [1] 3 4 3 2 1 4

或者三个...

b3 <- c(1, 2, 3)

sapply(seq_along(b1), function(x) a <<- a[- which(a == x)[1]])[[3]]
# [1] 4 3 2 1 4

Answer 4

我不认为以下是最好的解决方案（vecsets 方法让我觉得最好），但@Aaron 关于可能使用 Rcpp 的评论让我觉得很有趣。这是我第一次使用那个包。如果不出意外，我能够在不到 20 分钟的时间内获得工作代码这一事实强调了他的观点，即 Rcpp 使它相对容易：

library(Rcpp)
cppFunction('
  NumericVector difference(NumericVector xs, NumericVector ys){
    int m = xs.size();
    int n = ys.size();
    float flag = 1 + abs(max(xs)) + abs(max(ys)); //occurs in neither xs nor ys
    NumericVector zs = clone(xs);
    for(int i = 0; i < n; i++){
      double y = ys[i];
      int j = 0;
      while(j < m && zs[j]!= y) j++;
      if(j < m) zs[j] = flag;
    }
    int count = 0;
    for(int k = 0; k < m; k++){
      if(zs[k] < flag) count++;
    }
    NumericVector ws(count);
    int k = 0;
    for(int j = 0; j < m; j++){
      if(zs[j] < flag){
        ws[k] = zs[j];
        k++;
      }
    }
    return ws;
  }
')

在你获得这个之后：

> a = c(1,1,2,2,3,3,4,4)
> b = c(1,2,1)
> difference(a,b)
[1] 2 3 3 4 4

由于这是我第一次尝试编写此类代码，我确信它可以通过多种方式进行改进。

Answer 5

语法顺序有点令人沮丧，但 Reduce 和 which 只使用 Base R。

Reduce(b, a) a[-which(a==b)[1]], a, b)