refactor: standardize on scipy.special.comb over binom

Replace the last two uses of scipy.special.binom with comb for
consistency with the rest of the codebase.

Co-Authored-By: Claude Opus 4.6 <noreply@anthropic.com>

Author: leotrs

xgi/algorithms/simpliciality.py

@@ -472,11 +472,11 @@ def _missing_subfaces(t, face, min_size=1):
 
 
 def _max_number_of_subfaces(min_size, max_size):
-    from scipy.special import binom
+    from scipy.special import comb
 
     d = 2**max_size - 2  # subtract 2 for the face itself and the empty set
     for i in range(1, min_size):
-        d -= binom(max_size, i)
+        d -= comb(max_size, i, exact=True)
     return int(d)
 
 

xgi/utils/tensor.py

@@ -74,10 +74,10 @@ def banerjee_coeff(size, max_size):
     Sinan Aksoy, Ilya Amburg, Stephen Young,
     https://doi.org/10.1137/23M1584472
     """
-    from scipy.special import binom
+    from scipy.special import comb
 
     return sum(
-        ((-1) ** j) * binom(size, j) * (size - j) ** max_size for j in range(size + 1)
+        ((-1) ** j) * comb(size, j) * (size - j) ** max_size for j in range(size + 1)
     )