library(pacman)
p_load(rmdformats,readr,readxl,ggplot2,plotly,DT,xfun,gridExtra,leaflet,TSstudio, dplyr, psych, GGally, corrplot)#Regresión lineal múltiple
##Importar datos
#Datos de concentración y movilidad
aire <- read_excel("Concentracion.xlsx")
CM <- read_excel("Concentracion_Mov.xlsx")
datatable(aire)##Análisis de relación entre variables
datos <- as.data.frame(aire)
round(cor(x = datos, method="pearson"), 3)## PM10 Reactivacion_Comercial Supermercado_Farmacia
## PM10 1.000 0.281 0.280
## Reactivacion_Comercial 0.281 1.000 0.875
## Supermercado_Farmacia 0.280 0.875 1.000
## Parques_Centros 0.065 0.814 0.655
## Estaciones_Transito 0.270 0.944 0.872
## Lugares_Trabajo 0.100 0.581 0.492
## Residencia -0.238 -0.842 -0.729
## Parques_Centros Estaciones_Transito Lugares_Trabajo
## PM10 0.065 0.270 0.100
## Reactivacion_Comercial 0.814 0.944 0.581
## Supermercado_Farmacia 0.655 0.872 0.492
## Parques_Centros 1.000 0.770 0.371
## Estaciones_Transito 0.770 1.000 0.474
## Lugares_Trabajo 0.371 0.474 1.000
## Residencia -0.617 -0.758 -0.853
## Residencia
## PM10 -0.238
## Reactivacion_Comercial -0.842
## Supermercado_Farmacia -0.729
## Parques_Centros -0.617
## Estaciones_Transito -0.758
## Lugares_Trabajo -0.853
## Residencia 1.000