code: def central_likelihood(self, axis):
          
          result = self.cumulative_density(axis)
          result.histogram = 1 - 2 * np.abs(result.histogram - 0.5)
          return result
docs: Returns new histogram with all values replaced by their central likelihoods along axis.