如何在嵌套数据框中提取t检验结果的p值生成sig变量

假设你已经通过嵌套分组和t.test()得到了包含t_test列的dataNested数据框，以下是几种提取p.value生成sig列的方法：

方法1：使用tidyverse工具（dplyr + purrr） #

如果t_test列的每个元素是t.test()返回的htest对象，或是broom::tidy()整理后的结果表，可以用map_dbl()直接提取：

library(dplyr)
library(purrr)

# 提取p.value生成sig列
dataNested <- dataNested %>%
  mutate(sig = map_dbl(t_test, ~ .x$p.value))

# 更简洁的写法（当t_test元素是带p.value的结构时）
dataNested <- dataNested %>%
  mutate(sig = map_dbl(t_test, "p.value"))

方法2：基础R实现 #

如果不想加载tidyverse包，用基础R的lapply()结合unlist()也能实现：

dataNested$sig <- unlist(lapply(dataNested$t_test, function(test_result) test_result$p.value))

说明 #

• 上述代码会遍历t_test列的每个检验结果对象，提取其中的p.value数值，作为新的sig列添加到dataNested中。
• 如果你的t_test列是通过broom::tidy()处理过的结构化结果（每个元素是包含p.value列的小数据框），两种方法都能正常工作。

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