base.lmadd <-
lm(temp ~ location + treatment, data = Ice1930 %>% filter(phase == "b"))
anova(base.lmadd)
plot(base.lmadd, w = c(5, 2))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.