0

过去已经问过我关于在 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

这就是我想要的：

    Groups, Treatments and means
a        2.1     51.17547 
ab       4.1     50.7529 
abc      3.1     47.36229 
 bcd     1.1     45.81229 
  cd     5.1     44.55313 
   de    4.0     41.81757 
    ef   2.0     38.79482 
    ef   1.0     36.91257 
     f   3.0     36.34383 
     f   5.0     35.69507

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

在此先感谢，耶勒

4

1 回答 1

0

首先按aov功能构建例如 anova 模型

my_model <- stats::aov(data$value ~ data$treatment) 
summary(my_model)

使用 Tukey HSD 进行事后测试，使用 anova 模型作为输入数据，treatment原始数据集中的列是分组因子。

test_result <- agricolae::HSD.test(my_model, "my_model$treatment", group = TRUE, console = TRUE)

并在控制台中检查结果。除其他外，还有一张包含您询问的概述的表格。（在我的例子中并不重要）

  Groups, Treatments and means
 a   2008    6.755 
 a   2013    6.147 
 a   1975    4.653 
 a   1997    3.622

可以可视化治疗结果。

agricolae::bar.group(test_result$groups, ylim = c(0,10), density = 10, border = "blue")

Tukey HSD 的可视化

于 2017-07-31T14:13:57.057 回答