Diferenças

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

--- roteiro:meta_chuva [2012/05/28 12:30] – adalardo
+++ roteiro:meta_chuva [2024/01/09 18:18] (atual) – edição externa 127.0.0.1
@@ Linha 2: / Linha 2: @@
 {{:roteiro:chuvaprop.jpg?150 | Ênio e Beto andando na chuva http://desenhos.kids.sapo.pt/a-chuva.htm}}
-Neste exercício((roteiro produzido pelo monitor Marcel Vaz a partir de material do curso de [[http://ecologia.ib.usp.br/bie5786|Ecologia de Populações]] do nosso programa de pós-graduação em ecologia)) , iremos simular a dinâmica de uma metapopulação onde as probabilidades de extinção e de colonização de manchas são constantes (modelo de chuva de propágulos). Neste modelo há sempre uma fonte constante de imigrantes que podem colonizar qualquer mancha vazia. Relembrando, a variação da fração de manchas ocupadas no tempo é descrito pela seguinte fórmula geral:
+Neste exercício, iremos simular a dinâmica de uma metapopulação onde as probabilidades de extinção e de colonização de manchas são constantes (modelo de chuva de propágulos). Neste modelo há sempre uma fonte constante de imigrantes que podem colonizar qualquer mancha vazia. Relembrando, a variação da fração de manchas ocupadas no tempo é descrito pela seguinte fórmula geral:
 $$\frac{df}{dt}=I - E $$
@@ Linha 13: / Linha 13: @@
-onde $p_i$ é a probabilidade de imigração ou colonização, $ p_e$ é a probabilidade de extinção e $ f $ é a fração de manchas ocupadas (número de manchas ocupadas / número total de manchas).
+onde $p_i$ é a taxa de imigração ou colonização, $ p_e$ é a taxa de extinção e $ f $ é a fração de manchas ocupadas (número de manchas ocupadas / número total de manchas).
 Em primeiro lugar, vamos estabelecer a probabilidade de colonização de manchas vazias (pi), a probabilidade de extinção em manchas ocupadas (pe) e a fração inicial de manchas ocupadas (fi) como 30%, 15% e 40%, respectivamente.
@@ Linha 47: / Linha 47: @@
 paisag[,,1]<-sample(c(ocor,nocor))
 paisag[,,1]
+image(paisag[,,1], col=c("white", "green"))
 </code>
 Agora sim, estamos prontos para a simulação.
@@ Linha 66: / Linha 67: @@
 Se você observar atentamente, algumas manchas que estavam ocupadas (1) ficaram vazias (0):
 <code>
-paisag[,,1] # paisagem no tempo inicial
+par(mfrow=c(1,2))
-paisag[,,2] # paisagem depois de um passo no tempo
+image(paisag[,,1], col=c("white", "green"))
+image(paisag[,,2], col=c("white", "green"))
 </code>
 Estes foram eventos de extinção local. Note que o número de manchas ocupadas diminuiu.
@@ Linha 124: / Linha 127: @@
 Agora, para ficar mais bonito ainda, vamos criar um gráfico com esses resultados:
 <code>
+par(mfrow=c(1,1))
 plot(1:10,resultado$f,type="l",xlab="Tempo",ylab="Fração de manchas ocupadas",
 ylim=c(0,1),main="Dinâmica de ocupação de manchas",font.lab=2,lwd=2)
@@ Linha 130: / Linha 134: @@
 Como nossa filosofia de vida é (ou pelo menos deveria ser) melhorar sempre, vamos acrescentar a esse gráfico uma informação muito importante: a fração de manchas ocupadas no equilíbrio (F). O F é calculado da seguinte forma:
-<m14>F=p_i/{p_i+p_e}</m>
+$$F=\frac{p_i}{p_i+p_e}$$
 Portanto...
@@ Linha 180: / Linha 184: @@
 </box>
+<box 80% red |Exercício>
 Agora é com vocês! Tentem responder às seguintes perguntas:
-  * Quando você rodar várias mais de uma vez a função com os mesmos parâmetros o resultado é o mesmo? Por quê?((Você deve entender isto para responder várias das questões seguintes.))
+  * O que acontece você rodar mais de uma vez a função com os mesmos parâmetros ? Por quê?((Você deve entender isto para responder várias das questões seguintes.))
-  * O que acontece quando o número de manchas é muito grande? E quando é muito pequeno?
+  * O que acontece com a variação na fração de manchas ocupadas quando:
+      *  o número de manchas é muito grande? E quando é muito pequeno?
+      * e se as probabilidades de extinção e colonização são muito altas?
+      * e se forem muito baixas?
   * E quando a fração inicial (fi) é muito diferente da fração em equilíbrio (F)?
   * Quais condições levam à extinção da população na paisagem neste modelo?
-  * O que acontece quando as probabilidades de extinção e de colonização são muito altas?
-  * E quando elas são muito baixas?
+</box>
 == DICA ==
@@ Linha 202: / Linha 209: @@
 </code>
-----
@@ Linha 230: / Linha 236: @@
 anima(simula2)
 </code>
-E aí, gostaram? :-D