Calidadairetemporal

Calidad de aire

library(dplyr)
Sys.setenv(TZ = "UTC") # don't use local time zone
ruta<-"C:/Users/cguer/Documents/Claudia/Midropbox/Investigacion y escritos/karamanis/data/"
infile <- "calairetorreon.CSV"
nomarchi<-paste0(ruta,infile)
caliairehora<-read.csv(nomarchi,stringsAsFactors = FALSE)
caliairehora$t<-as.POSIXct(caliairehora$Fecha,format = "%d/%m/%Y")
caliairehora$PM10<-as.numeric(caliairehora$PM10)
caliairehora$PM2.5<-as.numeric(caliairehora$PM2.5)
caliairehora$O3<-as.numeric(caliairehora$O3)
caliairehora$NO2<-as.numeric(caliairehora$NO2)
caliairehora$SO2<-as.numeric(caliairehora$SO2)
caliairehora$CO<-as.numeric(caliairehora$CO)

caliairehoradia  <-caliairehora%>% 
  filter(t == ymd("2025-03-17"))

       #aes(x=as.POSIXct(caliairehora$Fecha,format = "%d/%m/%Y %H:%M"),

dia<-caliairehoradia$t[1]
a<-ggplot(caliairehoradia,
 #         aes(x=caliairehoradia$Fecha,
           aes(x=as.POSIXct(caliairehoradia$Fecha,format = "%d/%m/%Y %H:%M"),

           y=caliairehoradia$PM10))+
       
     geom_point()+ geom_line() + 
      theme(text = element_text(size=10),
        axis.text.x = element_text(angle=90, hjust=1)) +
        ggtitle(paste0("PM10 por hora ",dia)) +
          xlab("hora") + ylab("PM10 u/m3")

a

b<-ggplot()+
   geom_line(data=caliairehoradia, aes(x=as.POSIXct(caliairehoradia$Fecha,format = "%d/%m/%Y %H:%M"),y= caliairehoradia$PM10, colour = "PM10"))+
   geom_line(data=caliairehoradia, aes(x=as.POSIXct(caliairehoradia$Fecha,format = "%d/%m/%Y %H:%M"),y= caliairehoradia$PM2.5, colour ="PM2.5"))+
   
   theme(text = element_text(size=10),
        axis.text.x = element_text(angle=90, hjust=1)) +
        ggtitle(paste0("PM10 y PM2.5 por hora ",dia)) +
          xlab("hora") + ylab(" u/m3")

b

c<-ggplot()+
   geom_line(data=caliairehoradia, aes(x=as.POSIXct(caliairehoradia$Fecha,format = "%d/%m/%Y %H:%M"),y= caliairehoradia$O3, colour = "O3"))+
   geom_line(data=caliairehoradia, aes(x=as.POSIXct(caliairehoradia$Fecha,format = "%d/%m/%Y %H:%M"),y= caliairehoradia$NO2, colour ="NO2"))+
   geom_line(data=caliairehoradia, aes(x=as.POSIXct(caliairehoradia$Fecha,format = "%d/%m/%Y %H:%M"),y= caliairehoradia$SO2, colour ="SO2"))+
   geom_line(data=caliairehoradia, aes(x=as.POSIXct(caliairehoradia$Fecha,format = "%d/%m/%Y %H:%M"),y= caliairehoradia$CO, colour ="CO"))+

   theme(text = element_text(size=10),
        axis.text.x = element_text(angle=90, hjust=1)) +
        ggtitle(paste0("PM10 y PM2.5 por hora ",dia)) +
          xlab("hora") + ylab(" u/m3")

c