Replace agg(harmonic_mean) with custom function

[Pennycook]

Calling agg(harmonic_mean) results in a KeyError when attempting to aggregate over a group with a single member.

Related issues

N/A

Proposed changes

• Wrap statistics.harmonic_mean to ensure that it works when passed to agg.