TALLER_VAR_ALEATORIAS2

1. 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 de fallas en un semestre (150 días) y calcular la media y desviación estándar.

set.seed(777)

fallas <- rpois(150, 3)     
mean(fallas)        

## [1] 2.973333

sd(fallas)                

## [1] 1.749196

barplot(table(fallas),
main = "Fallas por día (Poisson λ=3)",
xlab = "Número de fallas",
ylab = "Frecuencia")

2. 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

set.seed(777)

vidas <- rexp(1000, rate = 1/500) 
mean(vidas > 700)                 

## [1] 0.254

hist(vidas, breaks = 30,
main = "Vida útil de 1000 componentes",
xlab = "Horas")

3. En una línea de ensamblaje, la probabilidad de que un producto sea defectuoso es del 5%. Simular 100 lotes de 50 productos y calcular el número promedio de productos defectuosos por lote.

set.seed(777)

defectuosos <- rbinom(100, 50, 0.05)  
mean(defectuosos)

## [1] 2.53

barplot(table(defectuosos),
main = "Defectuosos por lote (Binomial)",
xlab = "Número de defectuosos",
ylab = "Frecuencia")

4. La demanda diaria de energía (en MW) sigue una distribución normal con media de 100 MW y desviación estándar de 15 MW. Simular la demanda de un año (365 días) y calcular la probabilidad de que un día supere los 130 MW. y realizar el histograma.

set.seed(777)

demanda <- rnorm(365, mean = 100, sd = 15)
mean(demanda > 130)                          

## [1] 0.03013699

hist(demanda, breaks = 30,
main = "Demanda diaria de energía",
xlab = "Demanda (MW)")

5. 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.

1. Generar 1000 tiempos de vida del capacitor aplicando el método de la transformada inversa.
2. Estimar la media y la varianza de los tiempos generados y compararlas con los valores teóricos.
3. Gracar el histograma de los tiempos de vida simulados junto con la densidad teórica de la distribución exponencial.
4. Calcular la probabilidad de que un capacitor dure menos de 940 horas usando la simulación.

set.seed(777)

# PUNTO A
u <- runif(1000)             
beta <- 1000
tiempos <- -beta * log(u)     

# PUNTO B
# Valores teóricos: media=1000, varianza=1000000

mean(tiempos)   # Media simulada

## [1] 1044.722

var(tiempos)    # Varianza simulada

## [1] 1075728

# PUNTO C
hist(tiempos, breaks = 30, freq = FALSE,
     main = "Tiempos de vida de capacitores",
     xlab = "Horas")

# PUNTO D
#Probabilidad de duracion <940hrs
mean(tiempos < 940)

## [1] 0.592