Importar la base de datos

bd<-read.csv("/Users/Karen/Downloads/rentadebicis_2.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

Observaciones

#1. Porque los dias llegan hasta el 19 y no hasta el 31?
#1 Que significan los numeros en las estaciones? R. 1 es primavera, 2 es verano, 3 otono, 4 invierno
plot(bd$temperatura,bd$rentas_totales,main="Influencia de la temperatura sobre las rentas totales",xlab="Temperatura",ylab="Cantidad")

Herramienta el generador de valor de datos

Paso 1. Definir el area del negocio que buscamos impactar o mejorar y su KPI

El departamento de operaciones en el indicador de rentas mensuales de bicicletas

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 especificos

Elaborar un modelo predicitivo de rentas mensuales de bicis

Paso 4. Reunir los datos requeridos

Elaborar una base de datos con la variable dependiente (rentas de bicis) y las variables independientes

Paso 5. Plan de ejecución

Operaciones elaborará el plan de abastecimiento de bicicletas por mes

Sistemas asegurará la captura del markdown en las bases de datos

Generar regresión 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
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:12, año=2013, sensacion_termica=24, humedad=62,velocidad_del_viento=13)
predict(regresion,datos_nuevos)
##        1        2        3        4        5        6        7        8 
## 279.1478 286.7215 294.2952 301.8690 309.4427 317.0164 324.5901 332.1638 
##        9       10       11       12 
## 339.7375 347.3112 354.8849 362.4587

Conclusiones

A través de esta herramienta, creamos una regresión a partir de un conjunto de datos, que mostrara la influencia de diversos factores sobre las rentas mensuales de bicicletas. Al evaluar las regresiones lineales realizadas, se demostró que las rentas de bicis dependen en gran medida de la hora, el mes, la sensación térmica, la humedad, y la velocidad del viento.

Siendo de esta forma, se construyó un modelo predictivo que permitiera visualizar la renta de bicis a lo largo de un año, en este caso, 2013. Y con los resultados arrojados, el departamento de operaciones podrá elaborar el plan de abastecimiento de bicicletas por mes.

El programa R ayudó a dar una predicción de resultados futuros con base en datos históricos, para precisar un panorama de toma de decisiones.

LS0tCnRpdGxlOiA8c3BhbiBzdHlsZT0iQ29sb3I6SG90UGluayI+IFJlbnRhIGRlIGJpY2lzCmF1dGhvcjogIkthcmVuIE5lcmkgLSBBMDE2NTc5NzQiCmRhdGU6ICIyMDIyLTA5LTA2IgpvdXRwdXQ6IAogaHRtbF9kb2N1bWVudDoKICAgIHRvYzogdHJ1ZQogICAgdG9jX2Zsb2F0OiB0cnVlCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlIAotLS0KIVtdKGh0dHBzOi8vd3d3Lm1vdW50YWluYmlrZS5lcy91cGxvYWRzL3MxLzMyLzczLzg1LzkvNWVjZmNkNGIwZGU2OTQwYjI1MzQ5M2Y1LTEwLWJpY2lzLWNhbWJpYXJvbi1oaXN0b3JpYS1tb3VudGFpbi1iaWtlLmpwZWcpe3dpZHRoPSc2MDBweCd9CgojIyA8c3BhbiBzdHlsZT0iQ29sb3I6SG90UGluayI+IEltcG9ydGFyIGxhIGJhc2UgZGUgZGF0b3MKYGBge3J9CmJkPC1yZWFkLmNzdigiL1VzZXJzL0thcmVuL0Rvd25sb2Fkcy9yZW50YWRlYmljaXNfMi5jc3YiKQpgYGAKCiMjIDxzcGFuIHN0eWxlPSJDb2xvcjpIb3RQaW5rIj4gRW50ZW5kZXIgbGEgYmFzZSBkZSBkYXRvcyAKYGBge3J9CnJlc3VtZW4gPC0gc3VtbWFyeShiZCkKcmVzdW1lbgpgYGAKCiMjIDxzcGFuIHN0eWxlPSJDb2xvcjpIb3RQaW5rIj4gT2JzZXJ2YWNpb25lcwpgYGB7cn0KIzEuIFBvcnF1ZSBsb3MgZGlhcyBsbGVnYW4gaGFzdGEgZWwgMTkgeSBubyBoYXN0YSBlbCAzMT8KIzEgUXVlIHNpZ25pZmljYW4gbG9zIG51bWVyb3MgZW4gbGFzIGVzdGFjaW9uZXM/IFIuIDEgZXMgcHJpbWF2ZXJhLCAyIGVzIHZlcmFubywgMyBvdG9ubywgNCBpbnZpZXJubwpwbG90KGJkJHRlbXBlcmF0dXJhLGJkJHJlbnRhc190b3RhbGVzLG1haW49IkluZmx1ZW5jaWEgZGUgbGEgdGVtcGVyYXR1cmEgc29icmUgbGFzIHJlbnRhcyB0b3RhbGVzIix4bGFiPSJUZW1wZXJhdHVyYSIseWxhYj0iQ2FudGlkYWQiKQpgYGAKCiMjIDxzcGFuIHN0eWxlPSJDb2xvcjpIb3RQaW5rIj4gSGVycmFtaWVudGEgZWwgZ2VuZXJhZG9yIGRlIHZhbG9yIGRlIGRhdG9zIAoKKipQYXNvIDEuIERlZmluaXIgZWwgYXJlYSBkZWwgbmVnb2NpbyBxdWUgYnVzY2Ftb3MgaW1wYWN0YXIgbyBtZWpvcmFyIHkgc3UgS1BJKioKCkVsIGRlcGFydGFtZW50byBkZSBvcGVyYWNpb25lcyBlbiBlbCBpbmRpY2Fkb3IgZGUgcmVudGFzIG1lbnN1YWxlcyBkZSBiaWNpY2xldGFzCgoqKlBhc28gMi4gU2VsZWNjaW9uYXIgcGxhbnRpbGxhICgtcykgcGFyYSBjcmVhciB2YWxvciBhIHBhcnRpciBkZSBsb3MgZGF0b3MgZGUgbG9zIGNsaWVudGVzKioKCioqVmlzacOzbioqIC8gU2VnbWVudGFjacOzbiAvIFBlcnNvbmFsaXphY2nDs24gLyBDb250ZXh0dWFsaXphY2nDs24gCgoqKlBhc28gMy4gR2VuZXJhciBpZGVhcyBvIGNvbmNlcHRvcyBlc3BlY2lmaWNvcyoqCgpFbGFib3JhciB1biBtb2RlbG8gcHJlZGljaXRpdm8gZGUgcmVudGFzIG1lbnN1YWxlcyBkZSBiaWNpcwoKKipQYXNvIDQuIFJldW5pciBsb3MgZGF0b3MgcmVxdWVyaWRvcyAqKgoKRWxhYm9yYXIgdW5hIGJhc2UgZGUgZGF0b3MgY29uIGxhIHZhcmlhYmxlIGRlcGVuZGllbnRlIChyZW50YXMgZGUgYmljaXMpIHkgbGFzIHZhcmlhYmxlcyBpbmRlcGVuZGllbnRlcyAKCioqUGFzbyA1LiBQbGFuIGRlIGVqZWN1Y2nDs24gKioKCk9wZXJhY2lvbmVzIGVsYWJvcmFyw6EgZWwgcGxhbiBkZSBhYmFzdGVjaW1pZW50byBkZSBiaWNpY2xldGFzIHBvciBtZXMKClNpc3RlbWFzIGFzZWd1cmFyw6EgbGEgY2FwdHVyYSBkZWwgbWFya2Rvd24gZW4gbGFzIGJhc2VzIGRlIGRhdG9zIAoKCiMjIDxzcGFuIHN0eWxlPSJDb2xvcjpIb3RQaW5rIj4gR2VuZXJhciByZWdyZXNpw7NuIGxpbmVhbCAKYGBge3J9CnJlZ3Jlc2lvbiA8LSBsbShyZW50YXNfdG90YWxlc35ob3JhICsgZGlhICsgbWVzICsgYcOxbyArIGVzdGFjaW9uICsgZGlhX2RlX2xhX3NlbWFuYSArIGFzdWV0byArIHRlbXBlcmF0dXJhICsgc2Vuc2FjaW9uX3Rlcm1pY2EgKyBodW1lZGFkICsgdmVsb2NpZGFkX2RlbF92aWVudG8sZGF0YT1iZCkKc3VtbWFyeSAocmVncmVzaW9uKQpyZWdyZXNpb24gPC0gbG0ocmVudGFzX3RvdGFsZXN+aG9yYSArIG1lcyArIGHDsW8gKyBzZW5zYWNpb25fdGVybWljYSArIGh1bWVkYWQgKyB2ZWxvY2lkYWRfZGVsX3ZpZW50byxkYXRhPWJkKQpzdW1tYXJ5IChyZWdyZXNpb24pCmBgYAoKIyMgPHNwYW4gc3R5bGU9IkNvbG9yOkhvdFBpbmsiPiBDb25zdHJ1aXIgbW9kZWxvIHByZWRpY3Rpdm8KYGBge3J9CmRhdG9zX251ZXZvcyA8LSBkYXRhLmZyYW1lKGhvcmE9MTIsIG1lcz0xOjEyLCBhw7FvPTIwMTMsIHNlbnNhY2lvbl90ZXJtaWNhPTI0LCBodW1lZGFkPTYyLHZlbG9jaWRhZF9kZWxfdmllbnRvPTEzKQpwcmVkaWN0KHJlZ3Jlc2lvbixkYXRvc19udWV2b3MpCmBgYAoKIyMgPHNwYW4gc3R5bGU9IkNvbG9yOkhvdFBpbmsiPiBDb25jbHVzaW9uZXMKQSB0cmF2w6lzIGRlIGVzdGEgaGVycmFtaWVudGEsIGNyZWFtb3MgdW5hIHJlZ3Jlc2nDs24gYSBwYXJ0aXIgZGUgdW4gY29uanVudG8gZGUgZGF0b3MsIHF1ZSBtb3N0cmFyYSBsYSBpbmZsdWVuY2lhIGRlIGRpdmVyc29zIGZhY3RvcmVzIHNvYnJlIGxhcyByZW50YXMgbWVuc3VhbGVzIGRlIGJpY2ljbGV0YXMuIEFsIGV2YWx1YXIgbGFzIHJlZ3Jlc2lvbmVzIGxpbmVhbGVzIHJlYWxpemFkYXMsIHNlIGRlbW9zdHLDsyBxdWUgbGFzIHJlbnRhcyBkZSBiaWNpcyBkZXBlbmRlbiBlbiBncmFuIG1lZGlkYSBkZSBsYSAqKmhvcmEsIGVsIG1lcywgbGEgc2Vuc2FjacOzbiB0w6lybWljYSwgbGEgaHVtZWRhZCwgeSBsYSB2ZWxvY2lkYWQgZGVsIHZpZW50byoqLiAKClNpZW5kbyBkZSBlc3RhIGZvcm1hLCBzZSBjb25zdHJ1ecOzIHVuIG1vZGVsbyBwcmVkaWN0aXZvIHF1ZSBwZXJtaXRpZXJhIHZpc3VhbGl6YXIgbGEgcmVudGEgZGUgYmljaXMgYSBsbyBsYXJnbyBkZSB1biBhw7FvLCBlbiBlc3RlIGNhc28sIDIwMTMuIFkgY29uIGxvcyByZXN1bHRhZG9zIGFycm9qYWRvcywgZWwgZGVwYXJ0YW1lbnRvIGRlIG9wZXJhY2lvbmVzIHBvZHLDoSBlbGFib3JhciBlbCAqKnBsYW4gZGUgYWJhc3RlY2ltaWVudG8gZGUgYmljaWNsZXRhcyBwb3IgbWVzKiouIAoKRWwgcHJvZ3JhbWEgUiBheXVkw7MgYSBkYXIgdW5hIHByZWRpY2Npw7NuIGRlIHJlc3VsdGFkb3MgZnV0dXJvcyBjb24gYmFzZSBlbiBkYXRvcyBoaXN0w7NyaWNvcywgcGFyYSBwcmVjaXNhciB1biBwYW5vcmFtYSBkZSB0b21hIGRlIGRlY2lzaW9uZXMuIAo=