Using the R package openair

M. Salmon

2019-01-31

In this vignette I want to show how easy and useful it is to use the openair package with Open AQ data.

library("ropenaq")
library("openair")
library("dplyr")

Load data for the Netherlands

measurementsNL <- aq_measurements(location="Amsterdam-Einsteinweg", parameter="pm25",
                                       date_from = as.character(Sys.Date() - 100),
                                       date_to = as.character(Sys.Date()))

# filter negative values
# and rename columns for compatibilities with openair
# I do not drop the old columns though
measurementsNL <- dplyr::mutate(measurementsNL, date=dateLocal,
                                pm25=value) %>%
  filter(value>=0)

# for now openair functions do not work with tbl_df objects
# on MY computer
measurementsNL <- as.data.frame(measurementsNL)

Timeplot

# useful timeplot
timePlot(mydata = measurementsNL, pollutant = "pm25")

Calendarplot

# cool calendar plot 
if(as.numeric(format(Sys.Date(), "%m")) < 2){
  calendarPlot(mydata=measurementsNL, pollutant = "pm25",
             year = as.numeric(format(Sys.Date(), "%Y")) - 1)
}else{
  calendarPlot(mydata=measurementsNL, pollutant = "pm25",
             year = as.numeric(format(Sys.Date(), "%Y")))
}