Laboratorio #2 - Dilema del Vendedor de Periodico

Cesar Tinoco Alvarez - 13003387

Creamos vectores con cada una de las probabilidades segun el dia

set.seed(25)
tipo_dia <- c("Excelente", "Bueno", "Malo")
prob_tipo <- c(0.35, 0.45, 0.2)
prob_demanda <- c(40, 50, 60, 70, 80, 90, 100)
prob_excel <- c(0.03, 0.05, 0.15, 0.2, 0.35, 0.15, 0.07)
prob_bueno <- c(0.1, 0.18, 0.4, 0.2, 0.08, 0.04, 0)
prob_malo <- c(0.44, 0.22, 0.16, 0.12, 0.06, 0, 0)
diatipo <- numeric(20)
dem<- numeric(20)
resultado<-NULL
Result_Final<-NULL

Se crea un ciclo para recorrer cada una de las combinaciones posibles

for  (x in c(40,50,60,70,80,90,100)){
  earn <- 0
for (d in 1:20) {
    diatipo[d]<- sample(tipo_dia, size=1, replace = TRUE, prob=prob_tipo) 
    if(diatipo[d]=="Excelente"){dem[d] <- sample(prob_demanda, size=1, replace = TRUE, prob=prob_excel)}
    else if(diatipo[d]=="Bueno"){dem[d] <- sample(prob_demanda, size=1, replace = TRUE, prob=prob_bueno)}
    else if(diatipo[d]=="Malo"){dem[d] <- sample(prob_demanda, size=1, replace = TRUE, prob=prob_malo)}
   
    if(dem[d]<x)
    { 
      ganancias <- (dem[d]*0.5)- (x*0.33)  + ((x-dem[d])*0.05)
    } 
    else if (dem[d]==x)
    {
      ganancias <- (dem[d]*0.5)-(x*0.33)
    }
    else if (dem[d]>x){
      ganancias <- (x*(0.5-0.33))  - ((dem[d]-x)*(0.5-0.33))
    }
    earn<-earn+ganancias 
    nuevo <- data.frame (dia=d,Compra=x,Clientes=dem[d],Ganancia=ganancias)
    resultado <- rbind(resultado, nuevo) 
  } 
  earn <- earn/ 20
  a <- data.frame (Compra=x,Ganancia=earn)
  Result_Final <- rbind(Result_Final, a) 
} 

Se imprime el resultado

print(resultado)

Se muestra el resultado final

print(Result_Final)

Se muestran las graficas donde relaciona la ganancia y la compra

plot(resultado)

plot(Result_Final)

boxplot(formula = Ganancia ~ Compra, main= "Ganancia x Compra", data =  resultado,xlab = "Demanda", ylab = "Ganancia")