d <- read.table("Ex20.1.txt", header=T)

cor(d)     #相関係数行列
plot(d)    #2変量間の散布図



# lm は　線形モデル

model1 <- lm(Y ~ X1 + X2 + X3 + X4, data =d)
summary(model1)



step(model1)   #step関数によるモデル選択




library(MuMIn)
options(na.action = "na.fail") 　　#

dredge(model1, rank="AIC")