1. CORRE PRIMERO ESTE CHUNK

Calculo presupuesto Ejemplo-1

Simulacion de MONTECARLO

Establecer el numero de simulaciones

n <- 947 # SIMULACIONES

Estimaciones aleatorias de los costos de las 6 Actividades

set.seed(123)
actA <- runif(n, min = 10e3, max = 20e3)
hist(actA)

summary(actA)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   10005   12535   14923   14986   17504   19994

Estimacion de las otras 5 actividades

set.seed(123)
actB <- runif(n, min = 15e3, max = 15e3)
actC <- runif(n, min = 7.5e3, max = 12e3)
actD <- runif(n, min = 4.8e3, max = 6.2e3)
actE <- runif(n, min = 20e3, max = 25e3)
actF <- runif(n, min = 5e3, max = 7e3)

Calculo de Montecarlo

Presupuesto1 <- actA + actB + actC + actD + actE  + actF
hist(Presupuesto1, breaks = 50, main = "Project Cost Distribution",
     xlab = "Total Cost", col = "skyblue", border = "white")

Histograma con curva normal superpuesta

hist(Presupuesto1, freq = FALSE, main = "Histograma con curva normal",
     xlab = "Presupuesto Total", col = "skyblue", border = "white", breaks = 50)

# Curva normarl estandar
curve(dnorm(x, mean = mean(Presupuesto1), sd = sd(Presupuesto1)), 
      col = "red", lwd = 2, add = TRUE)

Q-Q plot (gráfico cuantil-cuantil)

qqnorm(Presupuesto1, main = "Gráfico Q-Q Normal")
qqline(Presupuesto1, col = "red", lwd = 2)

# Restaurar configuración de gráficos
par(mfrow = c(2, 1))

# Gráfico de densidad
plot(density(Presupuesto1), main = "Gráfico de densidad", 
     xlab = "Presupuesto1 Total", ylab = "Densidad")
curve(dnorm(x, mean = mean(Presupuesto1), sd = sd(Presupuesto1)), 
      col = "red", lwd = 2, add = TRUE)

# Boxplot para ver simetría y outliers
boxplot(Presupuesto1, main = "Presupuesto1", horizontal = TRUE, col = "skyblue")