simulate.true: Simulate a birth-shift-death-process

In jasonxu90/bdsem: Analysis of discretely observed Birth-Death-Shift processes

Description

This function simulates one realization of a birth-shift-death CTMC with per-particle birth, death, and shift rates lambda, nu, and mu.

Usage

  simulate.true(t.end, lam, v, mu, initNum)

Arguments

t.end	A number giving the length of time for the simulation
lam	Per-particle birth rate
v	Per-particle shift rate
mu	Per-particle death rate
initNum	An integer, initial total number of particles

Details

The process begins with initNum particles, each assigned a random location indexed by a number between 1 and 100,000. Simulation occurs until time t.end. The process returns an ordered pair giving the number of initial locations (indices) still occupied, followed by the number of new indices present, by t.end.

Value

A pair of integers giving the number of initial indices still present followed by the number of new indices present in the population. Otherwise, returns '999' as an error code

Examples

simulate.true(2,.2,.12,.15,10)

jasonxu90/bdsem documentation built on May 18, 2019, 5:54 p.m.