The weekly mortality data recently published by the Human Mortality Database can be used to explore seasonality in mortality rates. library(dplyr) library(tidyr) library(ggplot2) library(tsibble) library(feasts) Download the data We will first grab the latest data, using similar code to what I used in my recent post on “excess deaths”. However, this time we will keep the mortality rates rather than the numbers of deaths. mrates <- readr::read_csv("https://www.mortality.org/Public/STMF/Outputs/stmf.csv", skip = 1) %>% janitor::clean_names() %>% select(country_code:sex, …