讀取各車站資料，將屬同縣市的車站資料合併，取得每日各縣市臺鐵進出站人數。

library(dplyr)
fp <- system.file('railway-data', '2005-2018_carry.RDS', package = 'railwayViz')
carry <- readRDS(fp) %>%
  group_by(date, county) %>%
  summarise(num_in = sum(num_in), num_out = sum(num_out)) %>%
  ungroup()

將資料合併至臺灣縣市界的地理資料

fp <- system.file('railway-data', 'twcounty-shp.RDS', package = 'railwayViz')
taiwan <- readRDS(fp)

計算指標：除夕、平日以及都市化指標

carry_nye_wdays <- compare_nye2wdays(carry, 2018, groupby = 'county')
carry_nye_wdays <- to_sf(carry_nye_wdays, taiwan)