R语言对0×0矩阵调用apply触发FUN(numeric(0))的原因及规避方法

这个现象不属于代码逻辑错误，是R内置apply底层C实现的边界处理疏漏。
异常触发条件非常特殊：只有当传入的矩阵/数组两个维度同时为0（即dim属性为c(0,0)的空矩阵）时，按指定MARGIN遍历才会出现多余的FUN调用。
其余常规空矩阵场景，比如0行3列、3行0列，按对应边执行apply时，会正确识别可迭代对象数量为0，不会触发FUN执行；但0x0矩阵场景下，底层拆分遍历切片的逻辑没有正确计算可迭代元素总数为0，反而会生成一个内容为numeric(0)的空切片，强行传入FUN执行一次，这就是观察到f被触发输出空向量的直接原因。

无法通过修改apply的输入参数解决。该问题属于底层实现的边界疏漏，apply的所有公开参数（包括simplify配置、通过...传入的额外参数）都没有提供控制“空维度下是否执行试探性调用”的逻辑，无论如何传参，只要输入0x0矩阵搭配对应MARGIN，都会触发那次多余的FUN调用。

零散的if判断跳过apply的思路本身没有问题，但可以优化成通用可复用的逻辑，避免重复写判断：

• 封装带前置校验的安全apply函数：在进入原apply执行流程前，先校验目标遍历维度的长度，只要长度为0就直接返回和原apply逻辑结构一致的空结果，完全绕开底层存在边界问题的C实现路径。参考实现：

safe_apply <- function(X, MARGIN, FUN, ...) {
  # 获取目标遍历维度的长度
  target_dim_len <- dim(X)[MARGIN]
  # 维度长度为0时直接返回空结果，不进入apply主流程
  if (is.null(target_dim_len) || target_dim_len == 0) {
    return(X[0])
  }
  apply(X, MARGIN, FUN, ...)
}

• 如果不方便批量替换原有apply的调用点，也可以在每个自定义FUN的入口增加空值判断，遇到长度为0的入参时直接返回空值，避免后续业务逻辑因为入参为空崩溃。但这种方式维护成本更高——每个用到的FUN都要单独加判断，可靠性远不如统一的前置维度校验。

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