* [[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}}