\begin{Schunk}
\begin{Sinput}
 m1 <- lm(EX ~ MET + GROW + YOUNG + OLD + WEST, data=dat)
 EXStar <- resid(m1)
 m2 <- lm(ECAB ~ MET + GROW + YOUNG + OLD + WEST, data=dat)
 ECABStar <- resid(m2)
 m3 <- lm(EXStar ~ ECABStar)
 hopePcor <- summary(m3)
 hopePcor$r.squared
\end{Sinput}
\begin{Soutput}
[1] 0.2097133
\end{Soutput}
\end{Schunk}
