#Practica en R ## Simulación del lanzamiento de dos dados para estimar la suma más probable

dado1=round(runif(10,1,6),0)
dado1
##  [1] 1 1 1 1 2 5 2 3 5 2
dado2=round(runif(10,1,6),0)
dado2
##  [1] 4 4 5 2 4 3 2 2 2 4
suma=dado1+dado2
suma
##  [1] 5 5 6 3 6 8 4 5 7 6
table(suma)
## suma
## 3 4 5 6 7 8 
## 1 1 3 3 1 1
barplot(table(suma))

#simulación 100 replicas

dado1=round(runif(100,1,6),0)
dado1
##   [1] 4 1 3 4 5 5 3 5 6 4 4 1 6 1 2 3 5 2 5 4 2 5 1 2 3 1 3 1 5 4 4 6 4 3 5 4 3
##  [38] 5 1 2 4 5 3 4 1 2 4 2 6 3 2 5 4 1 4 3 4 3 4 6 4 5 3 5 1 4 4 2 5 1 2 1 1 3
##  [75] 2 4 3 4 2 5 6 5 5 2 3 5 2 3 4 5 5 4 3 2 5 4 3 2 2 5
dado2=round(runif(100,1,6),0)
dado2
##   [1] 3 2 5 3 1 6 4 5 3 6 4 3 4 2 6 5 5 3 2 5 5 4 3 1 3 6 4 4 5 3 3 5 6 6 2 1 4
##  [38] 6 2 1 5 1 2 3 5 3 4 5 6 3 6 2 3 5 2 2 2 2 2 5 3 2 2 4 4 6 2 5 4 4 4 3 5 2
##  [75] 2 5 1 2 2 2 4 2 1 5 2 2 4 4 5 5 4 1 5 6 6 2 6 1 4 6
suma=dado1+dado2
suma
##   [1]  7  3  8  7  6 11  7 10  9 10  8  4 10  3  8  8 10  5  7  9  7  9  4  3  6
##  [26]  7  7  5 10  7  7 11 10  9  7  5  7 11  3  3  9  6  5  7  6  5  8  7 12  6
##  [51]  8  7  7  6  6  5  6  5  6 11  7  7  5  9  5 10  6  7  9  5  6  4  6  5  4
##  [76]  9  4  6  4  7 10  7  6  7  5  7  6  7  9 10  9  5  8  8 11  6  9  3  6 11
table(suma)
## suma
##  3  4  5  6  7  8  9 10 11 12 
##  6  6 13 17 23  8 11  9  6  1
barplot(table(suma))

#simulación 1000 replicas

dado1=round(runif(1000,1,6),0)
dado1
##    [1] 2 1 5 4 3 4 3 5 5 3 3 6 5 3 4 5 2 5 6 3 5 6 2 5 6 6 1 3 3 5 6 3 2 2 5 4 6
##   [38] 6 5 6 6 3 3 3 5 3 3 4 2 3 4 4 2 1 3 2 6 4 3 4 4 2 3 1 5 4 4 4 5 5 6 5 6 3
##   [75] 4 5 3 1 3 3 1 4 3 1 4 3 4 3 3 3 5 4 5 2 3 3 4 5 3 1 5 4 2 5 1 5 3 3 5 1 6
##  [112] 4 1 1 5 1 4 2 3 2 3 3 1 3 3 5 4 2 2 4 4 2 3 3 5 3 3 5 2 3 3 5 5 3 3 2 1 3
##  [149] 5 2 6 2 3 5 6 6 3 4 4 5 1 4 2 1 3 3 4 5 5 1 6 3 3 4 5 4 5 2 1 4 2 6 5 1 2
##  [186] 4 4 4 1 1 2 1 2 5 4 6 2 3 2 4 4 5 5 4 4 4 5 6 4 1 5 5 5 2 2 2 6 4 2 6 4 5
##  [223] 3 4 6 2 5 2 6 1 5 2 3 6 5 2 4 6 5 4 4 1 3 1 2 2 3 5 3 5 2 4 3 3 5 1 3 6 2
##  [260] 2 4 3 3 2 2 2 6 6 4 5 5 1 4 2 4 3 2 4 4 4 3 4 5 3 4 5 5 3 5 3 6 2 4 3 4 2
##  [297] 3 5 1 3 3 5 4 5 5 3 3 5 1 6 3 5 4 1 5 6 4 2 3 4 4 4 5 6 5 1 3 4 5 1 4 4 4
##  [334] 2 4 5 5 2 5 6 2 2 4 3 3 1 3 2 3 1 2 4 5 4 2 5 4 1 3 3 1 1 3 4 4 1 2 6 1 5
##  [371] 5 5 4 5 6 2 4 1 2 3 5 3 6 5 2 5 4 5 5 4 6 1 6 6 3 2 1 4 5 6 4 5 2 6 6 4 2
##  [408] 2 4 2 3 1 4 4 4 1 2 3 5 5 4 4 4 5 3 3 2 2 6 2 4 2 5 5 1 4 6 3 5 3 6 5 2 3
##  [445] 4 5 5 3 5 3 2 3 4 2 4 5 5 5 4 3 5 1 4 2 3 5 2 2 2 3 1 3 5 3 5 5 5 5 2 2 3
##  [482] 4 2 4 4 6 3 4 2 2 4 5 3 2 6 1 2 2 1 2 1 4 4 5 3 4 5 2 2 4 6 3 4 1 2 3 2 3
##  [519] 2 1 3 5 4 3 1 1 5 5 6 4 2 4 3 3 5 5 2 3 2 3 2 1 5 2 3 3 6 6 3 4 3 4 2 6 1
##  [556] 5 5 4 4 4 2 1 1 3 2 5 6 1 1 3 6 4 6 1 3 4 1 6 2 3 5 3 6 3 5 4 5 6 5 4 2 4
##  [593] 2 1 1 4 6 3 3 5 1 4 5 3 1 3 2 3 5 1 3 4 2 2 5 4 3 5 6 4 4 3 1 5 2 2 4 2 4
##  [630] 5 5 3 4 6 6 4 2 1 3 5 1 4 5 5 4 5 5 2 3 2 2 3 1 4 6 5 4 3 2 5 1 3 3 5 1 2
##  [667] 3 5 3 5 1 2 3 4 5 5 1 5 6 5 4 2 2 5 3 1 5 3 3 5 3 6 4 3 4 1 2 6 1 6 3 4 2
##  [704] 4 5 3 2 4 2 2 3 4 4 3 5 5 2 3 2 4 1 2 5 1 3 3 5 5 6 3 2 5 4 4 5 3 5 4 5 5
##  [741] 5 3 4 5 3 5 4 3 2 2 2 6 5 5 5 1 3 5 5 1 4 2 4 1 6 3 5 2 3 2 5 3 3 5 6 4 1
##  [778] 2 2 6 5 4 3 3 2 6 2 6 6 2 5 2 4 3 4 6 2 4 5 2 4 4 5 3 5 5 3 3 2 4 6 2 3 3
##  [815] 4 5 3 3 3 1 3 3 4 3 6 3 5 6 2 4 2 5 5 4 3 2 3 3 3 3 3 5 2 2 4 3 3 3 3 3 3
##  [852] 4 3 4 5 3 2 3 6 4 5 2 6 4 2 4 4 5 1 3 2 4 5 4 1 4 1 2 5 1 4 5 2 4 3 6 2 2
##  [889] 3 2 4 2 3 3 5 5 3 2 3 4 2 2 5 5 5 4 4 3 2 3 6 5 4 3 2 5 3 2 3 3 4 5 2 4 5
##  [926] 2 5 3 2 3 1 5 2 4 4 1 2 3 4 4 4 1 3 2 1 3 2 3 5 1 2 3 2 2 5 3 2 1 5 4 5 3
##  [963] 5 4 2 4 6 5 3 1 4 3 2 5 3 6 5 4 4 3 6 6 5 2 3 4 3 2 5 3 2 3 4 6 6 5 5 4 5
## [1000] 3
dado2=round(runif(1000,1,6),0)
dado2
##    [1] 3 3 4 2 4 6 3 5 3 2 4 5 4 2 4 3 5 3 6 5 5 4 5 4 4 1 2 5 3 3 2 2 3 2 5 1 5
##   [38] 5 4 5 3 4 4 2 2 4 4 2 3 5 4 5 5 3 4 5 2 4 2 3 2 4 2 5 4 4 5 5 2 6 1 4 3 4
##   [75] 5 2 3 5 3 6 3 2 3 2 3 4 5 4 5 3 5 5 4 6 5 6 1 1 2 3 1 5 3 4 1 2 5 2 5 5 2
##  [112] 4 2 2 6 1 5 6 3 4 3 4 3 3 4 5 4 6 5 4 5 3 5 4 5 5 2 4 3 2 2 6 5 2 5 2 5 2
##  [149] 4 1 2 5 5 3 2 3 1 1 2 4 4 6 5 5 1 5 5 5 3 6 3 1 6 2 4 5 5 6 3 2 3 5 2 2 3
##  [186] 5 3 5 6 2 1 3 5 5 4 6 3 5 4 5 3 3 3 3 3 3 5 5 1 3 3 3 1 6 3 4 5 5 6 5 3 3
##  [223] 5 4 4 4 4 4 5 5 5 5 4 4 5 6 2 2 5 3 6 4 5 4 6 3 2 2 2 2 5 4 6 4 2 4 1 3 1
##  [260] 4 4 5 3 6 4 6 6 6 2 2 1 5 5 6 6 2 5 5 2 2 1 5 2 3 2 6 2 4 3 2 3 2 3 5 4 2
##  [297] 4 3 5 2 2 6 3 2 6 6 1 1 4 6 2 3 4 5 3 2 1 5 2 1 4 1 2 4 6 5 1 3 3 5 2 3 4
##  [334] 4 6 2 1 5 4 6 3 4 2 4 4 5 1 3 1 6 5 4 4 4 4 4 5 4 5 1 5 4 6 6 6 2 4 2 5 2
##  [371] 3 3 5 3 2 2 4 2 5 2 2 3 3 5 5 2 4 5 4 3 3 5 6 5 3 3 2 2 2 3 2 1 3 4 3 4 4
##  [408] 2 3 4 2 4 5 2 1 3 2 4 5 5 5 3 3 2 6 2 4 2 4 2 3 2 5 2 4 2 4 5 5 3 3 5 6 6
##  [445] 6 3 3 5 5 3 5 2 3 6 3 5 6 4 4 3 5 5 2 3 5 5 2 1 5 5 5 6 3 5 4 5 4 1 4 3 4
##  [482] 2 6 2 5 4 1 2 3 3 6 2 6 2 2 3 6 5 2 6 6 3 5 5 2 6 5 5 3 6 6 4 5 4 2 4 4 5
##  [519] 5 3 1 1 1 1 6 3 5 4 5 3 3 1 5 4 2 4 4 2 4 3 5 6 2 5 4 3 4 4 5 4 3 3 5 3 6
##  [556] 1 3 5 5 5 4 3 2 4 1 3 4 3 4 4 1 6 2 4 3 5 4 5 4 5 3 6 3 2 4 3 4 2 3 4 2 6
##  [593] 2 5 3 4 2 5 6 4 4 4 6 2 2 2 4 5 4 2 2 3 5 5 4 4 3 2 5 1 2 3 4 6 6 5 5 5 1
##  [630] 1 3 2 3 3 5 5 5 1 5 4 6 1 3 6 1 3 5 2 2 3 3 4 6 3 2 2 4 6 4 2 3 6 4 2 2 5
##  [667] 1 2 5 5 4 3 3 1 6 5 2 5 2 5 6 5 5 3 2 3 3 2 1 3 3 3 2 6 4 3 3 5 3 3 5 5 5
##  [704] 5 5 5 2 3 4 4 5 1 5 6 2 2 2 4 4 6 4 5 4 3 3 3 5 4 2 5 3 4 5 3 1 5 4 2 4 4
##  [741] 2 5 4 6 2 3 4 2 3 4 4 5 3 6 2 2 3 3 2 5 3 4 3 4 6 5 2 2 5 4 6 6 5 2 5 3 5
##  [778] 5 5 3 3 4 2 2 1 1 5 2 4 2 4 2 6 3 2 3 4 4 4 2 4 5 3 4 2 3 3 6 6 2 2 5 3 3
##  [815] 4 6 3 5 3 4 3 3 2 3 3 1 2 2 4 5 5 3 3 1 5 5 1 6 4 4 6 3 6 2 6 4 1 5 5 2 4
##  [852] 5 6 1 5 6 2 5 5 3 5 2 4 4 5 3 6 3 3 3 1 3 3 5 2 6 5 5 4 5 1 3 3 3 1 6 4 5
##  [889] 2 4 3 6 4 4 5 5 2 3 5 5 2 1 4 4 4 3 1 5 2 2 6 3 4 3 4 3 2 3 3 3 3 2 5 4 4
##  [926] 2 2 4 4 3 4 3 5 3 4 3 2 4 6 5 4 4 1 5 4 5 5 5 5 4 2 4 4 4 4 1 2 5 6 5 3 4
##  [963] 3 2 6 4 6 1 5 3 2 4 6 4 2 2 5 4 4 1 1 4 4 2 5 2 3 2 2 5 4 5 1 4 3 4 3 5 3
## [1000] 5
suma=dado1+dado2
suma
##    [1]  5  4  9  6  7 10  6 10  8  5  7 11  9  5  8  8  7  8 12  8 10 10  7  9
##   [25] 10  7  3  8  6  8  8  5  5  4 10  5 11 11  9 11  9  7  7  5  7  7  7  6
##   [49]  5  8  8  9  7  4  7  7  8  8  5  7  6  6  5  6  9  8  9  9  7 11  7  9
##   [73]  9  7  9  7  6  6  6  9  4  6  6  3  7  7  9  7  8  6 10  9  9  8  8  9
##   [97]  5  6  5  4  6  9  5  9  2  7  8  5 10  6  8  8  3  3 11  2  9  8  6  6
##  [121]  6  7  4  6  7 10  8  8  7  8  9  5  8  7 10  8  5  9  5  5  5 11 10  5
##  [145]  8  4  6  5  9  3  8  7  8  8  8  9  4  5  6  9  5 10  7  6  4  8  9 10
##  [169]  8  7  9  4  9  6  9  9 10  8  4  6  5 11  7  3  5  9  7  9  7  3  3  4
##  [193]  7 10  8 12  5  8  6  9  7  8  8  7  7  7 10 11  5  4  8  8  6  8  5  6
##  [217] 11  9  8 11  7  8  8  8 10  6  9  6 11  6 10  7  7 10 10  8  6  8 10  7
##  [241] 10  5  8  5  8  5  5  7  5  7  7  8  9  7  7  5  4  9  3  6  8  8  6  8
##  [265]  6  8 12 12  6  7  6  6  9  8 10  5  7  9  6  6  4  9  7  6  6 11  7  7
##  [289]  8  5  9  4  7  8  8  4  7  8  6  5  5 11  7  7 11  9  4  6  5 12  5  8
##  [313]  8  6  8  8  5  7  5  5  8  5  7 10 11  6  4  7  8  6  6  7  8  6 10  7
##  [337]  6  7  9 12  5  6  6  7  7  6  4  5  4  7  7  8  9  8  6  9  9  5  8  4
##  [361]  6  5  9 10 10  3  6  8  6  7  8  8  9  8  8  4  8  3  7  5  7  6  9 10
##  [385]  7  7  8 10  9  7  9  6 12 11  6  5  3  6  7  9  6  6  5 10  9  8  6  4
##  [409]  7  6  5  5  9  6  5  4  4  7 10 10  9  7  7  7  9  5  6  4 10  4  7  4
##  [433] 10  7  5  6 10  8 10  6  9 10  8  9 10  8  8  8 10  6  7  5  7  8  7 10
##  [457] 11  9  8  6 10  6  6  5  8 10  4  3  7  8  6  9  8  8  9 10  9  6  6  5
##  [481]  7  6  8  6  9 10  4  6  5  5 10  7  9  4  8  4  8  7  3  8  7  7  9 10
##  [505]  5 10 10  7  5 10 12  7  9  5  4  7  6  8  7  4  4  6  5  4  7  4 10  9
##  [529] 11  7  5  5  8  7  7  9  6  5  6  6  7  7  7  7  7  6 10 10  8  8  6  7
##  [553]  7  9  7  6  8  9  9  9  6  4  3  7  3  8 10  4  5  7  7 10  8  5  6  9
##  [577]  5 11  6  8  8  9  9  5  9  7  9  8  8  8  4 10  4  6  4  8  8  8  9  9
##  [601]  5  8 11  5  3  5  6  8  9  3  5  7  7  7  9  8  6  7 11  5  6  6  5 11
##  [625]  8  7  9  7  5  6  8  5  7  9 11  9  7  2  8  9  7  5  8 11  5  8 10  4
##  [649]  5  5  5  7  7  7  8  7  8  9  6  7  4  9  7  7  3  7  4  7  8 10  5  5
##  [673]  6  5 11 10  3 10  8 10 10  7  7  8  5  4  8  5  4  8  6  9  6  9  8  4
##  [697]  5 11  4  9  8  9  7  9 10  8  4  7  6  6  8  5  9  9  7  7  4  7  6 10
##  [721]  5  7  9  4  6  6 10  9  8  8  5  9  9  7  6  8  9  6  9  9  7  8  8 11
##  [745]  5  8  8  5  5  6  6 11  8 11  7  3  6  8  7  6  7  6  7  5 12  8  7  4
##  [769]  8  6 11  9  8  7 11  7  6  7  7  9  8  8  5  5  3  7  7  8 10  4  9  4
##  [793] 10  6  6  9  6  8  9  4  8  9  8  7  7  8  6  9  8  6  8  7  6  6  8 11
##  [817]  6  8  6  5  6  6  6  6  9  4  7  8  6  9  7  8  8  5  8  7  4  9  7  7
##  [841]  9  8  8  4 10  7  4  8  8  5  7  9  9  5 10  9  4  8 11  7 10  4 10  8
##  [865]  7  7 10  8  4  6  3  7  8  9  3 10  6  7  9  6  5  8  5  7  4 12  6  7
##  [889]  5  6  7  8  7  7 10 10  5  5  8  9  4  3  9  9  9  7  5  8  4  5 12  8
##  [913]  8  6  6  8  5  5  6  6  7  7  7  8  9  4  7  7  6  6  5  8  7  7  8  4
##  [937]  4  7 10  9  8  5  4  7  5  8  7  8 10  5  4  7  6  6  9  4  4  6 11  9
##  [961]  8  7  8  6  8  8 12  6  8  4  6  7  8  9  5  8 10  8  8  4  7 10  9  4
##  [985]  8  6  6  4  7  8  6  8  5 10  9  9  8  9  8  8
table(suma)
## suma
##   2   3   4   5   6   7   8   9  10  11  12 
##   3  25  78 119 148 183 188 130  79  35  12
barplot(table(suma))

dado1=round(runif(1000000,1,6),0)
#dado1
dado2=round(runif(1000000,1,6),0)
#dado2
suma=dado1+dado2
#suma
table(suma)
## suma
##      2      3      4      5      6      7      8      9     10     11     12 
##  10111  39910  80101 120242 160416 179230 159462 120280  79967  40297   9984
barplot(table(suma))

prop.table(table(suma))
## suma
##        2        3        4        5        6        7        8        9 
## 0.010111 0.039910 0.080101 0.120242 0.160416 0.179230 0.159462 0.120280 
##       10       11       12 
## 0.079967 0.040297 0.009984
barplot(prop.table(table(suma)))