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 64f0800f10c1b7440e5360660f33e267260b273c..6b6b6f90e3956e1219d43b0b2e6a1fb253287745 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 @@ -207,9 +207,9 @@ 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[:1]: for trend_test_class in [GevLocationChangePointTest, GevScaleChangePointTest, - GevShapeChangePointTest][:1]: + GevShapeChangePointTest][2:]: years = [None][:] for last_starting_year in years: for days in [1, 3][1:]: diff --git a/experiment/trend_analysis/univariate_test/abstract_gev_change_point_test.py b/experiment/trend_analysis/univariate_test/abstract_gev_change_point_test.py index 2524d32c8ff69da1ddec530226780f2be69dfc9a..93a49a124a241e3d1219355a11a3bd363b04b4cf 100644 --- a/experiment/trend_analysis/univariate_test/abstract_gev_change_point_test.py +++ b/experiment/trend_analysis/univariate_test/abstract_gev_change_point_test.py @@ -62,17 +62,23 @@ class AbstractGevChangePointTest(AbstractUnivariateTest): @property def non_stationary_nllh(self): if self.crashed: - return -np.inf + return np.nan else: return self.non_stationary_estimator.result_from_fit.nllh @property def stationary_deviance(self): - return self.stationary_estimator.result_from_fit.deviance + if self.crashed: + return np.nan + else: + return self.stationary_estimator.result_from_fit.deviance @property def non_stationary_deviance(self): - return self.non_stationary_estimator.result_from_fit.deviance + if self.crashed: + return np.nan + else: + return self.non_stationary_estimator.result_from_fit.deviance @property def is_significant(self) -> bool: diff --git a/extreme_estimator/extreme_models/utils.py b/extreme_estimator/extreme_models/utils.py index 0ba41b4ae625b00f0d76a7cb7e91d55848e1ff2b..c2fe6ea0c62534f09001e462712ce4fc2560d955 100644 --- a/extreme_estimator/extreme_models/utils.py +++ b/extreme_estimator/extreme_models/utils.py @@ -87,6 +87,7 @@ def safe_run_r_estimator(function, data=None, use_start=False, threshold_max_abs run_successful = False res = None f = io.StringIO() + # Warning print will not work in this part with redirect_stdout(f): while not run_successful: current_parameter = parameters.copy()