R Notebook

Problema del Vendedor de Periódicos

Información del problema

tiposDia <- c("Excelente", "Bueno", "Malo")
probabilidadTipo <- c(0.35, 0.45, 0.2)
demandas <- c(40, 50, 60, 70, 80, 90, 100)
probabilidadesExcelente <- c(0.03, 0.05, 0.15, 0.2, 0.35, 0.15, 0.07)
probabilidadesBueno <- c(0.1, 0.18, 0.4, 0.2, 0.08, 0.04, 0)
probabilidadesMalo <- c(0.44, 0.22, 0.16, 0.12, 0.06, 0, 0)
set.seed(25)
diatipo <- numeric(20)
demanda<- numeric(20)
resultado<-NULL
resultadoFinal<-NULL

##Ciclo por pedido
for (pedido in c(40,50,60,70,80,90,100)){
  utilidad <- 0
  ##Ciclo por día
  for (dia in 1:20) {
    ##Se toma el tipo de venta por día como una variable que cambia día con día.
    diatipo[dia]<- sample(tiposDia, size=1, replace = TRUE, prob=probabilidadTipo)
    
    if(diatipo[dia]=="Excelente"){
      demanda[dia] <- sample(demandas, size=1, replace = TRUE, prob=probabilidadesExcelente)
    }
    else if(diatipo[dia]=="Bueno"){
      demanda[dia] <- sample(demandas, size=1, replace = TRUE, prob=probabilidadesBueno)
    }
    else if(diatipo[dia]=="Malo"){
      demanda[dia] <- sample(demandas, size=1, replace = TRUE, prob=probabilidadesMalo)
    }
    if(demanda[dia]<pedido)
    { 
      ganancias <- (demanda[dia]*0.5)- (pedido*0.33)  + ((pedido-demanda[dia])*0.05)
    } 
    else if (demanda[dia]==pedido)
    {
      ganancias <- (demanda[dia]*0.5)-(pedido*0.33)
    }
    else if (demanda[dia]>pedido){
      ganancias <- (pedido*(0.5-0.33))-((demanda[dia]-pedido)*(0.5-0.33))
    }
    utilidad<-utilidad+ganancias
    nuevo <- data.frame(Dia=dia,Compra=pedido,Clientes=demanda[dia],Utilidad=ganancias)
    ##Data frame con los resultados de la simulación
    resultado <- rbind(resultado, nuevo) 
  } 
  utilidad <- utilidad/ 20
  promedioGananciasPedido <- data.frame (Compra=pedido,Utilidad=utilidad)
  ##Data frame con el promedio de ganancias por cada pedido
  resultadoFinal <- rbind(resultadoFinal, promedioGananciasPedido) 
} 

View(print(resultado))

boxplot(formula = Utilidad ~ Compra, data =  resultado,main = "Utilidad x Demanda", 
        xlab = "Demanda", ylab = "Utilidad",
       col = c("orange", "yellow", "grey", "green", "blue", "red", "brown"))

Warning in gzfile(file, "wb") :
  cannot open compressed file 'C:/Users/Oscar Padilla/Desktop/Galileo/Postgrado en Investigacion de Operaciones/TRIMESTRE VII/Modelacion y Simulacion I/Laboratorio # 2/.Rproj.user/shared/notebooks/F4F4F36-Laboratorio # 2/1/F28451AFB616D64E/cke6ghesi9jss_t/2d174fb2996e452e824c8f853150203b.snapshot', probable reason 'No such file or directory'
Error in gzfile(file, "wb") : cannot open the connection