R语言如何实现矩阵列表批量合并为大矩阵，避免逐个调用do.call

你需要实现的是按行优先排列的分块矩阵拼接，以下是两种通用实现方案，均无需手动逐个编写拼接语句，适配任意规模的矩阵列表：

方案1：基础R通用函数（无额外依赖） #

核心思路是自动将小矩阵列表按分块数分组，每组对应大矩阵的一列分块，组内先按行拼接后再整体按列拼接：

# 通用分块矩阵拼接函数
# 参数说明：
# mat_list: 存储所有小矩阵的列表，要求所有小矩阵尺寸一致
# n_block: 最终大矩阵每行/每列的分块数量，例如9个小矩阵拼3行3列分块，n_block=3
block_bind <- function(mat_list, n_block) {
  # 输入校验：小矩阵总数必须等于分块数的平方
  stopifnot(length(mat_list) == n_block^2)
  # 按分块数拆分列表，每一组对应大矩阵的一列分块
  block_col_groups <- split(mat_list, ceiling(seq_along(mat_list) / n_block))
  # 每组内按行拼接，再将所有列拼接为最终大矩阵
  do.call(cbind, lapply(block_col_groups, function(group) do.call(rbind, group)))
}

你的示例调用代码：

# 示例数据
mat_list <- list(matrix(1, 2, 2), matrix(2, 2, 2), matrix(3, 2, 2), 
                 matrix(4, 2, 2), matrix(5, 2, 2), matrix(6, 2, 2), 
                 matrix(7, 2, 2), matrix(8, 2, 2), matrix(9, 2, 2))
# 调用函数，3行3列分块
result <- block_bind(mat_list, n_block = 3)

输出结果和你预期的6×6矩阵完全一致。

方案2：高维数组实现（性能更优，适合大规模数据） #

如果你的小矩阵数量很大，推荐用数组维度重排的方法，执行效率远高于嵌套的do.call调用：

n_sub <- nrow(mat_list[[1]]) # 单个小矩阵的行数/列数，示例中为2
n_block <- 3 # 分块行列数
# 将所有小矩阵转为4维数组，维度依次为：小矩阵行、小矩阵列、分块行、分块列
arr <- array(unlist(mat_list), dim = c(n_sub, n_sub, n_block, n_block))
# 调整维度顺序后转为最终大矩阵
result <- matrix(aperm(arr, c(1, 3, 2, 4)), nrow = n_sub * n_block, ncol = n_sub * n_block)

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