a sequence of progressively more complex models:
R code:
o = lm(y ~ diet + drug + diet:drug, data = data1) anova(o)