U1A7

Análisis estadístico de calidad del aire del estado de Sonora

Contaminación

Importar

Importar paquetes

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

Importar datos

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

Transformar

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

#Tabla
datatable(datos)

#Grafico

ggplot(data=calidad) +
  geom_point(aes(Fecha, O3,  colour ="O3" )) +
  geom_point(aes(Fecha, NO2, colour = "NO2")) +
  xlab("Tiempo") +
  ylab("Concentración") + 
  labs (colour = "Estados") +
  ggtitle("Calidad del aire para Junio en Hermosillo (ERNO, UNAM)")

Analísis de correlación

Grafico de dispeción

names(No2O3)

## [1] "NO2" "O3"

#grafico pairs
pairs(No2O3)

Grado de correlación lineal

• para cuantificar el grado de correlacion lineal calculamos la matris de correlacion

cor(No2O3)

##           NO2        O3
## NO2 1.0000000 0.1462665
## O3  0.1462665 1.0000000

Recata de minimos cuadrados

regresion <- lm(O3 ~ NO2, data = datos)
summary(regresion)

## 
## Call:
## lm(formula = O3 ~ NO2, data = datos)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -31.2458  -7.2727  -0.1263   7.0646  30.6130 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  29.1888     0.5405  54.002  < 2e-16 ***
## NO2          11.6975     3.0564   3.827 0.000142 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 10.11 on 670 degrees of freedom
## Multiple R-squared:  0.02139,    Adjusted R-squared:  0.01993 
## F-statistic: 14.65 on 1 and 670 DF,  p-value: 0.0001418

Conclusión

No es posible visualizar de manera correcta la correlación que hay entre estos dos gases debido a que al parecer el sensor no esta registrando correctamente los datos, ya que hay valores negativos entre estos, pero se puede observar que se produce más cantidad de ozono (O3) en el día que en la noche, como ya sabemos, esto se debe a la circualación de la mayoria de transportes de la ciudad (automoviles, trailes, camiones, etc), ya que el dióxido de carbono (CO2) entra en contacto con el dióxido de nitrogeno (NO2) produciendo una reacción fotoquímica generando el ozono (03).