Тест Тьюки показывает только средние различия

Question

Я пытался использовать следующий код из предыдущего поста: Tukey Graphing Проблемы в R

SigOnly <- Tukey 
SigOnly$species <- SigOnly$species[SigOnly$species[,'p adj'] < .05, ]
plot(SigOnly)

Но получите следующую ошибку при попытке построения графика:

Ошибка в xy.coords(x, y, xlabel, ylabel, log): длины 'x' и 'y' различаются

Что я могу делать не так?

Answer 1

Поскольку вы не предоставляете воспроизводимый код с образцами данных, вот пример, использующий warpbreaks Пример данных:

fm1 <- aov(breaks ~ wool + tension, data = warpbreaks)
res <- TukeyHSD(fm1, "tension", conf.level = 0.95)

# Filter rows with padj < 0.05
res$tension <- res$tension[res$tension[, 4] < 0.05, ];

# Plot results
plot(res);