当尝试从序列中排除特定数字时,我注意到which(!(0:10 %in% 2:3))
返回的结果(删除 3 和 4 并添加 11(与setdiff(0:10,2:3)
或which(!(1:10 %in% 2:3))
不同。
which(!(1:10 %in% 2:3))
[1] 1 4 5 6 7 8 9 10
which(!(0:10 %in% 2:3))
[1] 1 2 5 6 7 8 9 10 11
setdiff(0:10,2:3)
[1] 0 1 4 5 6 7 8 9 10
这似乎是一个简单的逻辑问题,但我无法弄清楚这是由于什么原因?对于大序列,setdiff
是否与which(!())
一样快?
总结注释,which
给出了元素的位置,而不是元素本身,上面的例子不清楚,但这些例子很明显:
which(!(10:20 %in% 12:13))
[1] 1 2 5 6 7 8 9 10 11
0:10 %in% 2:3
[1] FALSE FALSE TRUE TRUE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
因此,选择元素的正确方法是:
myseq<-0:10
myseq[which(!(0:10 %in% 2:3))]
关于速度,直接从矢量中选择元素稍微快一些:
myseq<-0:1000000
> microbenchmark::microbenchmark(
+ wh=myseq[which(!(0:1000000 %in% 2:3))],
+ sd=setdiff(0:1000000,2:3),
+ seq=(0:1000000)[!0:1000000 %in% 2:3]
+ )
Unit: milliseconds
expr min lq mean median uq max neval cld
wh 18.62157 18.85489 25.17644 25.37830 26.89162 152.9487 100 a
sd 36.09655 42.83383 50.22088 44.16595 45.96227 178.2949 100 b
seq 17.51332 17.98346 25.00993 24.39265 25.91137 174.1875 100 a