Scheffe 사후검정

ANOVA 사후분석 중 하나인 scheffe 검증에 대해 공부하고 있습니다.

 

집단간 자료수가 동일하지 않은 경우에도 사용할 수 있다고 하니 꼭 알아두어야 할 검정방법인 듯 합니다.

 

agricolae 패키지의 샘플로 공부를 시작합니다.

sweetpotato 데이터의 건수가 매우 작아 aov 대신 비모수검정을 해야하는 것 아닌가 궁금합니다만 일단 지나갈께요. 

library(agricolae)

data(sweetpotato)
# 비모수 해야하는 것 아닌가?
model<-aov(yield~virus, data=sweetpotato)
summary(model)
comparison <- scheffe.test(model,"virus", group=TRUE,console=TRUE,
                           main="Yield of sweetpotato\nDealt with different virus")

# Old version scheffe.test()
df<-df.residual(model)
MSerror<-deviance(model)/df
Fc<-anova(model)["virus",4]
out <- with(sweetpotato,scheffe.test(yield, virus, df, MSerror, Fc))
print(out)

comparison 즉 테스트 결과를 확인해 보았는데, 눈에 쏙 들어오지가 않네요.

 

 

 

Wiki Scheffé's method

 

테스트 결과 해석을 위해 위키를 보았는데 첩첩산중이네요. 도무지 무슨 소린지...

 

출처:https://en.wikipedia.org/wiki/Scheff%C3%A9%27s_method

 

다음 블로그 내용이 도움이 될 듯합니다. 시간을 갖고 Wiki 내용과 비교하면서 이해를 좀 해보아야겠습니다.

분산분석에서는 둘 이상의 수준(factor levels) 평균들을 포함하는 비교를 하기 위해 각 수준들의 평균 μi 를 하나의 직선식으로 표현한 것을 대비(Contrasts)라고 하며 L로 표기하면, L의 추정량(estimated contrast L^)과 추정량의 분산(estimated variance of L^), 샤폐계수 T (scheffe's T statistics)는 아래와 같습니다.

출처: https://rfriend.tistory.com/134 [R, Python 분석과 프로그래밍의 친구 (by R Friend):티스토리]

 

 

ANOVA의 사후분석(Scheffe) 보는 방법

SPSS 에서 Scheffe 검증을 수행하고 결과를 해석하는 유튜브 영상이 이해하기가 쉬웠는데, R 결과와 연결이 아직 안되네요.

 

 

출처 : ANOVA의 사후분석(Scheffe) 보는 방법(https://youtu.be/W7Z0JEu62K8)