powered by
Harmonic mean times mean of weights
harmonic_times_mean_weight(data)
numeric vector of length nlevels(data$stratum.code)
nlevels(data$stratum.code)