source('../../Data/get_alcohol_data.R') DF = get_alcohol_data() m = lm(alcohol ~ tobacco, data=DF) print(summary(m)) plot(DF$tobacco, DF$alcohol, 'p', xlab='Tobacco', ylab='Alcohol') grid() abline(m)