R语言中自动化批量生成双向频数表的实现方案

针对这类“固定行变量集合、固定列变量集合，批量生成所有两两配对交叉频数表”的场景，不需要手动逐行写table()语句，完全可以通过组合apply族函数/循环实现自动化，和你之前用mapply()做单变量频数统计的逻辑一脉相承，没有额外的学习成本。

方案1：基础R实现（零额外依赖，和你现有使用习惯匹配） #

• 先定义需要参与交叉的变量集合：
  • 行变量（你示例里的cyl、cyl2、cyl3）：row_vars <- c("cyl", "cyl2", "cyl3")
  • 列变量（你示例里的vs、am）：col_vars <- c("vs", "am")
• 生成所有变量两两配对的组合：var_pairs <- expand.grid(row_vars, col_vars, stringsAsFactors = FALSE)
• 批量调用table()生成所有交叉表，结果存为列表方便后续读取：

  batch_cross_tables <- mapply(
    function(row_col, col_col) {
      # dnn参数自动给表设置维度名，避免分不清哪个变量对应哪个维度
      table(
        dat[[row_col]], 
        dat[[col_col]], 
        dnn = c(row_col, col_col),
        useNA = "ifany" # 保留缺失值统计，和手动运行结果保持一致
      )
    },
    var_pairs$Var1,
    var_pairs$Var2,
    SIMPLIFY = FALSE
  )
  

• 直接打印batch_cross_tables就能看到所有生成的交叉表，每个元素自动标注对应的变量名，不会出现手动运行时混淆表对应关系的问题。

方案2：tidyverse生态实现（更易梳理，适合大规模变量批量处理） #

如果平时常用dplyr/purrr系列的工具，可以用更结构化的方式存储结果，后续导出、筛选都更方便：

library(purrr)
library(dplyr)

batch_cross_tables <- expand.grid(
  row_var = c("cyl", "cyl2", "cyl3"),
  col_var = c("vs", "am"),
  stringsAsFactors = FALSE
) %>% 
  mutate(
    cross_table = map2(
      .x = row_var,
      .y = col_var,
      .f = ~table(dat[[.x]], dat[[.y]], dnn = c(.x, .y), useNA = "ifany")
    )
  )

• 这种方式生成的结果是结构化数据框，每一行对应一组变量交叉组合，旁边存储对应的频数表，比纯列表的可读性高很多。
• 如果需要把所有结果汇总成长格式方便后续透视、导出，可以追加一段转换代码：

  batch_long_result <- batch_cross_tables %>% 
    group_by(row_var, col_var) %>% 
    reframe(as.data.frame.table(cross_table[[1]]))
  

效率说明 #

后续如果需要增加新的行变量或列变量，只需要修改最开头row_vars、col_vars两个向量的取值即可，不需要重复编写大量table()语句，同时所有结果自动绑定对应变量名，完全解决手动运行效率低、结果难梳理的问题。

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