Objetivo.

Resolver cuestiones de casos de probabilidad en casos mediante la identificación de variables aleatorias, funciones de probabilidad, funciones acumuladas, media, varianza y desviación estándar de distribuciones de variables discretas; visualización gráfica relacionada con variables discretas.

Descripción.

Identificar casos relacionados con variables discretas para elaborar mediante programación R y markdown las variables discretas, las funciones de probabilidad de cada variable, la función acumulada, su visualización gráfica para su correcta implementación.

Se incluye en el caso, media, varianza y desviación estándar de distribuciones de variables discretas.

Los casos son identificados de la literatura relacionada con variables aleatorias discretas. Se deben elaborar tres ejercicios en este caso 15 encontrados en la literatura, se pueden apoyar de los mismos ejercicios del caso 14.

1. Cargar librerías.

library(ggplot2)
library(stringr)
library(stringi)
library(gtools)
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(knitr)

2. Identificar ejercicios de la literatura.

Para cada ejercicio realizar lo siguiente:

Describir y definir el contexo del ejercicio, poner la referencia (enlace) del mismo.

Elaborar su tabla de probabilidad que contenga los valores de la variable aleatoria, la función de probabilidad y su función acumulada, la gráfica de barra de los valores de las variables aleatoria y la gráfica lineal de la función acumulada.

En algunos ejercicios trae consigo el espacio muestral de preferencia incorporarlo en el documento

Identificar a diferencia del caso 14 o ¿Cuál es el valor de la media de la distribución? o ¿Cuál es el valor de la varianza de la distribución? o ¿Cuál es el valor de desviación de la distribución?.

Ejercicio 1.

Se venden 5000 billetes para una rifa a 1 euro cada uno. Existe un único premio de cierta cantidad, calcular los valores de las variables aleatorias y sus probabilidades para 0 para no gana y 1 para si gana cuando un comprador adquiere tres billetes. (Hero, n.d.)

Tabla de probabilidad.

discretas <- c(0,1)   
n <- 5000
casos <- c(4997,50)
probabilidades <- casos / n

acumulada <- cumsum(probabilidades)   

tabla <- data.frame(x=discretas, 
                    casos = casos,
                    f.prob.x = probabilidades,
                    F.acum.x = acumulada,
                    x.f.prob.x= (discretas*probabilidades))
kable(tabla)
x casos f.prob.x F.acum.x x.f.prob.x
0 4997 0.9994 0.9994 0.00
1 50 0.0100 1.0094 0.01

Valor esperado.

Se determina el valor esperado de acuerdo a la fórmula: μ=∑xP(x)

VE es el valor esperado.

VE <- sum(tabla$x * tabla$f.prob.x)

VE
## [1] 0.01

Varianza.

tabla <- cbind(tabla, 'VE' = VE, 'x-VE.cuad.f.prob.x' = (tabla$x - VE)^2 * tabla$f.prob.x)

kable(tabla, caption = "Tabla de probabilidad con valor esperado y columnas para varianza")
Tabla de probabilidad con valor esperado y columnas para varianza
x casos f.prob.x F.acum.x x.f.prob.x VE x-VE.cuad.f.prob.x
0 4997 0.9994 0.9994 0.00 0.01 0.0000999
1 50 0.0100 1.0094 0.01 0.01 0.0098010

var= varianza de la distribución.

varianza <- sum((tabla$x - VE)^2 * tabla$f.prob.x)
varianza
## [1] 0.00990094
desv.std <- sqrt(varianza)
desv.std
## [1] 0.09950347

Ejercicio 2.

Las ventas de automóviles de una empresa

Durante los últimos 300 días de operación, los datos de ventas muestran que hubo

Tabla de probabilidad o contingencia

discretas <- c(0:5)   
n <- 300
casos <- c(54, 117, 72, 42, 12, 3)
probabilidades <- casos / n

acumulada <- cumsum(probabilidades)   

tabla <- data.frame(x=discretas, 
                    casos = casos,
                    f.prob.x = probabilidades,
                    F.acum.x = acumulada,
                    x.f.prob.x= (discretas*probabilidades))
kable(tabla)
x casos f.prob.x F.acum.x x.f.prob.x
0 54 0.18 0.18 0.00
1 117 0.39 0.57 0.39
2 72 0.24 0.81 0.48
3 42 0.14 0.95 0.42
4 12 0.04 0.99 0.16
5 3 0.01 1.00 0.05

Valor esperado.

Se determina el valor esperado de acuerdo a la fórmula: μ=∑xP(x)

VE es el valor esperado

VE <- sum(tabla$x * tabla$f.prob.x)

VE
## [1] 1.5

Varianza.

tabla <- cbind(tabla, 'VE' = VE, 'x-VE.cuad.f.prob.x' = (tabla$x - VE)^2 * tabla$f.prob.x)

kable(tabla, caption = "Tabla de probabilidad con valor esperado y columnas para varianza")
Tabla de probabilidad con valor esperado y columnas para varianza
x casos f.prob.x F.acum.x x.f.prob.x VE x-VE.cuad.f.prob.x
0 54 0.18 0.18 0.00 1.5 0.4050
1 117 0.39 0.57 0.39 1.5 0.0975
2 72 0.24 0.81 0.48 1.5 0.0600
3 42 0.14 0.95 0.42 1.5 0.3150
4 12 0.04 0.99 0.16 1.5 0.2500
5 3 0.01 1.00 0.05 1.5 0.1225
varianza <- sum((tabla$x - VE)^2 * tabla$f.prob.x)
varianza
## [1] 1.25
desv.std <- sqrt(varianza)
desv.std
## [1] 1.118034

Ejercicio 3.

En Estados Unidos un porcentaje de los niños de cuarto grado no pueden leer un libro adecuado a su edad. La tabla siguiente muestra, de acuerdo con las edades de entre 6 y 14 años, el número de niños que tienen problemas de lectura. La mayoría de estos niños tienen problemas de lectura que debieron ser detectados y corregidos antes del tercer grado.(Anderson et al., 2008)

Tabla de probabilidad o Contingencia

discretas <- c(6:14)   

casos <- c(37369, 87436, 160840,239719,286719,306533,310787,302604,289168)
probabilidades <- casos / n

acumulada <- cumsum(probabilidades)   

tabla <- data.frame(x=discretas, 
                    casos = casos,
                    f.prob.x = probabilidades,
                    F.acum.x = acumulada,
                    x.f.prob.x= (discretas*probabilidades))
kable(tabla)
x casos f.prob.x F.acum.x x.f.prob.x
6 37369 124.5633 124.5633 747.380
7 87436 291.4533 416.0167 2040.173
8 160840 536.1333 952.1500 4289.067
9 239719 799.0633 1751.2133 7191.570
10 286719 955.7300 2706.9433 9557.300
11 306533 1021.7767 3728.7200 11239.543
12 310787 1035.9567 4764.6767 12431.480
13 302604 1008.6800 5773.3567 13112.840
14 289168 963.8933 6737.2500 13494.507

Valor esperado.

Se determina el valor esperado de acuerdo a la fórmula: μ=∑xP(x)

VE es el valor esperado.

VE <- sum(tabla$x * tabla$f.prob.x)

VE
## [1] 74103.86

Varianza.

tabla <- cbind(tabla, 'VE' = VE, 'x-VE.cuad.f.prob.x' = (tabla$x - VE)^2 * tabla$f.prob.x)

kable(tabla, caption = "Tabla de probabilidad con valor esperado y columnas para varianza")
Tabla de probabilidad con valor esperado y columnas para varianza
x casos f.prob.x F.acum.x x.f.prob.x VE x-VE.cuad.f.prob.x
6 37369 124.5633 124.5633 747.380 74103.86 6.839141e+11
7 87436 291.4533 416.0167 2040.173 74103.86 1.600179e+12
8 160840 536.1333 952.1500 4289.067 74103.86 2.943477e+12
9 239719 799.0633 1751.2133 7191.570 74103.86 4.386896e+12
10 286719 955.7300 2706.9433 9557.300 74103.86 5.246862e+12
11 306533 1021.7767 3728.7200 11239.543 74103.86 5.609300e+12
12 310787 1035.9567 4764.6767 12431.480 74103.86 5.686992e+12
13 302604 1008.6800 5773.3567 13112.840 74103.86 5.537104e+12
14 289168 963.8933 6737.2500 13494.507 74103.86 5.291107e+12
varianza <- sum((tabla$x - VE)^2 * tabla$f.prob.x)
varianza
## [1] 3.698583e+13
desv.std <- sqrt(varianza)
desv.std
## [1] 6081598

Interpretación del caso.

En este caso número 15 se habla sobre varianza, media y la desviación de distribución de variables discretas como en el caso anterior y se utilizaron los mismos ejercicios, a continuación haré una interpretación sobre cada problema:

Ejercicio 1: nos referimos a los 5000 billetes que se venden para una rifa de 1 euro cada uno y una persona compra tres billetes. En este ejercicio vemos que la variable aleatoria encontramos que es la variable con el nombre de casos ya que son las posibles probabilidades que son 4997 y 50, la dicha variable aleatoria toma esos dos resultados, el espacio muestral lo tomamos como casos ya que casos son todos los elementos posibles que son 4997 y 50, solo hay 2 elementos en el espacio muestral, el numero de casos es pequeño ya que hay 2 casos para este ejercicio y las probabilidades mas altas de cada variable aleatoria es para el 4997 es: 99.94% y para el 50 es del 1%, ya después encontramos el valor esperado fue del 0.01 ya que se multiplica los x casos con las probabilidades, la varianza salió 0.0000999 y en el otro caso salió 0.0098010 y finalmente la varianza de la distribución que fue 0.00990094 y la desviación estándar que fue del 0.09959347.

Ejercicio 2: trata de la empresa de ventas que toman los datos de venta de los últimos 300 días y en 54 días no se vendió ningún auto, en 117 días nomas 1 auto, 72 días 2 autos, 42 días 3 autos, 12 días 4 autos y en 3 días 5 automóviles. La variable aleatoria es la de casos que tiene los datos siguientes: 54,117,72,12,3, la variable aleatoria tiene 6 casos, el espacio muestral son todos los elementos y la que tiene todos los elementos es la de casos, hay 6 elementos en el espacio muestral, las probabilidades mas altas fueron: 54 son: 18%, en el 117 es: 39%, en el 72 es: 24%, en el 42 es: 14%, en el 12 es: .04% y en el 3 es del .01%, encontramos el valor esperado que es del 1.5, la varianza de cada caso es de: 1 es: 0.4050, 2 es: 0.0975, 3 es: 0.0600, 4 es: 0.3150, 5 es: 0.2500 y en el 6 es: 0.1225, la varianza de distribución es el 1.25, y la desviación estándar es del 1.118.

Ejercicio 3: es de porcentaje de los niños de cuarto grado que no pueden leer en estados unidos, adecuados a su edad con las edades de entre 6 y 14 años. La variable aleatoria tiene las probabilidades finitas o posibles que es la de casos, dicha variable aleatoria tiene 9 casos, el espacio muestral es en donde se encuentran todos los elementos y al igual que la variable aleatoria tiene los mismos, se encuentran los 9 elementos del espacio muestral, las probabilidades mas altas para cada elemento es del: 1 son: .018%, en el 2 es: .043%, en el 3 es: .079%, en el 4 es: 11%, en el 5 es: 14% y en el caso 6 es del 15%, en el 7 es del: 15%, en el 8 es del 14% y en el 9 es del 14%. El valor esperado es de 10.99, la varianza de cada caso es de: 1 es: 0.4620, 2 es: 0.6918, 3 es: 0.7157, 4 es: 0.4740, 5 es: 0.1416, 6 es: 0.0000001, 7 es: 0.1540, 8 es: 0.5993, 9 es: 1.2883, la varianza de distribución es 4.5271 y la desviación estándar es de: 2.1276.