library(ggplot2)
library(stringr) # String
library(stringi) # String
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)
options(scipen = 999) # Notación normal
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:
discretas <- 1:4
casos <- c(30, 40, 50, 60)
probabilidades <- c(0.4, 0.2, 0.1, 0.3)
acumulada <- cumsum(probabilidades)
p.30 <- 0.4
p.40 <- 0.2
p.50 <- 0.1
p.60 <- 0.3
tabla.funcion <- data.frame(x=discretas,
casos = casos,
p.x = probabilidades,
f.x = acumulada)
tabla.funcion
## x casos p.x f.x
## 1 1 30 0.4 0.4
## 2 2 40 0.2 0.6
## 3 3 50 0.1 0.7
## 4 4 60 0.3 1.0
VE <- sum(tabla.funcion$x * tabla.funcion$p.x)
VE
## [1] 2.3
tabla <- cbind(tabla.funcion, 'VE' = VE, 'x-VE.cuad.f.prob.x' = (tabla.funcion$x - VE)^2 * tabla.funcion$p.x)
kable(tabla, caption = "Tabla de probabilidad con valor esperado y columnas para varianza")
| x | casos | p.x | f.x | VE | x-VE.cuad.f.prob.x |
|---|---|---|---|---|---|
| 1 | 30 | 0.4 | 0.4 | 2.3 | 0.676 |
| 2 | 40 | 0.2 | 0.6 | 2.3 | 0.018 |
| 3 | 50 | 0.1 | 0.7 | 2.3 | 0.049 |
| 4 | 60 | 0.3 | 1.0 | 2.3 | 0.867 |
varianza <- sum((tabla.funcion$x - VE)^2 * tabla.funcion$p.x)
varianza
## [1] 1.61
desv.std <- sqrt(varianza)
desv.std
## [1] 1.268858
discretas <- 1:3
casos <- c(0,1,2) # Numero de caras posibles
probabilidades <- c(0.25, 0.5, 0.25)
acumulada <- cumsum(probabilidades)
p.0 <- 0.25
p.1 <- 0.5
p.2 <- 0.25
tabla.funcion <- data.frame(x=discretas,
casos = casos,
p.x = probabilidades,
f.x = acumulada)
tabla.funcion
## x casos p.x f.x
## 1 1 0 0.25 0.25
## 2 2 1 0.50 0.75
## 3 3 2 0.25 1.00
VE <- sum(tabla.funcion$x * tabla.funcion$p.x)
tabla <- cbind(tabla.funcion, 'VE' = VE, 'x-VE.cuad.f.prob.x' = (tabla.funcion$x - VE)^2 * tabla.funcion$p.x)
kable(tabla, caption = "Tabla de probabilidad con valor esperado y columnas para varianza")
| x | casos | p.x | f.x | VE | x-VE.cuad.f.prob.x |
|---|---|---|---|---|---|
| 1 | 0 | 0.25 | 0.25 | 2 | 0.25 |
| 2 | 1 | 0.50 | 0.75 | 2 | 0.00 |
| 3 | 2 | 0.25 | 1.00 | 2 | 0.25 |
varianza <- sum((tabla.funcion$x - VE)^2 * tabla.funcion$p.x)
varianza
## [1] 0.5
desv.std <- sqrt(varianza)
desv.std
## [1] 0.7071068
casos <- 1:6 # Numero de puntos posibles #en este caso tambien cuentan como la variable discreta
probabilidades <- c((1/6), (1/6), (1/6), (1/6), (1/6), (1/6))
acumulada <- cumsum(probabilidades)
p.1 <- 1/6
p.2 <- 1/6
p.3 <- 1/6
p.4 <- 1/6
p.5 <- 1/6
p.6 <- 1/6
tabla.funcion <- data.frame(x=casos,
p.x = probabilidades,
f.x = acumulada)
tabla.funcion
## x p.x f.x
## 1 1 0.1666667 0.1666667
## 2 2 0.1666667 0.3333333
## 3 3 0.1666667 0.5000000
## 4 4 0.1666667 0.6666667
## 5 5 0.1666667 0.8333333
## 6 6 0.1666667 1.0000000
VE <- sum(tabla.funcion$x * tabla.funcion$p.x)
tabla <- cbind(tabla.funcion, 'VE' = VE, 'x-VE.cuad.f.prob.x' = (tabla.funcion$x - VE)^2 * tabla.funcion$p.x)
kable(tabla, caption = "Tabla de probabilidad con valor esperado y columnas para varianza")
| x | p.x | f.x | VE | x-VE.cuad.f.prob.x |
|---|---|---|---|---|
| 1 | 0.1666667 | 0.1666667 | 3.5 | 1.0416667 |
| 2 | 0.1666667 | 0.3333333 | 3.5 | 0.3750000 |
| 3 | 0.1666667 | 0.5000000 | 3.5 | 0.0416667 |
| 4 | 0.1666667 | 0.6666667 | 3.5 | 0.0416667 |
| 5 | 0.1666667 | 0.8333333 | 3.5 | 0.3750000 |
| 6 | 0.1666667 | 1.0000000 | 3.5 | 1.0416667 |
varianza <- sum((tabla.funcion$x - VE)^2 * tabla.funcion$p.x)
varianza
## [1] 2.916667
desv.std <- sqrt(varianza)
desv.std
## [1] 1.707825