Introducción

Este taller práctico tiene como objetivo principal poner en práctica las técnicas de simulación de variables aleatorias mediante el uso del entorno estadístico R Studio. La simulación constituye una herramienta esencial para representar fenómenos aleatorios de alta complejidad y obtener estimaciones consistentes de sus parámetros, siendo ampliamente utilizada en ingeniería, finanzas y gestión de riesgos.

A lo largo de los ejercicios se abordarán las distribuciones de probabilidad más representativas para la modelación de fenómenos discretos y continuos: la distribución de Poisson (adecuada para eventos poco frecuentes en un intervalo), la distribución exponencial (empleada en tiempos de espera o vida útil), la distribución binomial (para conteos de éxitos y fracasos) y la distribución normal (aplicable a variables continuas con comportamiento centrado). Dichos modelos se contextualizarán en escenarios operativos relevantes como la producción industrial, la confiabilidad de dispositivos y la variación de la demanda energética.

Proceso de Fallas en un Sistema de Producción
Un sistema de producción tiene fallas según un proceso de Poisson con una tasa de 3 fallas por día. Simular el número total de fallas en un semestre (150 días) y calcular la media y desviación estándar.

Análisis Teórico: La Distribución de Poisson es ideal para modelar el número de ocurrencias de eventos discretos en un intervalo fijo de tiempo, como las fallas del sistema

set.seed(123)
fallas <- rpois(150, lambda = 3) # Simulación de 150 días con media de 3 fallas
mean_fallas <- mean(fallas) # Cálculo de la media
sd_fallas <- sd(fallas) # Cálculo de la desviación estándar

mean_fallas

## [1] 3

sd_fallas

## [1] 1.658818

Análisis
El código proporcionado simula el número de fallas diarias durante un semestre (150 días), asumiendo un proceso Poisson con una tasa media (λ) de 3 fallas por día.

Vida útil de un Componente Electrónico El tiempo de servicio (expresado en horas) de un elemento electrónico obedece a una distribución exponencial con una expectativa de 500 horas. Se requiere generar 1000 muestras de estos elementos y calcular la probabilidad estimada de que un componente individual permanezca funcional por más de 700 horas.

siendo utilizada frecuentemente para describir la duración, o el intervalo de tiempo, hasta la ocurrencia de un evento (como una avería o fin de vida

set.seed(123) # Semilla para obtener siempre el mismo resultado
vida <- rexp(1000, rate = 1/500) # Simulación de 1000 tiempos de vida con media 500
prob_mas700 <- mean(vida > 700) # Estimar probabilidad de durar más de 700 horas

prob_mas700

## [1] 0.255

Análisis
La probabilidad estimada por simulación de que un componente dure más de 700 horas es del 25.2%. Este valor se acerca significativamente al valor teórico (P(X>700)=e −0.002×700 ≈0.2466), lo que valida la simulación como una buena aproximación del modelo Exponencial.

Productos defectuosos en una Línea de Emsablaje
La vida útil (en horas) de un componente electrónico sigue una distribución exponencial con un promedio de 500 horas. Simular 1000 componentes y estimar la probabilidad de que un componente dure más de 700 horas.

La Distribución Binomial permite modelar el número de éxitos (defectuosos) en un número de fijo de ensayos (productos).

set.seed(123) # Semilla
defectuosos <- rbinom(100, size = 50, prob = 0.05) # 100 lotes de 50 productos
mean_defectuosos <- mean(defectuosos) # Promedio de defectuosos por lote

mean_defectuosos

## [1] 2.48

Análisis
El promedio de defectuosos por lote se acerca al valor esperado 50 x 0.05 = 2.5. Esto evidencia la consistencia del modelo binomial para este tipo de procesos industriales.

Demanda Diaria de Energía
La solicitud energética diaria (medida en MW) se ajusta a una distribución normal con una media (μ) de 100 MW y una desviación estándar (σ) de 15 MW. Se pide simular la demanda para el periodo de un año (365 días), calcular la probabilidad de que la demanda en un día cualquiera exceda los 130 MW, y generar el histograma correspondiente.

La Distribución Normal modela adecuadamente variables continuas con fluctuaciones alrededor de una media.

set.seed(123) # Semilla
demanda <- rnorm(365, mean = 100, sd = 15) # Simulación de 365 días
prob_mas130 <- mean(demanda > 130) # Probabilidad de superar 130 MW

prob_mas130

## [1] 0.03013699

# Creación del data frame para graficar
df_demanda <- data.frame(demanda)

# Histograma + Curva de Densidad
ggplot(df_demanda, aes(x = demanda)) +
  geom_histogram(aes(y = ..density..), bins = 30, fill = "black", color = "red", alpha = 0.7) +
  geom_density(color = "yellow", size = 1) +
  geom_vline(xintercept = 130, linetype = "dashed", color = "blue") +
  labs(title = "Histograma de la Demanda Diaria (MW)",
       x = "Demanda (MW)", y = "Densidad")

Análisis
El histograma confirma la forma de campana característica de la normal. La probabilidad de superar los 130 MW es baja, lo cual concuerda con que dicho valor está más de 2 desviaciones estándar por encima de la media.

Tiempo de Vida de Capacitadores Electrónicos Una empresa de manufactura electrónica quiere simular el tiempo de vida (en horas) de un nuevo modelo de capacitor. Basado en datos históricos, se ha determinado que el tiempo de vida sigue una distribución exponencial con parámetro β = 1000 horas, que representa el tiempo medio de vida de los capacitores.

a) Generar 1000 tiempos de vida del capacitor aplicando el método de la transformada inversa. X = −β⋅ln(1−U), U∼U(0,1)

set.seed(123) # Semilla
n <- 1000 # Número de Simulaciones
beta <- 1000 # Parámetro de la Distribución Exponencial
U <- runif(n)
vida_cap <- -beta * log(1 - U)

b) Estimar la media y la varianza de los tiempos generados y compararlas con los valores teóricos.

Valores Teóricos
E|X| = β = 1000, Var[X] = β^2 = 1000000

media_sim <- mean(vida_cap) # Media simulada
var_sim <- var(vida_cap) # Varianza simulada

media_sim

## [1] 986.1544

var_sim

## [1] 954966.2

Los valores simulados se aproximan a los teóricos, validando la simulación.

c) Graficar el histograma de los tiempos de vida simulados junto con la densidad teórica de la distribución exponencial.

df_vida <- data.frame(vida_cap)
ggplot(df_vida, aes(x = vida_cap)) +
  geom_histogram(aes(y = ..density..), bins = 30, fill = "red", color = "black", alpha = 0.7) +
  stat_function(fun = dexp, args = list(rate = 1/beta), color = "blue", size = 1) +
  labs(title = "Tiempo de Vida de Capacitores",
       x = "Tiempo de vida (horas)", y = "Densidad")

Los Valores simulados se aproximan a los teóricos, validando la simulación.

d) Calcular la probabilidad de que un capacitor dure menos de 940 horas usando la simulación.

prob_menos940 <- mean(vida_cap < 940) # Probabilidad Simulada
prob_menos940

## [1] 0.611

El resultado simulado concuerda con el cálculo teórico P(X < 940) = 1 -e ^-940/1000. Esto permite estimar la confiabilidad de los capacitadores en la práctica.

Variables Aleatorias 2

Jose anaya

2025

Introducción

Conclusión