AgricolaeパッケージからHSD.test（）を使用するときのグループの注文方法は？

Question

Tukey HSDテストがすでに過去に尋ねられた後に均質なサブセットを表示することについて私の質問がありますが（ディスカッションの最後には参照）、決して答えられませんでした。私は今同じ問題を抱えています。この機能がこのようなグループを表示する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 

は、私はすでに機能の後ろにスクリプトをチェックしますが、私の知識はかなり低いので、私はcouldn何も見つからない、少なくとも真っ直ぐなものは何もない。誰かが私に解決策を教えてくれますか？私は答えに非常に満足しています。事前に

おかげで、 Jelle

Answer 1

最初のビルド例えばaov機能

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

ポストホック検定TukeyのHSDとによるANOVAモデルは、元のデータセット内の列treatment因子をグループ化され、入力データとしてANOVAモデルを使用します。

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