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:29] – 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 |
- | - Copie o arquivo para um diretório em seu computador de trabalho | + | == 1. Carregue as funções no R == |
- | - Abra o R neste mesmo diretório e digite " | + | |
- | ==== Questões ==== | + | Abra o R, copie o código abaixo e cole na linha de comando do R: |
- | | + | < |
- | - Calcule | + | ## Simula crescimento populacional a tempo discreto com estocasticidade demografica e ambiental |
- | - Calcule | + | ## Parametros: |
- | - Faça um histograma | + | ## N0: tamanho inicial da populacao |
- | - Faça uma análise | + | ## 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 | ||
+ | muskox <- function(N0, anos, sobr, R.mu, R.sd){ | ||
+ | ## Um vetor sem valores para guardar os resultados | ||
+ | N <- rep(NA, anos) | ||
+ | ## primeiro elemento do vetor é o tamanho populacional inicial | ||
+ | N[1] <- N0 | ||
+ | ## O loop | ||
+ | for(i in 2:anos){ | ||
+ | ## Sorteio | ||
+ | R <- rlnorm(1, meanlog = R.mu, sdlog = R.sd) | ||
+ | ## Fecundidade no ano (diferenca da taxa de crescimento e sobrevivência) | ||
+ | ## Para os casos em que R < pS: | ||
+ | ## decidimos que não houve reprodução e portanto a taxa de sobrevivencia é a propria taxa de crescimento | ||
+ | if(R < sobr){ | ||
+ | pS <- R | ||
+ | pF <- 0} | ||
+ | if(R >= sobr){ | ||
+ | pS <- sobr | ||
+ | pF <- R - sobr} | ||
+ | ## Se a populacao nao se extinguiu o crescimento prossegue | ||
+ | if(N[i-1]> | ||
+ | ## Numero | ||
+ | Nv <- rbinom(1, size = N[i-1], 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[i] <- Nv + Nf | ||
+ | } | ||
+ | ## Se a populacao se extinguiu prossegue com zero | ||
+ | | ||
+ | N[i] <- 0 | ||
+ | } | ||
+ | return(N) | ||
+ | } | ||
+ | |||
+ | ## Repete | ||
+ | ## Argumentos adicionais: | ||
+ | ## nrep: numero de repeticoes | ||
+ | repete <- function(N0, | ||
+ | ## Uma matrix para guardar os resultados | ||
+ | results <- matrix( nrow = anos, ncol = nrep) | ||
+ | ## repete as simulacoes | ||
+ | for(i in 1:nrep) | ||
+ | results[,i] <- muskox(N0, anos, sobr, R.mu, R.sd) | ||
+ | return(results) | ||
+ | } | ||
+ | |||
+ | </ | ||
+ | |||
+ | == 2. Use a função para executar as simulações == | ||
+ | |||
+ | O código que você acabou de excutar | ||
+ | |||
+ | * '' | ||
+ | * '' | ||
+ | * '' | ||
+ | * '' | ||
+ | | ||
+ | * '' | ||
+ | |||
+ | 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.1507523348.txt.gz · Última modificação: 2024/01/12 10:39 (edição externa)