"harmonic.mean" <- 
function(x) { 1/(mean(1/x)) }