We propose a numerical, algorithm for statistical investigation of the
grouping of galaxies, via the construction of a hierarchy of subsyste
ms in a given system. This algorithm uses the position and velocity in
formation for the individual galaxies.