# 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