Refactor clustering based on averages expression within devStage

The plyr way isn’t shorter, but seems better to model. I suspect we should really be exploiting reshape2 here.
Delete unnecessary par() fussing