Calculate harmonic mean
harmonic(data)
data
Data.
numeric vector of length nlevels(data$stratum.code)
nlevels(data$stratum.code)