harmonic mean code snippets in r