Elaborar una simulación de extracción de monedas de un frasco para encontrar el espacio muestral y con ello determinar probabilidades.
El caso tiene sustento del Ejercicio 4.5. del libro de Mendengall (2006) de cuatro monedas. Un frasco contiene cuatro monedas de diferente denominación, se seleccionan al azar tres monedas del frasco y se determinan probabilidades. (Mendenhall, Beaver, and Beaver 2006)
El caso tiene describe la simulación de extracción de monedas de un frasco. En el frasco hay monedas de $1, $2, $5, $10, $20 y $50 pesos.
Se pide a una persona que no sabe de las características de las monedas (dimensiones, colores, texturas, entre otras cosas), digamos un niño o un extranjero seleccione (extraer) aleatoriamente una o tres monedas en tres experimentos diferentes posibles.
En un primer experimento, , extrae una moneda y se requiere conocer cuál es la probabilidad de que esa moneda sea de alguna denominación en particular.
En un segundo experimento con el estado inicial de todas las monedas en el frasco, es decir, deposita la moneda y revuelve el frasco y hace otro experimento, extrae tres monedas, y se debe encontrar la probabilidad de que la suma de las tres monedas sea igual o mayor que $50 pesos.
En un tercer experimento se extraen tres monedas de nuevo y se requiere conocer la probabilidad de que la suma de las tres monedas sea igual o mayor a $70 pesos.
El proceso del caso sería de la siguiente manera:
En un frasco contiene seis monedas de las donominaciones siguientes:
1 moneda de $1 peso
1 moneda de $2 pesos
1 moneda de $5 pesos
1 moneda de $10 pesos
1 moneda $20 pesos
1 moneda de $50 pesos
Haga una lista de los eventos simples en S
El primer experimento es simular extraer una moneda del frasco y determine probabilidades
Realizar permutaciones para grupos de tres monedas y determinar el nuevo espacio muestral.
El segundo experimento consiste en simular extraer tres monedas de las seis posibles.
Determinar la probabilidad para que la suma de las tres monedas sea igual o superior a $50 pesos
En un tercer experimento consiste en simular de nuevo extraer tres monedas de las seis posibles.
Determinar la probabilidad para que la sume de las tres monedas sea igual o superior a $70 pesos
Se solicita que la interpretación responda a las siguientes preguntas:
¿Cuál es la probabilidad de que la selección de una moneda contenga la moneda de $50 pesos?
¿Cuál es la probabilidad de que la suma total de extraer tres monedas sea igual o mayor a $50 pesos?
¿Cuál es la probabilidad de que la suma total de extraer tres monedas sea igual o superior a $70 pesos?
Se define el concepto de probabilidad y sus requerimiento básicos:
La probabilidad asignada a cada resultado experimental debe estar entre 0 y 1. Si denota con Ei el i-ésimo resultado experimental y con P(Ei) su probabilidad,
entonces exprese este requerimiento como
0≤P(Ei)≤1
La suma de las probabilidades de los resultados experimentales debe ser igual a 1 o al 100 en términos porcentuales. [@anderson_estadistica_2008].
Para resultados experimentales n se suman todas las probabilidades:
P(E1)+P(E2)+P(E3)+…P(En)=1
[@anderson_estadistica_2008].
La manera más sencilla de encontrar probabilidades es mediante el cálculo de frecuencia relativa. Si un experimento se realiza n veces, entonces la frecuencia relativa de un suceso particular es la cantidad de sucesos de ese tipo entre el valor de todos los posibles.[@mendenhall_introduccion_2006].
frecuencia.relativa=frecuencian
entonces
Prob(Ei)=frecuencia.relativa
Se cargan las librerías necesarias
Se define la lista de los eventos simples en S.
S=1,2,5,10,20,50
Simular la extracción de una moneda y determinar la probabilidad de que la selección contenga la moneda de $ 50 pesos
Simular la extracción de tres monedas haciendo permutaciones y determinar la probabilidad de que la suma de las tres monedas sea igual o mayor que $50 pesos.
Simular la extracción de tres monedas y con el espacio muestral de las permutaciones, determinar la probabilidad de que la suma total de extraer tres monedas sea igual o superior a $70 pesos?
Se carga la librería dplyr para filtrar registros con la función filter(), seleccionar con la función select() y generar nuevas variables con la función mutate() de un conjunto de datos o data.frame.
Se carga la librería gtools para generar permutaciones y combinaciones
Se carga la librería knitr para imágenes y tablas
library(dplyr) library(gtools)
library(knitr)Se define la lista de los eventos simples en un espacio muestral identificado como S. Se utiliza la variable tipo vector llamado S.
S={1,2,5,10,20,50
S <- c(1, 2, 5, 10, 20, 50)¿Cuál es la probabilidad de que la selección contenga la moneda de 50 pesos?
length() determina el número de elemenos (cuantos) de un vector S Los [] determina la posición del valor comparado which() identifica cuales elementos corresponden a una expresión S==50 La variable cuantas es sinónimo de frecuencia con respeto al total de elementos.
n <- length(S)
cuantas <- length(S[which(S==50)])
prob <- cuantas / n
paste("La probabilidad de que la selección contenga la moneda de $50 pesos es: ",cuantas ,"/",n, " o sea: ", prob, " o ", round(prob * 100,2), "%")¿Cuál es la probabilidad de que la suma total de extraer tres monedas sea igual o superior a 50 pesos?
Valores iniciales de n y r n es el total elementos del espacio muetral S, de los eentos simples. El valor de n=6 r es de cuantos en cuantos o grupos para las permutaciones, de r=3, dado que se pide extraer tres monedas
n <- length(S)
r <- 3 Hacer las permutaciones posibles y determinar el nuevo espacio muestral. La variable que se usa para el nuevo espacio muestal es S.permuta
S.permuta <- permutations(n,r,S)
S.permutaComprobar el número de puntos muesrales de permutaciones cuando n=6 y r=3 conforme a la fórmula de permutaciones vista en el caso 7. S=Pr(nr)=n!/(n−r)!
factorial(n) / factorial(n-r)Transformar el resultado de las permutaciones a un tipo de estructura data.frame para su mejor trato.
S.permuta <- data.frame(S.permuta) # Se convierte a data.framePoner nombres a las columnas como moneda1, moneda2, y moneda3 o más práctico names(S.permuta) <- c(‘m1,’ ‘m2,’ ‘m3’).
Con ello se debe observar como cambian los nombres de las columnas.
Con la función str(S.permuta) se describe la estructura del conjunto de datos S.permuta
names(S.permuta) <- c('m1', 'm2', 'm3')
str(S.permuta)Con las permutaciones sacar la suma de cada renglón
Se construye la columna suma mediante la función mutate()
S.permuta <- mutate(S.permuta, suma = m1 + m2 + m3) # Genera nueva columna
S.permutaN <- nrow(S.permuta)
paste("Número de casos posibles de permutaciones", N)Al extraer tres monedas, se determinan el conjunto de opciones posibles para que suma sea mayor o igual a 50 con la función filter().
Se determinan ¿cuáles y cuántas?
cuales <- filter(S.permuta, suma >= 50)
cuales
cuantas <- nrow(cuales)Se responde a las pregunta de ¿cuál es la probabilidad de que la suma de extraer tres monedas sea igual o superior a $50 pesos.
paste("En ", cuantas, " ocasiones la suma es mayor o igual 50 de ", N, " posibles ")paste("La probabilidad de que la suma sea mayor o igual a $50 es del ", cuantas , "/", N, " o sea", cuantas / N, " que significa el:", round(cuantas / N * 100, 2), " % ")¿Cuál es la probabilidad de que al extraer tres monedas, la suma sea mayor o igual a $70 pesos
Se utiliza el mismo espacio muestral de grupos de tres dado que la extracción de monedas es tres de seis posibles.
Se determinan las nuevas cantidades ¿cuáles monedas pudieran ser y cuántas la suma es igual o superior a $70?
cuales <- filter(S.permuta, suma >= 70)
head(cuales); tail(cuales)cuantas <- nrow(cuales)Se responde a las pregunta de ¿cuál es la probabilidad de que la suma de extraer tres monedas sea igual o superior a $70 pesos
paste("En ", cuantas, " ocasiones la suma es mayor o igual 70 de ", N, " posibles ")paste("La probabilidad de que la suma sea mayor o igual a $70 es del ", cuantas , "/", N, " o sea", cuantas / N, " que significa el:", round(cuantas / N * 100, 2), " % ")
Si en el frasco se tienen dos moneda de $100 pesos, dos de a $50, tres a $20, cuatro de a $10, cinco de $5, diez de a $2 y veinte monedas de $1. ¿cómo cambian las probabilidades?
S={100,100,50,50,20,20,20,10,10,10,10,5,5,5,5,5,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}
S <- c(100,100,50,50,20,20,20,10,10,10,10,5,5,5,5,5,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)Se debe hacer permutaciones de r=3
n <- length(S)
r <- 3 ¿Existe la posibilidad de que la suma de extraer tres de todo el nuevo espacio muestral habiendo hecho permutaciones de tres sea menor que $5?.
¿Cuál es la probabilidad de que al extraer tres monedas de todo el espacio muestral habiendo hecho permutaciones, la suma esté entre $40 y $60 pesos?.
¿Cuál es la probabilidad de que al extraer tres monedas de todo el espacio muestral habiendo hecho permutaciones, la suma esté entre $150 y $200 pesos?.
Se utilizará una función llamada rep() para generar números repetidos, con ello, se simula la cantidad de monedas iniciales que existen en el frasco.
S <- c(rep(100,2), rep(50, 2), rep(20,3), rep(10,4), rep(5,5), rep(2, 10), rep(1,20))
SSe identifica en n la cantidad de elementos del espacio muestral.
Se identifica r como la cantidad de cinco monedas.
n <- length(S)
nr <- 3
rDe acuerdo a la fórmula en donde n=46y r = 3 se determina cuántas posibilidades se deben generar en el experimento de extraer cinco monedas de un total de 46. La variable cuantas tiene el resultado
cuantas <- factorial(n) / factorial(n-r)
cuantasSe generan permutaciones en grupos de tres dado que se quiere experimentar con extraer tres monedas
S.permuta <- permutations(n=n, r=r, S, set = FALSE)
head(S.permuta)tail(S.permuta)Transformar el resultado de las permutaciones a un tipo de estructura data.frame para su mejor trato.
S.permuta <- data.frame(S.permuta) # Se convierte a data.framePoner nombres a las columnas como moneda1, moneda2, y moneda3 o más práctico names(S.permuta) <- c(‘m1,’ ‘m2,’ ‘m3’).
Con ello se debe observar como cambian los nombres de las columnas.
Con la función str(S.permuta) se describe la estructura del conjunto de datos S.permuta
names(S.permuta) <- c('m1', 'm2', 'm3')
str(S.permuta)
Con las permutaciones sacar la suma de cada renglón
S.permuta <- mutate(S.permuta, suma = m1 + m2 + m3) # Genera nueva columna
head(S.permuta)N <- nrow(S.permuta)
paste("Número de casos posibles de permutaciones", N)
Al extraer tres monedas, se determinan el conjunto de opciones posibles para que suma sea mayor o igual a las preguntas del experimento 4 con la función filter().
Se determinan ¿cuáles y cuántas?
¿Existe la posibilidad de que la suma de extraer tres de todo el nuevo espacio muestral habiendo hecho permutaciones de tres sea menor que $5?.
cuales <- filter(S.permuta, suma < 5)
str(cuales)cuantas <- nrow(cuales)¿Cuál es la probabilidad de que al extraer tres monedas de todo el espacio muestral habiendo hecho permutaciones, la suma esté entre $40 y $60 pesos?.
cuales <- filter(S.permuta, suma >= 40 & suma <=60)
str(cuales)cuantas <- nrow(cuales)La probabilidad
prob <- cuantas / N
prob <- round(prob * 100, 2)
probpaste ("La probabilidad es del ", prob, "%")
¿Cuál es la probabilidad de que al extraer tres monedas de todo el espacio muestral habiendo hecho permutaciones, la suma esté entre $150 y $200 pesos?.
cuales <- filter(S.permuta, suma >= 150 & suma <=200)
str(cuales)cuantas <- nrow(cuales)
Para encontrar la probabilidad se utiliza nuevamente la frecuencia relativa
prob <- cuantas / N
prob <- round(prob * 100, 2)
probpaste ("La probabilidad es del ", prob, "%")
Las preguntas del 1 al 3 se refieren a los experimentos uno al tres respectivamente.
¿Cuál es la probabilidad de que la selección de una moneda contenga la moneda de $50 pesos?. Utilizar experimento 1. 16.67%
¿Cuál es la probabilidad de que la suma total de extraer tres monedas sea igual o mayor a $50 pesos?. Utilizar experimento 2. 50%
¿Cuál es la probabilidad de que la suma total de extraer tres monedas sea igual o superior a $70 pesos?. Utilizar experimento 3.
20%
Realizar el cuarto experimento y responder:
Las siguientes preguntas son del experimento 4.
¿Existe la posibilidad de que la suma de extraer tres monedas de todo el nuevo espacio muestral habiendo hecho permutaciones de tres sea menor que $5?. Si
¿Cuál es la probabilidad de que al extraer cinco monedas de todo el espacio muestral habiendo hecho permutaciones, la suma esté entre $40 y $60 pesos?. 8.74%
¿Cuál es la probabilidad de que al extraer cinco monedas de todo el espacio muestral habiendo hecho permutaciones, la suma esté entre $150 y $200 pesos?. 1.12 %
Mendenhall, William, Robert J. Beaver, and Barbara M. Beaver. 2006. Introducción a La Probabilidad y Estadística. 13a Edición.