\begin{Schunk}
\begin{Sinput}
 dat$xlog <- log(dat$x)
 dat$ylog <- log(dat$y)
 m1 <- lm(ylog ~ xlog, data = dat)
 head(model.matrix(m1))
\end{Sinput}
\end{Schunk}
