R语言学习单因素协方差分析.doc

R语言学习-单因素协方差分析
单因素协方差分析
单因素协方差分析扩展了单因素方差分析，包含一个或多个定量的协变量。示例数据来源于multcomp包中的litter数据集。怀孕小&组，每个小组接受不同剂量的药物处理。产下的幼崽的体重均值为R语言学习-单因素协方差分析
单因素协方差分析
单因素协方差分析扩展了单因素方差分析，包含一个或多个定量的协变量。示例数据来源于multcomp包中的litter数据集。怀孕小&组，每个小组接受不同剂量的药物处理。产下的幼崽的体重均值为因变量，怀孕时间为协变量
library(multcomp)
table(litter$dose)
0550500
20191817
#可以看出每种剂量产下的幼崽数不同
aggregate(litter$weight,by=list(litter$dose),FUN=mean)





#显示各组体重均值
fit<-aov(weight~gesttime+dose,data=litter)#协变量放在主要因子前面
summary(fit)
DfSumSqMeanSqFvaluePr(>F)
**
*

:0‘***'‘**'‘*'‘.'‘'1
从单因素协方差的F检验可以看出怀孕时间与幼崽体重有关
结果可视化
HH包中的ancova()函数可以绘制因变量，协变量和因子之间的关系图
latticegridlatticeExtraRColorBrewergridExtra
library(HH)载入需要的程辑包载入需要的程辑包载入需要的程辑包载入需要的程辑包载入需要的程辑包
ancova(weight~gesttime+dose,data=litter)AnalysisofVarianceTable
Response:weight
DfSumSqMean