Ejercicio 1: Lanza un dado justo de 6 caras. ¿Cuál es la probabilidad de obtener un número par?
total_caras <- 6
caras_pares <- 3
probabilidad_par <- caras_pares / total_caras
probabilidad_par## [1] 0.5Ejercicio 3: Una urna contiene 10 bolas numeradas del 1 al 10. ¿Cuál es la probabilidad de extraer una bola con un número menor o igual a 4?
total_bolas <- 10
#bolas con número menor o igual a 4
bolas_menor_igual_4 <- 4
# Probabilidad
probabilidad_bola_menor_igual_4 <- bolas_menor_igual_4 / total_bolas
probabilidad_bola_menor_igual_4## [1] 0.4Ejercicio Actuarial 1: Un asegurador ofrece una póliza con un premio que varía uniformemente entre 200 y 500 unidades monetarias. Calcula la probabilidad de que un cliente pague un premio superior a 400 unidades.
#Limites
limite_inferior <- 200
limite_superior <- 500
# Premios arriba de 400
premio_superior_400 <- 400
# Probabilidad de que el premio sea superior a 400
probabilidad_premio_superior_400 <- (limite_superior - premio_superior_400) / (limite_superior - limite_inferior)
probabilidad_premio_superior_400## [1] 0.3333333Ejercicio 1: Se lanza una moneda justa. ¿Cuál es la probabilidad de obtener cara?
# Probabilidad de obtener cara en una moneda justa
p_cara <- 0.5
p_cara## [1] 0.5Ejercicio 2: Un test de diagnóstico tiene una tasa de éxito del 85%. ¿Cuál es la probabilidad de que el test resulte exitoso?
# Probabilidad de éxito en el test de diagnóstico
p_exito <- 0.85
p_exito## [1] 0.85Ejercicio Actuarial 1: En una póliza de seguro, la probabilidad de un evento catastrófico es de 0.1. Si el evento es modelado por una distribución de Bernoulli, ¿cuál es la probabilidad de que el evento ocurra en una póliza dada?
# Probabilidad de que ocurra un evento catastrófico
p_evento_catastrofico <- 0.1
p_evento_catastrofico## [1] 0.1Ejercicio Actuarial 2: Supón que un asegurador ofrece un seguro con una probabilidad de siniestro de 0.3. Calcula la probabilidad de que no haya siniestro en un año dado.
# Probabilidad de siniestro
p_siniestro <- 0.3
p_no_siniestro <- 1 - p_siniestro
p_no_siniestro## [1] 0.7Para resolver estos ejercicios utilizando la distribución geométrica en R, utilizaremos la función dgeom() para calcular la probabilidad de que el primer éxito ocurra en un intento específico. La distribución geométrica modela el número de fracasos antes de obtener el primer éxito en una serie de ensayos de Bernoulli.
La fórmula para la probabilidad geométrica es:
\[ P(X = k) = (1 - p)^{k - 1} \cdot p \]
Donde:
\(p\) es la probabilidad de éxito,
\(k\) es el número de intentos (incluyendo el primer éxito).
Ejercicio 1: ¿Cuál es la probabilidad de obtener el primer seis en el tercer lanzamiento de un dado justo?
# Probabilidad de éxito (obtener un 6)
p_seis <- 1/6
# Probabilidad de que el primer 6 ocurra en el tercer intento
prob_primer_seis_tercer_intento <- dgeom(2, prob = p_seis)
prob_primer_seis_tercer_intento## [1] 0.1157407Ejercicio 2: Si el tiempo hasta que se presenta un defecto en una producción sigue una distribución geométrica con p=0.2, ¿cuál es la probabilidad de que el primer defecto ocurra en el cuarto intento?
# Probabilidad de éxito (ocurrencia de un defecto)
p_defecto <- 0.2
# Probabilidad de que el primer defecto ocurra en el cuarto intento
prob_primer_defecto_cuarto_intento <- dgeom(3, prob = p_defecto)
prob_primer_defecto_cuarto_intento## [1] 0.1024Ejercicio 3: En un proceso de pruebas, la probabilidad de éxito en cada prueba es de 0.4. ¿Cuál es la probabilidad de necesitar exactamente 5 pruebas para tener el primer éxito?
# Probabilidad de éxito en una prueba
p_exito_prueba <- 0.4
# Probabilidad de que el primer éxito ocurra en el quinto intento
prob_exito_quinto_intento <- dgeom(4, prob = p_exito_prueba)
prob_exito_quinto_intento## [1] 0.05184Ejercicio Actuarial 1: Un asegurador evalúa el número de años hasta que se presenta el primer siniestro en una póliza de seguro. Si el tiempo hasta el primer siniestro sigue una distribución geométrica con p=0.2, calcula la probabilidad de que el primer siniestro ocurra en el tercer año.
# Probabilidad de siniestro por año
p_siniestro <- 0.2
# Probabilidad de que el primer siniestro ocurra en el tercer año
prob_siniestro_tercer_ano <- dgeom(2, prob = p_siniestro)
prob_siniestro_tercer_ano## [1] 0.128Ejercicio Actuarial 2: En el análisis de riesgo, el tiempo hasta que un cliente realiza una reclamación se modela con una distribución geométrica con p=0.4 ¿Cuál es la esperanza y la varianza del tiempo esperado hasta la primera reclamación?
Para una distribución geométrica con parámetro \(p\), la esperanza y la varianza están dadas por:
\[ E(X) = \frac{1}{p} \]
\[ Var(X) = \frac{1 - p}{p^2} \]
Donde:
\(p\) es la probabilidad de éxito.
# Probabilidad de una reclamación
p_reclamacion <- 0.4
# Esperanza (tiempo hasta la primera reclamación)
esperanza_reclamacion <- 1 / p_reclamacion
varianza_reclamacion <- (1 - p_reclamacion) / (p_reclamacion^2)
esperanza_reclamacion## [1] 2.5varianza_reclamacion## [1] 3.75Para resolver estos ejercicios utilizando la distribución de Poisson en R, podemos emplear la función dpois() para calcular la probabilidad de observar un número específico de eventos en un intervalo dado. La distribución de Poisson es utilizada para modelar el número de eventos en un intervalo de tiempo dado, con una tasa promedio de ocurrencia 𝜆.
La fórmula para la probabilidad de Poisson es:
\[ P(X = k) = \frac{\lambda^k e^{-\lambda}}{k!} \]
Donde:
\(\lambda\) es la tasa de ocurrencia por intervalo,
\(k\) es el número de eventos.
Ejercicio 1: ¿Cuál es la probabilidad de que ocurran exactamente 3 eventos en un intervalo de tiempo si la tasa de eventos es 2 por intervalo?
# Parámetro de Poisson (tasa de eventos por intervalo)
lambda_1 <- 2
# Probabilidad de que ocurran exactamente 3 eventos
prob_3_eventos <- dpois(3, lambda = lambda_1)
prob_3_eventos## [1] 0.180447Ejercicio 2: En una línea de producción, el número de defectos en una hora sigue una distribución de Poisson con una tasa de λ=5 ¿Cuál es la probabilidad de observar exactamente 4 defectos en una hora?
# Parámetro de Poisson (tasa de defectos por hora)
lambda_2 <- 5
# Probabilidad de observar exactamente 4 defectos
prob_4_defectos <- dpois(4, lambda = lambda_2)
prob_4_defectos## [1] 0.1754674Ejercicio 3: Si el número de llamadas a un centro de atención telefónica sigue una distribución Poisson con λ=10 llamadas por hora, ¿cuál es la probabilidad de recibir menos de 8 llamadas en una hora?
# Parámetro de Poisson (tasa de llamadas por hora)
lambda <- 10
# Probabilidad de recibir menos de 8 llamadas (es decir, 7 o menos)
prob_menos_8_llamadas <- ppois(7, lambda = lambda)
prob_menos_8_llamadas## [1] 0.2202206Ejercicio Actuarial 1: El número de reclamaciones de seguros en una región sigue una distribución de Poisson con una tasa de λ=12 reclamaciones por año. ¿Cuál es la probabilidad de recibir exactamente 10 reclamaciones en un año?
# Parámetro de Poisson (tasa de reclamaciones por año)
lambda_3 <- 12
# Probabilidad de recibir exactamente 10 reclamaciones
prob_10_reclamaciones <- dpois(10, lambda = lambda_3)
prob_10_reclamaciones## [1] 0.1048373Ejercicio Actuarial 2: En un análisis de riesgo de una póliza, se estima que el número de siniestros en un año sigue una distribución de Poisson con λ=8. Calcula la probabilidad de que haya entre 5 y 7 siniestros en un año.
# Parámetro de Poisson (tasa de siniestros por año)
lambda_4 <- 8
# Probabilidad de que haya entre 5 y 7 siniestros
prob_5_7_siniestros <- sum(dpois(5:7, lambda = lambda_4))
prob_5_7_siniestros## [1] 0.3533284Para resolver estos ejercicios que se basan en la distribución binomial, podemos usar la función dbinom() para obtener la probabilidad de un número exacto de éxitos y pbinom() para probabilidades acumuladas.
La fórmula para la probabilidad en una distribución binomial es:
\[ P(X = k) = \binom{n}{k} p^k (1 - p)^{n - k} \]
Donde: - \(n\) es el número de pruebas, - \(k\) es el número de éxitos deseados, - \(p\) es la probabilidad de éxito en cada prueba, - \(\binom{n}{k}\) es el coeficiente binomial, calculado como \(\frac{n!}{k!(n - k)!}\).
Ejercicio 1: ¿Cuál es la probabilidad de obtener exactamente 3 caras en 5 lanzamientos de una moneda justa?
# Parámetros
n <- 5 # Número de lanzamientos
p <- 0.5 # Probabilidad de obtener cara
k <- 3 # Número de caras exactas
# Probabilidad de obtener exactamente 3 caras
prob_3_caras <- dbinom(k, size = n, prob = p)
prob_3_caras## [1] 0.3125Ejercicio 2: Un test tiene una probabilidad de éxito del 0.7. ¿Cuál es la probabilidad de obtener exactamente 4 éxitos en 6 pruebas?
# Parámetros
n <- 6 # Número de pruebas
p <- 0.7 # Probabilidad de éxito
k <- 4 # Número de éxitos exactos
# Probabilidad de obtener exactamente 4 éxitos
prob_4_exitos <- dbinom(k, size = n, prob = p)
prob_4_exitos## [1] 0.324135Ejercicio 3: Si un dado se lanza 4 veces, ¿cuál es la probabilidad de obtener al menos un seis?
# Parámetros
n <- 4 # Número de lanzamientos
p <- 1/6 # Probabilidad de obtener un seis
# Probabilidad de no obtener ningún seis en 4 lanzamientos
prob_ningun_seis <- dbinom(0, size = n, prob = p)
# Probabilidad de obtener al menos un seis
prob_al_menos_un_seis <- 1 - prob_ningun_seis
prob_al_menos_un_seis## [1] 0.5177469Ejercicio Actuarial 1: Un modelo de riesgo en seguros usa una distribución binomial para representar el número de siniestros en un grupo de 15 pólizas, con una probabilidad de siniestro de 0.2 por póliza. Calcula la probabilidad de tener exactamente 4 siniestros.
# Parámetros
n <- 15 # Número de pólizas
p <- 0.2 # Probabilidad de siniestro
k <- 4 # Número exacto de siniestros
# Probabilidad de tener exactamente 4 siniestros
prob_4_siniestros <- dbinom(k, size = n, prob = p)
prob_4_siniestros## [1] 0.1876042Ejercicio Actuarial 2: En un estudio actuarial, se estima que la probabilidad de un evento de pérdida en una póliza de seguro es p=0.3 y el número de pólizas es n=10. ¿Cuál es la probabilidad de que haya al menos 3 eventos de pérdida?
# Parámetros
n <- 10 # Número de pólizas
p <- 0.3 # Probabilidad de un evento de pérdida
k <- 3 # Número mínimo de eventos de pérdida
# Probabilidad de tener menos de 3 eventos de pérdida
prob_menos_de_3 <- pbinom(k - 1, size = n, prob = p)
# Probabilidad de tener al menos 3 eventos de pérdida
prob_al_menos_3 <- 1 - prob_menos_de_3
prob_al_menos_3## [1] 0.6172172La distribución binomial negativa en R implica utilizar la función dnbinom(), que calcula la probabilidad de obtener un número específico de fracasos antes de alcanzar el 𝑟-ésimo éxito
Ejercicio 1: ¿Cuál es la probabilidad de obtener el tercer éxito en el quinto intento, si la probabilidad de éxito en cada intento es 0.4?
# Parámetros
r <- 3 # Número de éxitos deseados
x <- 4 # Número de fracasos antes del tercer éxito (el quinto intento incluye el éxito)
p <- 0.4 # Probabilidad de éxito
# Probabilidad de obtener el tercer éxito en el quinto intento
probabilidad <- dnbinom(x, size = r - 1, prob = p)
probabilidad## [1] 0.10368Ejercicio 2: En un experimento, se necesita 7 intentos para obtener el segundo éxito. Si la probabilidad de éxito en cada intento es 0.5, ¿cuál es la probabilidad de que se necesiten exactamente 7 intentos?
# Parámetros
r <- 2 # Número de éxitos deseados
x <- 6 # Número de fracasos antes del segundo éxito (el séptimo intento incluye el éxito)
p <- 0.5 # Probabilidad de éxito
# Probabilidad de obtener el segundo éxito en el séptimo intento
probabilidad <- dnbinom(x, size = r - 1, prob = p)
probabilidad## [1] 0.0078125Ejercicio 3: Si la probabilidad de éxito en cada intento es 0.3 y se necesitan 4 éxitos, ¿cuál es la probabilidad de obtener el cuarto éxito en el décimo intento?
# Parámetros
r <- 4 # Número de éxitos deseados
x <- 9 # Número de fracasos antes del cuarto éxito (el décimo intento incluye el éxito)
p <- 0.3 # Probabilidad de éxito
# Probabilidad de obtener el cuarto éxito en el décimo intento
probabilidad <- dnbinom(x, size = r - 1, prob = p)
probabilidad## [1] 0.05992511Ejercicio Actuarial 1: Un asegurador modela el número de años hasta que un cliente presenta su tercer siniestro utilizando una distribución binomial negativa con r=3 y p=0.2. ¿Cuál es la probabilidad de que el tercer siniestro ocurra en el sexto año?
# Parámetros
r <- 3 # Número de éxitos deseados
x <- 5 # Número de fracasos antes del tercer siniestro (el sexto año incluye el éxito)
p <- 0.2 # Probabilidad de éxito
# Probabilidad de obtener el tercer siniestro en el sexto año
probabilidad <- dnbinom(x, size = r - 1, prob = p)
probabilidad## [1] 0.0786432Ejercicio Actuarial 2: En el análisis de riesgos, se modela el número de reclamaciones hasta que un cliente presenta su cuarto siniestro con una distribución binomial negativa con r=4 y p=0.25. ¿Cuál es la probabilidad de que se necesiten 10 reclamaciones para que ocurra el cuarto siniestro?
# Parámetros
r <- 4 # Número de éxitos deseados
x <- 9 # Número de fracasos antes del cuarto siniestro (el décimo intento incluye el éxito)
p <- 0.25 # Probabilidad de éxito
# Probabilidad de obtener el cuarto siniestro en el décimo intento
probabilidad <- dnbinom(x, size = r - 1, prob = p)
probabilidad## [1] 0.0645259Para resolver problemas basados en la distribución hipergeométrica en R, utilizamos la función dhyper(), que calcula la probabilidad de obtener un número específico de éxitos en una muestra de una población finita sin reemplazo.
Ejercicio 1: En una urna con 20 bolas rojas y 30 bolas verdes, se extraen 10 bolas sin reemplazo. ¿Cuál es la probabilidad de que exactamente 4 de las bolas extraídas sean rojas?
# Parámetros
m <- 20 # Número total de bolas rojas
n <- 30 # Número total de bolas verdes
k <- 10 # Número total de bolas extraídas
x <- 4 # Número de bolas rojas extraídas
# Probabilidad de obtener exactamente 4 bolas rojas
probabilidad <- dhyper(x, m, n, k)
probabilidad## [1] 0.2800586Ejercicio 2: Se seleccionan 5 cartas de una baraja de 52 cartas sin reemplazo. ¿Cuál es la probabilidad de que exactamente 2 de las cartas sean tréboles?
# Parámetros
m <- 13 # Número total de tréboles
n <- 39 # Número total de cartas que no son tréboles
k <- 5 # Número total de cartas seleccionadas
x <- 2 # Número de tréboles seleccionados
# Probabilidad de obtener exactamente 2 tréboles
probabilidad <- dhyper(x, m, n, k)
probabilidad## [1] 0.2742797Ejercicio Actuarial 1: Un asegurador selecciona al azar una muestra de 15 pólizas de una cartera de 100 pólizas, de las cuales 20 tienen un siniestro. Calcula la probabilidad de que exactamente 5 de las pólizas seleccionadas tengan un siniestro.
# Parámetros
m <- 20 # Número total de pólizas con siniestro
n <- 80 # Número total de pólizas sin siniestro
k <- 15 # Número total de pólizas seleccionadas
x <- 5 # Número de pólizas con siniestro seleccionadas
# Probabilidad de obtener exactamente 5 pólizas con siniestro
probabilidad <- dhyper(x, m, n, k)
probabilidad## [1] 0.1007633Para resolver problemas basados en la distribución uniforme continua en R, utilizamos la fórmula de la función de densidad uniforme. La fórmula para calcular la probabilidad en una distribución uniforme continua 𝑈(𝑎,𝑏) es:
La fórmula para calcular la probabilidad en una distribución uniforme continua \(U(a, b)\) es:
\[ P(x_1 \leq X \leq x_2) = \frac{x_2 - x_1}{b - a} \]
donde:
\(a\) y \(b\) son los límites inferior y superior del intervalo de la distribución uniforme,
\(x_1\) y \(x_2\) son los límites del intervalo en el cual queremos calcular la probabilidad.
Ejercicio 2: Una variable continua sigue una distribución uniforme en el intervalo [0,10]. ¿Cuál es la probabilidad de que el valor esté entre 4 y 6?
# Parámetros
a <- 0 # Límite inferior de la distribución
b <- 10 # Límite superior de la distribución
x1 <- 4 # Límite inferior del intervalo de interés
x2 <- 6 # Límite superior del intervalo de interés
# Probabilidad de que el valor esté entre 4 y 6
probabilidad <- (x2 - x1) / (b - a)
probabilidad## [1] 0.2Ejercicio 3: La duración de un proceso se distribuye uniformemente entre 10 y 20 minutos. ¿Cuál es la probabilidad de que el proceso dure más de 15 minutos?
# Parámetros
a <- 10 # Límite inferior de la distribución
b <- 20 # Límite superior de la distribución
x1 <- 15 # Límite inferior del intervalo de interés
# Probabilidad de que el proceso dure más de 15 minutos
probabilidad <- (b - x1) / (b - a)
probabilidad## [1] 0.5Ejercicio Actuarial 1: Un asegurador estima que el tiempo de resolución de reclamaciones sigue una distribución uniforme continua entre 5 y 15 días. ¿Cuál es la probabilidad de que la resolución de una reclamación tome entre 8 y 12 días?
# Parámetros
a <- 5 # Límite inferior de la distribución
b <- 15 # Límite superior de la distribución
x1 <- 8 # Límite inferior del intervalo de interés
x2 <- 12 # Límite superior del intervalo de interés
# Probabilidad de que la resolución de una reclamación tome entre 8 y 12 días
probabilidad <- (x2 - x1) / (b - a)
probabilidad## [1] 0.4Para resolver estos problemas basados en la distribución normal en R, utilizamos las funciones pnorm() para calcular las probabilidades acumuladas y dnorm() para obtener la densidad.
Ejercicio 2: La estatura de una población sigue una distribución normal con media 170 cm y desviación estándar 8 cm. ¿Cuál es la probabilidad de que una persona seleccionada al azar tenga una estatura entre 160 y 180 cm?
# Parámetros
media <- 170
desviacion_estandar <- 8
x1 <- 160
x2 <- 180
# Probabilidad de que la estatura esté entre 160 y 180 cm
probabilidad <- pnorm(x2, mean = media, sd = desviacion_estandar) - pnorm(x1, mean = media, sd = desviacion_estandar)
probabilidad## [1] 0.7887005Ejercicio 3: En un test estandarizado con media 100 y desviación estándar 15, ¿cuál es la probabilidad de obtener una puntuación superior a 120?
# Parámetros
media <- 100
desviacion_estandar <- 15
valor <- 120
# Probabilidad de obtener una puntuación superior a 120
probabilidad <- 1 - pnorm(valor, mean = media, sd = desviacion_estandar)
probabilidad## [1] 0.09121122Ejercicio Actuarial 1: Un asegurador utiliza una distribución normal para modelar la pérdida total en una póliza con una media de $100,000 y una desviación estándar de $20,000. ¿Cuál es la probabilidad de que la pérdida total sea menor de $80,000?
# Parámetros
media <- 100000
desviacion_estandar <- 20000
valor <- 80000
# Probabilidad de que la pérdida total sea menor de $80,000
probabilidad <- pnorm(valor, mean = media, sd = desviacion_estandar)
probabilidad## [1] 0.1586553Ejercicio Actuarial 2: La duración de vida útil de un producto sigue una distribución normal con media de 4 años y desviación estándar de 0.5 años. ¿Cuál es la probabilidad de que el producto dure entre 3.5 y 4.5 años?
# Parámetros
media <- 4
desviacion_estandar <- 0.5
x1 <- 3.5
x2 <- 4.5
# Probabilidad de que el producto dure entre 3.5 y 4.5 años
probabilidad <- pnorm(x2, mean = media, sd = desviacion_estandar) - pnorm(x1, mean = media, sd = desviacion_estandar)
probabilidad## [1] 0.6826895