n <- 10000; dati <- NULL; for (i in 1:n) dati[i]=sum(floor(runif(10)*2))
hist(dati, seq(-0.5, 10.5, 1), probability=TRUE, col="yellow")
abline(h=axTicks(2), col="blue",lty=3)
summary(dati)
#  Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
# 0.000   4.000   5.000   4.983   6.000  10.000
dev.new()
boxplot(dati,horizontal=TRUE, col="yellow", range=0)
# la visualizzazione sul boxplot del 5° e del 95° percentile:
points(quantile(dati,0.05),1,pch=20); points(quantile(dati,0.95),1,pch=20)