# NOT RUN {
data(icecream)
# Linear model, temperature as covariate
fit_ice <- lm(sales ~ temperature*location, data=icecream)
# Try to guess from coefficients where the sales were higher:
summary(fit_ice)
# What about now?
with(icecream, plot(temperature, sales, pch=19, col=location))
legend("topleft", levels(icecream$location), fill=palette())
# }
Run the code above in your browser using DataLab