fixed inconsistencies for routines sem() in stats and mstats. Both now have a ddof parameter. For backward compatability reasons, their default values however differ! Denominator sqrt(n-1) was replaced by sqrt(n), as this is the proper approach. However, needs to be checked, if sqrt(n-ddof) would actually be actually the right way.
BUGfix in routine mstats.tvar(): the normalization factor n/(n-1.) was missing. Unittests for the routine with given limits are still missing.
BUGFIX: modified mstats.linregress(): stats.linregress() and mstats.linregress() now return the same values. Inconsistencies have been resolved by explicitely calling stats.linregress() instead of doing separate calculations in mstats.
BUGFIX: modified mstats.skewtest(): stats.skewtest() and mstats.skewtest() did return different results for the p-value. In fact, the mstats approach seemed to be wrong. Instead of fixing the problem, I changed mstats.skewtest() in a way that it just acts as a wrapper for the stats.skewtest() function
BUGFIX: fixed error in the unittesting of routine test_obrientransform(). Note that results are transformed. This was however kept like it is to be backward compatible
BUGFIX (MAJOR): major bugfix in mstats.tsem(); returned results were inconsistent with the corresponding stats.tsem() function. The return value was incorrectly normalized. Instead of sd/sqrt(n), sd/n was returned. It is also ensured that variance is calculated with ddof=1.
editorial
BUGFIX: stats.describe() and mstats.describe() were inconsistent. Added new parameter ddof in both routines, so one can specify the degree of freedom when calling. Changes are backward compatible which results still in different default values for ddof for both routines
BUGFIX: small bugfix in mstats.trima(); in case that both, upper/lower limit were None, the condition was not handled properly. Now the array  is directly returned.