U1A7

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

Contaminación del aire

Importar

Importar paquetes

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

Importar datos

datos <- read.csv("calidadaire.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 Agosto en Hermosillo (ERNO, UNAM)")

Análisis de correlación

names(NO2O3)

## [1] "NO2" "O3"

# Gráfico pairs
pairs(NO2O3)

cor(NO2O3)

##             NO2          O3
## NO2  1.00000000 -0.02779472
## O3  -0.02779472  1.00000000

Recta de mínimos cuadrados

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

## Warning: In lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
##  extra argument 'dsta' will be disregarded

summary(regresion)

## 
## Call:
## lm(formula = O3 ~ NO2, dsta = 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

Redacción personal

Se puede ver la relación entre ambos gases con el comando “cor”, el cual muestra un resultado negativo, en donde se asume que los datos son erróneos. Sin embargo, se puede evaluar que la concentración de estos gases es mayor en día, puesto que transitan más carros que generan esto.