如何根据dat1取值规则将dat2对应值除以dat3对应列值？

这是一个典型的元素级匹配运算需求，我们可以利用R的矩阵索引特性高效实现，以下是具体步骤和代码：

核心思路 #

dat1的每个元素值对应dat3的列名，我们先把dat1的取值转化为dat3的列索引，提取出对应位置的dat3数值，最后让dat2和这些数值做逐元素除法即可。

步骤1：统一数据格式（可选） #

如果你的实际数据中dat3的列名是数值型（而非示例中的字符型"1""2""3"），先转成字符型避免索引错误：

colnames(dat3) <- as.character(colnames(dat3))

如果dat1/dat2是数据框，建议转为矩阵提升运算效率：

dat1 <- as.matrix(dat1)
dat2 <- as.matrix(dat2)
dat3 <- as.matrix(dat3)

步骤2：实现运算 #

直接用向量化索引提取dat3的对应值，再和dat2做除法：

# 提取dat3中匹配的除数：每行对应dat1的行，每列对应dat1元素值对应的dat3列
divisor_values <- dat3[cbind(seq(nrow(dat1)), match(as.character(dat1), colnames(dat3)))]
# 把除数转换成和dat2同维度的矩阵，保留列名
divisor_matrix <- matrix(divisor_values, nrow = nrow(dat1), dimnames = dimnames(dat2))
# 计算最终结果
result <- dat2 / divisor_matrix

验证示例数据 #

把你提供的示例数据代入，验证结果：

dat1 = cbind(a = c(2,1,3), b = c(2,2,3), c = c(2,2,2))
dat2 = cbind(a = c(20, 20, 20), b = c(30, 30, 30), c = c(50, 50, 50))
dat3 = cbind("1" = c(100,100,100), "2" = c(200, 200, 200), "3" = c(300, 300, 300))

# 执行运算
divisor_values <- dat3[cbind(seq(nrow(dat1)), match(as.character(dat1), colnames(dat3)))]
divisor_matrix <- matrix(divisor_values, nrow = nrow(dat1), dimnames = dimnames(dat2))
result <- dat2 / divisor_matrix

# 查看格式化后的结果
round(result, 3)

输出完全符合你的目标：

a     b     c
[1,] 0.100 0.150 0.25
[2,] 0.200 0.150 0.25
[3,] 0.067 0.067 0.25

简化写法 #

如果想少写几行，也可以把运算合并成一行：

result <- dat2 / matrix(dat3[cbind(seq(nrow(dat1)), match(as.character(dat1), colnames(dat3)))], 
                        nrow = nrow(dat1), 
                        dimnames = dimnames(dat2))

注意事项 #

• 确保dat1的所有取值都存在于dat3的列名中，否则match会返回NA，导致结果出现缺失值
• 若实际数据是大型矩阵，这种向量化方法比循环更高效

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