====== Pure Birth Process ======
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.

===Algoritmo descrito no cap.2 do Renshaw===
==Descrição==
==Código em R==
{{:laporta:renshaw.R|}}
==Resultados==

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

====Código em R==
{{:reservada:simu_discreta.r|}}
\\
{{:mandai:simu_discreta_min_max.r|}}
==Resultados==
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~~