#-- 第9章 差の分析 --#

# データの読み込み
data<- read.csv("9_anova.csv", header=TRUE, fileEncoding="CP932")
head(data)

# 二元配置分散分析
res<- aov(信用度~国籍*年代,data=data)
summary(res)

# 多重比較検定
TukeyHSD(aov(信用度~年代,data=data))