p.140 교호작용 분산분석 할때 던넷을 사용하려고 하는데 reference 가 바뀌질 않는군요.
attach(war)
tension=factor(tension,levels=c("L","M","H"))
> levels(tension)
[1] "L" "M" "H"
> tension=relevel(tension,ref="L")
ref 를 L로 바꿨습니다.
> library(multcomp)
> out = lm(breaks~wool*tension, data=war)
> dunnett = glht(out, linfct=mcp(tension="Dunnett"))
summary(dunnett)
Estimate Std. Error t value Pr(>|t|)
L - H == 0 20.0000 5.1573 3.878 0.000626 ***
M - H == 0 -0.5556 5.1573 -0.108 0.991515
여전히 (L-H) 로 나오는데 어떻게 해야 하나요?
첫댓글 tension=relevel(tension,ref="L") 대신 다음을 시도해보세요.
war$tension=relevel(war$tension,ref="L")