Mercurial > repos > bimib > marea_2
diff marea_2/flux_to_map.py @ 248:419c6bfad590 draft
Uploaded
author | luca_milaz |
---|---|
date | Sun, 04 Aug 2024 13:29:55 +0000 |
parents | 6c3e7c338b93 |
children | a2adb47e630d |
line wrap: on
line diff
--- a/marea_2/flux_to_map.py Sun Aug 04 12:56:43 2024 +0000 +++ b/marea_2/flux_to_map.py Sun Aug 04 13:29:55 2024 +0000 @@ -814,18 +814,34 @@ return { id : list(map(utils.Float("Dataset values, not an argument"), values)) for id, values in dataset.items() }, IDs def computeEnrichmentMedoids(metabMap :ET.ElementTree, class_pat :Dict[str, List[List[float]]], ids :List[str]) -> None: - + #class_pat 462 * cellule utils.logWarning( - str(len(class_pat['1'])), + ids, ARGS.out_log) + medians = {} + means = {} + for key, value in class_pat.items(): + values = np.array(class_pat[key]) + median = np.median(values, axis=1) + mean = np.mean(values, axis=1) + medians[key] = median + means[key] = means + + max_flux_medians = max(medians.values(), key=abs) + max_flux_means = max(means.values(), key=abs) + + for key, value in medians.items(): + medians[key] = medians[key] / max_flux_medians + for key, value in means.items(): + means[key] = means[key] / max_flux_means + + + + pass - #for key, value in class_pat.items(): - #max_flux = max(class_pat.values(), key=abs) - - - + ############################ MAIN ############################################# def main() -> None: """