Abstract:The Gillespie algorithm provides statistically exact methods for simulating stochastic dynamics modeled as interacting sequences of discrete events including systems of biochemical reactions or earthquake occurrences, networks of queuing processes or spiking neurons, and epidemic and opinion formation processes on social networks. Empirically, the inter-event times of various phenomena obey long-tailed distributions. The Gillespie algorithm and its variants either assume Poisson processes (i.e., exponentially distributed inter-event times), use particular functions for time courses of the event rate, or work for non-Poissonian renewal processes, including the case of long-tailed distributions of inter-event times, but at a high computational cost. In the present study, we propose an innovative Gillespie algorithm for renewal processes on the basis of the Laplace transform. The algorithm makes use of the fact that a class of point processes is represented as a mixture of Poisson processes with different event rates. The method is applicable to multivariate renewal processes whose survival function of inter-event times is completely monotone. It is an exact algorithm and works faster than a recently proposed Gillespie algorithm for general renewal processes, which is exact only in the limit of infinitely many processes. We also propose a method to generate sequences of event times with a tunable amount of positive correlation between inter-event times. We demonstrate our algorithm with exact simulations of epidemic processes on networks, finding that a realistic amount of positive correlation in inter-event times only slightly affects the epidemic dynamics.

Probabilistic sampling of finite renewal processes

Sampling Point Processes on Stable Unbounded Regions and Exam Simulation of Queues

Statistics of a Large Number of Renewals in Equilibrium and Non-Equilibrium Renewal Processes

Replicating a renewal process at random times

Several Properties of a Nonstandard Renewal Counting Process and Their Applications

APPLICATION OF RENEWAL REWARD PROCESSES IN HOMOGENEOUS DISCRETE MARKOV CHAIN

Non-equilibrium dynamics of stochastic point processes with refractoriness

A Gillespie Algorithm for Non-Markovian Stochastic Processes

A fractional generalization of the Poisson processes

Renewal and Regenerative Processes: short review

Delayed Renewal Process with Uncertain Random Inter-Arrival Times.

Fuzzy Random Renewal Process and Renewal Reward Process

A New Fractal Point Process for Modeling Self-Similar Traffic

Aging Renewal Point Processes and Exchangeability of Event Times

Renewal Processes Represented as Doubly Stochastic Poisson Processes

Distribution theory for dependent renewal–reward processes and their first-passage times using saddlepoint and residue expansions

Web Renewal Counting Processes and Their Applications in Insurance

On renewal theory for cluster processes

Stochastic Renewal Equation for the Waiting Time Statistics for the First Occurrence of a Specific Sequence of States Successively Visited by an Alternating Renewal Process

Sampling using Adaptive Regenerative Processes