# Calcoli svolti seguendo questo esempio
# morti (in migliaia) per classi di eta` in Italia nel 1881-1890
mor <- c(3818,343,303,398,360,384,495,1177,708)
interv <- c(0,5,10,20,30,40,50,60,75,100) # un dato in piu' di mor
datiN <- function(pop,interv) {
totale <- sum(pop); fr_perc <- pop/totale*100
rap <- 1e5/sum(pop); freq <- pop*rap; num_int <- length(freq)
dati <- c(seq(interv[1],interv[2],by=(interv[2]-interv[1])/freq[1]))
for (i in 2:num_int) {dati <- c(dati,seq(interv[i],interv[i+1],by=(interv[i+1]-interv[i])/freq[i]))};
dati}
dati1 <- datiN(mor,interv)
# morti (in centinaia) per classi di eta` in Italia nel 1951
mor <- c(729,35,77,132,134,285,457,1401,1569)
dati2 <- datiN(mor,interv)
# morti (in centinaia) per classi di eta` in Italia nel 1988
mor <- c(68,10,31,66,70,161,423,1516,2983)
dati3 <- datiN(mor,interv)
# morti italiani nel 2006 (ogni 10000)
mor <- c(43, 5, 6,16,23,26,28,37,56,88,141,225,350,
518,814,1245,1771,2027,1643,727,195, 17)
interv <- c(0,5,10,15,20,25,30,35,40,45,50,55, 60, 65,
70, 75, 80, 85, 90, 95,100,105,110)
dati4 <- datiN(mor,interv)
boxplot(dati1,dati2,dati3,dati4,horizontal=TRUE,range=0)
summary(dati1); summary(dati2); summary(dati3); summary(dati4)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.000 2.615 7.553 28.190 57.750 100.000
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.00 43.43 66.00 58.32 80.80 100.00
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.00 64.98 77.67 74.82 88.84 100.00
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.00 75.50 83.89 81.13 90.25 110.00