This example produces a two-way ANOVA analysis of the data in table 7.4 from Plane Answers to Complex Questions by Christensen. The calculations are done using native R functions, and duplicated via custom R functions that implement concepts from linear model theory. The relevant sections from Christensen are referenced in the comments.