apoio:exercicios
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anteriorRevisão anteriorPróxima revisão | Revisão anterior | ||
apoio:exercicios [2017/10/09 04:52] – [Coelhos na Austrália] prado | apoio:exercicios [2024/01/12 10:39] (atual) – edição externa 127.0.0.1 | ||
---|---|---|---|
Linha 1: | Linha 1: | ||
====== Desafios ====== | ====== Desafios ====== | ||
- | Terminou a lista de exercícios? | + | Alguns |
===== Coelhos na Austrália ===== | ===== Coelhos na Austrália ===== | ||
Linha 60: | Linha 60: | ||
{{https:// | {{https:// | ||
- | No fim da última era glacial | + | No fim da última era glacial |
Neste exercício, adaptado de Akçakaya et al (1999), vamos avaliar alguns índices de viabilidade da população inicial, considerando os riscos criados pela estocasticidade demográfica e ambiental. Para isso vamos usar os seguintes parâmetros, | Neste exercício, adaptado de Akçakaya et al (1999), vamos avaliar alguns índices de viabilidade da população inicial, considerando os riscos criados pela estocasticidade demográfica e ambiental. Para isso vamos usar os seguintes parâmetros, | ||
Linha 68: | Linha 68: | ||
* Desvio-padrão da taxa intríseca de crescimento: | * Desvio-padrão da taxa intríseca de crescimento: | ||
* Probabilidade de sobrevivência anual: 0,921 | * Probabilidade de sobrevivência anual: 0,921 | ||
+ | |||
+ | |||
+ | ==== Para entender as simulações passo a passo==== | ||
+ | |||
+ | Copie o código abaixo em um script no R e acompanhe a explicação dos professores: | ||
+ | |||
+ | < | ||
+ | ## Tamanhos populacionais (Spencer & Lensik 1970, Tab 1) | ||
+ | Nobs <- c(31, 50, 76, 49, 57, 65, 61, 76, 77, 90, 100, 116, 126, 143, 181, 206, 256, 293, 353, 406, 467, 514, 569, 651, 714) | ||
+ | ano <- c(1936, 1938, 1942, 1947, 1948:1968) | ||
+ | plot(ano, Nobs) | ||
+ | ## Calculo das taxas de crescimento | ||
+ | ## Taxa de crescimento instrínseco | ||
+ | rh <- diff(log(Nobs)) / diff(ano) | ||
+ | ## Taxa de crescimento discreto | ||
+ | Rh <- exp(rh) | ||
+ | ## Media e desvio-padrao dos logs | ||
+ | (R.mu <- mean(log(Rh))) | ||
+ | (R.sd <- sd(log(Rh))) | ||
+ | ## Histograma da taxas e logormal ajustada a elas | ||
+ | hist(Rh, prob=TRUE) | ||
+ | curve(dlnorm(x, | ||
+ | |||
+ | ## Acertando parametros da simulacao | ||
+ | ## Semente de numeros aleatorios | ||
+ | set.seed(4242) | ||
+ | ## Tamanho inicial da população | ||
+ | N <- 31 | ||
+ | ## Probabilidade de sobrevivência (Akçakaia et al) | ||
+ | sobr <- 0.921 | ||
+ | |||
+ | ## Aqui começa a simulação de um ano | ||
+ | ## Sorteia um valor de taxa de crescimento anual | ||
+ | (R <- rlnorm(1, meanlog = R.mu, sdlog = R.sd)) | ||
+ | ## Fecundidade: | ||
+ | (pS <- sobr) | ||
+ | (pF <- R - pS) | ||
+ | ## Numero de sobreviventes no ano seguinte | ||
+ | (Nv <- rbinom(1, size = N, prob = pS)) | ||
+ | ## Numero de filhotes | ||
+ | ## Por sobrevivente | ||
+ | (ninhadas <- rpois(n = Nv, lambda = pF)) | ||
+ | ## Total de filhotes | ||
+ | (Nf <- sum(ninhadas)) | ||
+ | ## Total da populacao: sobreviventes + filhotes | ||
+ | (N <- Nv + Nf) | ||
+ | </ | ||
==== Para fazer as simulações ==== | ==== Para fazer as simulações ==== | ||
- | Para responder as perguntas | + | Para responder as perguntas |
== 1. Carregue as funções no R == | == 1. Carregue as funções no R == | ||
Linha 138: | Linha 185: | ||
== 2. Use a função para executar as simulações == | == 2. Use a função para executar as simulações == | ||
- | A função '' | + | O código que você acabou de excutar no R cria as funções para realizar as simulações. Você vai precisar apenas da função '' |
- | * N0: tamanho inicial da populacao | + | * '' |
- | * anos: numero de passos de tempo a simular | + | * '' |
- | * sobr: taxa de sobrevivência a cada passo de tempo | + | * '' |
- | * R.mu: média do logarítimo da taxa de crescimento | + | * '' |
- | * R.sd: desvio-padrão do logarítimo da taxa de crescimento | + | * '' |
- | * nrep: número de repetições da simulação. | + | * '' |
Por exemplo, para executar quatro simulações da dinâmica dos primeiros cinco anos, com os parâmetros da população de bois almiscarados de Nunivak use o comando: | Por exemplo, para executar quatro simulações da dinâmica dos primeiros cinco anos, com os parâmetros da população de bois almiscarados de Nunivak use o comando: | ||
Linha 176: | Linha 223: | ||
Você pode então fazer os cálculos com este objeto no R, ou exportá-lo para uma planilha com o comando | Você pode então fazer os cálculos com este objeto no R, ou exportá-lo para uma planilha com o comando | ||
- | < | + | < |
+ | | ||
==== Questões ==== | ==== Questões ==== | ||
- | - Faça um histograma da distribuição | + | - Faça um gráfico |
- | - Calcule a probabilidade de extinção da população | + | |
- | - Calcule a probabilidade | + | |
- | - Faça um histograma de tamanhos populacionais mínimos nos doze primeiros anos após a introdução. | + | - Faça uma análise de sensibilidade |
- | - Faça uma análise de sensibilidade | + | |
+ | ==== Você é rápido com a teoria de biogeografia de ilhas?==== | ||
+ | * [https:// | ||
====Referências==== | ====Referências==== | ||
* Akçakaya, H. R., Burgman, M. A. & Guinzburg, L. R. 1999. Applied Population Ecology. Sinauer Associates, Sunderland, Massachusetts. | * Akçakaya, H. R., Burgman, M. A. & Guinzburg, L. R. 1999. Applied Population Ecology. Sinauer Associates, Sunderland, Massachusetts. | ||
* Spencer, D.L. & Lensink, C.J. 1970. The muskox of Nunivak island, Alaska. The Journal of Wildlife Management, 1-15. | * Spencer, D.L. & Lensink, C.J. 1970. The muskox of Nunivak island, Alaska. The Journal of Wildlife Management, 1-15. |
apoio/exercicios.1507524751.txt.gz · Última modificação: 2024/01/12 10:39 (edição externa)