library(tidyverse)
library(dplyr)
library(gtools)
library(knitr)S <- c(1,2,5,10,20)n <- length(S)
cuantas <- length(S[which(S==10)])
prob <- cuantas / n
paste("La probabilidad de que la selección contenga la moneda de 10 pesos es: 1 /",n, " o sea: ", prob, " o ", prob * 100, "%")## [1] "La probabilidad de que la selección contenga la moneda de 10 pesos es: 1 / 5 o sea: 0.2 o 20 %"n <- length(S)
r <- 4permuta <- permutations(n,r,S)
permuta## [,1] [,2] [,3] [,4]
## [1,] 1 2 5 10
## [2,] 1 2 5 20
## [3,] 1 2 10 5
## [4,] 1 2 10 20
## [5,] 1 2 20 5
## [6,] 1 2 20 10
## [7,] 1 5 2 10
## [8,] 1 5 2 20
## [9,] 1 5 10 2
## [10,] 1 5 10 20
## [11,] 1 5 20 2
## [12,] 1 5 20 10
## [13,] 1 10 2 5
## [14,] 1 10 2 20
## [15,] 1 10 5 2
## [16,] 1 10 5 20
## [17,] 1 10 20 2
## [18,] 1 10 20 5
## [19,] 1 20 2 5
## [20,] 1 20 2 10
## [21,] 1 20 5 2
## [22,] 1 20 5 10
## [23,] 1 20 10 2
## [24,] 1 20 10 5
## [25,] 2 1 5 10
## [26,] 2 1 5 20
## [27,] 2 1 10 5
## [28,] 2 1 10 20
## [29,] 2 1 20 5
## [30,] 2 1 20 10
## [31,] 2 5 1 10
## [32,] 2 5 1 20
## [33,] 2 5 10 1
## [34,] 2 5 10 20
## [35,] 2 5 20 1
## [36,] 2 5 20 10
## [37,] 2 10 1 5
## [38,] 2 10 1 20
## [39,] 2 10 5 1
## [40,] 2 10 5 20
## [41,] 2 10 20 1
## [42,] 2 10 20 5
## [43,] 2 20 1 5
## [44,] 2 20 1 10
## [45,] 2 20 5 1
## [46,] 2 20 5 10
## [47,] 2 20 10 1
## [48,] 2 20 10 5
## [49,] 5 1 2 10
## [50,] 5 1 2 20
## [51,] 5 1 10 2
## [52,] 5 1 10 20
## [53,] 5 1 20 2
## [54,] 5 1 20 10
## [55,] 5 2 1 10
## [56,] 5 2 1 20
## [57,] 5 2 10 1
## [58,] 5 2 10 20
## [59,] 5 2 20 1
## [60,] 5 2 20 10
## [61,] 5 10 1 2
## [62,] 5 10 1 20
## [63,] 5 10 2 1
## [64,] 5 10 2 20
## [65,] 5 10 20 1
## [66,] 5 10 20 2
## [67,] 5 20 1 2
## [68,] 5 20 1 10
## [69,] 5 20 2 1
## [70,] 5 20 2 10
## [71,] 5 20 10 1
## [72,] 5 20 10 2
## [73,] 10 1 2 5
## [74,] 10 1 2 20
## [75,] 10 1 5 2
## [76,] 10 1 5 20
## [77,] 10 1 20 2
## [78,] 10 1 20 5
## [79,] 10 2 1 5
## [80,] 10 2 1 20
## [81,] 10 2 5 1
## [82,] 10 2 5 20
## [83,] 10 2 20 1
## [84,] 10 2 20 5
## [85,] 10 5 1 2
## [86,] 10 5 1 20
## [87,] 10 5 2 1
## [88,] 10 5 2 20
## [89,] 10 5 20 1
## [90,] 10 5 20 2
## [91,] 10 20 1 2
## [92,] 10 20 1 5
## [93,] 10 20 2 1
## [94,] 10 20 2 5
## [95,] 10 20 5 1
## [96,] 10 20 5 2
## [97,] 20 1 2 5
## [98,] 20 1 2 10
## [99,] 20 1 5 2
## [100,] 20 1 5 10
## [101,] 20 1 10 2
## [102,] 20 1 10 5
## [103,] 20 2 1 5
## [104,] 20 2 1 10
## [105,] 20 2 5 1
## [106,] 20 2 5 10
## [107,] 20 2 10 1
## [108,] 20 2 10 5
## [109,] 20 5 1 2
## [110,] 20 5 1 10
## [111,] 20 5 2 1
## [112,] 20 5 2 10
## [113,] 20 5 10 1
## [114,] 20 5 10 2
## [115,] 20 10 1 2
## [116,] 20 10 1 5
## [117,] 20 10 2 1
## [118,] 20 10 2 5
## [119,] 20 10 5 1
## [120,] 20 10 5 2permuta <- data.frame(permuta) # Se convierte a data.frame
names(permuta) <- c('muestraN1', 'muestraN2', 'muestraN3', 'muestraN4')
permuta <- mutate(permuta, suma = muestraN1 + muestraN2 + muestraN3 + muestraN4) # Genera nueva columna
permuta## muestraN1 muestraN2 muestraN3 muestraN4 suma
## 1 1 2 5 10 18
## 2 1 2 5 20 28
## 3 1 2 10 5 18
## 4 1 2 10 20 33
## 5 1 2 20 5 28
## 6 1 2 20 10 33
## 7 1 5 2 10 18
## 8 1 5 2 20 28
## 9 1 5 10 2 18
## 10 1 5 10 20 36
## 11 1 5 20 2 28
## 12 1 5 20 10 36
## 13 1 10 2 5 18
## 14 1 10 2 20 33
## 15 1 10 5 2 18
## 16 1 10 5 20 36
## 17 1 10 20 2 33
## 18 1 10 20 5 36
## 19 1 20 2 5 28
## 20 1 20 2 10 33
## 21 1 20 5 2 28
## 22 1 20 5 10 36
## 23 1 20 10 2 33
## 24 1 20 10 5 36
## 25 2 1 5 10 18
## 26 2 1 5 20 28
## 27 2 1 10 5 18
## 28 2 1 10 20 33
## 29 2 1 20 5 28
## 30 2 1 20 10 33
## 31 2 5 1 10 18
## 32 2 5 1 20 28
## 33 2 5 10 1 18
## 34 2 5 10 20 37
## 35 2 5 20 1 28
## 36 2 5 20 10 37
## 37 2 10 1 5 18
## 38 2 10 1 20 33
## 39 2 10 5 1 18
## 40 2 10 5 20 37
## 41 2 10 20 1 33
## 42 2 10 20 5 37
## 43 2 20 1 5 28
## 44 2 20 1 10 33
## 45 2 20 5 1 28
## 46 2 20 5 10 37
## 47 2 20 10 1 33
## 48 2 20 10 5 37
## 49 5 1 2 10 18
## 50 5 1 2 20 28
## 51 5 1 10 2 18
## 52 5 1 10 20 36
## 53 5 1 20 2 28
## 54 5 1 20 10 36
## 55 5 2 1 10 18
## 56 5 2 1 20 28
## 57 5 2 10 1 18
## 58 5 2 10 20 37
## 59 5 2 20 1 28
## 60 5 2 20 10 37
## 61 5 10 1 2 18
## 62 5 10 1 20 36
## 63 5 10 2 1 18
## 64 5 10 2 20 37
## 65 5 10 20 1 36
## 66 5 10 20 2 37
## 67 5 20 1 2 28
## 68 5 20 1 10 36
## 69 5 20 2 1 28
## 70 5 20 2 10 37
## 71 5 20 10 1 36
## 72 5 20 10 2 37
## 73 10 1 2 5 18
## 74 10 1 2 20 33
## 75 10 1 5 2 18
## 76 10 1 5 20 36
## 77 10 1 20 2 33
## 78 10 1 20 5 36
## 79 10 2 1 5 18
## 80 10 2 1 20 33
## 81 10 2 5 1 18
## 82 10 2 5 20 37
## 83 10 2 20 1 33
## 84 10 2 20 5 37
## 85 10 5 1 2 18
## 86 10 5 1 20 36
## 87 10 5 2 1 18
## 88 10 5 2 20 37
## 89 10 5 20 1 36
## 90 10 5 20 2 37
## 91 10 20 1 2 33
## 92 10 20 1 5 36
## 93 10 20 2 1 33
## 94 10 20 2 5 37
## 95 10 20 5 1 36
## 96 10 20 5 2 37
## 97 20 1 2 5 28
## 98 20 1 2 10 33
## 99 20 1 5 2 28
## 100 20 1 5 10 36
## 101 20 1 10 2 33
## 102 20 1 10 5 36
## 103 20 2 1 5 28
## 104 20 2 1 10 33
## 105 20 2 5 1 28
## 106 20 2 5 10 37
## 107 20 2 10 1 33
## 108 20 2 10 5 37
## 109 20 5 1 2 28
## 110 20 5 1 10 36
## 111 20 5 2 1 28
## 112 20 5 2 10 37
## 113 20 5 10 1 36
## 114 20 5 10 2 37
## 115 20 10 1 2 33
## 116 20 10 1 5 36
## 117 20 10 2 1 33
## 118 20 10 2 5 37
## 119 20 10 5 1 36
## 120 20 10 5 2 37paste("Número de casos posibles de permutaciones", nrow(permuta))## [1] "Número de casos posibles de permutaciones 120"cuales <- filter(permuta, suma >= 25)
cuales## muestraN1 muestraN2 muestraN3 muestraN4 suma
## 1 1 2 5 20 28
## 2 1 2 10 20 33
## 3 1 2 20 5 28
## 4 1 2 20 10 33
## 5 1 5 2 20 28
## 6 1 5 10 20 36
## 7 1 5 20 2 28
## 8 1 5 20 10 36
## 9 1 10 2 20 33
## 10 1 10 5 20 36
## 11 1 10 20 2 33
## 12 1 10 20 5 36
## 13 1 20 2 5 28
## 14 1 20 2 10 33
## 15 1 20 5 2 28
## 16 1 20 5 10 36
## 17 1 20 10 2 33
## 18 1 20 10 5 36
## 19 2 1 5 20 28
## 20 2 1 10 20 33
## 21 2 1 20 5 28
## 22 2 1 20 10 33
## 23 2 5 1 20 28
## 24 2 5 10 20 37
## 25 2 5 20 1 28
## 26 2 5 20 10 37
## 27 2 10 1 20 33
## 28 2 10 5 20 37
## 29 2 10 20 1 33
## 30 2 10 20 5 37
## 31 2 20 1 5 28
## 32 2 20 1 10 33
## 33 2 20 5 1 28
## 34 2 20 5 10 37
## 35 2 20 10 1 33
## 36 2 20 10 5 37
## 37 5 1 2 20 28
## 38 5 1 10 20 36
## 39 5 1 20 2 28
## 40 5 1 20 10 36
## 41 5 2 1 20 28
## 42 5 2 10 20 37
## 43 5 2 20 1 28
## 44 5 2 20 10 37
## 45 5 10 1 20 36
## 46 5 10 2 20 37
## 47 5 10 20 1 36
## 48 5 10 20 2 37
## 49 5 20 1 2 28
## 50 5 20 1 10 36
## 51 5 20 2 1 28
## 52 5 20 2 10 37
## 53 5 20 10 1 36
## 54 5 20 10 2 37
## 55 10 1 2 20 33
## 56 10 1 5 20 36
## 57 10 1 20 2 33
## 58 10 1 20 5 36
## 59 10 2 1 20 33
## 60 10 2 5 20 37
## 61 10 2 20 1 33
## 62 10 2 20 5 37
## 63 10 5 1 20 36
## 64 10 5 2 20 37
## 65 10 5 20 1 36
## 66 10 5 20 2 37
## 67 10 20 1 2 33
## 68 10 20 1 5 36
## 69 10 20 2 1 33
## 70 10 20 2 5 37
## 71 10 20 5 1 36
## 72 10 20 5 2 37
## 73 20 1 2 5 28
## 74 20 1 2 10 33
## 75 20 1 5 2 28
## 76 20 1 5 10 36
## 77 20 1 10 2 33
## 78 20 1 10 5 36
## 79 20 2 1 5 28
## 80 20 2 1 10 33
## 81 20 2 5 1 28
## 82 20 2 5 10 37
## 83 20 2 10 1 33
## 84 20 2 10 5 37
## 85 20 5 1 2 28
## 86 20 5 1 10 36
## 87 20 5 2 1 28
## 88 20 5 2 10 37
## 89 20 5 10 1 36
## 90 20 5 10 2 37
## 91 20 10 1 2 33
## 92 20 10 1 5 36
## 93 20 10 2 1 33
## 94 20 10 2 5 37
## 95 20 10 5 1 36
## 96 20 10 5 2 37cuantas <- nrow(cuales)
paste("En ", cuantas, " ocasiones la suma es 25 o mas de ", nrow(permuta), " posibles ")## [1] "En 96 ocasiones la suma es 25 o mas de 120 posibles "paste("La probabilidad de que la suma sea 25 o mas es del ", cuantas , "/", nrow(permuta), " o sea", cuantas / nrow(permuta), " que significa el:", cuantas / nrow(permuta) * 100, " % ")## [1] "La probabilidad de que la suma sea 25 o mas es del 96 / 120 o sea 0.8 que significa el: 80 % "