diff --git a/extreme_trend/abstract_gev_trend_test.py b/extreme_trend/abstract_gev_trend_test.py index bf235c3129b5d44fb35e046892d2503339cffbe4..c50de9988708c695597fbe0328a924fa846c9e51 100644 --- a/extreme_trend/abstract_gev_trend_test.py +++ b/extreme_trend/abstract_gev_trend_test.py @@ -1,4 +1,5 @@ from math import ceil, floor +import numpy.testing as npt import matplotlib.pyplot as plt import numpy as np @@ -113,7 +114,9 @@ class AbstractGevTrendTest(object): @property def aic(self): aic = 2 * self.total_number_of_parameters_for_unconstrained_model + self.unconstrained_model_deviance - assert np.equal(self.unconstrained_estimator.result_from_model_fit.aic, aic) + assert np.equal(self.total_number_of_parameters_for_unconstrained_model, self.unconstrained_estimator.margin_model.nb_params) + npt.assert_almost_equal(self.unconstrained_estimator.result_from_model_fit.aic, aic, decimal=12) + npt.assert_almost_equal(self.unconstrained_estimator.aic(), aic, decimal=12) return aic @property