Simulación del lanzamiento de dos dados para estimar la suma más
probable.
Para 10 réplicas
dado1 <- sample(1:6,10,replace = T)
dado1
## [1] 5 5 4 1 1 4 2 1 1 1
dado2 <- sample(1:6,10,replace = T)
dado2
## [1] 4 2 5 5 6 4 6 4 1 3
suma <- dado1+dado2
suma
## [1] 9 7 9 6 7 8 8 5 2 4
tabla_suma <- table(suma)
tabla_suma
## suma
## 2 4 5 6 7 8 9
## 1 1 1 1 2 2 2
barplot(tabla_suma)
Para 100 réplicas
dado1 <- sample(1:6,100,replace = T)
dado1
## [1] 2 1 2 1 2 3 1 5 3 3 5 4 2 2 1 3 3 6 5 3 5 2 3 3 1 2 3 4 5 1 6 4 4 6 5 5 3
## [38] 6 5 4 1 6 5 3 5 1 4 3 6 2 2 1 2 1 5 3 6 6 1 3 2 5 1 2 2 4 2 5 3 1 5 1 6 6
## [75] 2 4 6 6 2 6 2 6 5 1 4 3 1 2 4 3 6 2 5 2 2 5 6 3 3 1
dado2 <- sample(1:6,100,replace = T)
dado2
## [1] 3 4 6 1 2 6 3 6 4 2 1 1 5 5 5 1 6 4 3 1 6 2 3 5 6 2 3 1 6 1 6 1 2 2 3 4 5
## [38] 6 5 3 4 4 3 2 4 6 4 2 3 2 4 4 5 2 1 2 3 3 1 2 1 6 5 2 1 5 6 1 2 5 6 4 4 5
## [75] 6 1 1 1 5 2 2 2 2 4 1 6 3 5 6 1 2 1 2 4 3 2 4 5 5 4
suma <- dado1+dado2
suma
## [1] 5 5 8 2 4 9 4 11 7 5 6 5 7 7 6 4 9 10 8 4 11 4 6 8 7
## [26] 4 6 5 11 2 12 5 6 8 8 9 8 12 10 7 5 10 8 5 9 7 8 5 9 4
## [51] 6 5 7 3 6 5 9 9 2 5 3 11 6 4 3 9 8 6 5 6 11 5 10 11 8
## [76] 5 7 7 7 8 4 8 7 5 5 9 4 7 10 4 8 3 7 6 5 7 10 8 8 5
tabla_suma <- table(suma)
tabla_suma
## suma
## 2 3 4 5 6 7 8 9 10 11 12
## 3 4 11 19 11 14 15 9 6 6 2
barplot(tabla_suma)
Para 1000 réplicas
dado1 <- sample(1:6,1000,replace = T)
dado2 <- sample(1:6,1000,replace = T)
suma <- dado1+dado2
tabla_suma <- table(suma)
tabla_suma
## suma
## 2 3 4 5 6 7 8 9 10 11 12
## 28 66 86 103 151 169 143 100 63 64 27
prop.table(tabla_suma)
## suma
## 2 3 4 5 6 7 8 9 10 11 12
## 0.028 0.066 0.086 0.103 0.151 0.169 0.143 0.100 0.063 0.064 0.027
barplot(tabla_suma)
Para 1000000 réplicas
dado1 <- sample(1:6,1000000,replace = T)
dado2 <- sample(1:6,1000000,replace = T)
suma <- dado1+dado2
tabla_suma <- table(suma)
tabla_suma
## suma
## 2 3 4 5 6 7 8 9 10 11 12
## 27971 55614 83564 110834 138692 167480 138648 110407 83411 55543 27836
prop.table(tabla_suma)
## suma
## 2 3 4 5 6 7 8 9
## 0.027971 0.055614 0.083564 0.110834 0.138692 0.167480 0.138648 0.110407
## 10 11 12
## 0.083411 0.055543 0.027836
barplot(tabla_suma)
