检查矩阵行是否等于 R 中的矢量（矢量化）

我很惊讶这个问题没有被问到，也许答案会澄清原因。我想将矩阵的行与向量进行比较，并返回该行是否==处处的向量。请参阅下面的示例。我想要一个矢量化解决方案，没有应用函数，因为矩阵太大而无法缓慢循环。假设也有很多行，所以我想避免重复向量。

set.seed(1)

M = matrix(rpois(50,5),5,10)

v = c(3 ,   2 ,   7  ,  7 ,   4   , 4   , 7  ,  4  ,  5, 6)

M
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]    4    8    3    5    9    4    5    6    7     7
[2,]    4    9    3    6    3    1    5    7    6     1
[3,]    5    6    6   11    6    4    5    2    7     5
[4,]    8    6    4    4    3    8    3    6    5     6
[5,]    3    2    7    7    4    4    7    4    5     6

输出应该是

FALSE  FALSE FALSE FALSE TRUE

一种可能性是

rowSums(M == v[col(M)]) == ncol(M)
## [1] FALSE FALSE FALSE FALSE  TRUE

或者类似地

rowSums(M == rep(v, each = nrow(M))) == ncol(M)
## [1] FALSE FALSE FALSE FALSE  TRUE

Or

colSums(t(M) == v) == ncol(M)
## [1] FALSE FALSE FALSE FALSE  TRUE

v[col(M)]只是一个较短的版本rep(v, each = nrow(M))它创建一个大小相同的向量M（矩阵只是一个向量，尝试c(M)）然后使用每个元素与其对应的元素进行比较==。幸运的是==是一个泛型函数，它有一个array方法（参见methods("Ops") and is.array(M)）这允许我们运行rowSums (or colSums）以确保我们的匹配数量为ncol(M)