Ice3 <-
Ice2 %>% filter(time == 1930) %>%
spread(location, temp)
anova(lm(surface - intramuscular ~ treatment,
data = Ice3 %>% filter(phase == "t")))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.