多亚组t检验实现求助：基于R语言处理分层实验数据

用你已经在用的rstatix包就能完成需求，直接按性别、年龄分层后，批量对三个指标做t检验，自动保留分组标注，代码如下：

library(rstatix)
library(dplyr)

# 执行分层t检验：按sex、age分组，每个组内做control vs drug的t检验
subgroup_tests <- big.df %>%
  group_by(sex, age) %>%
  t_test(c(height, weight, speed) ~ trial) %>%
  # 可选：筛选并整理结果列，保留关键信息
  select(sex, age, variable, p, p.signif) %>%
  arrange(variable, sex, age)

# 查看结果
subgroup_tests

结果说明 #

运行后会得到一个数据框，包含以下关键列：

• sex：性别（m/f）
• age：年龄组（young/old）
• variable：检测指标（height/weight/speed）
• p：t检验的p值
• p.signif：显著性标记（比如*代表p<0.05，**代表p<0.01等）

正好对应你需要的12个分组对比结果（4个亚组 × 3个指标），所有p值的分组信息清晰标注。如果需要查看完整的t检验统计量（如t值、自由度），删掉代码里的select(...)行即可。

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