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:58] 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 181: Linha 185:
  
 <box 80% red |Exercício> <box 80% red |Exercício>
-Agora é com vocês! Responda às seguintes perguntas baseado em simulações+Agora é com vocês! Tentem responder às seguintes perguntas: 
-  * O que acontece quando roda o modelo com os mesmos parâmetros mais de uma vez? 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>  </box> 
 == DICA == == DICA ==
Linha 230: 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