fn save_genders(
    clusters: Vec<i32>,
    genders: HashMap<i32, ClusterStats>,
    outf: &Path
) -> Result<()>