From e9409ea33ce56ccb366583b11d36ba96e9ae2b47 Mon Sep 17 00:00:00 2001
From: Le Roux Erwan <erwan.le-roux@irstea.fr>
Date: Tue, 16 Jun 2020 16:55:52 +0200
Subject: [PATCH] [contrasting] add aic test in abstract_gev_trend_test.py
---
extreme_trend/abstract_gev_trend_test.py | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/extreme_trend/abstract_gev_trend_test.py b/extreme_trend/abstract_gev_trend_test.py
index bf235c31..c50de998 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
--
GitLab