Optimal change point detection and localization in sparse dynamic networks

From MaRDI portal



Abstract: We study the problem of change point localization in dynamic networks models. We assume that we observe a sequence of independent adjacency matrices of the same size, each corresponding to a realization of an unknown inhomogeneous Bernoulli model. The underlying distribution of the adjacency matrices are piecewise constant, and may change over a subset of the time points, called change points. We are concerned with recovering the unknown number and positions of the change points. In our model setting we allow for all the model parameters to change with the total number of time points, including the network size, the minimal spacing between consecutive change points, the magnitude of the smallest change and the degree of sparsity of the networks. We first identify a region of impossibility in the space of the model parameters such that no change point estimator is provably consistent if the data are generated according to parameters falling in that region. We propose a computationally-simple algorithm for network change point localization, called Network Binary Segmentation, that relies on weighted averages of the adjacency matrices. We show that Network Binary Segmentation is consistent over a range of the model parameters that nearly cover the complement of the impossibility region, thus demonstrating the existence of a phase transition for the problem at hand. Next, we devise a more sophisticated algorithm based on singular value thresholding, called Local Refinement, that delivers more accurate estimates of the change point locations. Under appropriate conditions, Local Refinement guarantees a minimax optimal rate for network change point localization while remaining computationally feasible.


The study refers to a dynamic network model consisting of a sequence (A(t))\(^T_{t=1}\) of \(n\times n\) adjacency matrices of independent inhomogeneous Bernoulli networks with means (\(\theta\)(t))\(^T_{t=1}\) that satisfy a set of conditions, summarized as Assumption~1. The model is described by the parameters: \(\Delta\), \(k_0\), \(n\) and \(\rho\). When \(T\) grows unbounded, the parameters will change as functions of \(T\) and the number \(K\) of change points will change, too. One is concerned with the problem of estimating the number and location of the change points, based on one observation \((A(1),\dots, A(T))\) of adjacency matrices satisfying Assumption~1. The feasibility of a consistent estimation of the change point location is studied in the second section. One develops an algorithm called network binary segmentation, NBS, which is able to localize multiple change points consistently. In the third section, one shows that under stronger conditions on model and scaling, a two-step procedure can achieve a minimax optimal localization rate. The procedure applies first NBS and then a local refinement procedure, the LR algorithm. The authors prove that the LR procedure indeed improves localization rates, it is nearly minimax optimal. New assumptions are introduced. Simulation results are presented in the fourth section. Proofs of the main results are presented in the Appendix part. Some details, discussion and technical elements can be found in the supplementary material to this article, \url{doi:10.1214/20-AOS1953SUPP;.pdf}.




Cited in
(48)


Describes a project that uses

Uses Software






This page was built for publication: Optimal change point detection and localization in sparse dynamic networks

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q97726)