---

œƒf[ƒ^“Η‚έ‚±‚έ
> mm <- read.table("Box5_R.tab")@@¦¨ƒf[ƒ^ƒtƒ@ƒCƒ‹FBox5_R.tabCBox5_R.data
> mm
REP V N DATA
1 1 V1 N1 2373
2 1 V1 N2 4076
3 1 V1 N3 7254
4 1 V2 N1 4007
5 1 V2 N2 5630
6 1 V2 N3 7053
7 1 V3 N1 2620
8 1 V3 N2 4676
9 1 V3 N3 7666
10 1 V4 N1 2726
11 1 V4 N2 4838
12 1 V4 N3 6881
13 1 V5 N1 4447
14 1 V5 N2 5549
15 1 V5 N3 6880
16 1 V6 N1 2572
17 1 V6 N2 3896
18 1 V6 N3 1556
19 2 V1 N1 3958
20 2 V1 N2 6431
21 2 V1 N3 6808
22 2 V2 N1 5795
23 2 V2 N2 7334
24 2 V2 N3 8284
25 2 V3 N1 4508
26 2 V3 N2 6672
27 2 V3 N3 7328
28 2 V4 N1 5630
29 2 V4 N2 7007
30 2 V4 N3 7735
31 2 V5 N1 3276
32 2 V5 N2 5340
33 2 V5 N3 5080
34 2 V6 N1 3724
35 2 V6 N2 2822
36 2 V6 N3 2706
37 3 V1 N1 4384
38 3 V1 N2 4889
39 3 V1 N3 8582
40 3 V2 N1 5001
41 3 V2 N2 7177
42 3 V2 N3 6297
43 3 V3 N1 5621
44 3 V3 N2 7019
45 3 V3 N3 8611
46 3 V4 N1 3821
47 3 V4 N2 4816
48 3 V4 N3 6667
49 3 V5 N1 4582
50 3 V5 N2 6011
51 3 V5 N3 6076
52 3 V6 N1 3326
53 3 V6 N2 4425
54 3 V6 N3 3214

> summary(mm)
REP V N DATA
Min. :1 V1:9 N1:18 Min. :1556
1st Qu.:1 V2:9 N2:18 1st Qu.:3970
Median :2 V3:9 N3:18 Median :5210
Mean :2 V4:9 Mean :5290
3rd Qu.:3 V5:9 3rd Qu.:6862
Max. :3 V6:9 Max. :8611

œˆφŽqŽw’θ
> mm$REP <- factor(mm$REP)
> mm$V <- factor(mm$V)
> mm$N <- factor(mm$N)

œBartlettŒŸ’θ
> bartlett.test(mm$DATA~mm$REP)

Bartlett test for homogeneity of variances

data: mm$DATA by mm$REP
Bartlett's K-squared = 0.4748, df = 2, p-value = 0.7887

> bartlett.test(mm$DATA~mm$V)

Bartlett test for homogeneity of variances

data: mm$DATA by mm$V
Bartlett's K-squared = 7.6346, df = 5, p-value = 0.1776

> bartlett.test(mm$DATA~mm$N)

Bartlett test for homogeneity of variances

data: mm$DATA by mm$N
Bartlett's K-squared = 7.0649, df = 2, p-value = 0.02923

œ‚Q—vˆφ Strip-plot •ͺŽU•ͺΝ
> fm <- aov(DATA~REP+V*N+Error(REP/V), data=mm) ¦…•½—vˆφ
> summary(fm)

Error: REP
Df Sum Sq Mean Sq
REP 2 9220962 4610481

Error: REP:V
Df Sum Sq Mean Sq F value Pr(>F)
V 5 57100201 11420040 7.6528 0.003372 **
Residuals 10 14922619 1492262
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
N 2 50676061 25338031 54.2579 1.245e-09 ***
V:N 10 23877979 2387798 5.1131 0.0005112 ***
Residuals 24 11207825 466993
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

> fm <- aov(DATA~REP+V*N+Error(REP/N), data=mm) ¦‚’Ό—vˆφ
> summary(fm)

Error: REP
Df Sum Sq Mean Sq
REP 2 9220962 4610481

Error: REP:N
Df Sum Sq Mean Sq F value Pr(>F)
N 2 50676061 25338031 34.069 0.003075 **
Residuals 4 2974908 743727
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
V 5 57100201 11420040 14.7956 2.454e-07 ***
V:N 10 23877979 2387798 3.0936 0.008028 **
Residuals 30 23155536 771851
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

> fm <- aov(DATA~REP+V*N+Error(REP/N*V), data=mm) ¦ŒπŒέμ—p—vˆφ
> summary(fm)

Error: REP
Df Sum Sq Mean Sq
REP 2 9220962 4610481

Error: V
Df Sum Sq Mean Sq
V 5 57100201 11420040

Error: REP:N
Df Sum Sq Mean Sq F value Pr(>F)
N 2 50676061 25338031 34.069 0.003075 **
Residuals 4 2974908 743727
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

Error: REP:V
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 10 14922619 1492262

Error: REP:N:V
Df Sum Sq Mean Sq F value Pr(>F)
V:N 10 23877979 2387798 5.8006 0.0004271 ***
Residuals 20 8232917 411646
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

œ•ͺŽU•ͺΝ•\i—v–ρj
Error: REP
Df Sum Sq Mean Sq
REP 2 9220962 4610481

Error: REP:V
Df Sum Sq Mean Sq F value Pr(>F)
V 5 57100201 11420040 7.6528 0.003372 **
Residuals 10 14922619 1492262

Error: REP:N
Df Sum Sq Mean Sq F value Pr(>F)
N 2 50676061 25338031 34.069 0.003075 **
Residuals 4 2974908 743727

Error: REP:N:V
Df Sum Sq Mean Sq F value Pr(>F)
V:N 10 23877979 2387798 5.8006 0.0004271 ***
Residuals 20 8232917 411646
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

---