I was teaching on a field course and when the students started analysing their data in R, one of them noticed that if they switched around the independent variables in an lm() they got different results with different methods of computing analysis of variance tables. I wanted to investigate it more, and this is the resulting R Markdown report that I wrote:
The report can be found here .