\begin{Schunk}
\begin{Sinput}
 Nsampsize <- 100
 b0 <- 3; b1 <- 0.25; stde <- 10; mult <- 0.05
 x <- rnorm(Nsampsize, m=50, sd=10)
 meanx <- mean(x)
 minx <- min(x)
 e <- rnorm(Nsampsize)
 y <- b0 + b1 * x + mult*(x-minx)*stde*e
 mod <- lm(y~x)
 plot(y~x, main="Heteroskedasticity")
 abline(mod)
\end{Sinput}
\end{Schunk}
