docs: def zinb_ll(data, P, R, Z):
      """
      Returns the zero-inflated negative binomial log-likelihood of the data.
      """
      lls = nb_ll(data, P, R)
      clusters = P.shape[1]
      for c in range(clusters):
          pass
      return lls