Implementação em R do código do Renshaw e outras alternativas, para simular um processo estocástico de nascimentos apenas, e verificar as previsões analíticas do modelo.
Vejam código da função ssa.d do pacote Gillespie SSA e o artigo de apresentação do pacote em Pineda-Kirsch (2008)
Não consigo colocar os gráficos aqui no wiki, mas é só rodar os scripts que sai lá. Como esperado, parece que quando aumenta n0 , as escadinhas vão se tornando cada vez mais uma curva que se aproxima da curva esperada, e por isso uma curva de crescimento exponencial (sem estocasticidade). Já quando a taxa de crescimento (lambda) é cada vez maior, as escadinhas também vão se tornando curvas, mas não se aproximam da curva esperada pelo cálculo da esperança.
~~DISCUSSION~~