Имея набор данных «data» и используя пакет openair, я создал timePlot, используя следующий код:
a=timePlot(data, pollutant="pm2.5", type="site", avg.time="24 hour", plot.type="s", ref.y=list(h=37, lty=2, col="black", lwd=1.6), smooth=TRUE, ci=FALSE, ylab="", main="Concentración diaria promedio PM2.5 por estación de monitoreo", ylim=c(0,50))
Получив следующий результат:
Обратите внимание на опорную линию, которую я добавил при y = 37, используя функцию ref.y, которая представляет собой дневной предел средней концентрации для pm2.5. Однако мне нужно 50 мкг / м3 с 1 марта по 30 июня и 37 мкг / м3 с 1 июля. Это соответствует Руководству по качеству воздуха в моей стране.
В некотором роде это выглядит так:
Я очень признателен за любую помощь , С наилучшими пожеланиями.