def bic(self): """ returns the bayesian information criterion of the model. :return: bic value, float """ num_data = self.likelihoodModule.num_data num_param_nonlinear = self.param_class.num_param()[0] num_param_linear = self.param_class.num_param_linear() num_param = num_param_nonlinear + num_param_linear bic = analysis_util.bic_model(self.best_fit_likelihood, num_data,num_param) return bic
def test_bic_model(self): bic = analysis_util.bic_model(0, np.e, 1) assert bic == 1