Special case np.mean. Check for MaskesdColumn.
Added tests for reduceat aggregation
Mention the substitution of ufuncs with reduceat