\begin{Schunk}
\begin{Sinput}
 dat$xctr <- scale(dat$x, scale = FALSE) 
 m1 <- lm(y ~ xctr + I(xctr^2), data = dat)
\end{Sinput}
\end{Schunk}
