ELABORAZIONE DATI CON R
Fonte: ISTAT, DATI. "Dati ISTAT." (2021)
https://www.istat.it/it/archivio/258104;
Tavole_dati_Indicatori_agro_ambientali_AEIs_2010_19_20_05_2021
https://www.istat.it/it/archivio/258104;
Tavole_dati_Indicatori_agro_ambientali_AEIs_2010_19_20_05_2021
Per il seguente lavoro sono state utilizzate tavole Excel ricavate dall'archivio ISTAT (pubblicamente fruibili).
L'obbiettivo è stato quello di estrarre ed elaborare i dati per poi creare grafici riassuntivi.
I dati riguardano i quantitativi assoluti(in kg) di fungicidi e insetticidi utilizzati nel periodo 2011-2017 in Italia;
La prima parte è relativa alle quantità totali usate nell'intero paese (anno per anno), suddivisa in fungicidi e insetticidi.
La seconda, invece, riguarda la media degli utilizzi annuali, divisi regione per regione, nel periodo complessivo di riferimento.
La prima parte è relativa alle quantità totali usate nell'intero paese (anno per anno), suddivisa in fungicidi e insetticidi.
La seconda, invece, riguarda la media degli utilizzi annuali, divisi regione per regione, nel periodo complessivo di riferimento.
SCRIPT
#imposto directory
setwd('./Esercizi/utilizzo fungicidi')
getwd()
setwd('./Esercizi/utilizzo fungicidi')
getwd()
#carico file Excel(fonte: ISTAT)
dati = './Esercizi/utilizzo fungicidi/utilizzo fungicidi.xlsx'
install.packages("readxl")
library(readxl)
dati = read_excel(dati)
dati = read_excel(dati)
#apporto modifiche al dataframe per adeguarlo alle operazioni successive
dati = dati[1:33,]
dati = dati[c(4:25,33),c(1,8,10,14,16,20,22,26,28,32,34,38,40,44,46)]
dati = dati[c(4:25,33),c(1,8,10,14,16,20,22,26,28,32,34,38,40,44,46)]
nomi_colonne = c('regione','2011 fungicidi', '2011 insetticidi/acaricidi','2012 fungicidi','2012 insetticidi/acaricidi','2013 fungicidi','2013 insetticidi/acaricidi','2014 fungicidi','2014 insetticidi/acaricidi','2015 fungicidi','2015 insetticidi/acaricidi','2016 fungicidi','2016 insetticidi/acaricidi','2017 fungicidi','2017 insetticidi/acaricidi')
length(nomi_colonne)
colnames(dati)<-nomi_colonne
nomi_righe = dati$regione
rownames(dati)<-nomi_righe
dati$regione<-NULL
rownames(dati)<-nomi_righe
rownames(dati)<-nomi_righe
dati = as.data.frame(dati)
dati = sapply(x, as.integer)
#grafico fungicidi anno per anno
ascisse_grafici = c('2011','2012','2013','2014','2015','2016','2017')
barplot(dati[,c(1,3,5,7,9,11,13)],
main = 'utilizzo fungicidi totale Italia',
names.arg = ascisse_grafici,cex.names = 0.8,
ylab = 'kg prodotti',col = 'green')
main = 'utilizzo fungicidi totale Italia',
names.arg = ascisse_grafici,cex.names = 0.8,
ylab = 'kg prodotti',col = 'green')
#grafico insetticidi anno per anno
barplot(dati[,c(2,4,6,8,10,12,14)],
main = 'utilizzo insetticidi totale Italia',
names.arg = ascisse_grafici, cex.names = 0.8,
ylab = 'Kg prodotti', col = 'red')
main = 'utilizzo insetticidi totale Italia',
names.arg = ascisse_grafici, cex.names = 0.8,
ylab = 'Kg prodotti', col = 'red')
#GRAFICO QUANTITATIVI REGIONI: per ogni regione calcolo la media delle quantità utilizzate di p.a (fung+ inst/ac) nell'intero periodo (2011-2017)
#creo nuova tabella dati per poi calcolarci le medie
#creo nuova tabella dati per poi calcolarci le medie
regioni_tot_prodotti = data.frame('Piemonte' = x[1,], "Valle d'Aosta"= x[2,],'Liguria' = x[3,],
'Lombardia'= x[4,], 'Trentino Alto Adige'= x[5,],
'Veneto' =x[8,], 'Friuli Venzia Giulia'= x[9,],
'Emilia Romagna'= x[10,], 'Toscana' = x[11,], 'Umbria'= x[12,],
'Marche'= x[13,],'Lazio'= x[14,], 'Abruzzo'= x[15,],
'Molise'=x[16,], 'Campania'= x[17,], 'Puglia'= x[18,],
'Basilicata'= x[19,], 'Calabria'= x[20,], 'Sicilia'= x[21,],
'Sardegna'= x[22,])
#calcolo delle medie
'Lombardia'= x[4,], 'Trentino Alto Adige'= x[5,],
'Veneto' =x[8,], 'Friuli Venzia Giulia'= x[9,],
'Emilia Romagna'= x[10,], 'Toscana' = x[11,], 'Umbria'= x[12,],
'Marche'= x[13,],'Lazio'= x[14,], 'Abruzzo'= x[15,],
'Molise'=x[16,], 'Campania'= x[17,], 'Puglia'= x[18,],
'Basilicata'= x[19,], 'Calabria'= x[20,], 'Sicilia'= x[21,],
'Sardegna'= x[22,])
#calcolo delle medie
regioni_medie= sapply(regioni_tot_prodotti, mean)
regioni_medie = as.integer(regioni_medie)
#creazione grafico
nomi_righe_r = c('Piem', "V.d.A", 'Lig','Lomb','Tr.Al.Ad',
'Ven', 'F.V.G','E-Rom', 'Tosc','Umbr','Marc','Laz',
'Abr','Mol', 'Camp','Pug','Basil','Cal','Sic',
'Sard')
'Ven', 'F.V.G','E-Rom', 'Tosc','Umbr','Marc','Laz',
'Abr','Mol', 'Camp','Pug','Basil','Cal','Sic',
'Sard')
barplot(regioni_medie, main = 'Medie p.a per regione 2011-2017',
names.arg = nomi_righe_r, cex.names = 1,
cex.axis = 0.8, ylim = c(0,3000000),
las = 2, col = rainbow(20))
names.arg = nomi_righe_r, cex.names = 1,
cex.axis = 0.8, ylim = c(0,3000000),
las = 2, col = rainbow(20))
