Este documento presenta la simulación de un sistema de soporte
técnico, donde se utilizan distribuciones de probabilidad para modelar
llegadas de solicitudes, tiempos entre llegadas, duración de la
atención y escalamiento de tickets.
El objetivo es observar cómo estos modelos permiten comprender mejor la
dinámica del sistema y apoyar decisiones de gestión.
# 1. Llegadas de solicitudes
set.seed(321)
llegadas <- rpois(5000, lambda = 5)
# Estadísticas descriptivas
media_llegadas <- mean(llegadas)
var_llegadas <- var(llegadas)
media_llegadas
## [1] 5.0066
var_llegadas
## [1] 5.053967
# Distribución de frecuencias
tabla_llegadas <- table(llegadas)
barplot(tabla_llegadas,
main = "Distribución de llegadas por hora",
xlab = "Número de solicitudes",
ylab = "Frecuencia observada",
col = "skyblue4")
# Probabilidad empírica de menos de 3 solicitudes
prob_llegadas_bajas <- mean(llegadas < 3)
prob_llegadas_bajas
## [1] 0.123
# 2. Tiempos entre llegadas
set.seed(321)
tiempos <- rexp(5000, rate = 1/12)
# Media empírica
mean(tiempos)
## [1] 11.81172
# Histograma con curva teórica
hist(tiempos,
breaks = 50,
col = "darkseagreen",
border = "white",
main = "Distribución de tiempos entre llegadas",
xlab = "Tiempo en minutos",
freq = FALSE)
curve(dexp(x, rate = 1/12), add = TRUE, col = "black", lwd = 2)
# 3. Duración de la atención (Normal con media=30, sd=5)
set.seed(321)
atencion <- rnorm(5000, mean = 30, sd = 5)
# Resumen
mean(atencion)
## [1] 30.07293
sd(atencion)
## [1] 5.016883
# Histograma
hist(atencion,
breaks = 35,
col = "lightgoldenrod",
main = "Duración de la atención",
xlab = "Minutos",
freq = FALSE)
curve(dnorm(x, mean = 30, sd = 5), add = TRUE, col = "darkred", lwd = 2)
# Probabilidad de que la atención dure más de 35 minutos
prob_atencion_larga <- mean(atencion > 35)
prob_atencion_larga
## [1] 0.1608
# 4. Tickets escalados a segundo nivel
set.seed(321)
escalados <- rbinom(1000, size = 10, prob = 0.2)
# Media empírica
mean(escalados)
## [1] 2.05
# Tabla de frecuencias
tabla_escalados <- table(escalados)
# Gráfico de barras
barplot(tabla_escalados,
main = "Tickets escalados por hora",
xlab = "Número de tickets escalados",
ylab = "Frecuencia",
col = "orchid")