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

Question

Я пытался использовать следующий код из предыдущего поста: Проблемы с графиком Тьюки в 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);