stats.mood generalization and testing

Added a new function stats.mood_nd and some tests for it and stats.mood. stats.mood_nd has an axis along which the testing of ndimensional arrays can be performed along the specified axis.
Note: since stats.rankdata is not vectorized, use a loop for ranking... Probably this has to be vectorized in the future.