R 三向量配对t检验

所以我有3个向量，需要对它们进行两两t检验

这是我到目前为止的代码

> x1=c(6.5,15.5,15.2,10.8,3.7,11.2)
> x2=c(11.2,10.5,6.8,5.1,5.7)
> x3=c(5.9,6.8,4.1,10.4)
> pairwise.t.test(x1,x2,x3,p.adj="bonferroni")

---

它总是给我错误。我猜数据必须是其他形式的？

成对t.test需要一个分组因子或“g”。在这种情况下，可以使用

mget

将单个

向量

转换为命名的

列表

，然后使用

堆栈

将其重塑为两列data.frame，然后使用

成对.t.test

with(stack(mget(paste0("x", 1:3))), 
     pairwise.t.test(values, ind, p.adj = 'bonferroni'))
#   Pairwise comparisons using t tests with pooled SD 

#data:  values and ind 

#   x1   x2  
#x2 0.79 -   
#x3 0.44 1.00