User:Carl Boettiger/Notebook/Stochastic Population Dynamics/2010/06/10
From OpenWetWare
Stochastic Population Dynamics  Main project page Previous entry Next entry 
QuestionsCurrent population dynamics model macroscopic equations: Interpreting μ terms as death events, b as birth, and the a terms as transitions between age classes gives a natural form for the master equation. Divergent Solutions to the linear noise approximation?Breaking the age classes into k subclasses makes waiting time within an age class gamma instead of exponential, which better matches observation. Also allows sustained oscillations. The numerical solution to the linear noise approximation shows divergent variance dynamics. Is this a numerical artifact? Differences between twostep (correlated events) and two single stepsCan we distinguish between a process where an element leaves one population and enters another (i.e. a transitions from one age class to the next) and a model where death rate in one class is equal to the birth rate in another? The linear noise approximation for a twostep process obeys a different equation than the onestep: where T is a symmetric matrix with offdiagonal elements i,j τ_{ij}(p) + τji(p) and diagonal elements chosen so that rows sum to zero. This approximation is implemented in noise_approx.R How different can this be from simply implementing all transitions as singlesteps? So far the difference seems small in my examples:
