Back

A fast numerical integration scheme for clonal expansion processes on graphs

Paterson, C. G. B.; Gao, M.; Hellier, J.; Luebeck, G.; Wedge, D. C.; Bozic, I.

2025-01-15 cancer biology
10.1101/2025.01.13.632690 bioRxiv
Show abstract

Compound birth-death processes are widely used to model the age-incidence curves of many cancers [1]. There are efficient schemes for directly computing the relevant probability distributions in the context of linear multi-stage clonal expansion (MSCE) models [2]. However, these schemes have not been generalised to models on arbitrary graphs, forcing the use of either full stochastic simulations or mean-field approximations, which can become inaccurate at late times or old ages [3, 4]. Here, we present a numerical integration scheme for directly computing survival probabilities of a first-order birth-death process on an arbitrary directed graph, without the use of stochastic simulations. As a concrete application, we show that this new numerical method can be used to infer the parameters of an example graphical model from simulated data.

Published in PLOS Computational Biology (predicted rank #2) · training set

Matching journals

The top 4 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.