使用R语言apply函数计算单行矩阵行和报错的解决咨询

嘿，这个问题我之前也踩过坑！咱们先搞清楚为什么会报错，再给你几个好用的解决方案：

报错原因 #

当你的矩阵只有1行时，mat[,1:15]这个操作会被R自动简化成一个向量（R的默认行为，为了节省内存），而apply()函数要求输入的对象必须有明确的维度（也就是dim()属性要有值），向量是没有维度的，所以就会抛出dim(X) must have a positive length的错误。而当矩阵有2行及以上时，mat[,1:15]会保持矩阵结构，所以apply()能正常工作。

解决方案 #

方案1：用drop=FALSE强制保持矩阵结构 #

在提取矩阵列的时候加上drop=FALSE，就能阻止R把单行矩阵简化成向量，这样apply()就能正常识别维度了：

mat <- matrix(rnorm(15), 1, 15)
apply(mat[,1:15, drop=FALSE], 1, sum)

方案2：直接用sum()函数（最简单） #

因为单行矩阵本质上和向量的元素是一样的，直接用sum()就能得到行和，完全不需要apply()：

sum(mat)

方案3：推荐使用rowSums()函数（最专业） #

R专门提供了rowSums()函数来计算矩阵的行和，它不管矩阵是单行还是多行，都能直接输出正确结果，而且比apply()的效率更高（底层是优化过的C代码）：

rowSums(mat)

验证示例 #

你可以运行下面的代码测试三种方法：

# 创建单行矩阵
mat_single <- matrix(rnorm(15), 1, 15)
# 创建多行矩阵
mat_multi <- matrix(rnorm(30), 2, 15)

# 单行矩阵测试
apply(mat_single[,1:15, drop=FALSE], 1, sum)
sum(mat_single)
rowSums(mat_single)

# 多行矩阵测试（三种方法都适用）
apply(mat_multi[,1:15], 1, sum)
rowSums(mat_multi)

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