Ecosistema urbano e trasporti pubblici in Sicilia

Caricamento librerie:

library(dplyr)
library(ggplot2)

Parte 1 : Dati

La redazione del quotidiano Sole 24 Ore ha fornito dati sulla qualità della vita in Italia nel 2020 e 2021, scaricabili da questo link : https://github.com/IlSole24ORE

L’Inail ha fornito dati sulle province e regioni italiane, scaricabili da questo link : https://dati.inail.it/opendata/elements/Provincia

Parte 2 : Esplorazione dati

ECOSISTEMA URBANO:

Gli indicatori di Ecosistema Urbano sono 18, sono riferiti alle città capoluogo e vengono utilizzati per i seguenti grafici. Derivano tutti da dati originali raccolti da Legambiente ad eccezione della capacità di depurazione e disponibilità di verde urbano (ISTAT), tasso di motorizzazione e incidenti stradali (ACI e ACI-ISTAT), uso efficiente del suolo (elaborato da Legambiente su dati ISPRA). Questi indicatori sono : BIOSSIDO DI AZOTO (NO2), POLVERI SOTTILI (Pm10), OZONO (O3),CONSUMI IDRICI DOMESTICI, DISPERSIONE DELLA RETE, CAPACITÀ DI DEPURAZIONE DELLE ACQUE REFLUE, PRODUZIONE DI RIFIUTI URBANI, RACCOLTA DIFFERENZIATA DEI RIFIUTI URBANI, PASSEGGERI DEL TRASPORTO PUBBLICO URBANO, OFFERTA DI TRASPORTO PUBBLICO URBANO, TASSO DI MOTORIZZAZIONE AUTO, INCIDENTALITÀ STRADALE, PISTE CICLABILI EQUIVALENTI, ISOLE PEDONALI, VERDE URBANO, ALBERI IN AREA URBANA, ENERGIE RINNOVABILI, USO EFFICIENTE DEL SUOLO.

Caricamento dati:

qualita <- read.csv("20201214_QDV2020_001.csv")
province <- read.csv("Provincia.csv", sep = ";")

qualita$INDICATORE <- as.factor(qualita$INDICATORE)
qualita$UNITA..DI.MISURA <- as.factor(qualita$UNITA..DI.MISURA)
df <- province[,c(3,6)]
colnames(qualita)[1]<- "DescrProvincia"
qualita <- merge(qualita,df, by.x = "DescrProvincia")
qualita$DescrRegione <- as.factor(qualita$DescrRegione)

ind <- "Ecosistema urbano"  

unita <-as.character(unique(qualita$UNITA..DI.MISURA[qualita$INDICATORE==ind]))
  
g<-qualita %>%
  filter(INDICATORE==ind) %>%
  group_by(DescrRegione) %>%
  summarise(totale=mean(VALORE)) %>%
  mutate(DescrRegione=reorder(DescrRegione,totale))  %>%
  ggplot(aes(DescrRegione,totale))+
  geom_bar(stat="identity", fill="green", colour="orange")+
  coord_flip()+
 geom_text(aes(label=round(totale,2)), hjust=0, size=2)+
  ylim(0,1)+
  guides(fill="none")+
  ylab(unita) +
  xlab("Regione") +
  ggtitle(ind, subtitle = unita)

print(g)

g<-qualita %>%
  filter(INDICATORE==ind, DescrProvincia %in% c("Agrigento","Caltanissetta","Catania","Enna","Messina","Palermo","Ragusa", "Siracusa", "Trapani")) %>%
  mutate(DescrProvincia=reorder(DescrProvincia,VALORE))  %>%
  ggplot(aes(DescrProvincia,VALORE))+
  geom_bar(stat="identity", fill="orange", colour="green")+
  coord_flip()+
 geom_text(aes(label=round(VALORE,2)), hjust=0, size=2)+
  ylim(0,1)+
  guides(fill="none")+
  ylab(unita) +
  xlab("Provincia") +
  ggtitle(ind, subtitle = unita)

 print(g) 

TRASPORTI:

Caricamento dati:

qualita <- read.csv("20211213_QDV2021_001.csv")
province <- read.csv("Provincia.csv", sep = ";")

qualita$INDICATORE <- as.factor(qualita$INDICATORE)
qualita$UNITA..DI.MISURA <- as.factor(qualita$UNITA..DI.MISURA)
df <- province[,c(3,6)]
colnames(qualita)[1]<- "DescrProvincia"
qualita <- merge(qualita,df, by.x = "DescrProvincia")
qualita$DescrRegione <- as.factor(qualita$DescrRegione)

ind <- "Offerta del trasporto pubblico" 

unita <-as.character(unique(qualita$UNITA..DI.MISURA[qualita$INDICATORE==ind]))
  
g<-qualita %>%
  filter(INDICATORE==ind) %>%
  group_by(DescrRegione) %>%
  summarise(totale=mean(VALORE)) %>%
  mutate(DescrRegione=reorder(DescrRegione,totale))  %>%
  ggplot(aes(DescrRegione,totale))+
  geom_bar(stat="identity", fill="red", colour="white")+
  coord_flip()+
 geom_text(aes(label=round(totale,2)), hjust=0, size=2)+
  theme(axis.title.x = element_text(size = 7), plot.subtitle =element_text(size = 8) )+
  ylim(0,50)+
  guides(fill="none")+
  ylab(unita) +
  xlab("Regione") +
  ggtitle(ind, subtitle = unita)

print(g)

g<-qualita %>%
  filter(INDICATORE==ind, DescrProvincia %in% c("Agrigento","Caltanissetta","Catania","Enna","Messina","Palermo","Ragusa", "Siracusa", "Trapani")) %>%
  mutate(DescrProvincia=reorder(DescrProvincia,VALORE))  %>%
  ggplot(aes(DescrProvincia,VALORE))+
  geom_bar(stat="identity", fill="orange", colour="white")+
  coord_flip()+
 geom_text(aes(label=round(VALORE,2)), hjust=0, size=2)+
  theme(axis.title.x = element_text(size = 7), plot.subtitle =element_text(size = 8) )+
  ylim(0,50)+
  guides(fill="none")+
  ylab(unita) +
  xlab("Provincia") +
  ggtitle(ind, subtitle = unita)

 print(g)