问题标签 [tukey]

我正在处理一个不平衡的设计/样本，最初是从aov(). 我现在知道，对于我的 ANOVA 测试，我需要使用类型 III 平方和，这涉及使用拟合lm()而不是使用aov().

问题是使用lm(). 我所做的所有研究都表明simint在multcomp包中使用会起作用，但现在它已更新，该命令似乎不可用。它似乎也依赖于经过aov()计算。

基本上我为 R 找到的所有 Tukey HSD 测试都假设您aov()用于比较而不是lm(). 要获得不平衡设计所需的 III 型平方和，我必须使用：

我如何使用 Tukey HSD 测试与我的 mod 使用lm()？或者相反，使用 III 型计算我的 ANOVA 并且仍然能够运行 Tukey HSD 测试？

谢谢！

我有一个数据集，我想用它来比较物种和栖息地对家庭范围大小的影响——同时在物种和栖息地内使用 III 型错误和成对比较。
这是数据的一个子集：

目前我正在拆分三个物种的数据，然后为每个物种运行单独的 ANOVA，但我相信使用一个 ANOVA 同时询问物种和栖息地更有意义。这是我为一个物种运行的 ANOVA 示例：

aov() 似乎使用 I 类错误。. . 我认为不合适；另外，我相信 Tukey 的测试对于成对比较来说可能过于保守。有人可以帮助我采用一种方法，让我运行一个 ANOVA，考虑物种和栖息地对家庭范围的影响，具有 III 型错误，还允许对物种和栖息地进行不太保守的成对比较？

在阅读StackOverflow 中相关问题的答案时，我看到了 matplotlib中用于计算 wiskers 位置和检测异常值的代码：

现在，这else部分非常有意义 - 根据Tukey boxplots的规范，我们在上四分位数的 1.5 IQR 内找到了最高基准。确实，这是-下面max(wish_hi)是最大的数据条目。 Q3+1.5*IQR

or然而，我不明白的部分。if len(wisk_hi) == 0翻译为...

这个条件如何适用？Q3 是通过在中位数上拆分数据，然后取上半部分的中位数，然后在其上加上 1.5*IQR 找到的- 怎么可能没有低于这个值的数据？

如果这是关于一个空数据集，那么第二部分or也没有意义（因为 Q3 或 IQR 没有数据就没有意义）。

可能遗漏了一些明显的东西 - 帮助？

我有多个类似于以下的 csv 文件。

它实际上看起来像这样：CSV 表

在这种情况下，每个 .csv 文件都以化学参数命名，例如“logtPSA.csv”。如您所知，标题是“库”和“值”。我有 4 个不同的库：FDA、SMRNA、VD_SM 和 VD_MV。

我一直在对所有这些进行 ANOVA 测试和 TukeyHSD 测试，但是 R 并没有以将结果分成单元格的格式给出它，所以我可以复制并粘贴到 Excel 中，因此如果我要手动进行它会变得相当乏味。

我想知道是否有某种方法可以编写一个函数，一次遍历我的所有文件，执行测试，并将其全部输出到两个整洁的表中（1 个用于 ANOVA，1 个用于 TukeyHSD）。

我无法导出我的 TukeyHSD 结果，以便当我在 Excel 等中打开结果时它们在单元格中分开。我尝试使用 write.csv() 但它说：

如何以一种我可以将它们复制并粘贴到 Excel 工作表中的方式捕获我的 TukeyUSD 结果？

我正在分析参与者对某些问题的二进制答案的数据集。我正在使用该glm函数来测试如何Var * Base_con影响Dec. 拟合后，我试图比较“Var”因子如何影响每个“Base_con”因子水平内的结果。在这个小插曲之后，我尝试了以下（失败的）方法，我相信它可以被复制（如果它不起作用，请告诉我）：

最后一个命令的输出是：

实际上，错误消息报告的两个元素的列数和长度是不同的：

到目前为止，这就是我所能取得的进展。有任何想法吗？谢谢大家。

过去已经问过我关于在 Tukey HSD 测试后显示同质子集的问题（请参阅https://stats.stackexchange.com/questions/31547/how-to-obtain-the-results-of-a-tukey -hsd-post-hoc-test-in-a-table-showing-groupe，在讨论的最后）但从未得到回答。我现在遇到了同样的问题，我迫切需要一个解决方案，因为这个函数似乎是 R 显示此类组的唯一可用函数。如果不是，请告诉我。简而言之，这就是我得到的：

Groups, Treatments and means a 140095-001 36.79 b 150004-001 32.1 b 136936-021 31.97 bc 137219-004 31.39 bc 136673-017 31.27 bc 136963-009 30.79 bcd 147328-016 30.63 bcd 0147592-01 30.55 bcde 140094-001 30.02 cde 136730-007 29.7 de 136963-066 29.49 ef 136936-004 28.4 efg 147414-004 28.2 efg 137109-036 28.2 efg 136765-001 28.06 efg 140089-001 27.82 efg 137186-020 27.8 fg 136936-006 27.48 fgh 147350-014 27.43 gh 136992-001 27.36 gh 136730-015 27.18 ghi 0147785-01 27.08 ghi 0147691-01 26.98 ghi 136891-010 26.7 ghij 0147792-01 26.49 ghijk 136947-014 26.3 ghijkl 140097-001 25.8

这就是我想要的：

我已经检查了函数背后的脚本，但由于我的知识相当低，我找不到任何东西，至少没有什么直截了当的。有人可以帮我解决吗？我会很高兴得到一个答案。

在此先感谢，耶勒

在 R 中，我有四个数值向量。我想测试这四个向量的平均值是否相同，如果不是，哪个平均值更大，哪个平均值更小（基本上我想通过它们对四个向量进行排名）。

这个帖子很适合我：

http://brownmath.com/stat/anova1

我正在遵循本指南：

https://www.r-bloggers.com/analysis-of-variance-anova-for-multiple-comparisons/

但是，这篇文章停止测试四个平均值是否相等，但如果答案不是，我怎么知道 R 中哪个平均值更大，哪个平均值更小？

非常感谢你，

帖子中的代码是：

我有以下数据（dat）

我希望对上述数据集执行 TukeysHSD 成对测试。从测试结果来看，我想将显着性比较合并到图中（显示显着组之间的“*”或“**”符号）。

这是尝试的代码：

我不知道是否（我可以）以及如何将测试的结果（“重要性”）直接合并到图中，而不将每个结果保存为外部数组，然后将其调用到 ggplot2 中。

我有以下数据（dat）

我有以下数据（dat）

我希望对上述数据集执行 TukeysHSD 成对测试。

但是，每次我尝试运行它时，它都会继续运行并且不会产生输出。有关如何纠正此问题的任何建议？

我还想使用函数 TukeyHSD() 执行相同的测试。但是，当我尝试使用宽/长格式时，我遇到了一个错误，上面写着