\begin{Schunk}
\begin{Sinput}
 plot(EX  ~  MET, data=dat, xlab="Metropolitan Population (%)",ylab="Public Expenditures (Per Capita, 1960)")
 pemod3 <- lm(EX ~ poly(MET,2), data=dat)
 rmet <- range(dat$MET)
 smet <- seq( rmet[1], rmet[2], length.out=100)
 polypred <- predict(pemod3, newdata=data.frame(MET=smet))
 lines(smet, polypred)
\end{Sinput}
\end{Schunk}
