#行の中で、「#」以降は、コメントアウトで処理しません。

#====================================================

#表2-3, 2-4の分析

#Cドライブにv04dataというフォルダを作り、v04c2table2-1.csvを置いているという設定）

#データの読み込み

setwd("c:/v04data")

data1=read.table("v04c2table2-1.csv",header=T,sep=",")

#変数のfactor化（1,2でコード化された量的変数を質的変数に変換する）

data1$treatment <- factor(data1$treatment)

data1$aptitude <- factor(data1$aptitude)

#表2-3 母数の推定値

summary(lm(score~1+aptitude+treatment+aptitude*treatment,data=data1))

#表2-4 分散分析表

summary(aov(score~aptitude+treatment+aptitude:treatment,data=data1),intercept=TRUE)

#表2-7の分析

#Cドライブにv04dataというフォルダを作り、v04c2table2-5.csvを置いているという設定）

#データの読み込み

setwd("c:/v04data")

data2=read.table("v04c2table2-5.csv",header=T,sep=",")

#図2-5の描画

#講義形式だけのデータ

dataclass <- subset(data2, subset = data2$classroom.e.learning == 0)

#e-Learningだけのデータ

dataelearning <- subset(data2, subset = data2$classroom.e.learning == 1)

#軸の最大値

xmax <- max(data2$minute)

ymax <- max(data2$score)

#プロット

plot(c(0,xmax),c(0,ymax),xlab="利用時間（分）",ylab="成績向上点")

par(new=T)

plot(dataclass$minute,dataclass$score,xlim=c(0,xmax),ylim=c(0,ymax),ann=F,axes=F,pch=1)

par(new=T)

plot(dataelearning$minute,dataelearning$score,xlim=c(0,xmax),ylim=c(0,ymax),ann=F,axes=F,pch=2)

#表2-7 重回帰分析（交互作用項なし）

reg0 <- lm(score ~ minute + classroom.e.learning, data=data2)

summary(reg0)

confint(reg0,level=0.95)

#表2-9 重回帰分析（交互作用項あり）

#TypeIIIの平方を出すために，carパッケージのAnova関数を利用します

#事前にcarパッケージをインストールしてください

library(car)

options(contrasts = c("contr.sum", "contr.sum"))

reg1 <- lm(score ~ 1 + minute + classroom.e.learning + minute * classroom.e.learning, data=data2)

summary(reg1)

confint(reg1,level=0.95)

#表2-10 分散分析表（テキストと数値が異なりますが，表2-9と検定結果が一致します）

Anova(reg1,type="III")

#偏η^2の出力（事前にlsrパッケージをインストールしてください）

library(lsr)

etaSquared(reg1, type=3)

#ダミー変数を-1,1で処理すると，テキストの表2-10が得られます。

#library(car)

#テキストの出力にあわせるため，ダミー変数を-1,1にします

data3 <- data2

data3$classroom.e.learning <- ifelse(data2$classroom.e.learning==0,-1,1)

reg2 <- lm(score ~ minute + classroom.e.learning + minute * classroom.e.learning, data=data3)

Anova(reg2,type="III")

library(lsr)

etaSquared(reg2, type=3)

#====================================================

#2015.11.08 熊谷龍一・荘島宏二郎