R语言中使用t.test计算p-value报错的问题排查与解决

首先给出数据框结构：

mat <- structure(list(ids = c("id1", "id2", "id3", "id4", "id5", "id6", 
"id7", "id8", "id9", "id10", "id11", "id12", "id13", "id14", 
"id15", "id16", "id17", "id18", "id19", "id20"), Group = c("A", 
"A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", 
"C", "C", "C", "A", "A", "B"), number = c(2L, 2L, 3L, 2L, 2L, 
44L, 172L, 34L, 78L, 27L, 31L, 55L, 23L, 34L, 14L, 18L, 25L, 
2L, 2L, 12L)), class = "data.frame", row.names = c(NA, -20L))

错误原因解析 #

第一种写法报错：object 'Group' not found #

t.test()的参数逻辑错误：当传入data参数时，要么用公式语法，要么明确指定向量。直接写t.test(Group, number, data=mat)时，函数会把Group当作第一个数值型参数，但它是分类变量，同时函数无法正确从data中解析列（参数顺序不符合要求），导致找不到对象的报错。

第二种写法报错：grouping factor must have exactly 2 levels #

公式的变量顺序写反了，t检验的公式应该是数值变量 ~ 分组变量，写反后逻辑完全错误。更核心的问题是：你的Group列有A、B、C三个分组水平，而t检验仅支持两组之间的差异比较，这才是触发该报错的根本原因。

正确获取p-value的操作 #

因为t检验只能做两组比较，所以需要指定具体的两组进行检验，以下是具体实现：

1. 指定两组做检验（以A vs B为例） #

两种写法都可以：

# 公式写法（更简洁）
t.test(number ~ Group, data = mat, subset = Group %in% c("A", "B"))

# 提取向量的写法
t.test(mat$number[mat$Group == "A"], mat$number[mat$Group == "B"])

2. 其他组间比较 #

• A vs C：

t.test(number ~ Group, data = mat, subset = Group %in% c("A", "C"))

• B vs C：

t.test(number ~ Group, data = mat, subset = Group %in% c("B", "C"))

单独提取p-value #

如果只需要结果中的p值，可以通过$p.value提取：

# 示例：提取A和B组比较的p值
a_b_pval <- t.test(number ~ Group, data = mat, subset = Group %in% c("A", "B"))$p.value
print(a_b_pval)

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