x72a <- read.csv("Table722a.csv")
x72a

#変数名を一部変えました（エラー回避のため）
attach(x72a)

​

#関数　t.test(従属変数~独立変数,等分散の仮説,対応の有無)

t.test(Il~Sr,data=x72a,var.equal=T,paired=TRUE)

​

tapply(Il,Sr,mean)

tapply(Il,Sr,sd)

​

#参考）効果量を求めるため  分散、n を求める

tapply(Il,Sr,var)

tapply(Il,Sr,length)

​

d2a <- abs(17.0-16.5)/sqrt((12*7.78+12*10.6)/(12+12))
d2a

​