如何在R中检查矩阵各元素是否处于另外两个对应矩阵元素的区间内？

这是个非常典型的矩阵逐元素区间校验需求，尤其是处理大规模矩阵时，向量化操作绝对是最优解——R的底层用C/Fortran实现了向量化运算，比手动写循环快几个数量级，内存效率也高得多。

核心思路 #

对于每个位置(i,j)，我们需要判断M[i,j]是否落在L[i,j]和U[i,j]构成的区间内（注意：L不一定是区间下限，U也不一定是上限，所以要先取两者的最小/最大值作为区间边界），符合条件记1，不符合记0。

一行代码搞定（基础R实现） #

# 生成0-1矩阵
M_comparisons <- as.integer(M >= pmin(L, U) & M <= pmax(L, U))
# 恢复矩阵维度（as.integer会把矩阵转成向量，所以要重置维度）
dim(M_comparisons) <- dim(M)

代码拆解 #

• pmin(L, U)：逐元素取L和U的较小值，得到每个位置的区间下限
• pmax(L, U)：逐元素取L和U的较大值，得到每个位置的区间上限
• M >= pmin(...) & M <= pmax(...)：逐元素判断是否在区间内，生成一个逻辑矩阵（TRUE/FALSE）
• as.integer()：把逻辑值转成0/1（TRUE→1，FALSE→0）
• 最后用dim()重置维度，确保结果还是和原矩阵同维度的矩阵

完整示例（匹配你的测试数据） #

# 构造你提供的测试矩阵
M <- matrix(c(-0.1278982, 0.4600544, 1.3271033,
              -0.1079272, 1.2196851, -0.5240536,
              0.9548449, -0.9979668, -0.6636296),
            nrow = 3, byrow = TRUE)

L <- matrix(c(-0.49416022, 1.0992247, -1.2834469,
              0.35491793, 0.2460946, 1.1411529,
              0.01755317, 0.2469062, 0.2843848),
            nrow = 3, byrow = TRUE)

U <- matrix(c(-0.84996737, 1.0036554, -0.04747612,
              0.35396050, -0.4612527, -1.47341334,
              -0.07716839, -1.5921625, 1.62187501),
            nrow = 3, byrow = TRUE)

# 执行判断
M_comparisons <- as.integer(M >= pmin(L, U) & M <= pmax(L, U))
dim(M_comparisons) <- dim(M)

# 查看结果
M_comparisons

运行后输出完全符合你的预期：

[,1] [,2] [,3]
[1,]    0    0    0
[2,]    0    0    1
[3,]    0    1    0

大规模矩阵的优化说明 #

• 这个方案完全是向量化操作，没有任何R层面的循环，处理几万行几万列的矩阵也能保持高效——因为所有运算都在底层C代码中完成，避免了R循环的巨大开销。
• 如果矩阵大到超出内存，可以考虑用Matrix包的稀疏矩阵类型存储结果，进一步节省内存；但常规规模下，基础R的实现已经足够。

内容的提问来源于stack exchange，提问作者user321627