def METHOD_NAME(self, parameters, neq):
    """
    Calculates the Harmonic mean of a vector of parameters
    """
    weight = neq.astype(float) / np.sum(neq)
    if np.shape(parameters)[0] != np.shape(weight)[0]:
        raise ValueError('Parameter vector not same shape as weights')
    else:
        average_value = np.zeros(np.shape(parameters)[0], dtype=float)
        id0 = np.logical_not(np.isnan(parameters))
        average_value = 1. / np.sum(weight[id0] / parameters[id0])
    return average_value