Bayesian Gaussian Mixture Modeling for Symmetric Matrix Variate Data

2026-08-06

AI Summary

This study addresses the challenges in modeling individual activity networks arising from insufficient data granularity and the high-dimensional complexity of symmetric matrix-variate covariates. The authors propose STRUCTURED, a novel model that characterizes GPS-derived social activity matrices using a mixture of symmetric matrix-variate normal distributions to uncover associations between spatiotemporal co-occurrence patterns and demographic attributes. A key innovation lies in imposing exchangeability constraints to parameterize the column precision matrix as a polynomial function of the row precision matrix, substantially reducing parameter dimensionality. Variable-order Bayesian inference is achieved via reversible-jump Markov chain Monte Carlo. Empirical evaluation demonstrates that STRUCTURED-RJ outperforms baseline methods under data sparsity, while STRUCTURED-FP excels with larger samples. The analysis further identifies local crime environments and youth employment density as pivotal drivers of heterogeneity in activity networks.

Abstract

Statistical inference on individual activity networks has been a historically difficult task due to the lack of available data at the appropriate granularity and the complexity of modeling individual mobility patterns. The recent availability of GPS data from individual devices, combined with highly detailed demographic information, suggests that one of these challenges can now be addressed. We introduce a new model which we call the Symmetric Matrix-Variate Normal Mixture Model (STRUCTURED) to estimate how demographic traits influence changes in human activity networks, using sociomatrices that capture the probabilistic spatial overlap between individuals over time. We exploit the commutativity constraint inherent in the symmetric matrix-variate normal distribution to parameterize the column precision matrix as a polynomial of the row precision matrix, reducing the effective parameter space by an order of magnitude. We develop two variants of STRUCTURED: STRUCTURED-FP, which estimates the full polynomial, and STRUCTURED-RJ, which uses reversible-jump MCMC to select a reduced-order parameterization. Simulation studies demonstrate that STRUCTURED-RJ outperforms existing methods in sparse-data regimes, whereas STRUCTURED-FP is preferred when sample sizes are large. We apply the model to GPS-derived sociomatrices of 293 individuals in King County, WA, finding that local crime environments and youth employment density are the dominant demographic factors explaining variation in weekly activity overlap patterns.

Problem

• activity networks
• symmetric matrix-variate data
• demographic traits
• human mobility
• statistical inference

Innovation

• Symmetric Matrix-Variate Normal
• Mixture Model
• Reversible-Jump MCMC
• Activity Networks
• Parameter Reduction