data(airquality)
attach(airquality)
lmodel <- lm(Ozone ~ Solar.R + Wind)
lboot <- lm.boot(lmodel, R = 100)
## Get sampling distributions for coefficients
s <- samples(lboot, "coef")
## Histogram for the intercept
hist(s[1,])
Run the code above in your browser using DataLab