From aaffc3b4a427984fbab2fa86ea7d48282b2aee31 Mon Sep 17 00:00:00 2001
From: Le Roux Erwan <erwan.le-roux@irstea.fr>
Date: Tue, 11 Jun 2019 14:15:00 +0200
Subject: [PATCH] [HYPERCUBE] display best starting year on mean log lik plots
---
.../altitude_hypercube_visualizer.py | 3 +++
.../main_hypercube_visualization.py | 15 +++++++++++----
2 files changed, 14 insertions(+), 4 deletions(-)
diff --git a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py
index cd477e12..31efe3ca 100644
--- a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py
+++ b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py
@@ -134,6 +134,9 @@ class AltitudeHypercubeVisualizer(AbstractHypercubeVisualizer):
for ax_idx, (ax, (ylabel, serie)) in enumerate(zip(axes_remaining, ylabel_to_series.items())):
xlabel_values = list(serie.index)
values = list(serie.values)
+ argmax_idx = np.argmax(values)
+ best_year = xlabel_values[argmax_idx]
+ plot_title += '{}'.format(best_year)
ax.plot(xlabel_values, values)
ax.set_ylabel(ylabel)
ax.set_title(plot_title)
diff --git a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/main_hypercube_visualization.py b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/main_hypercube_visualization.py
index 6b6b6f90..8908e616 100644
--- a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/main_hypercube_visualization.py
+++ b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/main_hypercube_visualization.py
@@ -117,10 +117,17 @@ def fast_altitude_without_trend_type_hypercube_extended():
trend_test_class=trend_test_class,
nb_data_reduced_for_speed=nb_data_reduced_for_speed,
last_starting_year=last_starting_year)
- # visualizer.visualize_year_trend_test()
visualizer.vsualize_year_trend_by_regions_and_altitudes()
+ # visualizer = AltitudeHypercubeVisualizerWithoutTrendType(altitude_to_visualizer,
+ # save_to_file=save_to_file,
+ # trend_test_class=trend_test_class,
+ # nb_data_reduced_for_speed=nb_data_reduced_for_speed,
+ # last_starting_year=last_starting_year)
+ # visualizer.visualize_year_trend_test()
+
+
def fast_altitude_year_hypercube_extended():
save_to_file = False
only_first_one = False
@@ -207,10 +214,10 @@ def full_altitude_year_hypercube__without_trend_extended():
only_first_one = False
nb_data_reduced_for_speed = False
altitudes = ALL_ALTITUDES[3:-6]
- for study_class in SCM_STUDIES[:1]:
+ for study_class in SCM_STUDIES[:]:
for trend_test_class in [GevLocationChangePointTest, GevScaleChangePointTest,
- GevShapeChangePointTest][2:]:
- years = [None][:]
+ GevShapeChangePointTest][:1]:
+ years = [2007, None][:]
for last_starting_year in years:
for days in [1, 3][1:]:
visualizers = [
--
GitLab