Reducing Spreading Processes on Networks to Markov Population Models
Pith reviewed 2026-05-25 14:15 UTC · model grok-4.3
The pith
Lumping via node partitioning reduces any epidemic model on a network to a Markov population model.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Lumping can be used to reduce any epidemic model to a Markov Population Model. A novel lumping scheme based on a partitioning of the nodes is proposed. By imposing different types of counting abstractions, coarse-grained Markov models with a natural MPM representation are obtained that approximate the original systems. This transfer makes the rich pool of approximation techniques developed for MPMs available for the computational analysis of complex networks' dynamics.
What carries the argument
A lumping scheme that partitions the nodes of the network and then imposes counting abstractions to produce a Markov population model representation of the process.
If this is right
- The reduced models approximate the original network dynamics with controllable error.
- Existing approximation techniques for Markov population models become directly applicable to network spreading processes.
- Accuracy of the approximation trades off against the size of the lumped state space and the type of counting abstraction used.
- Computational analysis of spreading phenomena on large networks becomes feasible through these smaller models.
Where Pith is reading between the lines
- The same partitioning approach might apply to non-epidemic stochastic processes on networks.
- Structure-aware or adaptive partitions could be tested to improve accuracy on heterogeneous graphs.
- Moment-closure methods already used for MPMs could be combined with the lumped models for further scaling.
Load-bearing premise
The chosen node partitioning and counting abstraction must preserve enough of the original interaction structure so that the resulting MPM still approximates the true dynamics on the full network.
What would settle it
A concrete network and epidemic process in which trajectories from the reduced Markov population model diverge substantially from direct stochastic simulations of the original full-state process.
Figures
read the original abstract
Stochastic processes on complex networks, where each node is in one of several compartments, and neighboring nodes interact with each other, can be used to describe a variety of real-world spreading phenomena. However, computational analysis of such processes is hindered by the enormous size of their underlying state space. In this work, we demonstrate that lumping can be used to reduce any epidemic model to a Markov Population Model (MPM). Therefore, we propose a novel lumping scheme based on a partitioning of the nodes. By imposing different types of counting abstractions, we obtain coarse-grained Markov models with a natural MPM representation that approximate the original systems. This makes it possible to transfer the rich pool of approximation techniques developed for MPMs to the computational analysis of complex networks' dynamics. We present numerical examples to investigate the relationship between the accuracy of the MPMs, the size of the lumped state space, and the type of counting abstraction.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that lumping via node partitioning and various counting abstractions can reduce any stochastic epidemic process on a network to a Markov Population Model (MPM) that approximates the original dynamics, thereby allowing transfer of MPM analysis techniques to network spreading models. Numerical examples are used to explore how approximation accuracy relates to lumped state-space size and abstraction type.
Significance. If the reduction is accompanied by explicit lumpability conditions and error controls, the approach would enable scalable analysis of large-network epidemic models by connecting them to existing MPM approximation methods. The numerical investigation of accuracy versus granularity is a useful empirical contribution.
major comments (2)
- [Abstract] Abstract: the assertion that the scheme reduces 'any epidemic model' to a closed MPM via arbitrary node partitioning plus counting abstraction is not supported by a lumpability condition; in general networks the infection rate for a susceptible node in one block depends on its specific neighbors in infected blocks, so two micro-states with identical block counts can have different total rates unless the partition is equitable.
- Proposed lumping scheme (as described after the abstract): no formal statement is given of the conditions under which the aggregated process remains Markov with transition rates depending only on compartment counts per block, nor are error bounds or convergence results supplied for the approximation under different counting abstractions.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for highlighting important points regarding the formal presentation of the lumping scheme. We address each major comment below.
read point-by-point responses
-
Referee: [Abstract] Abstract: the assertion that the scheme reduces 'any epidemic model' to a closed MPM via arbitrary node partitioning plus counting abstraction is not supported by a lumpability condition; in general networks the infection rate for a susceptible node in one block depends on its specific neighbors in infected blocks, so two micro-states with identical block counts can have different total rates unless the partition is equitable.
Authors: We agree that exact lumpability holds only for special partitions (e.g., equitable ones). The manuscript already qualifies the reduction as approximate ('approximate the original systems'), but the abstract could be clearer on this point. We will revise the abstract to state explicitly that the resulting MPM is an approximation for arbitrary partitions and add a brief remark on lumpability conditions in Section 3. revision: yes
-
Referee: [—] Proposed lumping scheme (as described after the abstract): no formal statement is given of the conditions under which the aggregated process remains Markov with transition rates depending only on compartment counts per block, nor are error bounds or convergence results supplied for the approximation under different counting abstractions.
Authors: The core contribution is the construction of approximating MPMs via node partitioning and counting abstractions, supported by numerical examples relating accuracy to state-space size. We will add a formal paragraph stating the (standard) lumpability condition under which the aggregated process is exactly Markov. Theoretical error bounds and convergence rates are not derived in the present work; the paper relies on empirical evaluation. We will expand the discussion to note this limitation and identify it as future work. revision: partial
Circularity Check
No significant circularity in the proposed lumping reduction
full rationale
The paper presents a methodological reduction: node partitioning plus counting abstractions are used to obtain approximate coarse-grained Markov models that admit an MPM representation. The abstract and description frame this as an approximation technique that transfers existing MPM analysis tools, without any self-definitional equations, fitted parameters renamed as predictions, or load-bearing self-citations. The derivation relies on standard lumping applied to network epidemic processes and is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we demonstrate that lumping can be used to reduce any epidemic model to a Markov Population Model (MPM). Therefore, we propose a novel lumping scheme based on a partitioning of the nodes. By imposing different types of counting abstractions...
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
q(y,y′) = 1/|L⁻¹(y)| ∑_{x∈L⁻¹(y)} ∑_{x′∈L⁻¹(y′)} q(x,x′)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
G. E. Allen and C. Dytham. An efficient method for stochastic simulation of biological populations in continuous time. Biosystems, 98(1):37–42, 2009
work page 2009
-
[2]
A.-L. Barab´ asi.Network science. Cambridge university press, 2016
work page 2016
-
[3]
L. Bortolussi. Hybrid behaviour of Markov population models. Information and Computation, 247:37–86, 2016
work page 2016
-
[4]
L. Bortolussi, J. Hillston, D. Latella, and M. Massink. Continuous approximation of collective system behaviour: A tutorial. Performance Evaluation, 70(5):317–349, 2013
work page 2013
- [5]
-
[6]
Y. Cao, D. T. Gillespie, and L. R. Petzold. Efficient step size selection for the tau-leaping simulation method. The Journal of chemical physics , 124(4)
-
[7]
L. Cardelli, M. Tribastone, M. Tschaikowski, and A. Vandin. Erode: a tool for the evaluation and reduction of ordinary differential equations. In International Conference on Tools and Algorithms for the Construction and Analysis of Systems , pages 310–328. Springer, 2017
work page 2017
-
[8]
W. Cota and S. C. Ferreira. Optimized gillespie algorithms for the simulation of markovian epidemic processes on large and heterogeneous networks. Computer Physics Communications, 219:303–312, 2017
work page 2017
-
[9]
K. Devriendt and P. Van Mieghem. Unified mean-field framework for susceptible- infected-susceptible epidemics on networks, based on graph partitioning and the isoperimetric inequality. Physical Review E, 96(5):052314, 2017
work page 2017
-
[10]
C. Gan, X. Yang, W. Liu, Q. Zhu, and X. Zhang. Propagation of computer virus under human intervention: a dynamical model. Discrete Dynamics in Nature and Society, 2012, 2012
work page 2012
-
[11]
M. Girvan and M. E. Newman. Community structure in social and biological networks. Proceedings of the national academy of sciences, 99(12):7821–7826, 2002
work page 2002
-
[12]
J. P. Gleeson. High-accuracy approximation of binary-state dynamics on networks. Physical Review Letters, 107(6):068701, 2011
work page 2011
-
[13]
J. P. Gleeson. Binary-state dynamics on complex networks: Pair approximation and beyond. Physical Review X, 3(2):021004, 2013
work page 2013
-
[14]
A. Goltsev, F. De Abreu, S. Dorogovtsev, and J. Mendes. Stochastic cellular automata model of neural networks. Physical Review E, 81(6):061921, 2010
work page 2010
-
[15]
J. Goutsias and G. Jenkinson. Markovian dynamics on complex reaction networks. Physics Reports, 529(2):199–264, 2013
work page 2013
-
[16]
P. Goyal and E. Ferrara. Graph embedding techniques, applications, and perfor- mance: A survey. Knowledge-Based Systems, 151:78–94, 2018
work page 2018
-
[17]
R. Grima. An effective rate equation approach to reaction kinetics in small volumes: Theory and application to biochemical reactions in nonequilibrium steady-state conditions. The Journal of Chemical Physics , 133(3):035101, 2010
work page 2010
-
[18]
R. Grima. A study of the accuracy of moment-closure approximations for stochastic chemical kinetics. The Journal of chemical physics , 136(15):04B616, 2012
work page 2012
-
[19]
G. Großmann, C. Kyriakopoulos, L. Bortolussi, and V. Wolf. Lumping the ap- proximate master equation for multistate processes on complex networks. In Qest, pages 157–172. Springer, 2018
work page 2018
-
[20]
Rejection-Based Simulation of Stochastic Spreading Processes on Complex Networks
G. Großmann and V. Wolf. Rejection-based simulation of stochastic spreading processes on complex networks. arXiv preprint arXiv:1812.10845 , 2018. Reducing Spreading Processes on Networks to Markov Population Models 17
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[21]
A. Hagberg, P. Swart, and D. S Chult. Exploring network structure, dynamics, and function using networkx. Technical report, Los Alamos National Lab.(LANL), Los Alamos, NM (United States), 2008
work page 2008
-
[22]
T. A. Henzinger, M. Mateescu, and V. Wolf. Sliding window abstraction for infinite markov chains. In International Conference on Computer Aided Verification, pages 337–352. Springer, 2009
work page 2009
-
[23]
P. Holme. Shadows of the susceptible-infectious-susceptible immortality transition in small networks. Physical Review E, 92(1):012804, 2015
work page 2015
-
[24]
M. J. Keeling and P. Rohani. Modeling infectious diseases in humans and animals . Princeton University Press, 2011
work page 2011
-
[25]
W. R. KhudaBukhsh, A. Auddy, Y. Disser, and H. Koeppl. Approximate lumpabil- ity for markovian agent-based models using local symmetries. arXiv:1804.00910
work page internal anchor Pith review Pith/arXiv arXiv
-
[26]
I. Z. Kiss, J. C. Miller, and P. L. Simon. Mathematics of epidemics on networks: from exact to approximate models. Forthcoming in Springer TAM series , 2016
work page 2016
- [27]
-
[28]
C. Kyriakopoulos, G. Grossmann, V. Wolf, and L. Bortolussi. Lumping of degree- based mean-field and pair-approximation equations for multistate contact pro- cesses. Physical Review E, 97(1):012301, 2018
work page 2018
- [29]
-
[30]
M. L´ opez-Garc´ ıa. Stochastic descriptors in an sir epidemic model for heterogeneous individuals in small networks. Mathematical biosciences, 271:42–61, 2016
work page 2016
-
[31]
M. Mateescu, V. Wolf, F. Didier, and T. Henzinger. Fast adaptive uniformisation of the chemical master equation. IET systems biology , 4(6):441–452, 2010
work page 2010
-
[32]
R. M. May and N. Arinaminpathy. Systemic risk: the dynamics of model banking systems. Journal of the Royal Society Interface , 7(46):823–838, 2009
work page 2009
-
[33]
M. Moslonka-Lefebvre, M. Pautasso, and M. J. Jeger. Disease spread in small- size directed networks: epidemic threshold, correlation between links to and from nodes, and clustering. Journal of theoretical biology, 260(3):402–411, 2009
work page 2009
-
[34]
T. W. Ng, G. Turinici, and A. Danchin. A double epidemic model for the sars propagation. BMC Infectious Diseases, 3(1):19, 2003
work page 2003
-
[35]
M. Pautasso, M. Moslonka-Lefebvre, and M. J. Jeger. The number of links to and from the starting node as a predictor of epidemic size in small-size directed networks. Ecological Complexity, 7(4):424–432, 2010
work page 2010
-
[36]
M. Porter and J. Gleeson. Dynamical systems on networks: A tutorial , volume 4. Springer, 2016
work page 2016
-
[37]
H. S. Rodrigues. Application of sir epidemiological model: new trends. arXiv:1611.02565, 2016
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[38]
H. S. Rodrigues, M. T. T. Monteiro, and D. F. Torres. Dynamics of dengue epi- demics when using optimal control. Mathematical and Computer Modelling , 52(9- 10):1667–1673, 2010
work page 2010
-
[39]
D. Schnoerr, G. Sanguinetti, and R. Grima. Approximation and inference meth- ods for stochastic biochemical kinetics - a tutorial review. Journal of Physics A , 51:169501, 2018
work page 2018
-
[40]
P. L. Simon, M. Taylor, and I. Z. Kiss. Exact epidemic models on graphs using graph-automorphism driven lumping. Journal of mathematical biology , 62(4):479– 508, 2011
work page 2011
-
[41]
A. Singh and J. P. Hespanha. Stochastic hybrid systems for studying biochemical processes. Royal Society A, 368(1930):4995–5011, 2010. 18 G. Großmann and L. Bortolussi
work page 1930
-
[42]
M. Soltani, C. A. Vargas-Garcia, and A. Singh. Conditional moment closure schemes for studying stochastic dynamics of genetic circuits. IEEE transactions on biomedical circuits and systems , 9(4):518–526, 2015
work page 2015
-
[43]
G. St-Onge, J.-G. Young, L. H´ ebert-Dufresne, and L. J. Dub´ e. Efficient sampling of spreading processes on complex networks using a composition and rejection algorithm. Computer Physics Communications , 2019
work page 2019
-
[44]
N. G. Van Kampen. Stochastic processes in physics and chemistry , volume 1. Elsevier, 1992
work page 1992
-
[45]
P. Van Mieghem, J. Omic, and R. Kooij. Virus spread in networks. IEEE/ACM Transactions on Networking, 17(1):1–14, 2009
work page 2009
-
[46]
J. A. Ward and J. Evans. A general model of dynamics on networks with graph automorphism lumping. In International Workshop on Complex Networks and their Applications, pages 445–456. Springer, 2018
work page 2018
-
[47]
D. J. Watts and P. S. Dodds. Influentials, networks, and public opinion formation. Journal of consumer research, 34(4):441–458, 2007
work page 2007
- [48]
-
[49]
X. Wei, N. C. Valler, B. A. Prakash, I. Neamtiu, M. Faloutsos, and C. Faloutsos. Competing memes propagation on networks: A network science perspective. IEEE Journal on Selected Areas in Communications , 31(6):1049–1060, 2013
work page 2013
-
[50]
L. Zhao, H. Cui, X. Qiu, X. Wang, and J. Wang. Sir rumor spreading model in the new media age. Physica A, 392(4):995–1003, 2013
work page 2013
-
[51]
L. Zhao, J. Wang, Y. Chen, Q. Wang, J. Cheng, and H. Cui. Sihr rumor spreading model in social networks. Physica A, 391(7):2444–2453, 2012. Reducing Spreading Processes on Networks to Markov Population Models 19 Appendix 8 Direct Construction of MPMs Here, we prosper a way of directly deriving the lumped MPMs from the contact network without building the ...
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.