Add new R script for NOx statistics