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


# lm は　線形モデル

model4 <- lm(Y ~ X + I(X^2) + I(X^3) + I(X^4), data =d)   # 4次の多項式回帰モデル
＃model4 <- lm(Y ~ poly(X, degree = 4, raw = TRUE) , data=d)



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

dredge(model4, rank="AIC")



P <- predict(model4)
plot(d$Y ~ d$X)
reorder <- order(d$X)
lines(d$X[reorder], P[reorder])