docs: def log_likelihood_pairwise(data, params):
      """Compute the log-likelihood of model parameters."""
      loglik = 0
      for winner, loser in data:
          loglik -= np.logaddexp(0, -(params[winner] - params[loser]))
      return loglik