问题标签 [tukey]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
10327 浏览

r - 将 Tukey 的重要字母添加到箱线图中

我正在尝试在 yaxis 上创建一个带有计数(MedMean)的 ggplot 和在 xaxis 上的各种独立样本(Site_Name)的箱线图。

我想将 Tukey 的重要字母添加到框中。

谢谢

0 投票
1 回答
17131 浏览

r - TukeyHSD.aov(my_data) 中的错误:拟合模型中没有因素

我已成功上传数据,在我的数据集上运行 shapiro、bartletts 和单向方差分析,但是无论我尝试什么,我似乎都无法运行 TukeyHSD 而不收到一些错误消息,例如上面的错误消息,这是我输入的内容我错过了什么?

任何帮助都会很棒!

0 投票
0 回答
254 浏览

python - 基于字母表示的所有成对比较的 Python 算法

在经过一些成对测试后,Python 中是否存在任何东西来获取字母,例如在 R 包 multcompView 中有来自 Piepho 2004 的 multcompLetter:https ://www.jstor.org/stable/1391186 ?

基本上,我有很多来自测试的成对布尔值,并且想基于 Python 构建我的字母。

编辑。

所以我想生成第二个图中的字母:https ://www.r-graph-gallery.com/84-tukey-test/

这是我测试后的数据示例:

0 投票
1 回答
717 浏览

r - 为单向混合效应方差分析的 Tukey 检验获取错误的 p 值

我正在尝试复制用于单向混合效应方差分析的 Tukey 测试的教科书示例(来自 Statistics,William L. Hays p 581-583),但我使用 lme 和 glht 获得的 p 值没有意义

该研究重复测量了四种治疗方法和 10 名受试者

数据

该模型

这个 F 值与 Hay 的 ( ) 足够接近F = 43.41,我很确定我的模型很好。

图基测试

这与书本不一致。Hay's 只报告了一些比较,并给出了 HSD 和均值,而不是 p,但对于 a3-a1 的对比,发现与 an 没有显着差异,HSD = 4.02并且mean = 3.9根据我的计算有p.value = 1-ptukey(3.9*sqrt(10/8.27),4,9)=0.05719563.

R 输出也没有意义,因为应该控制多重比较的 Tukey 检验 p 值远小于配对 t 检验的 p 值(p=0.0142,使用t.test(c(11, 14, 15, 17, 15, 14, 9, 17, 10, 12),c(5, 10, 6, 13, 6, 13, 7, 10, 14, 11), paired=TRUE).

知道我做错了什么以及如何在 R 中正确执行 Tukey 测试吗?

0 投票
1 回答
1161 浏览

r - 如何在数据帧列表上执行 Tukey HSD.test()?

我想对数据框列表执行 Tukey Post Hoc 测试。作为结果,我希望有字母表明哪些组彼此之间存在显着差异。agricolae 包的 HSD.test() 执行此操作,但我无法弄清楚如何一次将其应用于多个变量。这将为我节省大量时间,因为我的数据集包含很多变量。

这是我的数据的一部分:

我知道如何对一个变量进行 anova 和 HSD.test() 。

此外,我知道如何使用 sapply() 函数将方差分析应用于我的数据集中的所有变量:

但我不知道如何对这些结果执行 posthoc HSD.test() 。我试过这个:

我得到 NULL 作为结果,所以出了点问题,但我不知道是什么。我还尝试了 Tukey 事后测试的另一个功能,即:TukeyHSD()。

这适用于我的数据集,但这个函数没有给出我真正想要的字母。有人知道我如何对 HSD.test() 做同样的事情,以便获得这些信件吗?谢谢!

0 投票
0 回答
1074 浏览

python-3.x - 如何在 python 中计算 Multicompare Tukey HSD?

我正在尝试使用包含 5 个值列表(它是列表列表)的列表来计算多重比较 Tukey 测试。我正在阅读一些numpy.recarray似乎适合该主题的文档,但我真的不知道它是如何工作的。假设我的列表如下:

输出应该是这样的,需要比较 3 个元素:

此外,有没有办法在手段上绘制多重比较?提前致谢

0 投票
1 回答
351 浏览

glm - GLM 上带有 lsmeans Tukey 测试的“对比度错误”消息

我定义了一个广义线性模型如下:

我有以下重要的互动

我正在尝试进行事后测试(Tukey)来比较使用lsmeans包的 ParticleTypeSize 的交互。但是,一旦我继续,我就会收到以下消息:

我通过应用检查了 ParticleTypeSize 是否是一个有效因素:

我很困惑,不确定如何纠正这个错误信息。任何帮助将非常感激!

0 投票
2 回答
1006 浏览

python - 根据显着性标记来自 Tukey 测试结果的组

pairwise_tukeyhsd我有一个来自 python的 Tukey 测试表statsmodels.stats.multicomp

我有这张桌子pandas df。我想(用字母)标记表示统计关系的组(101-106)。对于这个特定示例,所需的结果是:(我不介意结果是 df、列表、字典)

如您所见,所有具有相似字母的组具有相同的均值(拒绝列 = False),具有不同字母的组(拒绝列 = True)具有不同的均值。例如,组 101 的平均值等于所有其他组的平均值,因为组 101 具有字母 ab,而所有其他组具有abab。另一方面,组 106 只有字母b,这表明它与除了只有字母a的组 102之外的所有组相似。

我找不到一个自动 python 解决方案。我看到R 有一个名为 的包,multcompLetters在 python 中有类似的东西吗?

0 投票
1 回答
269 浏览

r - Tukey 检验仅绘制显着均值差异

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

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

xy.coords(x, y, xlabel, ylabel, log) 中的错误:“x”和“y”长度不同

我可能做错了什么?

0 投票
1 回答
1405 浏览

r - 如何在 agricolae 包中进行事后测试,Tukey?

我做了这个

但我很感兴趣如何在 R 中的 TW ANOVA 之后执行事后测试。在agricolae包中?

在平均分离之后,我们如何用组号/字母绘制每个阶段的图表?

阴谋