在看R in action,方差分析里面的多重比较
> library(multcomp)
> contrast <- rbind("no drug vs. drug" = c(3, -1, -1, -1))
> summary(glht(fit, linfct=mcp(dose=contrast)))
看不太懂,这个linfct应该就是线性方程的意思,但是mcp是个啥呢,帮助里也没细说后面的(dose=contrast)为什么要这样写呢?感觉整个glht函数都怪模糊的有没有大神指点一下。