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

Importar

Importar paquetes

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

Importar datos

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

## Warning: Removed 25 rows containing missing values (geom_point).

## Warning: Removed 25 rows containing missing values (geom_point).

Análisis de correlación

Gráfico de disperión

names(no2o3)

## [1] "NO2" "O3"

#Gráfico pairs
pairs(no2o3)

Grado de correlación lineal

Para cuantificar el grado de correlación lineal calculamos la matriz de correlación

cor(no2o3)

##     NO2 O3
## NO2   1 NA
## O3   NA  1

Recta de mínimos cuadrados

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

## 
## Call:
## lm(formula = O3 ~ NO2, data = datos)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -18.491 -10.208  -3.825   9.116  32.998 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  15.1885     0.4887  31.077  < 2e-16 ***
## NO2         -10.1285     1.4380  -7.043 4.52e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 12.72 on 695 degrees of freedom
##   (25 observations deleted due to missingness)
## Multiple R-squared:  0.06663,    Adjusted R-squared:  0.06528 
## F-statistic: 49.61 on 1 and 695 DF,  p-value: 4.518e-12

Redaccion

En el ejercicio realizado elaboramos graficas para poder observar la relación entre los gases de O3 y NO2 en el mes de Enero. Pudimos observar que con poco NO2 se genera mucho O3 este dato lo pudimos agarrar en la primera grafica. En la segunda grafica analizamos su correlación. Mas adelante hicimos un summary de la regresión.

U1A7

Jose Ibarra

18/9/2020