r语言算巢式设计方差分析_逼迫学习5：用R语言做单因素方差分析及多重比较

MyPlan5

时间：6.10-6.16周内完成

SPSS在方差分析的应用已经做得非常好了，绝大多数的方差分析问题均可通过SPSS“点菜单”的方式得以解决，R语言在统计和可视化方面有自己的特色，我们不妨来对比着学习。

完成情况

选用R语言自带案例数据集PlantGrowth，研究两个处理和一个对照组对植物产量的影响，每组10例共3记录，主要考察处理对提高植物产量有无影响。

数据构成：因变量weight，因子变量group，三个水平依次为ctrl、trt1、trt2。

一、数据正态性检验

用夏皮罗-威尔克检验3组数据是否服从正态分布，样本容量3~5000均可，这一点比SPSS说得更明确。

data

shapiro.test(data[1:10])

shapiro.test(data[11:20])

shapiro.test(data[21:30])

ctrl、trt1、trt2三个分组的shapiro.test检验概率p依次为0.7475、0.4519、0.5643，均大于0.05，原假设成立(H0：假设数据服从正态分布)，表明3组数据均来自正态分布总体。

二、方差齐次检验

用bartlett.test检验3个分组数据方差是否一致。

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

直接看p-value = 0.2371>0.05，原假设成立(H0：假设3组数据方差相等)，表明3组数据的方差齐次。<