A4U1

setwd("~/Pye 1pm")

library(pacman)
p_load("datasets", "DT", "fdth", "prettydoc", "xfun", "ggplot2", "beanplot")

Ejercicio #4: Mediciones de la calidad del aire en Nueva York

Mediciones diarias de la calidad del aire en Nueva York de Mayo a Septiembre 1973.

##Conociendo los datos ## A través de una tabla interactiva

datatable(airquality )

##A través de una gráfica, usando ggplot2

xvol <- airquality 

ggplot(xvol, aes(Temp, Wind ) ) + geom_point(aes(col=Month), size=3) 

Medidas exploratorias de tendencia central

summary(xvol)

##      Ozone           Solar.R           Wind             Temp      
##  Min.   :  1.00   Min.   :  7.0   Min.   : 1.700   Min.   :56.00  
##  1st Qu.: 18.00   1st Qu.:115.8   1st Qu.: 7.400   1st Qu.:72.00  
##  Median : 31.50   Median :205.0   Median : 9.700   Median :79.00  
##  Mean   : 42.13   Mean   :185.9   Mean   : 9.958   Mean   :77.88  
##  3rd Qu.: 63.25   3rd Qu.:258.8   3rd Qu.:11.500   3rd Qu.:85.00  
##  Max.   :168.00   Max.   :334.0   Max.   :20.700   Max.   :97.00  
##  NA's   :37       NA's   :7                                       
##      Month            Day      
##  Min.   :5.000   Min.   : 1.0  
##  1st Qu.:6.000   1st Qu.: 8.0  
##  Median :7.000   Median :16.0  
##  Mean   :6.993   Mean   :15.8  
##  3rd Qu.:8.000   3rd Qu.:23.0  
##  Max.   :9.000   Max.   :31.0  
## 

Visualización

Histograma

hist(airquality$Temp,
     col="red",
     main= "Histograma de relacion entre las distintas temperaturas",
     xlab= "Temperatura",
     ylab= "Frecuencia absoluta")

Diferentes maneras de realizar histogramas

par(mfrow=c(1,2))
hist(xvol$Wind,xlim=c(0,8),ylim=c(0,40))
hist(xvol$Month,xlim=c(0,8),ylim=c(0,40))

## Gráfico de dispersión (scatterplot)

Crear un gráfico de dispersión

En este caso usaremos la función “plot”

plot ( xvol$Temp, xvol$Wind,
      col = "red",
      Main = "Gráfico de dispersión",
      xlab = "Temperatura",
      ylab = "Vientos",
      pch = 19
)

## Warning in plot.window(...): "Main" is not a graphical parameter

## Warning in plot.xy(xy, type, ...): "Main" is not a graphical parameter

## Warning in axis(side = side, at = at, labels = labels, ...): "Main" is not a
## graphical parameter

## Warning in axis(side = side, at = at, labels = labels, ...): "Main" is not a
## graphical parameter

## Warning in box(...): "Main" is not a graphical parameter

## Warning in title(...): "Main" is not a graphical parameter

Gráfico comparativo de “beanplot”

xxvol <- xvol
xxvol$Ozone <- NULL
xxvol$Solar.R <- NULL
xxvol$Day <- NULL
beanplot (xxvol, main="Mediciones de distribuición de vientos, temperatura y meses")

## Temperatura respecto a los vientos utilizando un gráfico de caja y bigotes

boxplot(Wind~Temp,
        data=xxvol,
        main="Temperaturas con los vientos",
        xlab="Viento ",
        ylab="Temperatura",
        col= "red",
        border="blue")

## Interpretando un gráfico de caja y bigote

#Extracción de subconjuntos

xt <- subset(xvol, Temp == “Temperatura”) xw <- subset(xvol, Wind == “Viento”)

xfun::embed_file("A4U1.Rmd")

Download A4U1.Rmd