This is a piece of code I implemented in 2004, which was supposed to be part of an R-package in multivariate testing (to be named, rather creatively, mvttests). Time has flown, I haven't still got around to implementing the said package, but people keep asking me for the sphericity.test function, so here it is, for posterity: sphericity.test More