対照区
施設名 1区 2区 3区 4区 5区
A 11.8 15.3 17.6 16.4 15.1
B 12.1 10.5 18.3 16.1 14.6
C 16.3 15.4 18.1 17.2 16.2
試験区
施設名 1区 2区 3区 4区 5区
A 22.1 21 19.5 22.3 23.1
B 18.3 18.3 20.1 23 15.9
C 20.1 20.2 18.1 16.3 22.8
まず次のようにデータを変換し、data.txt などに保存します(変換は R でやらず EXCEL などでやったほうが簡単です)。
F1 | F2 | Y |
taishouku | A | 11.8 |
taishouku | B | 12.1 |
taishouku | C | 16.3 |
taishouku | A | 15.3 |
taishouku | B | 10.5 |
taishouku | C | 15.4 |
taishouku | A | 17.6 |
taishouku | B | 18.3 |
taishouku | C | 18.1 |
taishouku | A | 16.4 |
taishouku | B | 16.1 |
taishouku | C | 17.2 |
taishouku | A | 15.1 |
taishouku | B | 14.6 |
taishouku | C | 16.2 |
shikenku | A | 22.1 |
shikenku | B | 18.3 |
shikenku | C | 20.1 |
shikenku | A | 21 |
shikenku | B | 18.3 |
shikenku | C | 20.2 |
shikenku | A | 19.5 |
shikenku | B | 20.1 |
shikenku | C | 18.1 |
shikenku | A | 22.3 |
shikenku | B | 23 |
shikenku | C | 16.3 |
shikenku | A | 23.1 |
shikenku | B | 15.9 |
shikenku | C | 22.8 |
あとは次のような R コマンドを実行すればOKです。
data <- read.table("data.txt", header=T) summary(aov(Y~F1+F2,data=data))
素晴らしい。ありがとうございます。ABC間の交互作用を求める入力過程もお願い致します。