I’ve made a new chart. Don’t look too closely though! And it looks much better during the WP edit.

To create this chart I created a new vector to remove the 2017 data. And now I’ve introduded geom_boxplot and facet functions:

I’ve created a vector from the new vector and decided to have month on the x-axis, and PM10 readings on the y-axis.

p <- ggplot(data=oxford_road_2018_2020, mapping = aes(x=month,y=PM10))

Now I want to generate the plot.

p + geom_boxplot(aes(group=day)) + facet_wrap(~year) + labs(title = “AQE PM10 by Year”)

There’s a mistake if you look more closely at fig 1. There’s more plots on the x-axis than there should be. I missed it the first time because I’d created a nice chart. Let’s try and fix that.

p + geom_boxplot(aes(group=month)) + facet_wrap(~year) + labs(title…