r - 最快的 R 实现

Question

就目前情况而言，这个问题不太适合我们的问答形式。我们希望答案得到事实、引用资料或专业知识的支持，但这个问题可能会引发辩论、争论、民意调查或扩展讨论。如果您觉得这个问题可以改进并可能重新开放， visit the help center以获得指导。

已关闭 10 年前。

我使用"%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 的工作。