x <- rnorm(100) # 100 random numbers from a normal(0,1) distribution y <- exp(x) + rnorm(100) # an exponential function with error result <- lsfit(x,y) # regress x on y and store the results ls.print(result) # print the regression results plot(x,y) # pretty obvious what this does abline(result) # add the regression line to the plot lines(lowess(x,y), col=2) # add a nonparametric regression line (a smoother) hist(result$residuals) # histogram of the residuals from the regression

Results