R语言中的方差分析

来源 | R友舍

一、单因子方差分析（one-way ANOVA）
1)建模：
我们采用multcomp包中的cholesterol数据集作为例子，其中response为响应变量，trt为预测变量，这个处理中有五种水平。从下面的箱形图中可观察到处理的不同水平对于响应变量的影响。再用aov函数建立单因子方差模型，从结果的P值可看到各组均值有显著不同。

aggregate(response, by=list(trt), FUN=mean)
bwplot(response~trt)
model=aov(response~trt)
summary(model)

2)多重比较：
方差分析只告诉我们这五组之间是不同的，但没有告诉我们哪两组之间有明显差别，此时需要使用TukeyHSD函数进行均值的多重比较分析，从结果中观察到有三个两两比较是不显著的。

(result=TukeyHSD(model))
plot(result)

3)假设检验：
方差分析需要一定的假设，即数据集应该符合正态和同方差，我们分别用下面的函数来进行检验，从P值观察到这两个假设是符合的。对于不符合假设的情况，我们就要用到非参数方法，例如Kruskal-Wallis秩和检验

shapiro.test(response)
bartlett.test(response~trt)

二、双因子方差分析（Two-way Factorial ANOVA）
我们用ToothGrowth数据集来举例双因子方差分析。其中supp和dose是预测变量，len是响应变量。我们仍然使用aov进行建模，然后使用HH包的绘图函数来展现双因子交互效果图

fit <- aov(len ~ supp*dose)
library(HH)
interaction2wt(len~supp*dose)

要注意在下面的情况下因子的先后顺序是有讲究的：

• 第一种情况是多因子非平衡情况下，此时重要的因子应该放在前面，

• 第二是在在有协变量情况下，此时协变量放在前面，然后是主因子和交互因子

—————————————

往期精彩：

• 美国“实体清单”再度升级，我国超级计算机今日起或面临“断供”！

• 我造的假我自己打，Adobe推出“反PS”
  

• 中国诚信全球垫底？讲讲《Science》现在的论文有多不靠谱