所用的数据是multcomp包中的数据,为不同治疗方法下针对胆固醇的疗效。
在画出上面这个图后,书上写:
从结果来看, 使用降低胆固醇的药物时, 一天四次5mg剂量比一天一次20mg剂量效果更佳,也优于候选药物drugD, 但药物DrugE比其他所有药物的疗法都更优。
请问这个是如何从上面的图看出来的?
我能看出来一天四次 比 一天一次更好,因为一天四次的那个均值高,但怎么就优于DrugD了?DrugD明明均值高于一天四次。
再者,一天四次和DrugE之间差异也不显著啊,不能比吧。
用的代码如下:
- library(multcomp)
- attach(cholesterol)
- table(trt)
- trt
- aggregate(response, by=list(trt), FUN=mean)
- aggregate(response, by=list(trt), FUN=sd)
- fit<-aov(response ~ trt)
- summary(fit)
- tuk<-glht(fit, linfct=mcp(trt="Tukey"))
- plot(cld(tuk, level=.05), col="lightgrey")


雷达卡




京公网安备 11010802022788号







