如何将含混合维度的嵌套表格展开为频率表数组？

需求说明 #

将一个嵌套表格（行标识为a/b，列标识为v1-v7）转换为一个包含7个频率表的多维数组：

• 每个频率表是2×2矩阵，行固定为0.1、0.3，列固定为a、b，单元格存储频次计数
• 输入表格中多数单元格为2维table对象，仅a行的v7列是单个数值

测试数据生成 #

先构造符合要求的嵌套测试数据：

set.seed(123)
# 生成随机2维频次表的工具函数
create_table <- function() {
  table(
    sample(c(0.1, 0.3), 10, replace = TRUE),
    sample(c("a", "b"), 10, replace = TRUE)
  )
}

# 构建嵌套表格：行a/b，列v1-v7，多数单元格为table，a行v7为数值
nested_df <- data.frame(
  row = c("a", "b"),
  v1 = I(list(create_table(), create_table())),
  v2 = I(list(create_table(), create_table())),
  v3 = I(list(create_table(), create_table())),
  v4 = I(list(create_table(), create_table())),
  v5 = I(list(create_table(), create_table())),
  v6 = I(list(create_table(), create_table())),
  v7 = I(list(5, create_table())) # a行v7为数值，b行v7为table
)

转换实现代码 #

# 标准化函数：将输入（数值或table）转换为统一结构的2×2矩阵
standardize_matrix <- function(x) {
  # 处理数值类型：生成默认矩阵，将数值填充到(0.1, a)位置，其余为0（可按需调整）
  if (is.numeric(x)) {
    mat <- matrix(0, nrow = 2, ncol = 2, 
                  dimnames = list(c("0.1", "0.3"), c("a", "b")))
    mat["0.1", "a"] <- x
    return(mat)
  }
  
  # 处理table类型：转换为矩阵并补全缺失的行列（缺失位置填0）
  mat <- as.matrix(x)
  target_rows <- c("0.1", "0.3")
  target_cols <- c("a", "b")
  
  standardized_mat <- matrix(0, nrow = 2, ncol = 2, 
                             dimnames = list(target_rows, target_cols))
  # 填充已有数据
  match_rows <- intersect(rownames(mat), target_rows)
  match_cols <- intersect(colnames(mat), target_cols)
  standardized_mat[match_rows, match_cols] <- mat[match_rows, match_cols]
  
  return(standardized_mat)
}

# 生成最终的3维数组（行×列×v列）
frequency_array <- array(
  data = 0,
  dim = c(2, 2, 7),
  dimnames = list(c("0.1", "0.3"), c("a", "b"), paste0("v", 1:7))
)

# 遍历每一列处理
for (col_num in 1:7) {
  col_name <- paste0("v", col_num)
  # 提取a行和b行的元素
  elem_a <- nested_df[nested_df$row == "a", col_name][[1]]
  elem_b <- nested_df[nested_df$row == "b", col_name][[1]]
  
  # 标准化为统一矩阵
  mat_a <- standardize_matrix(elem_a)
  mat_b <- standardize_matrix(elem_b)
  
  # 合并矩阵（这里以频次求和为例，可替换为其他逻辑如取最大值、单独存储等）
  combined_mat <- mat_a + mat_b
  
  # 存入数组对应位置
  frequency_array[, , col_num] <- combined_mat
}

# 查看结果示例（以v7为例）
frequency_array[, , "v7"]

关键说明 #

• standardize_matrix函数确保所有输入都转换为结构一致的矩阵，解决了数值和table混合的问题
• 合并逻辑可按需修改：如果不需要合并a/b行的频次，可将数组维度调整为c(2,2,2,7)（行×列×a/b行×v列），分别存储每行的频次表
• 最终数组的每个切片（frequency_array[, , "vX"]）对应vX列的合并后频率表

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