\begin{Schunk}
\begin{Sinput}
 set.seed(12323)
 x <- rnorm(100, m=50,s=10)
 y <-  2 + 3 *x + rnorm(100,m=0, s=20)
 y2 <- 2 + 3 * x + rnorm(100, m=10, s=80)
 mod1 <- lm(y ~ x)
 mod2 <- lm(y2 ~ x)
 plot(x,y, main="", xlab="Indep. Var.", ylab="Dep. Var",ylim=1.2*range(y2))
 abline(2,3,lty=4, col="red")
 abline(mod1)
 legend("bottomright",legend=c("True Relationship","Sample Estimate"), lty=c(4,1),col=c("red","black"))
\end{Sinput}
\end{Schunk}
