Examen Parcial 1

Examen Parcial, se utilizó una simulación de montecarlo

Para la producción de 10000 unidades

P1 <- 10000
simulacion1<-function(nsim=100){
  data<-data.frame(venta=0,
                   Ganancia=0)
  for(i in 1:nsim){
    venta <- (sample(demanda,1,prob = pdemanda))
    if (venta>= P1){
      Ganancia = (P1*(Precio-Costo)) - ((venta-P1)*(Costo)) + ((venta-P1)*Rescate)
    }
    if (venta < P1){
      Ganancia = (venta* (Precio - Costo)) - ((P1-venta)*(Costo)) + ((P1-venta)*Rescate)
    }
    data <- rbind(data, c(venta, Ganancia))
  }
  data<-data[2:nrow(data),]
  row.names(data)<-NULL
  return(data)
}  

Sim1<-simulacion1(10000)

hist(Sim1$Ganancia,xlab = "Ganancias ($)",main = "Simulación con Producción de 10000 UN")

Para la producción de 20000 unidades

P2 <- 20000
simulacion2<-function(nsim=100){
  data<-data.frame(venta=0,
                   Ganancia=0)
  for(i in 1:nsim){
    venta <- (sample(demanda,1,prob = pdemanda))
    if (venta>= P2){
      Ganancia = (P2*(Precio-Costo)) - ((venta-P1)*(Costo)) + ((venta-P2)*Rescate)
    }
    if (venta < P1){
      Ganancia = (venta* (Precio - Costo)) - ((P2-venta)*(Costo)) + ((P2-venta)*Rescate)
    }
    data <- rbind(data, c(venta, Ganancia))
  }
  data<-data[2:nrow(data),]
  row.names(data)<-NULL
  return(data)
}  

Sim2<-simulacion2(10000)

hist(Sim2$Ganancia,xlab = "Ganancias ($)",main = "Simulación con Producción de 20000 UN")

Para la producción de 40000 unidades

P3 <- 40000
simulacion3<-function(nsim=100){
  data<-data.frame(venta=0,
                   Ganancia=0)
  for(i in 1:nsim){
    venta <- (sample(demanda,1,prob = pdemanda))
    if (venta>= P3){
      Ganancia = (P3*(Precio-Costo)) - ((venta-P3)*(Costo)) + ((venta-P3)*Rescate)
    }
    if (venta < P3){
      Ganancia = (venta* (Precio - Costo)) - ((P3-venta)*(Costo)) + ((P3-venta)*Rescate)
    }
    data <- rbind(data, c(venta, Ganancia))
  }
  data<-data[2:nrow(data),]
  row.names(data)<-NULL
  return(data)
}  

Sim3<-simulacion3(10000)

hist(Sim3$Ganancia,xlab = "Ganancias ($)",main = "Simulación con Producción de 40000 UN")

Para la producción de 60000 unidades

P4 <- 60000
simulacion4<-function(nsim=100){
  data<-data.frame(venta=0,
                   Ganancia=0)
  for(i in 1:nsim){
    venta <- (sample(demanda,1,prob = pdemanda))
    if (venta>= P4){
      Ganancia = (P4*(Precio-Costo)) - ((venta-P4)*(Costo)) + ((venta-P4)*Rescate)
    }
    if (venta < P4){
      Ganancia = (venta* (Precio - Costo)) - ((P4-venta)*(Costo)) + ((P4-venta)*Rescate)
    }
    data <- rbind(data, c(venta, Ganancia))
  }
  data<-data[2:nrow(data),]
  row.names(data)<-NULL
  return(data)
}  

Sim4<-simulacion4(10000)

hist(Sim4$Ganancia,xlab = "Ganancias ($)",main = "Simulación con Producción de 60000 UN")

ResultFinales <- data.frame(Sim1$Ganancia,Sim2$Ganancia, Sim3$Ganancia,Sim4$Ganancia)
names(ResultFinales) <- c("Ganancias Sim1","Ganancias Sim2", "Ganancias Sim3","Ganancias Sim4")

Promedios <- data.frame(mean(Sim1$Ganancia), mean(Sim2$Ganancia), mean(Sim3$Ganancia), mean(Sim4$Ganancia))
names(Promedios) <- c("Ganancia con Producción 10000","Ganancia con Producción 20000", "Ganancia con Producción 40000","Ganancia con Producción 60000")


desvest <- data.frame(sd(Sim1$Ganancia), sd(Sim2$Ganancia), sd(Sim3$Ganancia), sd(Sim4$Ganancia))
names(desvest) <- c("Desviación  Producción 10000","Desviacióncon Producción 20000", "Desviacióncon Producción 40000","Desviación Producción 60000")

print(Promedios)

##   Ganancia con Producción 10000 Ganancia con Producción 20000
## 1                       -7598.8                       11922.4
##   Ganancia con Producción 40000 Ganancia con Producción 60000
## 1                       55760.4                       54129.8

print(desvest)

##   Desviación  Producción 10000 Desviacióncon Producción 20000
## 1                     22585.77                        21103.8
##   Desviacióncon Producción 40000 Desviación Producción 60000
## 1                       38815.17                    66165.61

La opción que da mayor ganancia y menos riesgo es la opción de 40000 tarjetas a producir, porque tiene un promedio de $55,760.40 de ganancia y una desviación de 38,815.17, es decir que en cualquier escenario existe una ganacia.

La opción que da una ganancia alta pero no la mayor y mayor riesgo es la opción con 60000 unidades de producción con una ganancia promedio de $54,129.80 y una desviación de 66,165.61, en este caso la desviación es mayor al promedio por lo que existe un riesgo muy alto de no tener ganancia algua.