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.

renshaw.r

Vejam código da função ssa.d do pacote Gillespie SSA e o artigo de apresentação do pacote em Pineda-Kirsch (2008)

simu_discreta.r
simu_discreta_min_max.r

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~~