Ejercicio 1 Modulo 2. BootStrap

Instalar paquetes y llamar Librerias

# install.packages("boot")
library(boot)
# Datos originales de la base de datos Hotel Mirage
ocupacion <- c(78, 82, 75, 80, 85, 88, 92, 90, 86, 84, 80, 83)
cat("Datos:", ocupacion, "\n\n")
## Datos: 78 82 75 80 85 88 92 90 86 84 80 83
#Función para calcular la media (para el bootstrap)
bootstrap_media <- function(data, indices) {
  return(mean(data[indices]))  
}

Preguntas

Ejercicio:

1-¿Cuál es el rango en el que podría esperarse que esté la ocupación hotelera para el siguiente mes

# Pregunta 1:
#Aplicar Bootstrap con 1000 repeticiones
set.seed(616)
b1 <- boot(data = ocupacion, statistic = bootstrap_media, R = 1000)
cat("Media:", round(mean(b1$t), 2), "% | 95% CI: [", round(boot.ci(b1, type = "perc")$percent[4], 2), ", ", round(boot.ci(b1, type = "perc")$percent[5], 2), "]\n\n")
## Media: 83.65 % | 95% CI: [ 81 ,  86.42 ]
# Respuesta: Media: 83.65 % | 95% CI: [ 81 ,  86.42 ]

2-¿Cómo cambia el intervalo si en vez de 1000 repeticiones hacemos 100,000?

# Pregunta 2: 
#Aplicar Bootstrap con Cambio en el intervalo 1,000 a 100,000. 
set.seed(616)
b2 <- boot(data = ocupacion, statistic = bootstrap_media, R = 100000)
cat("Media:", round(mean(b2$t), 2), "% | 95% CI: [", round(boot.ci(b2, type = "perc")$percent[4], 2), ", ", round(boot.ci(b2, type = "perc")$percent[5], 2), "]\n")
## Media: 83.58 % | 95% CI: [ 80.92 ,  86.25 ]
cat("Rango 100,000:", round(diff(range(b2$t)), 2), "%\n\n")
## Rango 100,000: 10.92 %
# Respuesta: Media: 83.58 % | 95% CI: [ 80.92 ,  86.25 ]
# Nuevo Rango 100,000: 10.92 %

3-El dueño del hotel se dio cuenta de que la ocupación hotelera para el mes de diciembre no fue de 83%, si no de 73%, ¿Cómo cambian los resultados (volver a hacer el ejercicio)?

# Pregunta 3: actualizacion de ultimo dato, de 83% a 73%
ocupacion2 <- ocupacion
ocupacion2[12] <- 73 # Aqui estamos modificando el numero o digito del apartado/mes de dicimebre= 12
b3 <- boot(data = ocupacion2, statistic = bootstrap_media, R = 1000)
cat("Media corregida:", round(mean(b3$t), 2), "% | 95% CI: [", round(boot.ci(b3, type = "perc")$percent[4], 2), ", ", round(boot.ci(b3, type = "perc")$percent[5], 2), "]\n")
## Media corregida: 82.78 % | 95% CI: [ 79.75 ,  86.08 ]
cat("Diferencia media:", round(mean(b1$t) - mean(b3$t), 2), "%\n\n")
## Diferencia media: 0.87 %

4-Si la ocupación cae por debajo del 80%, ¿qué acciones recomendarían para evitar que el hotel tenga pérdidas económicas?

Como futuro Licenciado en Inteligencia de Negocios, recomiendo que si la ocupación cae por debajo del 80%, se actúe en estas cuatro áreas:

1-En marketing, ofrecer paquetes especiales y descuentos atractivos.

2-En costos, reducir personal temporal y optimizar el consumo energético.

3-En servicios, organizar eventos corporativos y generar alianzas estratégicas.

4-En competencia, analizar precios del mercado y enfocarse en nichos desatendidos.

LS0tCnRpdGxlOiAiRWplcmNpY2lvIDEgTW9kdWxvIDIiCmF1dGhvcjogIlNlYmFzdGlhbiBFc3Bpbm96YSBBMDA4MzM3MDQiCmRhdGU6ICIyMDI1LTA4LTE4IgpvdXRwdXQ6IAogaHRtbF9kb2N1bWVudDogCiAgdG9jOiBUUlVFCiAgdG9jX2Zsb2F0OiBUUlVFCiAgY29kZV9kb3dubG9hZDogVFJVRQogIHRoZW1lOiB5ZXRpCi0tLQoKIVtdKGh0dHBzOi8vNjQubWVkaWEudHVtYmxyLmNvbS8zYmJlOGIzYzkwNWFiMjM5ZTIxNDJiYzc1NDQyYTQ2Zi8xMzQ1ZWNlY2ViZDA3NjAzLTQ3L3M1NDB4ODEwLzI1ZGZhYjExNzg5NTVkNmYxYmJkODczNTUzYWUzMGJkZmEwNjRiZDIuZ2lmdikKCgojIDxzcGFuIHN0eWxlPSdjb2xvcjpyZWQ7Jz4gIEVqZXJjaWNpbyAxIE1vZHVsbyAyLiBCb290U3RyYXAgIDwvc3Bhbj4KCgojIyA8c3BhbiBzdHlsZT0nY29sb3I6cmVkOyc+ICBJbnN0YWxhciBwYXF1ZXRlcyB5IGxsYW1hciBMaWJyZXJpYXMgIDwvc3Bhbj4KYGBge3J9CiMgaW5zdGFsbC5wYWNrYWdlcygiYm9vdCIpCmxpYnJhcnkoYm9vdCkKYGBgCgoKYGBge3J9CiMgRGF0b3Mgb3JpZ2luYWxlcyBkZSBsYSBiYXNlIGRlIGRhdG9zIEhvdGVsIE1pcmFnZQpvY3VwYWNpb24gPC0gYyg3OCwgODIsIDc1LCA4MCwgODUsIDg4LCA5MiwgOTAsIDg2LCA4NCwgODAsIDgzKQpjYXQoIkRhdG9zOiIsIG9jdXBhY2lvbiwgIlxuXG4iKQoKI0Z1bmNpw7NuIHBhcmEgY2FsY3VsYXIgbGEgbWVkaWEgKHBhcmEgZWwgYm9vdHN0cmFwKQpib290c3RyYXBfbWVkaWEgPC0gZnVuY3Rpb24oZGF0YSwgaW5kaWNlcykgewogIHJldHVybihtZWFuKGRhdGFbaW5kaWNlc10pKSAgCn0KCmBgYAoKIyA8c3BhbiBzdHlsZT0nY29sb3I6Z3JlZW47Jz4gIFByZWd1bnRhcyAgPC9zcGFuPgpFamVyY2ljaW86CgoxLcK/Q3XDoWwgZXMgZWwgcmFuZ28gZW4gZWwgcXVlIHBvZHLDrWEgZXNwZXJhcnNlIHF1ZSBlc3TDqSBsYSBvY3VwYWNpw7NuIGhvdGVsZXJhCnBhcmEgZWwgc2lndWllbnRlIG1lcwoKYGBge3J9CiMgUHJlZ3VudGEgMToKI0FwbGljYXIgQm9vdHN0cmFwIGNvbiAxMDAwIHJlcGV0aWNpb25lcwpzZXQuc2VlZCg2MTYpCmIxIDwtIGJvb3QoZGF0YSA9IG9jdXBhY2lvbiwgc3RhdGlzdGljID0gYm9vdHN0cmFwX21lZGlhLCBSID0gMTAwMCkKY2F0KCJNZWRpYToiLCByb3VuZChtZWFuKGIxJHQpLCAyKSwgIiUgfCA5NSUgQ0k6IFsiLCByb3VuZChib290LmNpKGIxLCB0eXBlID0gInBlcmMiKSRwZXJjZW50WzRdLCAyKSwgIiwgIiwgcm91bmQoYm9vdC5jaShiMSwgdHlwZSA9ICJwZXJjIikkcGVyY2VudFs1XSwgMiksICJdXG5cbiIpCgojIFJlc3B1ZXN0YTogTWVkaWE6IDgzLjY1ICUgfCA5NSUgQ0k6IFsgODEgLCAgODYuNDIgXQoKYGBgCgoyLcK/Q8OzbW8gY2FtYmlhIGVsIGludGVydmFsbyBzaSBlbiB2ZXogZGUgMTAwMCByZXBldGljaW9uZXMgaGFjZW1vcyAxMDAsMDAwPwpgYGB7cn0KIyBQcmVndW50YSAyOiAKI0FwbGljYXIgQm9vdHN0cmFwIGNvbiBDYW1iaW8gZW4gZWwgaW50ZXJ2YWxvIDEsMDAwIGEgMTAwLDAwMC4gCnNldC5zZWVkKDYxNikKYjIgPC0gYm9vdChkYXRhID0gb2N1cGFjaW9uLCBzdGF0aXN0aWMgPSBib290c3RyYXBfbWVkaWEsIFIgPSAxMDAwMDApCmNhdCgiTWVkaWE6Iiwgcm91bmQobWVhbihiMiR0KSwgMiksICIlIHwgOTUlIENJOiBbIiwgcm91bmQoYm9vdC5jaShiMiwgdHlwZSA9ICJwZXJjIikkcGVyY2VudFs0XSwgMiksICIsICIsIHJvdW5kKGJvb3QuY2koYjIsIHR5cGUgPSAicGVyYyIpJHBlcmNlbnRbNV0sIDIpLCAiXVxuIikKY2F0KCJSYW5nbyAxMDAsMDAwOiIsIHJvdW5kKGRpZmYocmFuZ2UoYjIkdCkpLCAyKSwgIiVcblxuIikKCiMgUmVzcHVlc3RhOiBNZWRpYTogODMuNTggJSB8IDk1JSBDSTogWyA4MC45MiAsICA4Ni4yNSBdCiMgTnVldm8gUmFuZ28gMTAwLDAwMDogMTAuOTIgJQoKYGBgCgozLUVsIGR1ZcOxbyBkZWwgaG90ZWwgc2UgZGlvIGN1ZW50YSBkZSBxdWUgbGEgb2N1cGFjacOzbiBob3RlbGVyYSBwYXJhIGVsIG1lcyBkZQpkaWNpZW1icmUgbm8gZnVlIGRlIDgzJSwgc2kgbm8gZGUgNzMlLCDCv0PDs21vIGNhbWJpYW4gbG9zIHJlc3VsdGFkb3MgKHZvbHZlcgphIGhhY2VyIGVsIGVqZXJjaWNpbyk/CgpgYGB7cn0KIyBQcmVndW50YSAzOiBhY3R1YWxpemFjaW9uIGRlIHVsdGltbyBkYXRvLCBkZSA4MyUgYSA3MyUKb2N1cGFjaW9uMiA8LSBvY3VwYWNpb24Kb2N1cGFjaW9uMlsxMl0gPC0gNzMgIyBBcXVpIGVzdGFtb3MgbW9kaWZpY2FuZG8gZWwgbnVtZXJvIG8gZGlnaXRvIGRlbCBhcGFydGFkby9tZXMgZGUgZGljaW1lYnJlPSAxMgpiMyA8LSBib290KGRhdGEgPSBvY3VwYWNpb24yLCBzdGF0aXN0aWMgPSBib290c3RyYXBfbWVkaWEsIFIgPSAxMDAwKQpjYXQoIk1lZGlhIGNvcnJlZ2lkYToiLCByb3VuZChtZWFuKGIzJHQpLCAyKSwgIiUgfCA5NSUgQ0k6IFsiLCByb3VuZChib290LmNpKGIzLCB0eXBlID0gInBlcmMiKSRwZXJjZW50WzRdLCAyKSwgIiwgIiwgcm91bmQoYm9vdC5jaShiMywgdHlwZSA9ICJwZXJjIikkcGVyY2VudFs1XSwgMiksICJdXG4iKQpjYXQoIkRpZmVyZW5jaWEgbWVkaWE6Iiwgcm91bmQobWVhbihiMSR0KSAtIG1lYW4oYjMkdCksIDIpLCAiJVxuXG4iKQpgYGAKCgo0LVNpIGxhIG9jdXBhY2nDs24gY2FlIHBvciBkZWJham8gZGVsIDgwJSwgwr9xdcOpIGFjY2lvbmVzIHJlY29tZW5kYXLDrWFuCnBhcmEgZXZpdGFyIHF1ZSBlbCBob3RlbCB0ZW5nYSBww6lyZGlkYXMgZWNvbsOzbWljYXM/CgpDb21vIGZ1dHVybyBMaWNlbmNpYWRvIGVuIEludGVsaWdlbmNpYSBkZSBOZWdvY2lvcywgcmVjb21pZW5kbyBxdWUgc2kgbGEgb2N1cGFjacOzbiBjYWUgcG9yIGRlYmFqbyBkZWwgODAlLCBzZSBhY3TDumUgZW4gZXN0YXMgY3VhdHJvIMOhcmVhczoKCjEtKipFbiBtYXJrZXRpbmcqKiwgb2ZyZWNlciBwYXF1ZXRlcyBlc3BlY2lhbGVzIHkgZGVzY3VlbnRvcyBhdHJhY3Rpdm9zLgoKMi0qKkVuIGNvc3RvcyoqLCByZWR1Y2lyIHBlcnNvbmFsIHRlbXBvcmFsIHkgb3B0aW1pemFyIGVsIGNvbnN1bW8gZW5lcmfDqXRpY28uCgozLSoqRW4gc2VydmljaW9zKiosIG9yZ2FuaXphciBldmVudG9zIGNvcnBvcmF0aXZvcyB5IGdlbmVyYXIgYWxpYW56YXMgZXN0cmF0w6lnaWNhcy4KCjQtKipFbiBjb21wZXRlbmNpYSoqLCBhbmFsaXphciBwcmVjaW9zIGRlbCBtZXJjYWRvIHkgZW5mb2NhcnNlIGVuIG5pY2hvcyBkZXNhdGVuZGlkb3MuCgoKIVtdKGh0dHBzOi8vMjUubWVkaWEudHVtYmxyLmNvbS82ZTg1NTk4OTRjNTRmYzE0MWFiZGM0MjkwZjA5OTY3ZS90dW1ibHJfbWhmMTB4RmpyZDFya2hweW9vMV80MDAuZ2lmKQoKIVtdKGh0dHBzOi8vbWVkaWEudGVub3IuY29tL1RHS2l1SWtsZExFQUFBQU0vdGlwLWhhdC1qYWNrLXNwYXJyb3cuZ2lmKQo=