Two-sample inference for high-dimensional Markov networks

Abstract

Markov networks are frequently used in sciences to represent conditional independence relationships underlying observed variables arising from a complex system. It is often of interest to understand how an underlying network differs between two conditions. In this paper, we develop methods for comparing a pair of high-dimensional Markov networks where we allow the number of observed variables to increase with the sample sizes. By taking the density ratio approach, we are able to learn the network difference directly and avoid estimating the individual graphs. Our methods are thus applicable even when the individual networks are dense as long as their difference is sparse. We prove finite-sample Gaussian approximation error bounds for the estimator we construct under significantly weaker assumptions than are typically required for model selection consistency. Furthermore, we propose bootstrap procedures for estimating quantiles of a max-type statistics based on our estimator, and show how they can be used to test the equality of two Markov networks or construct simultaneous confidence intervals. The performance of our methods is demonstrated through extensive simulations. The scientific usefulness is illustrated with an analysis of a new fMRI dataset.

arXiv comment81 pages, 18 figures, 10 tables

References in corpus (10)

• Confidence Intervals and Hypothesis Testing for High-Dimensional Regression
• Estimating divergence functionals and the likelihood ratio by convex risk minimization
• Least squares after model selection in high-dimensional sparse models
• Gaussian and bootstrap approximations for high-dimensional U-statistics and their applications
• Change-Point Estimation in High-Dimensional Markov Random Field Models
• Testing for Differences in Gaussian Graphical Models: Applications to Brain Connectivity
• Simultaneous Inference for Pairwise Graphical Models with Generalized Score Matching
• Post-Regularization Inference for Time-Varying Nonparanormal Graphical Models
• Joint Estimation and Inference for Data Integration Problems based on Multiple Multi-layered Gaussian Graphical Models
• Inter-Subject Analysis: Inferring Sparse Interactions with Dense Intra-Graphs

Cited by in corpus (3)

• Estimating Differential Latent Variable Graphical Models with Applications to Brain Connectivity
• Constrained High Dimensional Statistical Inference
• Posterior Ratio Estimation of Latent Variables