Diferenças

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

--- roteiro:meta_chuva [2012/05/28 13:06] – adalardo
+++ roteiro:meta_chuva [2024/01/09 18:18] (atual) – edição externa 127.0.0.1
@@ 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)