代码如下：
library(multcomp)
attach(cholesterol)
cholesterol
table(trt)
aggregate(response, by = list(trt), FUN=mean)
aggregate(response, by = list(trt), FUN=sd)
#方差分析，trt只能是因子变量
fit <- aov(response ~ trt)
summary(fit)
library(gplots)
plotmeans(response ~ trt, xlab="Treatment", ylab="Response",
main="mean plot")
虽然知道了组间有差异，但是并没有告诉哪个组和哪个组有差异，使用TukeyHSD(fit)
TukeyHSD(fit)
#画图查看差异
par(las=2)
par(mar=c(5,8,4,2))
plot(TukeyHSD(fit))
有无交互作用
原创文章（本站视频密码：66668888），作者：xujunzju，如若转载，请注明出处：https://zyicu.cn/?p=16596
微信扫一扫 
支付宝扫一扫 