R语言dcast特殊场景：将control组转为新列以开展多组t.test

嗨，我太懂这种“明明看起来简单却卡半天”的烦躁感了！咱们先把dcast的问题解决，再顺带着完成批量t检验的需求，一步一步来：

第一步：先明确问题根源 #

你之前用dcast没成功，大概率是因为没有指定行ID变量，导致dcast默认对数据做了聚合（比如取均值），而不是保留所有观测值。咱们先模拟一份和你类似的数据，再演示正确操作：

library(tidyverse)
library(data.table) # dcast是data.table包的函数，先加载它

# 模拟你的数据：group包含control和其他组，x是数值变量
set.seed(123)
df <- tibble(
  group = rep(c("control", "A", "B", "C"), each = 10), # 每个组10个观测
  x = rnorm(40, mean = c(5, 6, 4.5, 6.2), sd = 1)
)

第二步：用dcast正确转宽格式 #

先给每个组内的观测添加行号，这样dcast就知道如何把不同组的观测一一对应到同一行：

# 转换为data.table并添加组内行号
dt <- as.data.table(df)
dt[, row_id := seq_len(.N), by = group] # 每个组内生成1到10的行号

# 用dcast转宽：行号作为行，group作为列，x作为值
wide_df <- dcast(dt, row_id ~ group, value.var = "x")

# 转完后的数据结构：row_id + control + A + B + C列，每行是对应组的第n个观测
head(wide_df)

第三步：批量做t检验 #

现在wide_df里已经有了control列和其他组的列，咱们可以用purrr批量遍历每个非control列，和control做t检验：

# 提取所有需要检验的组名（排除control）
test_groups <- names(wide_df)[names(wide_df) != "control"]

# 批量执行t检验，结果存入列表
t_test_results <- map(test_groups, function(group_name) {
  t.test(wide_df[[group_name]], wide_df[["control"]])
})

# 给结果列表命名，方便查看
names(t_test_results) <- test_groups

# 查看任意组的结果，比如A组
t_test_results$A

额外技巧：不用转宽也能直接做 #

其实你完全可以跳过转宽的步骤，直接用dplyr+purrr批量处理，效率更高：

# 先提取对照组的数据
control_x <- df %>% filter(group == "control") %>% pull(x)

# 对每个非对照组单独做t检验
t_test_results2 <- df %>%
  filter(group != "control") %>%
  group_split(group) %>% # 按组拆分数据框
  set_names(map(., ~ .x$group[1])) %>% # 给每个拆分后的数据集命名
  map(function(group_data) {
    t.test(group_data$x, control_x)
  })

# 查看B组的结果
t_test_results2$B

这样不管是转宽还是直接处理，都能完成你要的“每个组和control做t检验”的需求啦！

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