Tree-structured Markov random fields with Poisson marginal distributions
Benjamin Côté, Hélène Cossette, Etienne Marceau
TL;DR
This work addresses modeling counts with tractable dependence by introducing a tree-structured multivariate Poisson MRF that preserves Poisson marginals with parameter $\lambda$ while encoding dependence along edges through $\alpha_e$. The authors provide a stochastic construction based on binomial thinning that yields explicit joint PMF and PGF, enabling exact sampling and efficient computation of the sum distribution, which is shown to be compound Poisson with a tree-structured secondary distribution. They establish positive dependence along edges, derive covariance along paths, and prove ordering results (supermodular and convex) that describe how increasing dependence affects the aggregate $M=\sum_v N_v$ under fixed trees, with no general order across different topologies. The framework supports fast FFT-based computations of the distribution of $M$ and the ordinary generating function of expected allocations (OGFEA), facilitating risk-sharing and allocation analyses for high-dimensional count data. The work also highlights the impact of tree topology on dependence and suggests directions for extending the approach to general graphs and topology-aware analyses of centrality and ordering.
Abstract
A new family of tree-structured Markov random fields for a vector of discrete counting random variables is introduced. According to the characteristics of the family, the marginal distributions of the Markov random fields are all Poisson with the same mean, and are untied from the strength or structure of their built-in dependence. This key feature is uncommon for Markov random fields and most convenient for applications purposes. The specific properties of this new family confer a straightforward sampling procedure and analytic expressions for the joint probability mass function and the joint probability generating function of the vector of counting random variables, thus granting computational methods that scale well to vectors of high dimension. We study the distribution of the sum of random variables constituting a Markov random field from the proposed family, analyze a random variable's individual contribution to that sum through expected allocations, and establish stochastic orderings to assess a wide understanding of their behavior.
