그룹화 된 쌍을 보여주는 표에서 Tukey HSD 사후 테스트 결과를 얻는 방법은 무엇입니까?

R과 함께 양방향 Anova 후 TukeyHSD 사후 테스트를 수행하여 중요한 차이로 그룹화 된 정렬 쌍을 포함하는 테이블을 얻고 싶습니다. (언어에 대해 미안하지만, 나는 여전히 통계에 익숙하지 않다.)

나는 이와 같은 것을 갖고 싶다 :

따라서 별이나 문자로 그룹화됩니다.

어떤 생각? 패키지 에서 함수 HSD.test()를 테스트 agricolae했지만 양방향 테이블을 처리하지 않는 것 같습니다.

이 agricolae::HSD.test함수는 정확히 그렇게하지만 상호 작용 항에 관심이 있음을 알려야 합니다. Stata 데이터 셋을 사용한 예는 다음과 같습니다.

library(foreign)
yield <- read.dta("http://www.stata-press.com/data/r12/yield.dta")
tx <- with(yield, interaction(fertilizer, irrigation))
amod <- aov(yield ~ tx, data=yield)
library(agricolae)
HSD.test(amod, "tx", group=TRUE)

결과는 다음과 같습니다.

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 

그것들은 다음 명령으로 얻을 수있는 것과 일치합니다 :

. webuse yield
. regress yield fertilizer##irrigation
. pwcompare fertilizer#irrigation, group mcompare(tukey)

-------------------------------------------------------
                      |                           Tukey
                      |     Margin   Std. Err.   Groups
----------------------+--------------------------------
fertilizer#irrigation |
                 1 0  |   36.91257   1.116571    AB    
                 1 1  |   45.81229   1.116571      CDE 
                 2 0  |   38.79482   1.116571    AB    
                 2 1  |   51.17547   1.116571         F
                 3 0  |   36.34383   1.116571    A     
                 3 1  |   47.36229   1.116571       DEF
                 4 0  |   41.81757   1.116571     BC   
                 4 1  |    50.7529   1.116571        EF
                 5 0  |   35.69507   1.116571    A     
                 5 1  |   44.55313   1.116571      CD  
-------------------------------------------------------
Note: Margins sharing a letter in the group label are
      not significantly different at the 5% level.

multcomp의 패키지는 ( '소형 편지 표시'를 참조 상징적 인 시각화를 제공합니다 : 비교 및 평가를 컴팩트 문자 표시에 대한 알고리즘을 의미 페어의 비교,이 표 형식으로 표시하지 않는 있지만 자세한 내용). 그러나 상자 그림을 사용하여 결과를 편리하게 표시 할 수있는 플로팅 방법이 있습니다. 프리젠 테이션 순서도 변경할 수 있으며 (옵션 decreasing=) 여러 비교를위한 더 많은 옵션이 있습니다. 이러한 기능을 확장하는 multcompView 패키지 도 있습니다.

다음은 동일한 예를 분석 한 것입니다 glht.

library(multcomp)
tuk <- glht(amod, linfct = mcp(tx = "Tukey"))
summary(tuk)          # standard display
tuk.cld <- cld(tuk)   # letter-based display
opar <- par(mai=c(1,1,1.5,1))
plot(tuk.cld)
par(opar)

동일한 문자를 공유하는 치료는 선택한 수준에서 크게 다르지 않습니다 (기본값, 5 %).

덧붙여서, 현재 R-Forge에서 호스팅되는 새로운 프로젝트가 있는데, 유망한 것으로 보입니다 : factorplot . 여기에는 라인 및 문자 기반 디스플레이와 모든 쌍별 비교에 대한 매트릭스 개요 (레벨 플롯을 통한)가 포함됩니다. 작업 논문은 여기에서 찾을 수 있습니다 : factorplot : GLM에서 단순 대비 표현 개선

이라는 함수입니다 TukeyHSD즉, 도움말 파일에 따라, 범위의 지정된 가족 현명한 확률 요인의 수준 수단의 차이점에 대한 신뢰 구간의 집합을 계산합니다. 간격은 Tukey의 "정직한 유의미한 차이"방법 인 스튜던트 화 된 범위 통계를 기반으로합니다. 이것이 당신이 원하는 것을합니까?

http://stat.ethz.ch/R-manual/R-patched/library/stats/html/TukeyHSD.html