User:Carl Boettiger/Notebook/Stochastic Population Dynamics/2010/05/26
Stochastic Population Dynamics | <html><img src="/images/9/94/Report.png" border="0" /></html> Main project page <html><img src="/images/c/c3/Resultset_previous.png" border="0" /></html>Previous entry<html> </html>Next entry<html><img src="/images/5/5c/Resultset_next.png" border="0" /></html> |
Rethinking Beetles Noise
and the larval fluctuations are essentially:
where [math]\displaystyle{ \beta_i }[/math] is the intrinsic noise of the age class. Hence a class i which propagates large noise to another class j has a large [math]\displaystyle{ \partial_i f_j }[/math]. If this term is a linear transition [math]\displaystyle{ \lambda X_i }[/math], then the same term appears in [math]\displaystyle{ f_i }[/math] and hence damps the noise [math]\displaystyle{ \sigma_i^2 }[/math] and cancels out. Hence noise must propagate into a class through nonlinear transition rates OR through an asymmetry in the transition (i.e. the c_1, c_2 large noise example in the generalized crowley).
Adding age delay to beetle dynamicsCurrent formulation has used exponential waiting times between stages. By subdividing the classes (increasing the system dimension size) and creating single jump within-state transitions, these become gamma-distributed waiting times. Adding ten steps to each phase and a little parameter fiddling introduces sustained oscillations. (Version-stable code).
Xo[1] = 100
and now we have noise in oscillatory, gamma-waiting model:
Misc / Code notes
|