Importar base de datos

#file.choose ()
bd <- read.csv("/Users/vanessaelizondo/Desktop/Tec/Semestre 7/CSV/rentadebicis.csv")

Entender la base de datos

resumen <- summary(bd)
resumen
##       hora            dia              mes              aÒo      
##  Min.   : 0.00   Min.   : 1.000   Min.   : 1.000   Min.   :2011  
##  1st Qu.: 6.00   1st Qu.: 5.000   1st Qu.: 4.000   1st Qu.:2011  
##  Median :12.00   Median :10.000   Median : 7.000   Median :2012  
##  Mean   :11.54   Mean   : 9.993   Mean   : 6.521   Mean   :2012  
##  3rd Qu.:18.00   3rd Qu.:15.000   3rd Qu.:10.000   3rd Qu.:2012  
##  Max.   :23.00   Max.   :19.000   Max.   :12.000   Max.   :2012  
##     estacion     dia_de_la_semana     asueto         temperatura   
##  Min.   :1.000   Min.   :1.000    Min.   :0.00000   Min.   : 0.82  
##  1st Qu.:2.000   1st Qu.:2.000    1st Qu.:0.00000   1st Qu.:13.94  
##  Median :3.000   Median :4.000    Median :0.00000   Median :20.50  
##  Mean   :2.507   Mean   :4.014    Mean   :0.02857   Mean   :20.23  
##  3rd Qu.:4.000   3rd Qu.:6.000    3rd Qu.:0.00000   3rd Qu.:26.24  
##  Max.   :4.000   Max.   :7.000    Max.   :1.00000   Max.   :41.00  
##  sensacion_termica    humedad       velocidad_del_viento
##  Min.   : 0.76     Min.   :  0.00   Min.   : 0.000      
##  1st Qu.:16.66     1st Qu.: 47.00   1st Qu.: 7.002      
##  Median :24.24     Median : 62.00   Median :12.998      
##  Mean   :23.66     Mean   : 61.89   Mean   :12.799      
##  3rd Qu.:31.06     3rd Qu.: 77.00   3rd Qu.:16.998      
##  Max.   :45.45     Max.   :100.00   Max.   :56.997      
##  rentas_de_no_registrados rentas_de_registrados rentas_totales 
##  Min.   :  0.00           Min.   :  0.0         Min.   :  1.0  
##  1st Qu.:  4.00           1st Qu.: 36.0         1st Qu.: 42.0  
##  Median : 17.00           Median :118.0         Median :145.0  
##  Mean   : 36.02           Mean   :155.6         Mean   :191.6  
##  3rd Qu.: 49.00           3rd Qu.:222.0         3rd Qu.:284.0  
##  Max.   :367.00           Max.   :886.0         Max.   :977.0

Generador de Valores

Pasos Importantes:.

Paso 1. Definir el área de negocios que deseamos impactar o mejorar y su KPI.

Impactar en el area de ventas, mejorando sus KPIs de conformidad y atracción de clientes.

Paso 2. Seleccionar plantilla (-s) para crear valor a partir de los datos de los clientes.

Visión / Segmentación / Personalización /Contextualización

Paso 3. Generar Ideas o conceptos específicos

Elaborar un modelo de ventas predictivas

Paso 4. Reunir los datos requeridos

Crear una nueva base de datos con la información esencial y las variables significativas que tienen mayor impacto y valor en la renta de las bicis.

Paso 5. Plan de ejecución

1. Mercadotecnia hace un plan de atracción de los clientes segmentados tomando en cuenta las variables con mayor importancia en la renta de bicis.

2. El area de Ventas conoce la información de impacto para el cliente.

3. Se lanza mayor publicidad en días con características relevantes para los clientes y la publicidad es atractiva al mencionar aquello que se encuentra como relevante para el cliente.

Observaciones

¿Porqué los dias llegan hasta el 19 y no el 31? ¿Qué significan las estaciones? R: 1 es primavera, 2 es verano y 3 invierno

plot(bd$temperatura, bd$rentas_totales, main = "Influencia de la Temperatura sobre las Rentas Totales", xlab = "Temperatura (·c)", ylab = "Rentas Totales")

Generar Regresion Lineal

regresion <- lm (rentas_totales ~ hora + dia + mes + aÒo + estacion + dia_de_la_semana + asueto + temperatura + sensacion_termica + humedad + velocidad_del_viento, data=bd)
summary (regresion)
## 
## Call:
## lm(formula = rentas_totales ~ hora + dia + mes + aÒo + estacion + 
##     dia_de_la_semana + asueto + temperatura + sensacion_termica + 
##     humedad + velocidad_del_viento, data = bd)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -305.52  -93.64  -27.70   61.85  649.10 
## 
## Coefficients:
##                        Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          -1.661e+05  5.496e+03 -30.217  < 2e-16 ***
## hora                  7.735e+00  2.070e-01  37.368  < 2e-16 ***
## dia                   3.844e-01  2.482e-01   1.549  0.12150    
## mes                   9.996e+00  1.682e+00   5.943 2.89e-09 ***
## aÒo                   8.258e+01  2.732e+00  30.225  < 2e-16 ***
## estacion             -7.774e+00  5.177e+00  -1.502  0.13324    
## dia_de_la_semana      4.393e-01  6.918e-01   0.635  0.52545    
## asueto               -4.864e+00  8.365e+00  -0.582  0.56089    
## temperatura           1.582e+00  1.038e+00   1.524  0.12752    
## sensacion_termica     4.748e+00  9.552e-01   4.971 6.76e-07 ***
## humedad              -2.115e+00  7.884e-02 -26.827  < 2e-16 ***
## velocidad_del_viento  5.582e-01  1.809e-01   3.086  0.00203 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 141.7 on 10874 degrees of freedom
## Multiple R-squared:  0.3891, Adjusted R-squared:  0.3885 
## F-statistic: 629.6 on 11 and 10874 DF,  p-value: < 2.2e-16

Evaluar y en caso necesario ajustar la regresion lineal

regresion <- lm (rentas_totales ~ hora + mes + aÒo + sensacion_termica + humedad + velocidad_del_viento, data=bd)
summary (regresion)
## 
## Call:
## lm(formula = rentas_totales ~ hora + mes + aÒo + sensacion_termica + 
##     humedad + velocidad_del_viento, data = bd)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -308.60  -93.85  -28.34   61.05  648.09 
## 
## Coefficients:
##                        Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          -1.662e+05  5.496e+03 -30.250  < 2e-16 ***
## hora                  7.734e+00  2.070e-01  37.364  < 2e-16 ***
## mes                   7.574e+00  4.207e-01  18.002  < 2e-16 ***
## aÒo                   8.266e+01  2.732e+00  30.258  < 2e-16 ***
## sensacion_termica     6.172e+00  1.689e-01  36.539  < 2e-16 ***
## humedad              -2.121e+00  7.858e-02 -26.988  < 2e-16 ***
## velocidad_del_viento  6.208e-01  1.771e-01   3.506 0.000457 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 141.7 on 10879 degrees of freedom
## Multiple R-squared:  0.3886, Adjusted R-squared:  0.3883 
## F-statistic:  1153 on 6 and 10879 DF,  p-value: < 2.2e-16

Construir Modelo Predictivo

datos_nuevos <- data.frame (hora=12, mes=1, aÒo=2013, sensacion_termica=24, humedad=62, velocidad_del_viento=13)
predict(regresion,datos_nuevos)
##        1 
## 279.1478

Conclusión

Se hizo una regresión lineal para conocer la revelancia de las variables y la relación que tienen entre ellas.

En este caso se buscaba ver el Total de Rentas de Bicis y con el uso de la regresión lineal pudimos identificar variables con mayor impacto en la información y la variable que buscabamos analizar. Hicimos una comparación entre la influencia de la Temperatura y la Sensación Térmica y nos dimos cuenta que tiene mayor impacto la Sensación Térmica con 4.74. Otras variables con mucha relevancia fueron el mes y el año para definir las Rentas Totales las cuales indican tener un máximo de 977 Rentas.

También realizamos un modelo predictivo para predecir situaciones o riesgos a futuro y pudimos jugar con las variables y reconocer su importancia, al final las variables con mayor impacto fueron seleccionadas, las cuales fueron: la hora, mes, año, sensación térmica, humedad y velocidad del viento. Estás para poder conocer cantidades y totales de las rentas de bicis. Y con la predicción se obtuvo el dato de que en enero del 2013, a las 12pm, con una sensación térmica de 24, humedad de 62 y velocidad de 13, se obtendría una renta de bicis totales de 279.

También con este modelo predictivo pudimos definir que el modelo tiene un R_squared de .38, por lo que se ve que tiene mucho margen de oportunidad y de mejora.

LS0tCnRpdGxlOiA8c3BhbiBzdHlsZSA9ICJjb2xvcjpncmVlbiI+KipSZW50YSBkZSBCaWNpcyoqCmF1dGhvcjogIlZhbmVzc2EgRWxpem9uZG8gLSBBMDA4Mjc2NzIiCmRhdGU6ICIyMDIyLTA5LTA2IgpvdXRwdXQ6IAogIGh0bWxfZG9jdW1lbnQ6CiAgICB0aGVtZTogZGFya2x5CiAgICBoaWdobGlnaHQ6IHplbmJ1cm4KICAgIHRvYzogdHJ1ZQogICAgdG9jX2Zsb2F0OiB0cnVlCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCi0tLQo8aW1nIHNyYz0gIi9Vc2Vycy92YW5lc3NhZWxpem9uZG8vRGVza3RvcC9UZWMvU2VtZXN0cmUgNy9FbnRyZWdhIDMuMS9mb3Rvcy9iaWNpcy5qcGVnIj4KCiMjIEltcG9ydGFyIGJhc2UgZGUgZGF0b3MKCmBgYHtyfQojZmlsZS5jaG9vc2UgKCkKYmQgPC0gcmVhZC5jc3YoIi9Vc2Vycy92YW5lc3NhZWxpem9uZG8vRGVza3RvcC9UZWMvU2VtZXN0cmUgNy9DU1YvcmVudGFkZWJpY2lzLmNzdiIpCgpgYGAKCiMjIEVudGVuZGVyIGxhIGJhc2UgZGUgZGF0b3MKYGBge3J9CnJlc3VtZW4gPC0gc3VtbWFyeShiZCkKcmVzdW1lbgpgYGAKCiMjIEdlbmVyYWRvciBkZSBWYWxvcmVzCjxzcGFuIHN0eWxlID0gImNvbG9yOmJsdWUiPlBhc29zIEltcG9ydGFudGVzOjwvc3Bhbj4uCgoqKlBhc28gMS4qKiBEZWZpbmlyIGVsIMOhcmVhIGRlIG5lZ29jaW9zIHF1ZSBkZXNlYW1vcyBpbXBhY3RhciBvIG1lam9yYXIgeSBzdSBLUEkuCgoqSW1wYWN0YXIgZW4gZWwgYXJlYSBkZSB2ZW50YXMsIG1lam9yYW5kbyBzdXMgS1BJcyBkZSBjb25mb3JtaWRhZCB5IGF0cmFjY2nDs24gZGUgY2xpZW50ZXMuKiAKCioqUGFzbyAyLioqIFNlbGVjY2lvbmFyIHBsYW50aWxsYSAoLXMpIHBhcmEgY3JlYXIgdmFsb3IgYSBwYXJ0aXIgZGUgbG9zIGRhdG9zIGRlIGxvcyBjbGllbnRlcy4KCiAgKlZpc2nDs24gLyAqKlNlZ21lbnRhY2nDs24qKiAvIFBlcnNvbmFsaXphY2nDs24gL0NvbnRleHR1YWxpemFjacOzbioKICAKKipQYXNvIDMuKiogR2VuZXJhciBJZGVhcyBvIGNvbmNlcHRvcyBlc3BlY8OtZmljb3MKCipFbGFib3JhciB1biBtb2RlbG8gZGUgdmVudGFzIHByZWRpY3RpdmFzKgoKKipQYXNvIDQuKiogUmV1bmlyIGxvcyBkYXRvcyByZXF1ZXJpZG9zCgoqQ3JlYXIgdW5hIG51ZXZhIGJhc2UgZGUgZGF0b3MgY29uIGxhIGluZm9ybWFjacOzbiBlc2VuY2lhbCB5IGxhcyB2YXJpYWJsZXMgc2lnbmlmaWNhdGl2YXMgcXVlIHRpZW5lbiBtYXlvciBpbXBhY3RvIHkgdmFsb3IgZW4gbGEgcmVudGEgZGUgbGFzIGJpY2lzLiogCgoqKlBhc28gNS4qKiBQbGFuIGRlIGVqZWN1Y2nDs24gCgoqMS4gTWVyY2Fkb3RlY25pYSBoYWNlIHVuIHBsYW4gZGUgYXRyYWNjacOzbiBkZSBsb3MgY2xpZW50ZXMgc2VnbWVudGFkb3MgdG9tYW5kbyBlbiBjdWVudGEgbGFzIHZhcmlhYmxlcyBjb24gbWF5b3IgaW1wb3J0YW5jaWEgZW4gbGEgcmVudGEgZGUgYmljaXMuKgoKKjIuIEVsIGFyZWEgZGUgVmVudGFzIGNvbm9jZSBsYSBpbmZvcm1hY2nDs24gZGUgaW1wYWN0byBwYXJhIGVsIGNsaWVudGUuKgoKKjMuIFNlIGxhbnphIG1heW9yIHB1YmxpY2lkYWQgZW4gZMOtYXMgY29uIGNhcmFjdGVyw61zdGljYXMgcmVsZXZhbnRlcyBwYXJhIGxvcyBjbGllbnRlcyB5IGxhIHB1YmxpY2lkYWQgZXMgYXRyYWN0aXZhIGFsIG1lbmNpb25hciBhcXVlbGxvIHF1ZSBzZSBlbmN1ZW50cmEgY29tbyByZWxldmFudGUgcGFyYSBlbCBjbGllbnRlLioKCgojIyBPYnNlcnZhY2lvbmVzCsK/UG9ycXXDqSBsb3MgZGlhcyBsbGVnYW4gaGFzdGEgZWwgMTkgeSBubyBlbCAzMT8Kwr9RdcOpIHNpZ25pZmljYW4gbGFzIGVzdGFjaW9uZXM/IApSOiAxIGVzIHByaW1hdmVyYSwgMiBlcyB2ZXJhbm8geSAzIGludmllcm5vCmBgYHtyfQpwbG90KGJkJHRlbXBlcmF0dXJhLCBiZCRyZW50YXNfdG90YWxlcywgbWFpbiA9ICJJbmZsdWVuY2lhIGRlIGxhIFRlbXBlcmF0dXJhIHNvYnJlIGxhcyBSZW50YXMgVG90YWxlcyIsIHhsYWIgPSAiVGVtcGVyYXR1cmEgKMK3YykiLCB5bGFiID0gIlJlbnRhcyBUb3RhbGVzIikKYGBgCgojIyBHZW5lcmFyIFJlZ3Jlc2lvbiBMaW5lYWwKYGBge3J9CnJlZ3Jlc2lvbiA8LSBsbSAocmVudGFzX3RvdGFsZXMgfiBob3JhICsgZGlhICsgbWVzICsgYcOSbyArIGVzdGFjaW9uICsgZGlhX2RlX2xhX3NlbWFuYSArIGFzdWV0byArIHRlbXBlcmF0dXJhICsgc2Vuc2FjaW9uX3Rlcm1pY2EgKyBodW1lZGFkICsgdmVsb2NpZGFkX2RlbF92aWVudG8sIGRhdGE9YmQpCnN1bW1hcnkgKHJlZ3Jlc2lvbikKYGBgCgojIyMgRXZhbHVhciB5IGVuIGNhc28gbmVjZXNhcmlvIGFqdXN0YXIgbGEgcmVncmVzaW9uIGxpbmVhbApgYGB7cn0KcmVncmVzaW9uIDwtIGxtIChyZW50YXNfdG90YWxlcyB+IGhvcmEgKyBtZXMgKyBhw5JvICsgc2Vuc2FjaW9uX3Rlcm1pY2EgKyBodW1lZGFkICsgdmVsb2NpZGFkX2RlbF92aWVudG8sIGRhdGE9YmQpCnN1bW1hcnkgKHJlZ3Jlc2lvbikKYGBgCgojIyBDb25zdHJ1aXIgTW9kZWxvIFByZWRpY3Rpdm8KYGBge3J9CmRhdG9zX251ZXZvcyA8LSBkYXRhLmZyYW1lIChob3JhPTEyLCBtZXM9MSwgYcOSbz0yMDEzLCBzZW5zYWNpb25fdGVybWljYT0yNCwgaHVtZWRhZD02MiwgdmVsb2NpZGFkX2RlbF92aWVudG89MTMpCnByZWRpY3QocmVncmVzaW9uLGRhdG9zX251ZXZvcykKYGBgCgojIyAqKkNvbmNsdXNpw7NuKioKU2UgaGl6byB1bmEgKipyZWdyZXNpw7NuIGxpbmVhbCoqIHBhcmEgY29ub2NlciBsYSByZXZlbGFuY2lhIGRlIGxhcyB2YXJpYWJsZXMgeSBsYSByZWxhY2nDs24gcXVlIHRpZW5lbiBlbnRyZSBlbGxhcy4gCgpFbiBlc3RlIGNhc28gc2UgYnVzY2FiYSB2ZXIgZWwgKlRvdGFsIGRlIFJlbnRhcyBkZSBCaWNpcyogeSBjb24gZWwgdXNvIGRlIGxhIHJlZ3Jlc2nDs24gbGluZWFsIHB1ZGltb3MgaWRlbnRpZmljYXIgdmFyaWFibGVzIGNvbiBtYXlvciBpbXBhY3RvIGVuIGxhIGluZm9ybWFjacOzbiB5IGxhIHZhcmlhYmxlIHF1ZSBidXNjYWJhbW9zIGFuYWxpemFyLiBIaWNpbW9zIHVuYSBjb21wYXJhY2nDs24gZW50cmUgbGEgaW5mbHVlbmNpYSBkZSBsYSBUZW1wZXJhdHVyYSB5IGxhIFNlbnNhY2nDs24gVMOpcm1pY2EgeSBub3MgZGltb3MgY3VlbnRhIHF1ZSB0aWVuZSBtYXlvciBpbXBhY3RvIGxhICpTZW5zYWNpw7NuIFTDqXJtaWNhKiBjb24gNC43NC4gT3RyYXMgdmFyaWFibGVzIGNvbiBtdWNoYSByZWxldmFuY2lhIGZ1ZXJvbiBlbCBtZXMgeSBlbCBhw7FvIHBhcmEgZGVmaW5pciBsYXMgUmVudGFzIFRvdGFsZXMgbGFzIGN1YWxlcyBpbmRpY2FuIHRlbmVyIHVuIG3DoXhpbW8gZGUgOTc3IFJlbnRhcy4gIAoKVGFtYmnDqW4gcmVhbGl6YW1vcyB1biAqKm1vZGVsbyBwcmVkaWN0aXZvKiogcGFyYSBwcmVkZWNpciBzaXR1YWNpb25lcyBvIHJpZXNnb3MgYSBmdXR1cm8geSBwdWRpbW9zIGp1Z2FyIGNvbiBsYXMgdmFyaWFibGVzIHkgcmVjb25vY2VyIHN1IGltcG9ydGFuY2lhLCBhbCBmaW5hbCBsYXMgdmFyaWFibGVzIGNvbiBtYXlvciBpbXBhY3RvIGZ1ZXJvbiBzZWxlY2Npb25hZGFzLCBsYXMgY3VhbGVzIGZ1ZXJvbjogKmxhIGhvcmEsIG1lcywgYcOxbywgc2Vuc2FjacOzbiB0w6lybWljYSwgaHVtZWRhZCB5IHZlbG9jaWRhZCBkZWwgdmllbnRvKi4gRXN0w6FzIHBhcmEgcG9kZXIgY29ub2NlciBjYW50aWRhZGVzIHkgdG90YWxlcyBkZSBsYXMgcmVudGFzIGRlIGJpY2lzLiBZIGNvbiBsYSBwcmVkaWNjacOzbiBzZSBvYnR1dm8gZWwgZGF0byBkZSBxdWUgZW4gZW5lcm8gZGVsIDIwMTMsIGEgbGFzIDEycG0sIGNvbiB1bmEgc2Vuc2FjacOzbiB0w6lybWljYSBkZSAyNCwgaHVtZWRhZCBkZSA2MiB5IHZlbG9jaWRhZCBkZSAxMywgc2Ugb2J0ZW5kcsOtYSB1bmEgcmVudGEgZGUgYmljaXMgdG90YWxlcyBkZSAqKjI3OSoqLiAKClRhbWJpw6luIGNvbiBlc3RlIG1vZGVsbyBwcmVkaWN0aXZvIHB1ZGltb3MgZGVmaW5pciBxdWUgZWwgbW9kZWxvIHRpZW5lIHVuICpSX3NxdWFyZWQqIGRlIC4zOCwgcG9yIGxvIHF1ZSBzZSB2ZSBxdWUgdGllbmUgbXVjaG8gbWFyZ2VuIGRlIG9wb3J0dW5pZGFkIHkgZGUgbWVqb3JhLiA=