在R语言中使用mapply实现矩阵列表的逐元素除法

首先，先理清楚你的场景：你已经构造了两个列表dfZ（包含两个10×10矩阵）和dfX（包含两个10×1列向量），想要批量对每个矩阵执行每一行元素除以对应行的列向量元素的操作，也就是手动实现的A1和A2的逻辑，但用mapply时得到了空列表，这是因为你的代码存在语法错误。

先看你给出的数据代码（顺便修正一个小笔误：10×10矩阵需要100个随机数，所以把rnorm(10)改成rnorm(100)）：

set.seed(12345)
df1z <- matrix(rnorm(100), nrow=10, ncol=10)
df2z <- df1z
dfZ <- list(df1z, df2z)

df1x <- matrix(rnorm(10), nrow=10, ncol=1)
df2x <- df1x
dfX <- list(df1x, df2x)

你手动实现的操作是：

A1 = t(t(df1z)/as.vector(df1x))
A2 = t(t(df2z)/as.vector(df2x))

你的错误代码问题所在 #

你写的mapply代码有两个语法错误：

1. return()的括号没有闭合，导致参数dfZ被错误地放在了return的括号内部；
2. 函数定义结束后缺少逗号，分隔函数和后续的参数。

正确的实现代码 #

修正后的mapply代码如下，运行后就能得到包含A1和A2的列表：

# 方案1：显式使用return
result <- mapply(function(a, b) {
  return(t(t(a)/as.vector(b)))
}, dfZ, dfX, SIMPLIFY = FALSE)

# 方案2：省略return，R函数默认返回最后一行结果（更简洁）
result <- mapply(function(a, b) t(t(a)/as.vector(b)), dfZ, dfX, SIMPLIFY = FALSE)

额外小优化 #

其实t(t(a)/as.vector(b))可以简化为a / as.vector(b)，因为R的广播规则会自动处理矩阵和向量的逐行除法：

result <- mapply(function(a, b) a / as.vector(b), dfZ, dfX, SIMPLIFY = FALSE)

这个写法和你手动实现的逻辑完全一致，结果也相同。

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