ryuyan 方差分析_R语言中实现方差分析及其可视化

方差分析(Analysis of Variance，简称ANOVA)，又称“变异数分析”，是R.A.Fisher发明的，用于两个及两个以上样本均数差别的显著性检验。 由于各种因素的影响，研究所得的数据呈现波动状。 造成波动的原因可分成两类，一是不可控的随机因素，另一是研究中施加的对结果形成影响的可控因素。

方差分析分类：

方差分析常用的为单因素方差分析和双因素方差分析，而多因素方差分析比较复杂用的不多，这里主要介绍前两种方差分析。

单因素方差分析

单因素方差分析比较的是分类因子定义的两个或多个组别中的因变量均值。与t检验区别：t检验适用于两列数据的均值比较。单因素方差分析适用于两列或更多列数据的均值比较。但对于两列数据的均值比较，单因素方差分析=等方差假设的双尾t检验。

以R中内置数据PlantGrowth集为例。不同地块( ‘ctrl’, ‘trt1’, and ‘trt2’)种植的植物的产量(weight)数据；解决的问题如下

1.不同地块产量是否存在差异

2.trt1与trt2是哪块差异更大

##正态性检验

my_data 

shapiro.test(my_data$weight)

#结果表明，p=0.8915 >0.05，接受原假设，说明数据在多个水平下都是正态分布的。

##方差齐性检验

bartlett.test(weight~group,data=my_data)

#结果表明，p=0.2371 >0.05，接受原假设，说明数据在不同水平下是等方