如何矢量化模数?

时间:2017-12-28 10:46:00

标签: r vectorization modulo

R中,我有一个带整数的向量,对于多个值,它应被检查为可分割而没有余数。

v <- c(3, 4, 10, 11, 12)
v %% 4 == 0 | v %% 10 == 0 #... this might include 100 checks or more

我正在寻找一个这样的矢量化(如v %% c(4, 10)),这可能会返回一个矩阵或类似的东西。

我编写了一个可以执行此操作的函数,但我想知道是否已经存在这样一个我可以重用的函数,而不是再次“重新发明”它。

`%V%`  <- function(x, y){
  ret <- matrix(nrow = length(x), ncol = length(y))
  for(i in 1:ncol(ret)){
    ret[, i] <- x %% y[i]
  }
  return(ret)
}
v %V% c(4, 10)

出于我最初的目的,我现在将运行apply(v %V% c(4, 10), 1, function(x) any(x == 0))或相应地更改我的功能。

1 个答案:

答案 0 :(得分:3)

您可以使用outer - 功能来实现您的目标:

outer(v, c(4,10), '%%') == 0

给出:

      [,1]  [,2]
[1,] FALSE FALSE
[2,]  TRUE FALSE
[3,] FALSE  TRUE
[4,] FALSE FALSE
[5,]  TRUE FALSE