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
Última revisãoAmbos lados da revisão seguinte
roteiro:meta_chuva [2012/05/28 12:53] adalardoroteiro:meta_chuva [2012/05/28 17:16] mortara
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. 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]<-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 66: 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 124: 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 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: 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 180: 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 229: Linha 236:
 anima(simula2) anima(simula2)
 </code> </code>
- 
roteiro/meta_chuva.txt · Última modificação: 2024/01/09 18:18 por 127.0.0.1
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