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

descarga.jpg

Importar

Importar paquetes

library(pacman)
p_load("readr", "tidyverse", "DT")
setwd("~/PYE1112ITSON")

Importar datos

datos <- read.csv("junio.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 Julio 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

Se puede observar que hay una gran cantidad de Ozono en el aire en este mes ya sea por los carros, las empesas, claro no se puede apreciar tan bien ya que el sensor no esta funcionando bien, como lo podemos observar ya que nos da valroes negativos,el ozono se debe a una reacción fotoquímica NO2 producido por los automoviles generando el ozono O3, como lo podemos observar para generar ozono ocupamos muy poco No2 y si reducimos un poco el No2 el ozono bajaria drasticamente esto se pude observar en la grafica de puntos. Los datos fueron representado en histogramas y tablas gracias a las bibliotecas ya utilizadas para la maxima comprención de los datos, se hizo como lo emos estado haciendo leyendo los datos,creandos data frames,y gracias a estos nos resulta mas facil analizar los datos de la calidad del aire en el mes de julio en Sonora.