[Spiegazione di come è stato ottenuto il grafico.  Con binom abbiamo memorizzato la distribuzione binomiale. Quindi basta che facciamo congiungere, per i da 0 a N=6, i punti (i,binom(i)). In plot metto n=N+1, type="p" in modi che tracci 7 punti; con pch=19 specifico la forma dei punti. Nel successivo plot type="h" sta per tipo "istogramma", in inglese "histogram", e comanda il tracciamento di linee vericali. Non devo dividere la media per il numero dei dati in quanto in binom ho già le densità di frequenza]

N <- 6; p <- 1/8
binom <- function(m) choose(N,m)*p^m*(1-p)^(N-m)
plot(binom, 0,N, n=N+1, type="p", pch=19)
abline(h=axTicks(2), col="blue",lty=3)
plot(binom,0,N,n=N+1,add=TRUE,type="h",col="blue",lwd=3)
plot(binom,0,N,n=N+1,add=TRUE,lty=3,col="blue",lwd=2)
x <- seq(0,N,1); media <- sum(x*binom(x)); media
  0.75
var <- sum((x-media)^2*binom(x)); var
  0.65625