Tukey测试仅绘制显着的平均差异

Question

我尝试使用上一篇文章中的以下代码：Tukey Graphing Problems in R

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

但是在尝试绘图时会出现以下错误：

  

xy.coords（x，y，xlabel，ylabel，log）中的错误：         ＆＃39; X＆＃39;并且＆＃39; y＆＃39;长度不同

我可能做错了什么？

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);