I am trying to do a trending analysis of air pollution data. My data is in the form (with more dates, some values of PM10 are NA):
structure(list(NetC = c("Cosenza Provincia", "Cosenza Provincia",
"Cosenza Provincia", "Cosenza Provincia", "Cosenza Provincia",
"Cosenza Provincia", "Cosenza Provincia", "Cosenza Provincia",
"Cosenza Provincia", "Reti Private", "Reti Private", "Reti Private",
"Reti Private", "Reti Private", "Reti Private"), ID = c("IT1938A",
"IT1938A", "IT1938A", "IT2086A", "IT2086A", "IT2086A", "IT2110A",
"IT2110A", "IT2110A", "IT1766A", "IT1766A", "IT1766A", "IT2090A",
"IT2090A", "IT2090A"), Stat = c("Citta dei Ragazzi", "Citta dei Ragazzi",
"Citta dei Ragazzi", "Rende", "Rende", "Rende", "Acri", "Acri",
"Acri", "Firmo", "Firmo", "Firmo", "Schiavonea", "Schiavonea",
"Schiavonea"), Data = c("1/1/2022", "1/2/2022", "1/3/2022", "1/1/2022",
"1/2/2022", "1/3/2022", "1/1/2022", "1/2/2022", "1/3/2022", "1/1/2022",
"1/2/2022", "1/3/2022", "1/1/2022", "1/2/2022", "1/3/2022"),
PM10 = c(13.29, 11.14, 9.08, 16.62, 12.98, 10.4, 16.2, 19.4,
15.7, 10.82, 12.29, 9.54, 24.54, 22.88, 27.33)), class = "data.frame", row.names = c(NA,
-15L))
I have tried using openair::TheilSen and openair::polarPlot, but it doesn't work (because I need to group the data for each station).
I actually need 5 graphs, but I can't group them.
p1 <- polarPlot(Acri, poll = "PM10", stati = "nwr",
main = "polarPlot of PM10 (Acri)")
p2 <-polarPlot(Citta, poll = "PM10", stati = "nwr",
main = "polarPlot of PM10 (Città dei Ragazzi)")
p3 <- polarPlot(Rende, poll = "PM10", stati = "nwr",
main = "polarPlot of PM10 (Rende)")
library(gridExtra)
grid.arrange(p1, p2, p3, ncol=1,top="polarPlot of PM10")
The output should look like:
