
Hello, thanks for answering. I hope this dataset helps.
I have tried the aqStats function:
aqStats(selectByDate(datao3, year = 2000:2017), pollutant = "o3", type = "code", percentile=c(5,25,50,75,90,95,99), data.thresh = 75, transpose=FALSE)
But I don't know how to manage the data thresh argument.