# Le probabilità con cui escono i vari numeri su un particolare
# dado truccato.
p <- c(0.095,0.19,0.125,0.125,0.18,0.285)
# Verifica se la somma è il 100%
sum(p)
# ottengo 1: OK
V <- function(x) ifelse( x < p[1],1,ifelse( x < sum(p[1:2]),2,
          ifelse( x < sum(p[1:3]),3,ifelse( x < sum(p[1:4]),4,
          ifelse( x < sum(p[1:5]),5,6)))))
V(runif(1))
# Questa è un'uscita del nostro dado.
# Nel modo seguente ottengo un istogramma sperimentale che approssima
# quello teorico (faccio 100 mila prove)
U <- V(runif(1e5)); hist(U, seq(0.5,6.5,1), probability=TRUE)
abline(h=axTicks(2),lty=3)