U1A7

Silvia Flores

18/9/2020

Análisis estadísticos de calidad del aire del estado de Sonora el mes de agosto de 2020

Calidad del aire

Importar

Importar paquetes

library(pacman)
p_load("readr", "tidyverse", "DT")

Importar datos

Se importan desde un archivo local .csv, datos horarios sobre la calidad del aire durante el mes de agosto, obtenidos de: https://www.ruoa.unam.mx/index.php?page=estaciones&st=erno&id=6

datos <- read.csv("calidadaire.csv")

Transformar

#Los datos están como characters y aquí los covertimos a los valores que necesitamos
Fecha = as.Date(datos$Fecha)
NO2 <- as.numeric(datos$NO2)
O3 <- as.numeric(datos$O3)
CALIDAD <- data.frame(Fecha, NO2, O3)
no2o3 <- data.frame(NO2,O3)

Visualizar

#Creamos tabla interactiva
datatable(datos)
#Grafico
ggplot(data=CALIDAD) +  #gráfica de puntos
  geom_point(aes(Fecha, O3, colour = "O3")) +  #Le damos un color al O3
  geom_point(aes(Fecha, NO2, colour = "NO2")) + #Le damos un color al NO2
  xlab("Tiempo") +
  ylab("Concetración") + 
  labs (colour = "Gases") +  
  ggtitle("Calidad del aire para Agosto en Hermosillo (ERNO, UNAM)") #Titulo de la gráfica

Análisis de correlación

Gráfico de dispersión

names(no2o3)
## [1] "NO2" "O3"
#Gráfico pairs, como se relacionas el no2 y el o3
pairs(no2o3)

Grado de correlación lineal

para cuantificar el grado de relación líneal, calculamos la matriz de correlación

cor(no2o3)
##             NO2          O3
## NO2  1.00000000 -0.02779472
## O3  -0.02779472  1.00000000

Recta de minimos cuadrados

regresion <- lm(O3 ~ NO2, data=datos  )
summary(regresion)
## 
## Call:
## lm(formula = O3 ~ NO2, data = datos)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -20.458  -7.333  -2.115   6.152  42.735 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  23.4710     0.5828  40.275   <2e-16 ***
## NO2          -3.0717     4.1934  -0.733    0.464    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9.94 on 694 degrees of freedom
## Multiple R-squared:  0.0007725,  Adjusted R-squared:  -0.0006673 
## F-statistic: 0.5366 on 1 and 694 DF,  p-value: 0.4641

Análisis estadísticos de calidad del aire del estado de Sonora el mes de Enero de 2020

Calidad del aire

Importar

Importar paquetes

library(pacman)
p_load("readr", "tidyverse", "DT")

Importar datos

Se importan desde un archivo local .csv, datos horarios sobre la calidad del aire durante el mes de enero, obtenidos de: https://www.ruoa.unam.mx/index.php?page=estaciones&st=erno&id=6

datos2 <- read.csv("calidadaireenero.csv")

Transformar

#Los datos están como characters y aquí los covertimos a los valores que necesitamos
Fecha <- as.Date(datos2$Fecha)
NO2 <- as.numeric(datos2$NO2)
O3 <- as.numeric(datos2$O3)
CALIDAD2 <- data.frame(Fecha, NO2, O3)
no2o3 <- data.frame(NO2,O3)

Visualizar

#Creamos tabla interactiva
datatable(datos2)
#Grafico
ggplot(data=CALIDAD2) +  #gráfica de puntos
  geom_point(aes(Fecha, O3, colour = "O3")) +  #Le damos un color al O3
  geom_point(aes(Fecha, NO2, colour = "NO2")) + #Le damos un color al NO2
  xlab("Tiempo") +
  ylab("Concetración") + 
  labs (colour = "Gases") +  
  ggtitle("Calidad del aire para Enero en Hermosillo (ERNO, UNAM)") #Titulo de la gráfica

Análisis de correlación

Gráfico de dispersión

names(no2o3)
## [1] "NO2" "O3"
#Gráfico pairs, como se relacionas el no2 y el o3
pairs(no2o3)

Grado de correlación lineal

para cuantificar el grado de relación líneal, calculamos la matriz de correlación

cor(no2o3)
##            NO2         O3
## NO2  1.0000000 -0.2586445
## O3  -0.2586445  1.0000000

Recta de minimos cuadrados

regresion <- lm(O3 ~ NO2, data=datos  )
summary(regresion)
## 
## Call:
## lm(formula = O3 ~ NO2, data = datos)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -20.458  -7.333  -2.115   6.152  42.735 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  23.4710     0.5828  40.275   <2e-16 ***
## NO2          -3.0717     4.1934  -0.733    0.464    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9.94 on 694 degrees of freedom
## Multiple R-squared:  0.0007725,  Adjusted R-squared:  -0.0006673 
## F-statistic: 0.5366 on 1 and 694 DF,  p-value: 0.4641

Conclusión

En este ejercicio vimos como fue cambiando la calidad del aire el mes de agosto y enero, se muestra que durante la noche una mejor calidad del aire debido a que hay una menor concentracion de gases contaminantes. En la gráfica se puede observar la relación que hay entre el NO2 y O3. Los datos no están relacionados porque nos está dando lecturas negativas. (No podemos hacer una buena lectura de los resultados ya que los datos con los que realizamos el ejercicio parecen ser incorrectos)