Instalar paquetes y llamar librerías
#install.packages("tidyverse")
library("tidyverse")
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.4 ✔ readr 2.1.5
## ✔ forcats 1.0.0 ✔ stringr 1.5.1
## ✔ ggplot2 3.5.2 ✔ tibble 3.3.0
## ✔ lubridate 1.9.4 ✔ tidyr 1.3.1
## ✔ purrr 1.1.0
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(readr)
Importar la base de datos
df<- read_csv("C:/Users/Salvador/Downloads/walmart.csv")
## Rows: 6435 Columns: 8
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (1): Date
## dbl (7): Store, Weekly_Sales, Holiday_Flag, Temperature, Fuel_Price, CPI, Un...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
Comprender la base de datos
summary(df)
## Store Date Weekly_Sales Holiday_Flag
## Min. : 1 Length:6435 Min. : 209986 Min. :0.00000
## 1st Qu.:12 Class :character 1st Qu.: 553350 1st Qu.:0.00000
## Median :23 Mode :character Median : 960746 Median :0.00000
## Mean :23 Mean :1046965 Mean :0.06993
## 3rd Qu.:34 3rd Qu.:1420159 3rd Qu.:0.00000
## Max. :45 Max. :3818687 Max. :1.00000
## Temperature Fuel_Price CPI Unemployment
## Min. : -2.06 Min. :2.472 Min. :126.1 Min. : 3.879
## 1st Qu.: 47.46 1st Qu.:2.933 1st Qu.:131.7 1st Qu.: 6.891
## Median : 62.67 Median :3.445 Median :182.6 Median : 7.874
## Mean : 60.66 Mean :3.359 Mean :171.6 Mean : 7.999
## 3rd Qu.: 74.94 3rd Qu.:3.735 3rd Qu.:212.7 3rd Qu.: 8.622
## Max. :100.14 Max. :4.468 Max. :227.2 Max. :14.313
str(df)
## spc_tbl_ [6,435 × 8] (S3: spec_tbl_df/tbl_df/tbl/data.frame)
## $ Store : num [1:6435] 1 1 1 1 1 1 1 1 1 1 ...
## $ Date : chr [1:6435] "05-02-2010" "12-02-2010" "19-02-2010" "26-02-2010" ...
## $ Weekly_Sales: num [1:6435] 1643691 1641957 1611968 1409728 1554807 ...
## $ Holiday_Flag: num [1:6435] 0 1 0 0 0 0 0 0 0 0 ...
## $ Temperature : num [1:6435] 42.3 38.5 39.9 46.6 46.5 ...
## $ Fuel_Price : num [1:6435] 2.57 2.55 2.51 2.56 2.62 ...
## $ CPI : num [1:6435] 211 211 211 211 211 ...
## $ Unemployment: num [1:6435] 8.11 8.11 8.11 8.11 8.11 ...
## - attr(*, "spec")=
## .. cols(
## .. Store = col_double(),
## .. Date = col_character(),
## .. Weekly_Sales = col_double(),
## .. Holiday_Flag = col_double(),
## .. Temperature = col_double(),
## .. Fuel_Price = col_double(),
## .. CPI = col_double(),
## .. Unemployment = col_double()
## .. )
## - attr(*, "problems")=<externalptr>
df$Date<- as.Date(df$Date,format="%d-%m-%Y")
str(df)
## spc_tbl_ [6,435 × 8] (S3: spec_tbl_df/tbl_df/tbl/data.frame)
## $ Store : num [1:6435] 1 1 1 1 1 1 1 1 1 1 ...
## $ Date : Date[1:6435], format: "2010-02-05" "2010-02-12" ...
## $ Weekly_Sales: num [1:6435] 1643691 1641957 1611968 1409728 1554807 ...
## $ Holiday_Flag: num [1:6435] 0 1 0 0 0 0 0 0 0 0 ...
## $ Temperature : num [1:6435] 42.3 38.5 39.9 46.6 46.5 ...
## $ Fuel_Price : num [1:6435] 2.57 2.55 2.51 2.56 2.62 ...
## $ CPI : num [1:6435] 211 211 211 211 211 ...
## $ Unemployment: num [1:6435] 8.11 8.11 8.11 8.11 8.11 ...
## - attr(*, "spec")=
## .. cols(
## .. Store = col_double(),
## .. Date = col_character(),
## .. Weekly_Sales = col_double(),
## .. Holiday_Flag = col_double(),
## .. Temperature = col_double(),
## .. Fuel_Price = col_double(),
## .. CPI = col_double(),
## .. Unemployment = col_double()
## .. )
## - attr(*, "problems")=<externalptr>
Agregar variables a la base de datos
df$Year<-format(df$Date, "%Y")
df$Year<-as.integer(df$Year)
df$Month<-format(df$Date, "%m")
df$Month<-as.integer(df$Month)
df$WeekYear<-format(df$Date, "%W")
df$WeekYear<-as.integer(df$WeekYear)
df$WeekDay<-format(df$Date, "%w")
df$WeekDay<-as.integer(df$WeekDay)
df$Day<-format(df$Date, "%d")
df$Day<-as.integer(df$Day)
Regresión lineal
regresion<- lm(Weekly_Sales~., data=df)
summary(regresion)
##
## Call:
## lm(formula = Weekly_Sales ~ ., data = df)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1094800 -382464 -42860 375406 2587123
##
## Coefficients: (2 not defined because of singularities)
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -2.384e+09 9.127e+09 -0.261 0.7940
## Store -1.538e+04 5.202e+02 -29.576 < 2e-16 ***
## Date -3.399e+03 1.266e+04 -0.268 0.7883
## Holiday_Flag 4.773e+04 2.706e+04 1.763 0.0779 .
## Temperature -1.817e+03 4.053e+02 -4.484 7.47e-06 ***
## Fuel_Price 6.124e+04 2.876e+04 2.130 0.0332 *
## CPI -2.109e+03 1.928e+02 -10.941 < 2e-16 ***
## Unemployment -2.209e+04 3.967e+03 -5.569 2.67e-08 ***
## Year 1.212e+06 4.633e+06 0.262 0.7937
## Month 1.177e+05 3.858e+05 0.305 0.7604
## WeekYear NA NA NA NA
## WeekDay NA NA NA NA
## Day 2.171e+03 1.269e+04 0.171 0.8642
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 520900 on 6424 degrees of freedom
## Multiple R-squared: 0.1495, Adjusted R-squared: 0.1482
## F-statistic: 113 on 10 and 6424 DF, p-value: < 2.2e-16
Ajustar la regresión
df_ajustada <-df %>% select(-Fuel_Price,-Date, -Year:-Date)
regresion_ajustada <- lm(Weekly_Sales~., data=df)
summary(regresion_ajustada)
##
## Call:
## lm(formula = Weekly_Sales ~ ., data = df)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1094800 -382464 -42860 375406 2587123
##
## Coefficients: (2 not defined because of singularities)
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -2.384e+09 9.127e+09 -0.261 0.7940
## Store -1.538e+04 5.202e+02 -29.576 < 2e-16 ***
## Date -3.399e+03 1.266e+04 -0.268 0.7883
## Holiday_Flag 4.773e+04 2.706e+04 1.763 0.0779 .
## Temperature -1.817e+03 4.053e+02 -4.484 7.47e-06 ***
## Fuel_Price 6.124e+04 2.876e+04 2.130 0.0332 *
## CPI -2.109e+03 1.928e+02 -10.941 < 2e-16 ***
## Unemployment -2.209e+04 3.967e+03 -5.569 2.67e-08 ***
## Year 1.212e+06 4.633e+06 0.262 0.7937
## Month 1.177e+05 3.858e+05 0.305 0.7604
## WeekYear NA NA NA NA
## WeekDay NA NA NA NA
## Day 2.171e+03 1.269e+04 0.171 0.8642
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 520900 on 6424 degrees of freedom
## Multiple R-squared: 0.1495, Adjusted R-squared: 0.1482
## F-statistic: 113 on 10 and 6424 DF, p-value: < 2.2e-16
LS0tDQp0aXRsZTogIldBTE1BUlQiDQphdXRob3I6ICJTYWx2YWRvciBOYXJ2YWV6IEEwMDU3MTg0OCINCmRhdGU6ICIyMDI1LTA5LTAxIg0Kb3V0cHV0OiANCiAgaHRtbF9kb2N1bWVudDoNCiAgICB0b2M6IFRSVUUNCiAgICB0b2NfZmxvYXQ6IFRSVUUNCiAgICBjb2RlX2Rvd25sb2FkOiBUUlVFDQogICAgdGhlbWU6IGJvb3RzdHJhcA0KLS0tDQoNCjxjZW50ZXI+DQohW10oaHR0cHM6Ly9pLmdpZmVyLmNvbS8zaGFnLmdpZikNCjwvY2VudGVyPg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcsOtYXMgPC9zcGFuPg0KYGBge3J9DQojaW5zdGFsbC5wYWNrYWdlcygidGlkeXZlcnNlIikNCmxpYnJhcnkoInRpZHl2ZXJzZSIpDQpsaWJyYXJ5KHJlYWRyKQ0KYGBgDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gSW1wb3J0YXIgbGEgYmFzZSBkZSBkYXRvcyA8L3NwYW4+DQpgYGB7cn0NCmRmPC0gcmVhZF9jc3YoIkM6L1VzZXJzL1NhbHZhZG9yL0Rvd25sb2Fkcy93YWxtYXJ0LmNzdiIpDQpgYGANCg0KIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBDb21wcmVuZGVyIGxhIGJhc2UgZGUgZGF0b3MgPC9zcGFuPg0KYGBge3J9DQpzdW1tYXJ5KGRmKQ0Kc3RyKGRmKQ0KZGYkRGF0ZTwtIGFzLkRhdGUoZGYkRGF0ZSxmb3JtYXQ9IiVkLSVtLSVZIikNCnN0cihkZikNCmBgYA0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IEFncmVnYXIgdmFyaWFibGVzIGEgbGEgYmFzZSBkZSBkYXRvcyA8L3NwYW4+DQpgYGB7cn0NCmRmJFllYXI8LWZvcm1hdChkZiREYXRlLCAiJVkiKQ0KZGYkWWVhcjwtYXMuaW50ZWdlcihkZiRZZWFyKQ0KDQpkZiRNb250aDwtZm9ybWF0KGRmJERhdGUsICIlbSIpDQpkZiRNb250aDwtYXMuaW50ZWdlcihkZiRNb250aCkNCg0KZGYkV2Vla1llYXI8LWZvcm1hdChkZiREYXRlLCAiJVciKQ0KZGYkV2Vla1llYXI8LWFzLmludGVnZXIoZGYkV2Vla1llYXIpDQoNCmRmJFdlZWtEYXk8LWZvcm1hdChkZiREYXRlLCAiJXciKQ0KZGYkV2Vla0RheTwtYXMuaW50ZWdlcihkZiRXZWVrRGF5KQ0KDQpkZiREYXk8LWZvcm1hdChkZiREYXRlLCAiJWQiKQ0KZGYkRGF5PC1hcy5pbnRlZ2VyKGRmJERheSkNCmBgYA0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFJlZ3Jlc2nDs24gbGluZWFsIDwvc3Bhbj4NCmBgYHtyfQ0KcmVncmVzaW9uPC0gbG0oV2Vla2x5X1NhbGVzfi4sIGRhdGE9ZGYpDQpzdW1tYXJ5KHJlZ3Jlc2lvbikNCmBgYA0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IEFqdXN0YXIgbGEgcmVncmVzacOzbiA8L3NwYW4+DQpgYGB7cn0NCmRmX2FqdXN0YWRhIDwtZGYgJT4lIHNlbGVjdCgtRnVlbF9QcmljZSwtRGF0ZSwgLVllYXI6LURhdGUpIA0KcmVncmVzaW9uX2FqdXN0YWRhIDwtIGxtKFdlZWtseV9TYWxlc34uLCBkYXRhPWRmKQ0Kc3VtbWFyeShyZWdyZXNpb25fYWp1c3RhZGEpDQpgYGANCg0KDQoNCg0KDQo=