Caricamento librerie

library(dplyr)
library(ggplot2)

Parte 1: Dati

La Protezione Civile fornisce quotidianamente i dati regionali e nazionali del Covid-19 in Italia da questo link : https://github.com/pcm-dpc/COVID-19

Caricamento dei dati:

regioni <- read.csv("dpc-covid19-ita-regioni.csv")
italia <- read.csv("dpc-covid19-ita-andamento-nazionale.csv")
province <- read.csv("dpc-covid19-ita-province.csv")

Parte 2: Esplorazione dati

df <-regioni %>%
  filter(denominazione_regione=="Sicilia")

n<-nrow(df)
df <- cbind(df,c=1:n)

v<-rep(NA,n)
v[2:n]<-diff(df$casi_testati)

df <- cbind(df,v=df$nuovi_positivi*100/v)
df %>%
  ggplot(aes(x=as.Date(data),y=terapia_intensiva, group=1))+
  geom_line(color="red")+
  #geom_text(aes(label=terapia_intensiva), vjust=0, size=2)+
  theme(axis.text.x = element_text(angle=45,hjust=1))+
  xlab("Giorni")+
  ggtitle("Terapie intensive per Covid-19 in Sicilia", subtitle = "al 7 novembre 2020")

df %>%
  filter(c>190) %>%
  ggplot(aes(x=as.Date(data),y=terapia_intensiva, group=1))+
  geom_line(color="red")+
  geom_text(aes(label=terapia_intensiva), vjust=0, size=2)+
  theme(axis.text.x = element_text(angle=45,hjust=1))+
  xlab("Giorni")+
  ggtitle("Terapie intensive per Covid-19 in Sicilia", subtitle = "al 7 novembre 2020")

df %>%
  ggplot(aes(x=as.Date(data),y=ricoverati_con_sintomi, group=1))+
  geom_line(color="red")+
  #geom_text(aes(label=ricoverati_con_sintomi), vjust=0, size=2)+
  theme(axis.text.x = element_text(angle=45,hjust=1))+
  xlab("Giorni")+
  ggtitle("Ricoverati con sintomi per Covid-19 in Sicilia", subtitle = "dal 24 febbraio al 7 novembre 2020")

df %>%
  filter(c>190) %>%
  ggplot(aes(x=as.Date(data),y=ricoverati_con_sintomi, group=1))+
  geom_line(color="red")+
  geom_text(aes(label=ricoverati_con_sintomi), vjust=0, size=2)+
  theme(axis.text.x = element_text(angle=45,hjust=1))+
  xlab("Giorni")+
  ggtitle("Ricoverati con sintomi per Covid-19 in Sicilia", subtitle = "al 7 novembre 2020")

df %>%
  filter(c>170) %>%
  ggplot(aes(x=as.Date(data),y=v, group=1))+
  geom_line(color="orange")+
  #geom_smooth(method = 'loess', formula= 'y ~ x') +
  geom_text(aes(label=round(v,2)), vjust=0, size=2)+
  theme(axis.text.x = element_text(angle=45,hjust=1))+
  xlab("Giorni")+
  ylab("Percentuale %")+
  ggtitle("Percentuale di nuovi positivi rispetto a nuovi tamponi di Covid-19 in Sicilia", subtitle = "da giugno al 7 novembre 2020, esclusi i tamponi di controllo")