Я пытаюсь понять, как работать с ANOVA и апостериорными тестами в R. Пока что я использовали aov () и TukeyHSD () для анализа моих данных. Пример:
uni2.anova <- aov(Sum_Uni ~ Micro, data= uni2)
uni2.anova
Call:
aov(formula = Sum_Uni ~ Micro, data = uni2)
Terms:
Micro Residuals
Sum of Squares 0.04917262 0.00602925
Deg. of Freedom 15 48
Residual standard error: 0.01120756
Estimated effects may be unbalanced
Моя проблема в том, что теперь у меня есть огромный список парных сравнений, но я ничего не могу с ним поделать:
TukeyHSD(uni2.anova)
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = Sum_Uni ~ Micro, data = uni2)
$Micro
diff lwr upr p adj
Act_Glu2-Act_Ala2 -0.0180017863 -0.046632157 0.0106285840 0.6448524
Ana_Ala2-Act_Ala2 -0.0250134285 -0.053643799 0.0036169417 0.1493629
NegI_Ala2-Act_Ala2 0.0702274527 0.041597082 0.0988578230 0.0000000
В этом наборе данных 40 строк ... В идеале, я бы хотел получить набор данных это выглядит примерно так:
Надеюсь, вы уловили суть. Пока я не нашел ничего сопоставимого в Интернете ... Я также попытался выбрать только значимые пары в файле, полученном в результате TukeyHSD, но файл не «подтверждает», что он состоит из строк и столбцов, что делает выбор невозможным .. .
Может быть, в моем подходе что-то не так?