（我手动微调了日期，以便在图中获得较好的分段效果。使用可视化来帮助你理解函数的行为，是一种非常强大且通用的技术。）

观察这个新变量如何影响一周中的其他日子是很有用的：

看起来不同学期之间存在显著差异，因此为每个学期分别拟合一个星期几效应是合理的。这改进了我们的模型，但没有达到我们期望的程度：

mod1 <- lm(n ~ wday, data = daily) mod2 <- lm(n ~ wday * term, data = daily) daily %>% gather_residuals(without_term = mod1, with_term = mod2) %>% ggplot(aes(date, resid, color = model)) + geom_line(alpha = 0.75)