R语言：列表内矩阵行按指定滞后阶数列和做除法的实现问题

我有一个包含100个3×51矩阵的列表，需要对每个矩阵执行以下操作：

• 第一行元素除以该矩阵列和的2阶滞后值
• 第二行元素除以该矩阵列和的3阶滞后值
• 第三行元素除以该矩阵列和的4阶滞后值
  用计算结果替换原矩阵对应行，得到新的列表数据集。

我已经掌握了向量列表的类似操作，代码如下：

Example_List2 <- lapply(Example_List1,function(x){x/lag(x,n=2)})

我的尝试代码存在问题，无法正确实现需求：

List2 <- List1
lapply(List2, `[`,1,) <- lapply(List1,function(x){lapply(x, `[`,1,)/lag(colSums(x),n=2)})
lapply(List2, `[`,2,) <- lapply(List1,function(x){lapply(x, `[`,2,)/lag(colSums(x),n=3)})
lapply(List2, `[`,3,) <- lapply(List1,function(x){lapply(x, `[`,3,)/lag(colSums(x),n=4)})

正确实现方法 #

问题出在嵌套lapply的错误使用，以及试图通过lapply直接赋值修改列表元素。正确的做法是在lapply内部对每个矩阵完成完整的行修改操作后返回新矩阵，最终生成新列表：

# 若使用dplyr的lag函数（避免base R中lag对非时间序列的特殊处理）
library(dplyr)

List2 <- lapply(List1, function(mat) {
  # 计算当前矩阵的列和
  col_totals <- colSums(mat)
  
  # 计算各阶滞后的列和，前n个位置会生成NA，可根据需求调整（比如用0填充或保留NA）
  lag_2 <- lag(col_totals, n = 2)
  lag_3 <- lag(col_totals, n = 3)
  lag_4 <- lag(col_totals, n = 4)
  
  # 替换对应行的元素
  mat[1, ] <- mat[1, ] / lag_2
  mat[2, ] <- mat[2, ] / lag_3
  mat[3, ] <- mat[3, ] / lag_4
  
  return(mat)
})

如果不想依赖dplyr，可以用base R手动实现滞后效果：

List2 <- lapply(List1, function(mat) {
  col_totals <- colSums(mat)
  # 手动构造滞后向量：前n个为NA，后面取前(length(col_totals)-n)个元素
  lag_2 <- c(rep(NA, 2), col_totals[1:(length(col_totals)-2)])
  lag_3 <- c(rep(NA, 3), col_totals[1:(length(col_totals)-3)])
  lag_4 <- c(rep(NA, 4), col_totals[1:(length(col_totals)-4)])
  
  mat[1, ] <- mat[1, ] / lag_2
  mat[2, ] <- mat[2, ] / lag_3
  mat[3, ] <- mat[3, ] / lag_4
  
  return(mat)
})

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