Diferenças

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

Link para esta página de comparações

Ambos lados da revisão anteriorRevisão anterior
Próxima revisão
Revisão anterior
roteiro:meta_chuva [2012/05/28 12:26] adalardoroteiro:meta_chuva [2024/01/09 18:18] (atual) – edição externa 127.0.0.1
Linha 1: Linha 1:
 ====== Chuva de Propágulos ====== ====== Chuva de Propágulos ======
-{{:roteiro:chuvaprop.jpg?150?150 | Ênio e Beto andando na chuva http://desenhos.kids.sapo.pt/a-chuva.htm}}+{{: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:
-<m14> +
-df/dt=I - E +
-</m>, onde I é a taxa de entrada de migrantes, e E a taxa de saída. A partir dele podemos definir um modelo simples para a dinâmica de ocupação de manchas que formam uma metapopulação:+
  
 +$$\frac{df}{dt}=I - E $$
  
-<m14> + onde $I$ é a taxa de entrada de migrantes, e $E$ a taxa de saída. A partir dele podemos definir um modelo simples para a dinâmica de ocupação de manchas que formam uma metapopulação:
-df/dt=p_i(1 - f)-p_e f +
-</m> +
  
-onde <m14>p_i</m> é a probabilidade de imigração ou colonização, <m14>p_e</m> é a probabilidade de extinção e <m14>f</m> é a fração de manchas ocupadas (número de manchas ocupadas / número total de manchas).+ 
 + 
 +$$\frac{df}{dt}=p_i(1 - f)-p_e f $$ 
 +  
 + 
 +onde $p_ié a taxa de imigração ou colonização, p_eé a taxa de extinção e é 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. 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 46: Linha 47:
 paisag[,,1]<-sample(c(ocor,nocor)) paisag[,,1]<-sample(c(ocor,nocor))
 paisag[,,1] paisag[,,1]
 +image(paisag[,,1], col=c("white", "green"))
 </code> </code>
 Agora sim, estamos prontos para a simulação.  Agora sim, estamos prontos para a simulação. 
Linha 65: Linha 67:
  
 Se você observar atentamente, algumas manchas que estavam ocupadas (1) ficaram vazias (0): Se você observar atentamente, algumas manchas que estavam ocupadas (1) ficaram vazias (0):
 +
 <code> <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> </code>
 Estes foram eventos de extinção local. Note que o número de manchas ocupadas diminuiu. Estes foram eventos de extinção local. Note que o número de manchas ocupadas diminuiu.
Linha 123: Linha 127:
 Agora, para ficar mais bonito ainda, vamos criar um gráfico com esses resultados: Agora, para ficar mais bonito ainda, vamos criar um gráfico com esses resultados:
 <code> <code>
 +par(mfrow=c(1,1))
 plot(1:10,resultado$f,type="l",xlab="Tempo",ylab="Fração de manchas ocupadas", 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) ylim=c(0,1),main="Dinâmica de ocupação de manchas",font.lab=2,lwd=2)
Linha 129: 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: 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... Portanto...
Linha 179: Linha 184:
 </box> </box>
  
 +<box 80% red |Exercício>
 Agora é com vocês! Tentem responder às seguintes perguntas: 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)?   * 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?   * 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 == == DICA ==
  
Linha 201: Linha 209:
 </code> </code>
    
----- 
  
  
Linha 229: Linha 236:
 anima(simula2) anima(simula2)
 </code> </code>
- 
-E aí, gostaram? :-D 
- 
roteiro/meta_chuva.1338207978.txt.gz · Última modificação: 2024/01/09 18:17 (edição externa)
www.chimeric.de Creative Commons License Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0