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 2021, scaricabili da questo link : https://github.com/IlSole24ORE/QDV2021

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

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)

Parte 2 : Esplorazione dati

indici <- c("Qualità della vita degli anziani","Qualità della vita dei bambini","Qualità della vita dei giovani","Qualità della vita delle donne" )

for (ind in indici) {
  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,650)+
  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,650)+
  guides(fill="none")+
  ylab(unita) +
  xlab("Provincia") +
  ggtitle(ind, subtitle = unita)

 print(g) 
 
}