# Ozone

plot(air$ozone, air$temperature)
lines(ksmooth(air$ozone, air$temperature, kernel = "triangle", bandwidth = 0.2))
lines(ksmooth(air$ozone, air$temperature, kernel = "triangle", bandwidth = 0.5))

ozone.lm <- lm(ozone ~ temperature, data = air)

plot(air$ozone, air$temperature)
lines(ozone.lm$fitted, air$temperature)

# time series

r.norm <- rnorm(100)
x.cor <- ts(r.norm[1:98] + r.norm[2:99] + r.norm[3:100])

par(mfrow = c(1,2))
tsplot(x.cor)
acf(x.cor)

# Faces

faces(t(cereal.attitude), labels = dimnames(cereal.attitude)[[2]], 

      ncol = 3)