aqStats function

Calculate summary statistics for air pollution data by year