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:35] – [Para fazer as simulações] 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 24: | Linha 24: | ||
Você pode também buscar outros dados para substituir estes ou complementá-los (p. ex. em bases como a [[https:// | Você pode também buscar outros dados para substituir estes ou complementá-los (p. ex. em bases como a [[https:// | ||
- | ==Fonte== | + | ====Fonte==== |
* Stodart E, Parer I (1988) Colonisation of Australia by the rabbit // | * Stodart E, Parer I (1988) Colonisation of Australia by the rabbit // | ||
Linha 52: | Linha 52: | ||
Consulte os dois artigos e faça sua avaliação: | Consulte os dois artigos e faça sua avaliação: | ||
- | ==Fontes== | + | ===Fontes=== |
* Sabath, M. D., Boughton, W. C., & Easteal, S. (1981). Expansion of the range of the introduced toad Bufo marinus in Australia from 1935 to 1974. Copeia, 676-680. | * Sabath, M. D., Boughton, W. C., & Easteal, S. (1981). Expansion of the range of the introduced toad Bufo marinus in Australia from 1935 to 1974. Copeia, 676-680. | ||
* Easteal, S., van Beurden, E. K., Floyd, R. B., & Sabath, M. D. (1985). Continuing geographical spread of Bufo marinus in Australia: range expansion between 1974 and 1980. Journal of Herpetology, | * Easteal, S., van Beurden, E. K., Floyd, R. B., & Sabath, M. D. (1985). Continuing geographical spread of Bufo marinus in Australia: range expansion between 1974 and 1980. Journal of Herpetology, | ||
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 == | ||
+ | |||
+ | Abra o R, copie o código abaixo e cole na linha de comando do R: | ||
- | - Abra o R, copie o código abaixo e cole na linha de comando do R: | ||
< | < | ||
## Simula crescimento populacional a tempo discreto com estocasticidade demografica e ambiental | ## Simula crescimento populacional a tempo discreto com estocasticidade demografica e ambiental | ||
Linha 132: | Linha 182: | ||
</ | </ | ||
- | - A 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. | ||
- | ==== Questões ==== | ||
- | - Faça um histograma | + | == 2. Use a função para executar as simulações == |
- | | + | |
- | - Calcule a probabilidade | + | 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 '' |
- | | + | |
- | - Faça uma análise | + | * '' |
+ | | ||
+ | * '' | ||
+ | * '' | ||
+ | | ||
+ | | ||
+ | |||
+ | Por exemplo, para executar quatro simulações da dinâmica dos primeiros | ||
+ | |||
+ | < | ||
+ | repete(N0 = 31, anos = 5, sobr = 0.921, | ||
+ | R.mu = 0.1222906, R.sd = 0.08011615, nrep=4) | ||
+ | </ | ||
+ | |||
+ | Que irá retornar | ||
+ | |||
+ | < | ||
+ | [,1] [,2] [,3] [,4] | ||
+ | [1,] | ||
+ | [2,] | ||
+ | [3,] | ||
+ | [4,] | ||
+ | [5,] | ||
+ | </ | ||
+ | |||
+ | == 3. Execute muitas simulações e guarde o resultado == | ||
+ | |||
+ | Para resolver os problemas propostos você terá que estimar probabilidades | ||
+ | |||
+ | < | ||
+ | sim1 <- repete(N0 = 31, anos = 5, sobr = 0.921, | ||
+ | R.mu = 0.1222906, R.sd = 0.08011615, nrep=4) | ||
+ | </ | ||
+ | |||
+ | Você pode então fazer os cálculos com este objeto no R, ou exportá-lo | ||
+ | |||
+ | < | ||
+ | |||
+ | ==== Questões ==== | ||
+ | - Faça um gráfico de probabilidade acumulada dos tamanhos populacionais mínimos previstos pelo modelo após cinco anos de introdução da população na ilha. Use este gráfico para estimar: | ||
+ | - A probabilidade de extinção da população. | ||
+ | - A probabilidade do tamanho populacional ter caído abaixo do número inicial no período. | ||
+ | - Faça uma análise de sensibilidade das duas probabilidades acima à variação dos parâmetros demográficos usados na simulação. Para isso, calcule a probabilidade de extinção em 12 anos variando os parâmetros em cinco a dez por cento por cento do valor original, para cima e para baixo. Faça variar um parâmetro por vez, mantendo os outros constantes. Use estes resultados para propor recomendações para um programa de reintrodução. | ||
+ | ==== Você é rápido com a teoria de biogeografia de ilhas?==== | ||
+ | * [https:// | ||
+ | ====Referências==== | ||
+ | * 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. |
apoio/exercicios.1507523732.txt.gz · Última modificação: 2024/01/12 10:39 (edição externa)