Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- apoio:exercicios [2017/10/10 14:52] – [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 ======
-Terminou a lista de exercícios? Aqui vão alguns problemas em aberto:
+Alguns problemas em aberto para trabalhar em sala:
 ===== Coelhos na Austrália =====
@@ Linha 72: / Linha 72: @@
 ==== Para entender as simulações passo a passo====
-Baixe {{:apoio:muskox_step_by_step.zip|este código}}, mude a extensão de "zip" para "R", abra no R e acompanhe a explicação dos professores.
+Copie o código abaixo em um script no R e acompanhe a explicação dos professores:
+<code>
+## 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, meanlog = R.mu, sdlog = R.sd), add=TRUE)
+## 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: diferenca entre a taxa de crescimento e a sobrevivencia
+(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)
+</code>
 ==== Para fazer as simulações ====
@@ Linha 190: / Linha 232: @@
   - 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://goo.gl/forms/ulPLygU9bU8YzjeQ2|teste aqui]]
 ====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.