attach(airquality) shapiro.test(Wind) shapiro.test(Ozone) layout(matrix(c(1,3,2,4),2,2)) hist(Wind) hist(Ozone) cOzone <- log(Ozone+10) shapiro.test(cOzone) plot(cOzone~Wind,yaxt="n",ylab="Ozone",ylim=c(0,log(210)), main="1973年のNYにおけるオゾン濃度\n(Ozone)と風速(Wind)との関係") yi <- 0:4*50 axis(2,log(yi+10),yi) res <- lm(cOzone~Wind) X <- data.frame(Wind=seq(min(Wind),max(Wind),length=20)) Y <- predict(res,X,interval="predict") matlines(X,Y,col=1,lty=c(1,2,2)) plot(exp(fitted.values(res))-10,residuals(res),main="残差プロット",xlab="Ozone (expected)") exp(predict(res,list(Wind=15),interval="confidence"))-10 exp(predict(res,list(Wind=25),interval="confidence"))-10 summary(res) detach(airquality)