- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我对方差分析进行了事后 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 轴刻度标记标签,但我不知道如何更改。
谢谢。
最佳答案
因此,由于他们如何为 TukeyHSD
类对象编写 plot()
方法,默认情况下您无法更改轴标签,此细节隐藏在 ?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 标签,则会显示模型中的默认标签。
可重现的例子:
fm1 <- aov(breaks ~ wool + tension, data = warpbreaks))
tuk <- TukeyHSD(fm1, "tension")
关于r - 如何在 R 中编辑 Tukey 测试图,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/60794019/
我正在尝试使用 yaxis 上的计数 (MedMean) 和 xaxis 上的各种独立样本 (Site_Name) 的 ggplot 创建箱线图。 ggplot(medianlist,aes(x=re
在实现对快速排序分区的改进时,我尝试使用 Tukey 的九分法来找到主元(几乎借用了 sedgewick 在 QuickX.java 中的实现) 每次打乱整数数组时,我的下面代码都会给出不同的结果。
我知道 FFT 实现是如何工作的 ( Cooley-Tuckey algorithm ),并且我知道有一个 CUFFT CUDA 库可以快速计算 1D 或 2D FFT,但我想知道在此过程中如何利用
我一直无法在任何标准包中找到这个函数,所以我写了下面的一个。然而,在将它扔向 Cheeseshop 之前,有人知道已经发布的版本吗?或者,请提出任何改进建议。谢谢。 def fivenum(v):
我对方差分析进行了事后 Tukey 检验,然后绘制了结果图。我似乎无法更改 x 轴标题或 y 轴标题。我收到此错误: Error in plot.default(c(xi[, "lwr"], xi[,
我正在尝试对死亡率数据进行 Tukey 检验,我想测试死亡率是否受铜含量(单向方差分析)以及铜和温度的组合(双向方差分析)的影响.这些是我的公式: lm2<-lm(Mortality~Cu) anov
我收集了 216 个人的数据。我测量了每个人体内相同的 7 种物质的浓度,用 Sub1:Sub7 表示。这些物质的浓度在不同地点的个体中可能不同。我对这些个体可以根据这些物质的浓度进行分组的细化程度感
问题:我想了解如何更改 R 中 Tukey 测试计算平均值并分配相应字母的样本顺序。下面是非常简单的示例。 我研究了iris数据,发现不同物种的萼片长度存在差异。这是箱线图: 我进行了方差分析测试,发
我和我的一些学生已经在许多地方寻找解决这个问题的方法,但没有成功,而且实际上已经花费了几个月的时间。我一直被称为我不想使用的 lme 命令。提供的输出不是我的同事或我自己使用了超过 15 年的输出。此
我有多个链接,例如 - http://example.com/link1/title-with-detail - http://example.com/link2/title-with-detail
我刚接触到 tucker,正在尝试在 pentaho 服务器(spring)中重写 url 查询参数。 我想做的是:重写 http://hostname:8080/pentaho/api/repos/
在阅读a related question in StackOverflow的答案时,我看到matplotlib中使用的代码to compute the wiskers positions and d
我正在分析用 Python 编写的 Cooley-Tukey 算法实现的复杂性(代码取自 here): def fft(x): N = len(x) print N, N//2 if N 10 的列
我正在尝试为 FFT 编写 Cooley Tukey 算法。现在,该算法运行良好,但仅适用于 2 个数字 - 没有别的。例如,我使用在线 FFT 计算,输入相同的数据并得到相同的结果。这是算法的代码:
试图找出如何使用 statsmodel 计算 Tukey 的 HSD。我可以使它工作并且结果看起来很棒但是有一个我看不到的均值差异图。一定是我在做的傻事。 它是对象 TukeyHSDResults 中
我已尝试查看 Python 的多个统计模块,但似乎找不到任何支持 one-way ANOVA 事后测试的模块。 最佳答案 单方面的方差分析可以像这样使用 from scipy import stats
我有一个包含多个变量的大型数据集。我需要进行双向方差分析,然后使用 Tukey HSD 进行事后成对多重比较。 我的前 25 个条目的数据头部是这样的: > head(my_data2, 25 )
我正在对我的数据运行事后 Tukey HSD,该数据有 10 个因子级别。该表很大,我希望在成对的表中向读者展示 p 值,将 45 行表留给附录。 这是一个示例数据集: set.seed(42) x
我正在尝试对从 posthoc Tukey 获得的临界值进行精细可视化。有一些good guidelines out there用于可视化成对比较,但我需要更精致的东西。我的想法是,我会有一个图,其中
我正在对我的数据运行事后 Tukey HSD,该数据有 10 个因子级别。该表很大,我希望在成对的表中向读者展示 p 值,将 45 行表留给附录。 这是一个示例数据集: set.seed(42) x
我是一名优秀的程序员,十分优秀!