操作R语言举办方差阐明
一、单因子方差分析（one-way ANOVA）
1)建模：
我们采用multcomp包中的cholesterol数据集作为例子，其中response为响应变量，trt为预测变量，这个处理中有五种水平。从下面的箱形图中可观察到处理的不同水平对于响应变量的影响。再用aov函数建立单因子方差模型，从结果的P值可看到各组均值有显著不同。
1 |
aggregate(response, by=list(trt), FUN=mean)
|
2 |
bwplot(response~trt)
|
3 |
model=aov(response~trt)
|
4 |
summary(model)
|
2)多重比较：
方差分析只告诉我们这五组之间是不同的，但没有告诉我们哪两组之间有明显差别，此时需要使用TukeyHSD函数进行均值的多重比较分析，从结果中观察到有三个两两比较是不显著的。
|