code: def reduced_chi2(self, model, error_map=0):
          
          chi2 = self.reduced_residuals(model, error_map)
          return np.sum(chi2**2) / self.num_data_evaluate()
docs: returns reduced chi2
          :param model:
          :param error_map:
          :return: