R로 사후분석하기 (R in post-hoc analysis)

ANOVA 분석을 할 때, 유의한 factor의 level 이 3개 이상 되는 경우, 각 level 간의 차이는 ANOVA가 설명해줄 수 없다.

이때 활용하게 되는 것이 사후분석 (post-hoc analysis)이다.


R에서는 사후분석을 쉽게도와주는 PostHocTest 라는 함수가 있다.

먼저 필요한 라이브러리를 설치하자

> install.packages("DescTools")

> library(DescTools)

그 후 PostHocTest 함수를 쓰면 된다. 단 대소문자에 주의하자

> PostHocTest( aov(formula = value ~ factor1 * factor2), method="lsd", conf.level="0.95")



자세한 파라미터 값 설정은 아래 페이지를 참고하면 된다.
http://finzi.psych.upenn.edu/R/library/DescTools/html/PostHocTest.html



--추가 '17.12.05
SNK를 쓰는 경우
agricolae packages에 있는 SNK.test 함수를 바로 사용하면 된다.

> install.packages("agricolae")
> library(agricolae)
> SNK.test(aov(Time ~ A*B*C, data=exp), "A", group = TRUE, console=TRUE)