Table of Contents
Fetching ...

On the Effect of Misspecifying the Embedding Dimension in Low-rank Network Models

Roddy Taing, Keith Levin

TL;DR

The paper addresses how adjacency spectral embedding (ASE) behaves when the embedding dimension is misspecified in latent-space network models, notably the random dot product graph (RDPG) and its binary and weighted variants. It develops a signal-plus-noise framework $\mathbf{A}=\mathbf{P}+\mathbf{E}$ with $\mathbf{P}=\rho_n\mathbf{X}\mathbf{X}^\top$ and proves a key delocalization result: all non-signal eigenvectors of $\mathbf{A}$ are delocalized under mild assumptions. The authors quantify misspecification effects: when the chosen dimension $d=r+k$ is too small ($k<0$), they obtain a lower bound on estimation error of order $\sqrt{\rho_n|k|}$, while for over-specification ($k>0$) they establish an upper bound showing consistency with a slower rate $\|\hat{\mathbf{X}}_{1:r+k}-\rho_n^{1/2}\mathbf{X}_{1:r+k}\|_{2,\infty} \lesssim \phi_n + \sqrt{k}\,r^2(\log n)^{5+6\gamma}/n^{1/4}$. The results extend to weighted and binary networks, with experiments demonstrating the predicted rates and showing consistency persists under over-specification. The work provides guidance for dimension selection in network embeddings and contributes a delocalization-toolkit for random-matrix-based spectral analysis in low-rank perturbations.

Abstract

As network data has become ubiquitous in the sciences, there has been growing interest in network models whose structure is driven by latent node-level variables in a (typically low-dimensional) latent geometric space. These "latent positions" are often estimated via embeddings, whereby the nodes of a network are mapped to points in Euclidean space so that "similar" nodes are mapped to nearby points. Under certain model assumptions, these embeddings are consistent estimates of the latent positions, but most such results require that the embedding dimension be chosen correctly, typically equal to the dimension of the latent space. Methods for estimating this correct embedding dimension have been studied extensive in recent years, but there has been little work to date characterizing the behavior of embeddings when this embedding dimension is misspecified. In this work, we provide theoretical descriptions of the effects of misspecifying the embedding dimension of the adjacency spectral embedding under the random dot product graph, a class of latent space network models that includes a number of widely-used network models as special cases, including the stochastic blockmodel. We consider both the case in which the dimension is chosen too small, where we prove estimation error lower-bounds, and the case where the dimension is chosen too large, where we show that consistency still holds, albeit at a slower rate than when the embedding dimension is chosen correctly.A range of synthetic data experiments support our theoretical results. Our main technical result, which may be of independent interest, is a generalization of earlier work in random matrix theory, showing that all non-signal eigenvectors of a low-rank matrix subject to additive noise are delocalized.

On the Effect of Misspecifying the Embedding Dimension in Low-rank Network Models

TL;DR

The paper addresses how adjacency spectral embedding (ASE) behaves when the embedding dimension is misspecified in latent-space network models, notably the random dot product graph (RDPG) and its binary and weighted variants. It develops a signal-plus-noise framework with and proves a key delocalization result: all non-signal eigenvectors of are delocalized under mild assumptions. The authors quantify misspecification effects: when the chosen dimension is too small (), they obtain a lower bound on estimation error of order , while for over-specification () they establish an upper bound showing consistency with a slower rate . The results extend to weighted and binary networks, with experiments demonstrating the predicted rates and showing consistency persists under over-specification. The work provides guidance for dimension selection in network embeddings and contributes a delocalization-toolkit for random-matrix-based spectral analysis in low-rank perturbations.

Abstract

As network data has become ubiquitous in the sciences, there has been growing interest in network models whose structure is driven by latent node-level variables in a (typically low-dimensional) latent geometric space. These "latent positions" are often estimated via embeddings, whereby the nodes of a network are mapped to points in Euclidean space so that "similar" nodes are mapped to nearby points. Under certain model assumptions, these embeddings are consistent estimates of the latent positions, but most such results require that the embedding dimension be chosen correctly, typically equal to the dimension of the latent space. Methods for estimating this correct embedding dimension have been studied extensive in recent years, but there has been little work to date characterizing the behavior of embeddings when this embedding dimension is misspecified. In this work, we provide theoretical descriptions of the effects of misspecifying the embedding dimension of the adjacency spectral embedding under the random dot product graph, a class of latent space network models that includes a number of widely-used network models as special cases, including the stochastic blockmodel. We consider both the case in which the dimension is chosen too small, where we prove estimation error lower-bounds, and the case where the dimension is chosen too large, where we show that consistency still holds, albeit at a slower rate than when the embedding dimension is chosen correctly.A range of synthetic data experiments support our theoretical results. Our main technical result, which may be of independent interest, is a generalization of earlier work in random matrix theory, showing that all non-signal eigenvectors of a low-rank matrix subject to additive noise are delocalized.
Paper Structure (16 sections, 33 theorems, 353 equations, 5 figures, 1 table)

This paper contains 16 sections, 33 theorems, 353 equations, 5 figures, 1 table.

Key Result

Lemma 2.1

Suppose that $\mathbf{A} = \rho_n\mathbf{X} \mathbf{X}^\top + \mathbf{E}$, and denote the $(r+k)$-dimensional ASE of $\mathbf{A}$ as $\mathbf{\hat{X}}_{1:r+k} = \mathbf{\hat{U}}_{1:r+k}|\mathbf{\hat{S}}|^{1/2}_{1:r+k}$. When $k<0$, On the other hand, when $k > 0$, suppose that there exists a sequence of $\mathbf{W}^* \in \mathbb{O}_r$ such that for some sequence $(\phi_n)_{n=1}^\infty$. Then the

Figures (5)

  • Figure 1: Estimation error of ASE in $(2,\infty)$-norm as a function of number of vertices $n$ for weighted networks as in Equation \ref{['eq:def:weightedDense']} under four different choices of noise models (a) normal (b) Laplace (c) exponential and (d) Poisson, under six choices of embedding dimension (colored lines). The correctly-specified embedding dimension ($r = 5, k= 0$; in gold) is display alongside five misspecified embedding dimensions ($r + k = 4, 6, 10, 20, 40$; orange, green, teal, purple and magenta, respectively). The black and gray dashed lines indicate, respectively, the convergence rate predicted for the correctly-specified setting and the setting where the embedding dimension is chosen too large. Error bars represent two standard errors of the mean.
  • Figure 2: Error in recovering the true latent positions in $(2,\infty)$-norm, as a function of embedding dimension for varying choices of network size $n$, indicated by line color. Both axes are on a logarithmic scale, and the dashed black line indicates the $\sqrt{k}$-like behavior predicted by our results. The black vertical line indicates the true embedding dimension, $\operatorname{rank} \mathbf{P} = 5$. Error bars indicate two standard errors of the mean.
  • Figure 3: Estimation error of ASE in $(2,\infty)$-norm as a function of network size for varying signal strength levels (indicated by line color) for six choices of embedding dimension. Subplots show behavior under correctly-specified embedding dimension ($r+k = r = 5$; top row, middle panel), as well as when the dimension is chosen too low ($r+k=4$; top row, left) and too high ($r + k = 6, 7, 10, 20$; top-row right panel and all three bottom row panels). The black (resp. gray) dashed line with a slope of $-1/2$ (resp. $-1/4$) indicates the predicted convergence rate of $n^{-1/2}$ (resp. $n^{-1/4}$) from Theorem \ref{['theorem:main_result:weighted']} for the correctly chosen (resp. incorrectly chosen) embedding dimension. Error bars represent two standard errors of the mean.
  • Figure 4: Average estimation error of ASE in $(2,\infty)$-norm as a function of number of vertices $n$ for different choices of embedding dimension $r+k$ (line colors) under the SBM described in Equation \ref{['eq:def:SBM']}. Both axes are on logarithmic scales, with dashed lines in black and grey indicating, respectively, the $n^{-1/2}$ and $n^{-1/4}$ rates predicted by our Conjecture \ref{['conj:maintext:deloc']} under the settings where the dimension is chosen, respectively, either correctly or too large. Error bars indicate two standard errors of the mean.
  • Figure 5: Estimation error of ASE in $(2,\infty)$-norm as a function of network size $n$ under the sparse binary model in Equation \ref{['eq:def:sparseDiri']} for varying sparsity $\rho_n = n^{-\gamma}$ ($\gamma=0,0.1,0.2,0.3,0.4,0.5$; indicated by line color) for varying choices of embedding dimension $r+k=4,5,6,7,10,20$. The $n^{-1/2}$ rate suggested by our theory for the well-specified case ($r+k=5$; top row, middle plot) is indicated by a black dashed line. The $n^{-1/4}$ rate suggest by our theory when the embedding dimension is too large ($r+k>5$; top-row right-hand plot and bottom row) is indicated by a grey dashed line. Errors bars indicate two standard errors of the mean.

Theorems & Definitions (68)

  • Definition 1: Adjacency Spectral Embedding
  • Definition 2: Random Dot Product Graph
  • Remark 1
  • Lemma 2.1: Misspecified Model Bounds
  • Theorem 2.2: lyzinski_community_2016, Theorem 5
  • Definition 3: GRDPG_original, Definition 2
  • Theorem 2.3: GRDPG_original, Theorem 1
  • Theorem 2.4: levin_recovering_nodate, Theorem 6
  • Remark 2
  • Theorem 3.1
  • ...and 58 more