r - 最快的 R 实现

Question

我使用"%within[]%" <- function(x,y){x>=y[1] & x<=y[2]} (意味着 x 在紧凑集 y 中)在 R 代码中很多，但我很确定它非常慢。你有更快的事吗？它需要适用于 > 的所有内容。已定义。

编辑:x可以是一个向量并且 y按升序排列的 2 个元素向量...

EDIT2:奇怪的是没有人(据我所知)编写了一个包 rOperator快速实现C运算符如 %w/i[]%, %w/i[[%, ...

编辑3:我意识到我的问题太笼统了，因为对 x,y 做出假设会修改任何结果，我认为我们应该关闭它，感谢您的输入。

Answer 1

"%within[]%" <- function(x,y){x>=y[1] & x<=y[2]}

x <- 1:10
y <- c(3,5)

x %within[]% y
"%within[]2%" <- function(x,y) findInterval(x,y,rightmost.closed=TRUE)==1
x %within[]2% y

library(microbenchmark)

microbenchmark(x %within[]% y,x %within[]2% y)

Unit: microseconds
             expr   min    lq median    uq    max
1  x %within[]% y 1.849 2.465 2.6185 2.773 11.395
2 x %within[]2% y 4.928 5.544 5.8520 6.160 37.265

x <- 1:1e6
microbenchmark(x %within[]% y,x %within[]2% y)

Unit: milliseconds
             expr      min       lq   median       uq      max
1  x %within[]% y 27.81535 29.60647 31.25193 56.68517 88.16961
2 x %within[]2% y 20.75496 23.07100 24.37369 43.15691 69.62122

这可能是 Rcpp 的工作。