Diferenças

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

--- roteiro:meta_inter [2012/05/28 12:39] – adalardo
+++ roteiro:meta_inter [2012/05/28 17:41] – adalardo
@@ Linha 4: / Linha 4: @@
-Podemos((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))  eliminar do modelo anterior o pressuposto de uma chuva de propágulos constante e fazer com que a colonização seja uma função do número de lugares ocupados. Em uma formulação simples desse modelo, a fonte de propágulos é unicamente interna (sistema fechado) e a probabilidade de colonização varia de forma linear à proporção de lugares ocupados.
+Podemos  eliminar do modelo anterior o pressuposto de uma chuva de propágulos constante e fazer com que a colonização seja uma função do número de lugares ocupados. Em uma formulação simples desse modelo, a fonte de propágulos é unicamente interna (sistema fechado) e a probabilidade de colonização varia de forma linear à proporção de lugares ocupados.
 Dessa forma, nosso modelo não terá mais uma probabilidade de colonização constante (**pi**), mas sim uma probabilidade de colonização dependente do número de manchas ocupadas:
@@ Linha 68: / Linha 68: @@
 <code>
-anima2=function(dados){
+meta.anima2=function(dados)
-	tf=dim(dados)[3]
+{
-	for(i in 1:tf){
+nsim=dim(dados)[3]
-	image(dados[,,i], main=("Ocupação de manchas"),col=c("white","red"),bty="n",xaxt='n',yaxt='n')
+ln=dim(dados)[1]
-	grid(dim(dados)[1],dim(dados)[2])
+cl=dim(dados)[2]
-	Sys.sleep(.2)
+image(0:ln, 0:cl, dados[,,1], col=c("white", "green") , breaks=c(0,0.99,5),main="Metapopulation Dynamics", sub=paste("Initial configuration from", nsim," simulations",  sep=""), xlab="", ylab="")
-	}
+grid(ln,cl)
+Sys.sleep(.5)
+conta12=dados[,,1]+ (2*dados[,,2])
+image(0:ln, 0:cl, conta12, col=c("white","red","lightgreen", "darkgreen") , breaks=c(0,0.9,1.9,2.9,3.9),main="Metapopulation Dynamics", sub=paste("red= extintion; light green= colonization; dark green = permanence \n maximum time = ", nsim, sep=""), xlab="", ylab="")
+	for(i in 3:nsim)
+	{
+	conta12=dados[,,(i-1)]+ (2*dados[,,i])
+	image(0:ln, 0:cl, conta12, col=c("white","red","lightgreen", "darkgreen") , breaks=c(0,0.9,1.9,2.9,3.9), xlab="", ylab="", add=TRUE)
+	Sys.sleep(.1)
 	}
+}
+################
 </code>
@@ Linha 83: / Linha 94: @@
 <code>
-anima2(dados=sim.int1)
-</code>
+meta.anima2(dados=sim.int1)
+</code>