Two possible errors in code

[bakuninpr]

To get the number of deaths in 2016 Line 45 of excess_est reads

mutate(Sep = Sep*(1/3)) %>% sum

But to get the deaths in 2016 from September 20 to December 31, 11 days of September should be considered, not 10. Thus the line should be

mutate(Sep = Sep*(11/30)) %>% sum

Also, all 102/365 should be 103/365. The general conclusion of the article will not change, but point estimates and interval bounds will be shifted down by 70.