stats.py fixes and test coverage

- Fixes test failure with older numpy.
- Fixes ``stats.tsem``, it was returning ``s/n`` instead of ``s/sqrt(n)`` (no tests is really bad, yet again).
- Provides some test coverage for all trimmed statistics.

There's still some inconsistencies in the signatures of the trimmed stats functions, that should be fixed sometime soon.