如何在R中编辑Tukey测试图

Question

我对ANOVA进行了事后Tukey测试，然后绘制了结果图。我似乎无法更改我的x轴标题或y轴标题。我收到此错误：

Error in plot.default(c(xi[, "lwr"], xi[, "upr"]), rep.int(yvals, 2L),  : 
  formal argument "xlab" matched by multiple actual arguments

这是我的相关代码：

tuk <- TukeyHSD(final)
plot(tuk,xlab="Differences in mean departure times", ylab="Comparisons")

我还需要更改y轴刻度线标签，但我不知道如何。

谢谢。

Answer 1

因此，由于他们是如何为plot()类对象编写TukeyHSD方法的，因此默认情况下您无法更改轴标签，此详细信息已埋在?TuketHSD手册页中。

但是您可以轻松地将确实允许您复制的副本合并在一起。首先用getAnywhere(plot.TukeyHSD)查找现有方法的代码。然后像这样修改它：

tuk_plot <- function (x, xlab, ylab, ylabels = NULL, ...) {
  for (i in seq_along(x)) {
    xi <- x[[i]][, -4L, drop = FALSE]
    yvals <- nrow(xi):1L
    dev.hold()
    on.exit(dev.flush())
    plot(c(xi[, "lwr"], xi[, "upr"]), rep.int(yvals, 2L), 
         type = "n", axes = FALSE, xlab = "", ylab = "", main = NULL, 
         ...)
    axis(1, ...)
    # change for custom axis labels
    if (is.null(ylabels)) ylabels <- dimnames(xi)[[1L]]

    axis(2, at = nrow(xi):1, labels = ylabels, 
         srt = 0, ...)
    abline(h = yvals, lty = 1, lwd = 0.5, col = "lightgray")
    abline(v = 0, lty = 2, lwd = 0.5, ...)
    segments(xi[, "lwr"], yvals, xi[, "upr"], yvals, ...)
    segments(as.vector(xi), rep.int(yvals - 0.1, 3L), as.vector(xi), 
             rep.int(yvals + 0.1, 3L), ...)
    title(main = paste0(format(100 * attr(x, "conf.level"), 
                               digits = 2L), "% family-wise confidence level\n"), 
          # change for custom axis titles
          xlab = xlab, ylab = ylab)

    box()
    dev.flush()
    on.exit()
  }
}

现在您可以调整x和y轴以及自定义y标签：

tuk_plot(tuk, "Hello X Axis", "Hello Y Axis", c("One", "Two", "Three"))

如果不提供y标签，则会显示模型中的默认y标签。

可复制的示例：

fm1 <- aov(breaks ~ wool + tension, data = warpbreaks))
tuk <- TukeyHSD(fm1, "tension")