
342
19
章 モデル構築
19.3.1
曜日
このプロットには、非常に強い曜日効果があるので、長期的傾向が掴みづらくなっています。まず
曜日ごとの便数の分布を調べることから始めます。
daily <- daily %>%
mutate(wday = wday(date, label = TRUE))
ggplot(daily, aes(wday, n)) +
geom_boxplot()
ほとんどの旅行が商用なので、週末の便数は少なくなっています。この効果は土曜日効果と呼ば
れます。月曜の会議のために日曜に出発することがありますが、土曜日は家族と家で過ごすことが多
いので、土曜出発は稀です。
この強いパターンを除去する1 つの方法がモデルを使うことです。まず、モデルを適合させて、そ
の予測を元のデータに重ねて表示します。
mod <- lm(n ~ wday, data = daily)
grid <- daily %>%
data_grid(wday) %>%
add_predictions(mod, "n")
ggplot(daily, aes(wday, n)) +
geom_boxplot() +
geom_point(data = grid, color = "red", size = 4)