如何在使用Agricolae软件包中的HSD.test（）时订购组？

我的问题是关于在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的功能。否则，请告诉我。总之，这是我得到：如何在使用Agricolae软件包中的HSD.test（）时订购组？

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

我已经检查了脚本功能落后，但因为我的知识是相当低我不能找不到任何东西，至少没有什么直接和明确的。任何人都可以帮我解决问题吗？我会很高兴回答。

由于提前，耶勒

来源

2016-11-16 Jelle

首先构建例如ANOVA模型通过aov功能

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")

来源

2017-07-31 14:13:57 tlask

如何在使用Agricolae软件包中的HSD.test（）时订购组？

回答

相关问题