* [[ecovirt:roteiro:den_ind:di_edrcmdr|{{:ecovirt:logorcmdr01.png?20|}}]]
* [[ecovirt:roteiro:den_ind:di_edr|{{:ecovirt:rlogo.png?20|}}]]
* [[ecovirt:roteiro:den_ind:di_edr_passo|{{:ecovirt:prompt.png?20|}}]]
====== Crescimento denso-independente com estocasticidade demográfica - Roteiro no R ======
{{section>ecovirt:roteiro:den_ind:di_ed_base#crescimento_denso-independente_com_estocasticidade_demografica}}
====A distribuição binomial ====
{{section>ecovirt:roteiro:avisos#roteiro_r}}
Vamos nos familiarizar com a ideia de distribuições de probabilidades, calculando valores da distribuição binomial. Abaixo vai o código no R que fará os gráficos de função de massa da distribuição binomial:
x <- 0:10 # para desenhar o eixo x do gráfico
plotDistr(x, dbinom(x, size=10, prob=0.5), xlab="Number of Successes",
ylab="Probability Mass",
main="Binomial Distribution: Binomial trials=10, Probability of success=0.5",
discrete=TRUE)
Neste código, $dbinom$ é a função que vai calcular a densidade de probabilidade de sucesso de cada número de tentativas. $size$ é o número de tentativas, no nosso caso $N_0$, e $prob$ é a probabilidade de sucesso a cada tentativa, no caso a probabilidade de cada indivíduo sobreviver.
Faça o gráfico de nosso exemplo:
- em $size$ coloque $2$
- em $prob$ coloque $0.5$
- execute a função
{{section>ecovirt:roteiro:den_ind:di_ed_base#a_distribuicao_binomial}}
==== Simulação no computador ====
{{section>ecovirt:roteiro:den_ind:di_ed_base#simulacao_no_computador}}
Vamos utilizar a função $estDem$ do pacote **EcoVirtual** no R.
{{section>ecovirt:roteiro:den_ind:di_ed_base#parametros}}
==== Distribuição dos tamanhos populacionais ====
Vamos inspecionar a distribuição dos tamanhos populacionais até o tempo $t=2$. Para isso temos que guardar os resultados em um objeto do R. Execute as simulações com os valores:
# guarde o resultado da função no objeto sim1
tmax = 2
nsim = 1000
N0 = 20
b = 0
d = 0.693
Execute a função, e os resultados estão gravados em uma lista de 1000 tabelas no R, que chamamos ''sim1''. Cada tabela tem os tempos em que a população perdeu um indivíduo e o tamanho da população a partir daquele momento, até o tempo máximo estipulado nas opções.
Para ver a primeira tabela, copie e execute o comando abaixo no R.
sim1[[1]]
Inspecione outras tabelas: execute de novo o comando, mudando o índice numérico que está entre os colchetes duplos, para qualquer valor entre um e mil.
{{section>ecovirt:roteiro:den_ind:di_ed_base#distribuicao_dos_tamanhos_populacionais}}
===== Nascimentos e mortes =====
{{section>ecovirt:roteiro:den_ind:di_ed_base#nascimentos_e_mortes}}
===== Para saber mais =====
{{section>ecovirt:roteiro:den_ind:di_ed_base#para_saber_mais}}
{{tag>R uma_população crescimento_exponencial tempo_discreto estocasticidade_demográfica}}