CasesDeathsOmikron

rm(list = ls())
library(tidyverse)
library(readr)
library(ggthemes)
library(lubridate)
library(zoo)
getwd()
data_folder <- file.path("")
url <- "https://covid.ourworldindata.org/data/owid-covid-data.csv"
name <- "owid-covid-data.csv"
download.file(url = url, destfile = paste0(data_folder,name))
setwd(data_folder)
Covidowid_covid_data <- read_csv(paste0(data_folder,name))
df<-Covidowid_covid_data
countries <- c(unique(df_owid$iso_code))
dfNLD <- df%>%
  dplyr::filter(iso_code == "NLD")%>%
  dplyr::select(date,iso_code,date,new_cases_per_million, new_deaths_per_million)%>%
  dplyr::mutate(cases_07da       = zoo::rollmean(new_cases_per_million, k = 7, fill = NA),
                deaths_07da      = zoo::rollmean(new_deaths_per_million, k = 7, fill = NA),
                deathdate_21plus = date - 21)
ggplot(dfNLD)+
  geom_line(aes(x=date, y=log(cases_07da), colour="New Cases"))+
  geom_line(aes(x=deathdate_21plus, y=log(deaths_07da), colour="New Deaths"))+
  scale_colour_manual(name="",
                      values=c('red', 'grey'), 
                      labels = c("Cases", "Deaths 21 days later"))+
  theme_bw()+
  theme(legend.position="bottom")+
  labs(x="Date", 
        y="New Cases and New Deaths (log Scale)", 
        title="New Cases vs New Deaths on a 7-day moving average log scale")