Hierarchical Clustering Algorithms on Poisson and Other Stationary Point Processes

This paper introduces a hierarchical clustering algorithm, the Clustroid Hierarchical Nearest Neighbor ($\mathrm{CHN}^2$), designed for datasets with a countably infinite number of data represented by points in the Euclidean space. The method builds clusters across successive levels by linking nearest-neighbor points or clusters using the clustroid distance. The properties of this algorithm make it suitable for very large datasets.
To evaluate its properties, we first apply the algorithm to the homogeneous Poisson point process, which serves as a natural null-hypothesis model with no intrinsic data aggregation. In this setting, the algorithm generates a random forest that is a deterministic factor of the Poisson point process, and hence unimodular. We prove that at every level, the level-$k$ graph has only finite connected components (a.s.) and derive bounds on their mean size. We also establish the existence of a limiting graph as the number of levels tends to infinity. In this limit, clusters are shown to be all infinite and one-ended, which induces a natural order within each component and supports a tree-like phylogenetic interpretation.
Beyond the Poisson case, we extend the analysis to a class of Cox and more general stationary point processes without second-order descending chains (introduced here), for which analogous results hold. Simulations show that comparing the Cox case with the Poisson baseline allows an efficient detection of aggregation, thereby linking the stochastic-geometric analysis to practical clustering tasks.