Блог

Итак, у меня есть 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")
  

Он просто продолжает выдавать мне ошибки. Я предполагаю, что данные должны быть в какой-то другой форме?

pairwise.t.test Требуется коэффициент группировки или «g». В этом случае отдельные vector s могут быть преобразованы в именованный list с mget помощью, затем преобразовать его в двухколоночный data.frame с stack помощью, а затем использовать pairwise.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