R语言中apply函数为何折叠零长度向量列表？如何规避？

问题重现 #

当使用apply处理全零矩阵时，返回结果的格式和存在非零元素的矩阵不一致：

M0 = matrix(0,10,3)
M1 = M0
M1[1,1] = 1 # M1仅第一列第一个元素为1，其余和M0一致

apply(M0,2,rep.int,x=1:10) # 返回 integer(0) —— 不符合预期
apply(M1,2,rep.int,x=1:10) # 返回 list(1,integer(0),integer(0)) —— 符合预期

预期M0的返回结果为list(integer(0),integer(0),integer(0))。

问题原因 #

apply函数内部存在自动简化逻辑：

• 当所有分组（这里是每一列）的返回结果都是长度为0的向量时，apply会直接将结果简化为一个单一的空向量（integer(0)）；
• 只要有一个分组返回非空结果，apply就无法完成向量化简化，只能保留列表结构，这也是M1返回列表的原因。

这种简化逻辑是为了在多数场景下返回更紧凑的结果，但在全空的边缘场景下会破坏格式一致性。

解决方法 #

方法1：给原函数套一层列表包装 #

让每个列的处理结果都被包裹为列表元素，这样apply会保留外层列表结构：

apply(M0, 2, function(col) list(rep.int(col, x=1:10)))
# 返回结果：[[1]] integer(0), [[2]] integer(0), [[3]] integer(0)

方法2：改用lapply处理 #

lapply始终返回列表结构，先将矩阵按列拆分后处理：

# 方法2.1：转成数据框后处理
lapply(as.data.frame(M0), rep.int, x=1:10)

# 方法2.2：用split按列拆分
lapply(split(M0, col(M0)), rep.int, x=1:10)

方法3：自定义apply包装器 #

如果需要频繁使用这种“强制返回列表”的逻辑，可以写一个专用包装器：

apply_list <- function(X, MARGIN, FUN, ...) {
  # 先让每个分组返回列表，再合并成最终列表
  result <- apply(X, MARGIN, function(x) list(FUN(x, ...)))
  # 解开内层的单元素列表，保持结构一致
  lapply(result, unlist, recursive = FALSE)
}

# 使用示例
apply_list(M0, 2, rep.int, x=1:10)
apply_list(M1, 2, rep.int, x=1:10)

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