r语言算巢式设计方差分析_R语言进行方差分析

方差分析在试验科学中有重要的地位，今天谈谈如何用

R

做方差分析

前提假设：独立、正态、方差齐次(各水平间)

例

子

：

x

0,21.2)

数据集用

5

个因子水平测量，问是否存在差异

光是这样是无法进行分析的，对数据

x

进行格式转化

b

b

得到结果如下(

gl

指定因子，

5

是水平，

4

是重复次数)

x a

1

25.6 1

2

22.2 1

3

28.0 1

4

29.8 1

5

24.4 2

6

30.0 2

7

29.0 2

8

27.5 2

9

25.0 3

10 27.7 3

11 23.0 3

12 32.2 3

13 28.8 4

14 28.0 4

15 31.5 4

16 25.9 4

17 20.6 5

18 21.2 5

19 22.0 5

20 21.2 5

在进行方差分析之前先对几条假设进行检验，由于随机抽取，假设总体满足独立、正态，考察方差齐

次性(用

bartlett

检验)

> bartlett.test(x~a,data=b)

Bartlett test of homogeneity of variances

data:

x by a

Bartlett's K-squared = 7.0966, df = 4, p-value = 0.1309

方差齐次性符合

下面进行方差分析

m1

summary(m1)

Df Sum Sq Mean Sq F value Pr(>F)