■ 分析前にやっておくこと

(1) Rのインストール
(2) 作業ディレクトリの確認と変更をしておく
(3) 作業ディレクトリにデータを入れておく
(4) 必要なパッケージをインストールしておく

 

---

■ 繰り返しのない二元配置分散分析（ANOVA君を使った方法）

データの読み込み

書籍使用データ（右クリックで「別名で保存」して作業ディレクトリに入れる）

dat <- read.csv("ch07-1.csv", header=TRUE, fileEncoding="CP932")
# fileEncodingはMacの場合入れる

もしくは，Excel を開いてデータ範囲をコピーしてから読み込む
（クリップボードからのデータ読み込み）

Macの場合

dat <- read.table(pipe("pbpaste"), header=TRUE)  
# header=TRUEで変数名を含む，FALSEにすると含まない

Windowsの場合

dat <- read.table("clipboard", header=TRUE) 
# header=TRUEで変数名を含む，FALSEにすると含まない

 

ANOVA君を使えるようにする

これを右クリックで「別名で保存」して作業ディレクトリに入れて
以下のように入力

source("anovakun_420.txt")

 

繰り返しのない二元配置分散分析の実行

anovakun(dat, "ABs", 2, 3, holm = T, peta = T)

“ABs”は被験者間を示し（sの左側が被験者間要因，右側が被験者内要因）
次の2と3で，要因Aが2水準，要因Bが3水準と指定
holm = T でHolmの方法による多重比較，peta = T で偏イータ2乗を出力
その他のオプションは， ANOVA君のサイトで確認

 

グラフ作図

dat$教室 などと書かなくていいようにattach

attach(dat)

要因を因子の型に変更

教室 <- factor(教室, labels=c("普通教室","CALL教室"))
クラスサイズ <- factor(クラスサイズ, labels=c("10人","20人","40人"))

par(family=””)でMacの場合，グラフに日本語表記をできるようにする
参照 http://oku.edu.mie-u.ac.jp/~okumura/stat/graphs.html

par(family="HiraKakuPro-W3")  #Macの場合

interaction.plot(クラスサイズ, 教室, 点数, type="b",  pch=c(1,2))

type=”b”は点と線のプロット，pch=c(1,2) でマーカーの形式指定
参照 http://cse.naro.affrc.go.jp/takezawa/r-tips/r/53.html

 

 

---

■ 繰り返しのある二元配置分散分析（ANOVA君を使った方法）

データの読み込み

書籍使用データ （右クリックで「別名で保存」して作業ディレクトリに入れる）

dat <- read.csv("ch07-2.csv", header=TRUE, fileEncoding="CP932")
# fileEncoding="CP932"はMacの場合入れる

もしくは，Excel を開いてデータ範囲をコピーしてから読み込む
（クリップボードからのデータ読み込み）

Macの場合

dat <- read.table(pipe("pbpaste"), header=TRUE)  
# header=TRUEで変数名を含む，FALSEにすると含まない

Windowsの場合

dat <- read.table(“clipboard", header=TRUE) 
# header=TRUEで変数名を含む，FALSEにすると含まない

 

ANOVA君を使えるようにする

これを右クリックで「別名で保存」して作業ディレクトリに入れて
以下のように入力

source("anovakun_420.txt")

 

繰り返しのある二元配置分散分析の実行

anovakun(dat, "AsB", 2, 3, mau=T, auto=T, holm = T, peta=T)

“AsB”のAは被験者間，Bは被験者内（sの左側が被験者間，右側が被験者内）の指定
次の2と3で，要因Aが2水準，要因Bが3水準と指定
mau = T とするとMauchlyの球面性検定
auto = T とすることで必要があればGreenhouse-Geisserの調整を行う
holm = T でHolmの方法による多重比較，peta = T で偏イータ2乗を出力
その他のオプションは， ANOVA君のサイトを確認

 

グラフ作図
スタック形式に変更
参照 http://133.100.216.71/R_analysis/basic_data_frame01.html

x <- stack(dat[2:4])
y <- data.frame(dat$指導法, x)

要因を因子の型に変更

y$dat.指導法 <- factor(y$dat.指導法, labels=c("処置群","対照群"))

名義変数の順序を指定（グラフ作図用)

y$ind <- factor(y$ind,levels=c("事前","事後","遅延"))

変数名を戻す

names(y) <- c("指導法", "点数", "時間")

y$指導法 などと書かなくていいようにattach

attach(y)

par(family=””)でMacの場合，グラフに日本語表記をできるようにする
参照 http://oku.edu.mie-u.ac.jp/~okumura/stat/graphs.html

par(family="HiraKakuPro-W3")

interaction.plot(時間, 指導法, 点数, type="b",  pch=c(1,2))

type=”b”は点と線のプロット，pch=c(1,2) でマーカーの形式指定
参照 http://cse.naro.affrc.go.jp/takezawa/r-tips/r/53.html

 

 

---

■ Rcommanderを使う場合

繰り返しのない二元配置分散分析をRコマンダーで行う方法は，
『Rによる心理・調査データ解析』緒賀 郷志 (著) に説明がある。

ただし，繰り返しのある二元配置分散分析は R コマンダーでは
実行できないため，上記に説明している，コンソールでコードを
書いて行うか，ANOVA4 on the Web を使うとよい。