---

●データ読みこみ
> mm <- read.table("Box1_R.tab")　　※→データファイルBox1_R.tab

［※行番号を含まず，見出し行のみを先頭行にもつ，より単純なデータファイルBox1_R.dataを入力データとする場合は「> mm <- read.table("Box1_R.data" header=T)」となる．→参照：データ入力方法］

> mm
TRT DATA
1 DM1 2537
2 DM1 2069
3 DM1 2104
4 DM1 1797
5 DM2 3366
6 DM2 2591
7 DM2 2211
8 DM2 2544
9 DDT 2536
10 DDT 2459
11 DDT 2827
12 DDT 2385
13 AZO 2387
14 AZO 2453
15 AZO 1556
16 AZO 2116
17 DB 1997
18 DB 1679
19 DB 1649
20 DB 1859
21 DK 1796
22 DK 1704
23 DK 1904
24 DK 1320
25 Con 1401
26 Con 1516
27 Con 1270
28 Con 1077

●因子指定
> mm$TRT <- factor(mm$TRT)

●Bartlett検定
> bartlett.test(mm$DATA~mm$TRT)

Bartlett test for homogeneity of variances

data: mm$DATA by mm$TRT
Bartlett's K-squared = 5.559, df = 6, p-value = 0.4744

●１要因完全無作為法の分散分析
> fm <- aov(DATA~TRT, data=mm)
> summary(fm)
Df Sum Sq Mean Sq F value Pr(>F)
TRT 6 5587175 931196 9.8255 3.329e-05 ***
Residuals 21 1990237 94773
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

●多重比較
> pairwise.t.test(mm$DATA, mm$TRT, p.adj="bonferroni")

Pairwise comparisons using t tests with pooled SD

data: mm$DATA and mm$TRT

AZO Con DB DDT DK DM1
Con 0.02596 - - - - -
DB 1.00000 0.81348 - - - -
DDT 1.00000 0.00026 0.04786 - - -
DK 1.00000 1.00000 1.00000 0.01364 - -
DM1 1.00000 0.02632 1.00000 1.00000 1.00000 -
DM2 0.41194 7e-05 0.01206 1.00000 0.00341 0.40684

P value adjustment method: bonferroni

---