R语言做单因素方差分析

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

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

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

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

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

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

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