1.单因素协方差分析
这里我们使用multcomp包中的litter数据集(74只怀孕小鼠被分成四组,自变量是每个小组分别接受0克,5克,50克和500克的药物剂量处理,因变量是产下的小鼠幼崽的体重,协变量是母鼠的孕期)来进行演示。如果只是看哪个效果更加明显可以直接计算平均值,如果要看差别的显著水平则需要用方差分析。
library(multcomp)#加载包
attach(litter)#用来加载数据,用attach函数就不需要用变量赋值操作
table(litter$dose)#统计分组情况
aggregate(weight,by=list(dose),FUN=mean)#分组统计平均数
fit<-aov(weight~gesttime+dose,data=litter)#进行方差分析,协变量在前,自变量在后
summary(fit)#总结结果

结果显示,药物对体重数据有明显影响,P值小于0.05;同时,幼崽体重与母鼠怀孕时间也有关系。
网友评论