use reduceat when aggregating to improve performance
fix mask bug
Special case np.mean. Check for MaskesdColumn.
Added tests for reduceat aggregation
Mention the substitution of ufuncs with reduceat
Update doc and CHANGES
Merge pull request #2625 from luzuku/group-aggregate

Use reduceat() in group aggregation where possible.  Closes #1546.