set.seed(123)
tiempos_atencion <- rnorm(20, mean = 15, sd = 5)
media_atencion <- mean(tiempos_atencion)
sd_atencion <- sd(tiempos_atencion)
error_atencion <- qt(0.975, df = 19) * sd_atencion / sqrt(20)
ic_atencion <- c(media_atencion - error_atencion, media_atencion + error_atencion)
print("Tiempos de atención (Normal):")
## [1] "Tiempos de atención (Normal):"
print(tiempos_atencion)
## [1] 12.197622 13.849113 22.793542 15.352542 15.646439 23.575325 17.304581
## [8] 8.674694 11.565736 12.771690 21.120409 16.799069 17.003857 15.553414
## [15] 12.220794 23.934566 17.489252 5.166914 18.506780 12.636043
tiempos_llegada <- rexp(20, rate = 1/3.5)
media_llegada <- mean(tiempos_llegada)
sd_llegada <- sd(tiempos_llegada)
error_llegada <- qt(0.975, df = 19) * sd_llegada / sqrt(20)
ic_llegada <- c(media_llegada - error_llegada, media_llegada + error_llegada)
print("Tiempos entre llegadas (Exponencial):")
## [1] "Tiempos entre llegadas (Exponencial):"
print(tiempos_llegada)
## [1] 5.3504309 4.7296638 4.7181557 4.0898514 5.6204832 5.2386000 5.4972839
## [8] 0.1111871 2.0924739 7.5874391 1.7731551 0.9084524 9.0891224 4.3015901
## [15] 2.7673862 2.2024803 4.3912435 2.0603962 3.9525151 1.4712768
resultados <- data.frame(
Distribucion = c("Normal (Atención)", "Exponencial (Llegadas)"),
Media = c(media_atencion, media_llegada),
Desviacion_Estandar = c(sd_atencion, sd_llegada),
IC_95_Inferior = c(ic_atencion[1], ic_llegada[1]),
IC_95_Superior = c(ic_atencion[2], ic_llegada[2])
)
print("Tabla resumen de resultados:")
## [1] "Tabla resumen de resultados:"
print(resultados)
## Distribucion Media Desviacion_Estandar IC_95_Inferior
## 1 Normal (Atención) 15.708119 4.863327 13.432012
## 2 Exponencial (Llegadas) 3.897659 2.250994 2.844162
## IC_95_Superior
## 1 17.984226
## 2 4.951157