Contexto

Hershey’s México empresa chocolatera con más de 100 años de historia a nivel global (fundación 1903) y establecida en México desde 1969, como una de confitería entre Hershey Food Corporation y Anderson Clayton & Co. S.A. formando Nacional de Dulces S.A. de C.V. en el Distrito Federal y después de 12 años cambiando sus instalaciones a El Salto, Jalisco e iniciando operaciones en el mes de febrero de 1981.

En 2022 Hershey´s México se consolidó como una de las empresas chocolateras del país, y con más fuerza en la venta de sus lechitas y ha sumado esfuerzos para que a pesar de la pandemia de la Covid-19 ofrezca el mejor y disponibilidad de productos a sus clientes.

Importar Librerías

Importar la base de datos

home <- read.xlsx("C:/Users/josea/OneDrive - Instituto Tecnologico y de Estudios Superiores de Monterrey/Escritorio/Tecnológico de Monterrey/LIT/Septimo Semestre/Econometria/SitProblema1/hersheysSales.xlsx")
home$Month <- as.Date(home$Month, origin = "1899-12-30")
describe(home)
## home 
## 
##  2  Variables      36  Observations
## --------------------------------------------------------------------------------
## Month 
##          n    missing   distinct       Info       Mean        Gmd        .05 
##         36          0         36          1 2018-06-16      375.1 2017-02-22 
##        .10        .25        .50        .75        .90        .95 
## 2017-04-16 2017-09-23 2018-06-16 2019-03-08 2019-08-16 2019-10-08 
## 
## lowest : 2017-01-01 2017-02-01 2017-03-01 2017-04-01 2017-05-01
## highest: 2019-08-01 2019-09-01 2019-10-01 2019-11-01 2019-12-01
## --------------------------------------------------------------------------------
## Sales 
##        n  missing distinct     Info     Mean      Gmd      .05      .10 
##       36        0       36        1    30415     3769    25781    26625 
##      .25      .50      .75      .90      .95 
##    27520    30475    32961    34562    34922 
## 
## lowest : 23740.11 25520.51 25868.43 26253.58 26996.11
## highest: 34324.12 34799.91 34846.17 35151.28 36133.07
## --------------------------------------------------------------------------------

1.- Utilizando modelos ARIMA (Box-Jenkins, ARMA, SARIMA) y los datos históricos de las ventas de leche saborizada ¿Cuál es el modelo que mejor se adapta a la serie?

#Modelo ARIMA
timeSeries <- ts(data=home$Sales, start = c(2017,1), end = c(2019,12),frequency=12)
timeSeries
##           Jan      Feb      Mar      Apr      May      Jun      Jul      Aug
## 2017 25520.51 23740.11 26253.58 25868.43 27072.87 27150.50 27067.10 28145.25
## 2018 28463.69 26996.11 29768.20 29292.51 29950.68 30099.17 30851.26 32271.76
## 2019 32496.44 31287.28 33376.02 32949.77 34004.11 33757.89 32927.30 34324.12
##           Sep      Oct      Nov      Dec
## 2017 27546.29 28400.37 27441.98 27852.47
## 2018 31940.74 32995.93 32197.12 31984.82
## 2019 35151.28 36133.07 34799.91 34846.17
modelo_arima <- auto.arima(timeSeries, D = 1)
pronostico_hersh <- forecast (modelo_arima, level = c(95), h = 5)
pronostico_hersh
##          Point Forecast    Lo 95    Hi 95
## Jan 2020       35498.90 34616.48 36381.32
## Feb 2020       34202.17 33155.29 35249.05
## Mar 2020       36703.01 35596.10 37809.92
## Apr 2020       36271.90 35141.44 37402.36
## May 2020       37121.98 35982.07 38261.90
plot(pronostico_hersh, main = "Pronóstico de venta de Leche Hershey's")

2.- ¿Qué modelo de regresión ofrece mejor exactitud predictiva?

multi_regression <- lm(Sales ~ Month, data = home)
summary(multi_regression)
## 
## Call:
## lm(formula = Sales ~ Month, data = home)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2089.63  -331.49    40.99   466.38  1536.91 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.432e+05  7.573e+03  -18.90   <2e-16 ***
## Month        9.807e+00  4.278e-01   22.92   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 811.1 on 34 degrees of freedom
## Multiple R-squared:  0.9392, Adjusted R-squared:  0.9374 
## F-statistic: 525.4 on 1 and 34 DF,  p-value: < 2.2e-16

De acuerdo con los resultados obtenidos en la creación de los modelos de regresión lineal, pudimos llegar a la conclusión de que el modelo que ofrece mejor exactitud predictiva es el modelo de regresión lineal múltiple. Para la regresión lineal abrimos el summary del modelo y pudimos observar que podia predecir el 93% de los casos de venta de Leche con transición. Mientras que para el otro modelo, abrimos el valor ar2, el cual es el equivalente para el modelo ARIMA y pudimos observar que regresaba un 60%,

3.-Según su mejor modelo ¿Cuál es la proyección de ventas en valor monetario para el siguiente a total de Hershey´s?

proyec_hersh <- forecast (modelo_arima, level = c(95), h = 5)
proyec_hersh
##          Point Forecast    Lo 95    Hi 95
## Jan 2020       35498.90 34616.48 36381.32
## Feb 2020       34202.17 33155.29 35249.05
## Mar 2020       36703.01 35596.10 37809.92
## Apr 2020       36271.90 35141.44 37402.36
## May 2020       37121.98 35982.07 38261.90

Se puede observar que la proyección de ventas en valor monetario para el siguiente año a total de Hershey’s son los mostrados en la tabla; por ejemplo, para enero del 2020 se proyecta un valor total de 35498.9 ventas.

4.- Considerando algunos imponderables establezca con los modelos construidos tres escenarios futuros ( proyecciones) A (escenario esperado) , B (escenario optimista) y C (escenario pesimista).

# Escenario Esperado
escen_esper <- forecast (modelo_arima, level = c(95), h = 5)
escen_esper
##          Point Forecast    Lo 95    Hi 95
## Jan 2020       35498.90 34616.48 36381.32
## Feb 2020       34202.17 33155.29 35249.05
## Mar 2020       36703.01 35596.10 37809.92
## Apr 2020       36271.90 35141.44 37402.36
## May 2020       37121.98 35982.07 38261.90
# Escenario Optimista
escen_opti <- escen_esper$upper
escen_opti
##           Jan      Feb      Mar      Apr      May
## 2020 36381.32 35249.05 37809.92 37402.36 38261.90
# Escenario Pesimista
escen_pesi <- escen_esper$lower
escen_pesi
##           Jan      Feb      Mar      Apr      May
## 2020 34616.48 33155.29 35596.10 35141.44 35982.07

Al generar el modelo de Auto-ARIMA, este por naturaleza está basado en probabilidades, no puede predecir el futuro. Sin embargo, nos ofrece un intervalo de confianza para los meses que estamos intentando predecir. Con esto nos referimos que hay un 95% de probabilidad de que el verdadero resultado de los meses que estamos intentando predecir se encontrará entre los valores mínimos y máximos que ofrece el modelo.

Estos valores son útiles para determinar los escenarios optimistas y pesimistas para una compañía, ya que es una estimación basada en datos con un 95% de confianza en que sucederán.”

5.- Con este análisis descriptivo, predictivo ¿Qué recomendaciones y medidas prescriptivas le puede dar a la compañía Hershey´s?

  1. Establecer un aumento de producción en los meses donde más ventas puede llegan a tener como lo son en Mayo 2020, considerar festividades que hay al rededor de estas ocasiones.

  2. Pensar en mercados y clientes específicos, por ejemplo temporada alta de regreso a clases en agosto, aumentar la producción y los esfuerzos de mercadotecnia para incentivar el consumo del producto en los jóvenes y niños que se conoce que son mercados meta clave considerando sus padres que son quienes pagan el producto.

  3. En meses de baja demanda, en vez de producir menos, puede que se hagan estrategias para el consumo y ser estrategas para incentivar su compra por ejemplo con exhibiciones, muestras en tienda, productos que se compran juntos o promociones.

LS0tDQp0aXRsZTogIkFjdGl2aWRhZCAyIg0KYXV0aG9yOiAiRXF1aXBvICM1OiBcbiBMdWPDrWEgRMOtYXogLSBBMDEyODM2NDggfCBZYW1pbGV0aCBVcmlhcyAtIEEwMDgyNjYyNyB8IEpvc8OpIEFudG9uaW8gR29uesOhbGV6IC0gQTAwODI0MTQyIg0KZGF0ZTogJzIwMjMtMDItMjInDQpvdXRwdXQ6IA0KICBodG1sX2RvY3VtZW50Og0KICAgIHRvYzogdHJ1ZQ0KICAgIHRvY19mbG9hdDogdHJ1ZQ0KICAgIGNvZGVfZG93bmxvYWQ6IHRydWUNCi0tLQ0KDQoNCg0KIyMgQ29udGV4dG8NCkhlcnNoZXnigJlzIE3DqXhpY28gZW1wcmVzYSBjaG9jb2xhdGVyYSBjb24gbcOhcyBkZSAxMDAgYcOxb3MgZGUgaGlzdG9yaWEgYSBuaXZlbCBnbG9iYWwgKGZ1bmRhY2nDs24gMTkwMykgeSBlc3RhYmxlY2lkYSBlbiBNw6l4aWNvIGRlc2RlIDE5NjksIGNvbW8gdW5hIGRlIGNvbmZpdGVyw61hIGVudHJlIEhlcnNoZXkgRm9vZCBDb3Jwb3JhdGlvbiB5IEFuZGVyc29uIENsYXl0b24gJiBDby4gUy5BLiBmb3JtYW5kbyBOYWNpb25hbCBkZSBEdWxjZXMgUy5BLiAgZGUgQy5WLiBlbiBlbCBEaXN0cml0byBGZWRlcmFsIHkgZGVzcHXDqXMgZGUgMTIgYcOxb3MgY2FtYmlhbmRvIHN1cyBpbnN0YWxhY2lvbmVzIGEgRWwgU2FsdG8sIEphbGlzY28gZSBpbmljaWFuZG8gb3BlcmFjaW9uZXMgZW4gZWwgbWVzIGRlIGZlYnJlcm8gZGUgMTk4MS4gIA0KDQpFbiAyMDIyIEhlcnNoZXnCtHMgTcOpeGljbyBzZSBjb25zb2xpZMOzIGNvbW8gdW5hIGRlIGxhcyBlbXByZXNhcyBjaG9jb2xhdGVyYXMgZGVsIHBhw61zLCB5IGNvbiBtw6FzIGZ1ZXJ6YSBlbiBsYSB2ZW50YSBkZSBzdXMgbGVjaGl0YXMgeSBoYSBzdW1hZG8gZXNmdWVyem9zIHBhcmEgcXVlIGEgcGVzYXIgZGUgbGEgcGFuZGVtaWEgZGUgbGEgQ292aWQtMTkgb2ZyZXpjYSBlbCBtZWpvciB5IGRpc3BvbmliaWxpZGFkIGRlIHByb2R1Y3RvcyBhIHN1cyBjbGllbnRlcy4NCg0KIyMgSW1wb3J0YXIgTGlicmVyw61hcw0KYGBge3IgaW5jbHVkZT1GQUxTRX0NCmxpYnJhcnkoSG1pc2MpDQpsaWJyYXJ5KGdncGxvdDIpDQpsaWJyYXJ5KGdsbW5ldCkNCmxpYnJhcnkob3Blbnhsc3gpDQpsaWJyYXJ5KHJlYWR4bCkNCmxpYnJhcnkoZ2dmb3J0aWZ5KQ0KbGlicmFyeSh0c2VyaWVzKQ0KbGlicmFyeShmb3JlY2FzdCkNCmBgYA0KDQojIyBJbXBvcnRhciBsYSBiYXNlIGRlIGRhdG9zDQpgYGB7cn0NCmhvbWUgPC0gcmVhZC54bHN4KCJDOi9Vc2Vycy9qb3NlYS9PbmVEcml2ZSAtIEluc3RpdHV0byBUZWNub2xvZ2ljbyB5IGRlIEVzdHVkaW9zIFN1cGVyaW9yZXMgZGUgTW9udGVycmV5L0VzY3JpdG9yaW8vVGVjbm9sw7NnaWNvIGRlIE1vbnRlcnJleS9MSVQvU2VwdGltbyBTZW1lc3RyZS9FY29ub21ldHJpYS9TaXRQcm9ibGVtYTEvaGVyc2hleXNTYWxlcy54bHN4IikNCmhvbWUkTW9udGggPC0gYXMuRGF0ZShob21lJE1vbnRoLCBvcmlnaW4gPSAiMTg5OS0xMi0zMCIpDQpkZXNjcmliZShob21lKQ0KYGBgDQoNCiMjIDEuLSBVdGlsaXphbmRvIG1vZGVsb3MgQVJJTUEgKEJveC1KZW5raW5zLCBBUk1BLCBTQVJJTUEpIHkgIGxvcyBkYXRvcyBoaXN0w7NyaWNvcyBkZSBsYXMgdmVudGFzIGRlIGxlY2hlIHNhYm9yaXphZGEgwr9DdcOhbCBlcyBlbCBtb2RlbG8gcXVlIG1lam9yIHNlIGFkYXB0YSBhIGxhIHNlcmllPw0KDQpgYGB7cn0NCiNNb2RlbG8gQVJJTUENCnRpbWVTZXJpZXMgPC0gdHMoZGF0YT1ob21lJFNhbGVzLCBzdGFydCA9IGMoMjAxNywxKSwgZW5kID0gYygyMDE5LDEyKSxmcmVxdWVuY3k9MTIpDQp0aW1lU2VyaWVzDQoNCm1vZGVsb19hcmltYSA8LSBhdXRvLmFyaW1hKHRpbWVTZXJpZXMsIEQgPSAxKQ0KcHJvbm9zdGljb19oZXJzaCA8LSBmb3JlY2FzdCAobW9kZWxvX2FyaW1hLCBsZXZlbCA9IGMoOTUpLCBoID0gNSkNCnByb25vc3RpY29faGVyc2gNCg0KcGxvdChwcm9ub3N0aWNvX2hlcnNoLCBtYWluID0gIlByb27Ds3N0aWNvIGRlIHZlbnRhIGRlIExlY2hlIEhlcnNoZXkncyIpDQpgYGANCg0KIyMgMi4tIMK/UXXDqSBtb2RlbG8gZGUgcmVncmVzacOzbiBvZnJlY2UgbWVqb3IgZXhhY3RpdHVkIHByZWRpY3RpdmE/DQoNCmBgYHtyfQ0KbXVsdGlfcmVncmVzc2lvbiA8LSBsbShTYWxlcyB+IE1vbnRoLCBkYXRhID0gaG9tZSkNCnN1bW1hcnkobXVsdGlfcmVncmVzc2lvbikNCmBgYA0KRGUgYWN1ZXJkbyBjb24gbG9zIHJlc3VsdGFkb3Mgb2J0ZW5pZG9zIGVuIGxhIGNyZWFjacOzbiBkZSBsb3MgbW9kZWxvcyBkZSByZWdyZXNpw7NuIGxpbmVhbCwgcHVkaW1vcyBsbGVnYXIgYSBsYSBjb25jbHVzacOzbiBkZSBxdWUgZWwgbW9kZWxvIHF1ZSBvZnJlY2UgbWVqb3IgZXhhY3RpdHVkIHByZWRpY3RpdmEgZXMgZWwgbW9kZWxvIGRlIHJlZ3Jlc2nDs24gbGluZWFsIG3Dumx0aXBsZS4NClBhcmEgbGEgcmVncmVzacOzbiBsaW5lYWwgYWJyaW1vcyBlbCBzdW1tYXJ5IGRlbCBtb2RlbG8geSBwdWRpbW9zIG9ic2VydmFyIHF1ZSBwb2RpYSBwcmVkZWNpciBlbCA5MyUgZGUgbG9zIGNhc29zIGRlIHZlbnRhIGRlIExlY2hlIGNvbiB0cmFuc2ljacOzbi4gTWllbnRyYXMgcXVlIHBhcmEgZWwgb3RybyBtb2RlbG8sIGFicmltb3MgZWwgdmFsb3IgYXIyLCBlbCBjdWFsIGVzIGVsIGVxdWl2YWxlbnRlIHBhcmEgZWwgbW9kZWxvIEFSSU1BIHkgcHVkaW1vcyBvYnNlcnZhciBxdWUgcmVncmVzYWJhIHVuIDYwJSwgDQoNCiMjIDMuLVNlZ8O6biBzdSBtZWpvciBtb2RlbG8gIMK/Q3XDoWwgZXMgbGEgcHJveWVjY2nDs24gZGUgdmVudGFzIGVuIHZhbG9yIG1vbmV0YXJpbyBwYXJhIGVsIHNpZ3VpZW50ZSBhIHRvdGFsIGRlIEhlcnNoZXnCtHM/DQpgYGB7cn0NCnByb3llY19oZXJzaCA8LSBmb3JlY2FzdCAobW9kZWxvX2FyaW1hLCBsZXZlbCA9IGMoOTUpLCBoID0gNSkNCnByb3llY19oZXJzaA0KYGBgDQpTZSBwdWVkZSBvYnNlcnZhciBxdWUgbGEgcHJveWVjY2nDs24gZGUgdmVudGFzIGVuIHZhbG9yIG1vbmV0YXJpbyBwYXJhIGVsIHNpZ3VpZW50ZSBhw7FvIGEgdG90YWwgZGUgSGVyc2hleSdzIHNvbiBsb3MgbW9zdHJhZG9zIGVuIGxhIHRhYmxhOyBwb3IgZWplbXBsbywgcGFyYSBlbmVybyBkZWwgMjAyMCBzZSBwcm95ZWN0YSB1biB2YWxvciB0b3RhbCBkZSAzNTQ5OC45IHZlbnRhcy4NCg0KIyMgNC4tIENvbnNpZGVyYW5kbyBhbGd1bm9zIGltcG9uZGVyYWJsZXMgZXN0YWJsZXpjYSBjb24gbG9zIG1vZGVsb3MgY29uc3RydWlkb3MgdHJlcyBlc2NlbmFyaW9zIGZ1dHVyb3MgKCBwcm95ZWNjaW9uZXMpICBBIChlc2NlbmFyaW8gZXNwZXJhZG8pICwgQiAoZXNjZW5hcmlvIG9wdGltaXN0YSkgIHkgQyAoZXNjZW5hcmlvIHBlc2ltaXN0YSkuDQpgYGB7cn0NCiMgRXNjZW5hcmlvIEVzcGVyYWRvDQplc2Nlbl9lc3BlciA8LSBmb3JlY2FzdCAobW9kZWxvX2FyaW1hLCBsZXZlbCA9IGMoOTUpLCBoID0gNSkNCmVzY2VuX2VzcGVyDQpgYGANCg0KYGBge3IgZWNobz1UUlVFfQ0KIyBFc2NlbmFyaW8gT3B0aW1pc3RhDQplc2Nlbl9vcHRpIDwtIGVzY2VuX2VzcGVyJHVwcGVyDQplc2Nlbl9vcHRpDQpgYGANCg0KYGBge3IgZWNobz1UUlVFfQ0KIyBFc2NlbmFyaW8gUGVzaW1pc3RhDQplc2Nlbl9wZXNpIDwtIGVzY2VuX2VzcGVyJGxvd2VyDQplc2Nlbl9wZXNpDQpgYGANCg0KQWwgZ2VuZXJhciBlbCBtb2RlbG8gZGUgQXV0by1BUklNQSwgZXN0ZSBwb3IgbmF0dXJhbGV6YSBlc3TDoSBiYXNhZG8gZW4gcHJvYmFiaWxpZGFkZXMsIG5vIHB1ZWRlIHByZWRlY2lyIGVsIGZ1dHVyby4gU2luIGVtYmFyZ28sIG5vcyBvZnJlY2UgdW4gaW50ZXJ2YWxvIGRlIGNvbmZpYW56YSBwYXJhIGxvcyBtZXNlcyBxdWUgZXN0YW1vcyBpbnRlbnRhbmRvIHByZWRlY2lyLiBDb24gZXN0byBub3MgcmVmZXJpbW9zIHF1ZSBoYXkgdW4gOTUlIGRlIHByb2JhYmlsaWRhZCBkZSBxdWUgZWwgdmVyZGFkZXJvIHJlc3VsdGFkbyBkZSBsb3MgbWVzZXMgcXVlIGVzdGFtb3MgaW50ZW50YW5kbyBwcmVkZWNpciBzZSBlbmNvbnRyYXLDoSBlbnRyZSBsb3MgdmFsb3JlcyBtw61uaW1vcyB5IG3DoXhpbW9zIHF1ZSBvZnJlY2UgZWwgbW9kZWxvLiANCg0KRXN0b3MgdmFsb3JlcyBzb24gw7p0aWxlcyBwYXJhIGRldGVybWluYXIgbG9zIGVzY2VuYXJpb3Mgb3B0aW1pc3RhcyB5IHBlc2ltaXN0YXMgcGFyYSB1bmEgY29tcGHDscOtYSwgeWEgcXVlIGVzIHVuYSBlc3RpbWFjacOzbiBiYXNhZGEgZW4gZGF0b3MgIGNvbiB1biA5NSUgZGUgY29uZmlhbnphIGVuIHF1ZSBzdWNlZGVyw6FuLiINCg0KIyMgNS4tIENvbiBlc3RlIGFuw6FsaXNpcyBkZXNjcmlwdGl2bywgIHByZWRpY3Rpdm8gwr9RdcOpIHJlY29tZW5kYWNpb25lcyB5IG1lZGlkYXMgcHJlc2NyaXB0aXZhcyBsZSBwdWVkZSBkYXIgYSBsYSBjb21wYcOxw61hIEhlcnNoZXnCtHM/IA0KDQoxLiBFc3RhYmxlY2VyIHVuIGF1bWVudG8gZGUgcHJvZHVjY2nDs24gZW4gbG9zIG1lc2VzIGRvbmRlIG3DoXMgdmVudGFzIHB1ZWRlIGxsZWdhbiBhIHRlbmVyIGNvbW8gbG8gc29uIGVuIE1heW8gMjAyMCwgY29uc2lkZXJhciBmZXN0aXZpZGFkZXMgcXVlIGhheSBhbCByZWRlZG9yIGRlIGVzdGFzIG9jYXNpb25lcy4NCg0KMi4gUGVuc2FyIGVuIG1lcmNhZG9zIHkgY2xpZW50ZXMgZXNwZWPDrWZpY29zLCBwb3IgZWplbXBsbyB0ZW1wb3JhZGEgYWx0YSBkZSByZWdyZXNvIGEgY2xhc2VzIGVuIGFnb3N0bywgYXVtZW50YXIgbGEgcHJvZHVjY2nDs24geSBsb3MgZXNmdWVyem9zIGRlIG1lcmNhZG90ZWNuaWEgcGFyYSBpbmNlbnRpdmFyIGVsIGNvbnN1bW8gZGVsIHByb2R1Y3RvIGVuIGxvcyBqw7N2ZW5lcyB5IG5pw7FvcyBxdWUgc2UgY29ub2NlIHF1ZSBzb24gbWVyY2Fkb3MgbWV0YSBjbGF2ZSBjb25zaWRlcmFuZG8gc3VzIHBhZHJlcyBxdWUgc29uIHF1aWVuZXMgcGFnYW4gZWwgcHJvZHVjdG8uDQoNCjMuIEVuIG1lc2VzIGRlIGJhamEgZGVtYW5kYSwgZW4gdmV6IGRlIHByb2R1Y2lyIG1lbm9zLCBwdWVkZSBxdWUgc2UgaGFnYW4gZXN0cmF0ZWdpYXMgcGFyYSBlbCBjb25zdW1vIHkgc2VyIGVzdHJhdGVnYXMgcGFyYSBpbmNlbnRpdmFyIHN1IGNvbXByYSBwb3IgZWplbXBsbyBjb24gZXhoaWJpY2lvbmVzLCBtdWVzdHJhcyBlbiB0aWVuZGEsIHByb2R1Y3RvcyBxdWUgc2UgY29tcHJhbiBqdW50b3MgbyBwcm9tb2Npb25lcy4=