R语言做单因素方差分析

R语言做单因素方差分析

方差分析是一种分析试验数据的有效的统计方法，它主要是根据观测数据的总变异按照变异原因的不同分解为因子效应和试验误差并对其作出数量分析，比较各种原因在总变异中所占的重要程度，作为进一步统计推断的依据。而R语言正好可以方便快速地解决这个问题，下面小编就简单叙述一下在R语言中如何进行单因素方差分析。

首先，我们打开R。从书上找一个例题进行演示。（例题来自《统计建模与R软件》薛毅版）

将例题中数据输入在R操作框内。注意，要以数据框形式输入。

接下来，我们调用aov()函数进行方差分析。

最后，我们用summary()函数提取方差分析的信息，回车即可得到方差分析的结果。

分析上述计算结果，Df表示自由度，Sum Sq 表示平方和，Mean Sq 表示均方，F value 是F值，Pr(>F)是p值，A即为因子A，Residuals 是残差。

但是我们注意到，这个结果并不完整，我们做方差分析时还需要计算总行和，所以我们这里调用一个自编函数，该函数可以将表中第一行与第二行求和，得到总行和的值。

运行后，得到的结果中出现了总行和值Total ，这就得到了一个完整的方差分析表了。就可以根据所得结果进行分析了。