R code and data (from Gorman & Johnson 2011)

You may also want to give Rbrul a spin.

[ ~ ]