docs: def objective(self, params):
          """Compute the negative penalized log-likelihood."""
          val = self._penalty * np.sum(params**2)
          for win, los in self._data:
              val += np.logaddexp(0, -(params[win] - params[los]))
          return val