\begin{Schunk}
\begin{Sinput}
 dat <- read.table("cystfibr.txt", header=T)
 plot(weight ~ height, data = dat) 
 mod1 <- lm(weight ~ height, data = dat)
 summary(mod1)
 abline(mod1)
\end{Sinput}
\end{Schunk}
