Structure Learning in Graphical Models from Indirect Observations

preprint2022arXivOpen access

Abstract

This paper considers learning of the graphical structure of a pp-dimensional random vector X∈RpX \in R^p using both parametric and non-parametric methods. Unlike the previous works which observe xx directly, we consider the indirect observation scenario in which samples yy are collected via a sensing matrix A∈Rd×pA \in R^{d\times p}, and corrupted with some additive noise ww, i.e, Y=AX+WY = AX + W. For the parametric method, we assume XX to be Gaussian, i.e., x∈Rp∼N(μ,Σ)x\in R^p\sim N(μ, Σ) and Σ∈Rp×pΣ\in R^{p\times p}. For the first time, we show that the correct graphical structure can be correctly recovered under the indefinite sensing system (d<pd < p) using insufficient samples (n<pn < p). In particular, we show that for the exact recovery, we require dimension d=Ω(p0.8)d = Ω(p^{0.8}) and sample number n=Ω(p0.8log⁡3p)n = Ω(p^{0.8}\log^3 p). For the nonparametric method, we assume a nonparanormal distribution for XX rather than Gaussian. Under mild conditions, we show that our graph-structure estimator can obtain the correct structure. We derive the minimum sample number nn and dimension dd as n≳(deg)4log⁡4nn\gtrsim (deg)^4 \log^4 n and d≳p+(deg⋅log⁡(d−p))β/4d \gtrsim p + (deg\cdot\log(d-p))^{β/4}, respectively, where deg is the maximum Markov blanket in the graphical model and β>0β> 0 is some fixed positive constant. Additionally, we obtain a non-asymptotic uniform bound on the estimation error of the CDF of XX from indirect observations with inexact knowledge of the noise distribution. To the best of our knowledge, this bound is derived for the first time and may serve as an independent interest. Numerical experiments on both real-world and synthetic data are provided confirm the theoretical results.

Explore this paper’s authors, topics and related work.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Reviews 0

Write a reviewWrite

No reviews yet.

Discussion 0

Add a commentWrite

No comments yet.