# Confronto tra l'età della popolazione italiana in vari anni con boxplot
# Elaboro i dati (già classificati) seguendo quanto indicato qui
# e metto man mano il file rielaborato in eta1871, eta1921, ...
pop1871 <- c(3097,2917,5106,8005,5337,1528,812)
eta1871 <- c(   0,   5,  10,  20,  40,  60, 70, 100)
pop1921 <- c(3549,4051,7930,11023,7273,2445,1503)
eta1921 <- c(   0,   5,  10,   20,  40,  60, 70, 100)
pop1971 <- c(4428,4618,8031,15109,12939,5216,3796)
eta1971 <- c(   0,   5,  10,   20,  40,  60, 70, 100)
pop2009 <- c(2845,2834,5766,15639,12938,6843,3781, 15)
eta2009 <- c(   0,   5,  10,  20,   40,  60,  70, 100, 105)
freq <- pop1871; interv <- eta1871
dev.new(); source("http://macosa.dima.unige.it/R/daticlas.txt")
title(1871); dat1871 <- XxXx
freq <- pop1921; interv <- eta1921
dev.new(); source("http://macosa.dima.unige.it/R/daticlas.txt")
title(1921); dat1921 <- XxXx
freq <- pop1971; interv <- eta1971
dev.new(); source("http://macosa.dima.unige.it/R/daticlas.txt")
title(1971); dat1971 <- XxXx
freq <- pop2009; interv <- eta2009
dev.new(); source("http://macosa.dima.unige.it/R/daticlas.txt")
title(2009); dat2009 <- XxXx

dev.new(); etichette <- c("1871","1921","1971","2009")
boxplot(dat1871,dat1921,dat1971,dat2009,horizontal=TRUE,range=0,names=etichette,col="yellow",las=1)
abline(v=axTicks(1),lty=3)