Using the openair package, I created a timeplot using the following code:
a=timePlot(data, pollutant="pm2.5", type="site", avg.time="24 hour", plot.type="s", smooth=TRUE, ci=FALSE, ylab="", main="Concentración diario promedio PM2.5 por estación de monitoreo", ylim=c(0,60))
(The resulting plot is in the image link down below as I can only post 1 image per post).
Notice I created a reference line at y=37 using the
ref.y=list(h=37) function, however, I need this value to be 50 from march 1st to june 30th and 37 from july 31 onwards, this is due to my country's air quality guidelines.
I want to create something like this:
I'd appreciate any help.