code: def join_complexes(*complexes):
      """Join a sequence of simplicial complexes.
  
      Returns the union of all the particular faces.
      """
      faces = {}
  
      nmax = np.array([len(c) for c in complexes]).max()
      for i in range(nmax):
          faces[i + 1] = set()
      for c in complexes:
          for i in range(nmax):
              if i + 1 in c:
                  faces[i + 1] = faces[i + 1].union(c[i + 1])
      return faces
  
file_path: nipy/nipy/algorithms/statistics/utils.py
imports: from itertools import combinations
  import numpy
  from scipy.stats import norm
label_desc: Use image node_id 6 for calling a global function with example usage: join_complexes() and returns: faces