Agricolae 패키지에서 HSD.test()를 사용할 때 그룹을 주문하는 방법은 무엇입니까?

Question

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 

내가 나오지 않았어 오히려 낮다 아무것도 발견하지 못했고, 최소한 똑바로 명확한 것은 아무것도 없습니다. 아무도 해결책을 내게 도울 수 있습니까? 나는 대답에 매우 만족할 것이다. 사전에

감사합니다, Jelle

Answer 1

먼저 빌드 예를 ANOVA 모델 aov 기능을 Tukey HSD와

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

사후 시험에 의해, 입력 데이터, 계수 그룹화되는 원래의 데이터 세트로서 열 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")