Changed the visualization to use the readings.aggregate file not the old readings.count.
Fixing CDH breakpoint at 18C for the time being.

Reduced the modelling section to use only
lm(...) rather than adding glm(...) for now.
I can find AIC simply by using AIC(...)
function.

Finally, created an area where I try to
make the most descriptive model possible
regardless of AIC or fixed effect meaning.