n <- 1e6; dadi <- trunc(runif(n, min=1, max=7))+trunc(runif(n, min=1, max=7)) n1=0; for(i in 1:n) if(dadi[i]>7) n1=n1+1; print(n1/n) #[1] 0.416937 # Per avere l'andamento delle uscite (10, 100, ... 10^6 uscite): for(k in 1:6) {tot=10^k; n1=0; for(i in 1:tot) if(dadi[i]>7) n1=n1+1; print(n1/tot)} #[1] 0.4 #[1] 0.45 #[1] 0.399 #[1] 0.4103 #[1] 0.4149 #[1] 0.416937 # # Per avere una (facile) rappresentazione grafica della distribuzione delle # uscite carico: source("http://macosa.dima.unige.it/r.R") # Quindi: noClassi=1; Istogramma(dadi, 1.5,12.5, 1) #[i lati della griglia distano il 1 %] # Frequenze e frequenze percentuali: #27801, 55540, 82972, 111333, 138900, 166620, 139240, 111049, 83185, 55674, 27686 #2.78,5.55,8.3,11.13,13.89,16.66,13.92,11.1,8.32,5.57,2.77 # Per altre statistiche usa il comando altrestat() altrestat() # Min. 1st Qu. Median Mean 3rd Qu. Max. # 2 5 7 7 9 12 # I pallini marroni sono 5° e 95° percentile # Il pallino rosso è la media