--- mod1:mat_apoio:meta_inter [2010/08/16 18:56] – paulo
+++ mod1:mat_apoio:meta_inter [2024/01/11 15:21] (atual) – edição externa 127.0.0.1
@@ Linha 22: / Linha 22: @@
 <code>
 meta.inter=function(tf,c,l,fi,i,pe){
-	paisag=array(0,dim=c(l,c,tf))
+         paisag=array(0,dim=c(l,c,tf))
-	paisag[,,1]=sample(c(rep(0,round(c*l-fi*c*l)),rep(1,round(fi*c*l))))
+         nmanchas=c*l
-	resultado=numeric()
+         paisag[,,1]=matrix(sample(c(rep(1,fi*nmanchas), rep(0,round((1-fi)*nmanchas)))), ncol=c)
-	for(t in 2:tf){
+         resultado=numeric()
-		 pi=i*sum(paisag[,,t-1])/(c*l)
+         for(t in 2:tf){
-	       paisag[,,t][paisag[,,(t-1)]==1]<-sample(c(0,1),sum(paisag[,,1]),replace=T,prob=c(pe,1-pe))
+           pi=i*sum(paisag[,,t-1])/(c*l)
-	       paisag[,,t][paisag[,,(t-1)]==0]<-sample(c(0,1),c*l-sum(paisag[,,1]),replace=T,prob=c(1-pi,pi))
+           paisag[,,t][paisag[,,(t-1)]==1]<-sample(c(0,1),sum(paisag[,,(t-1)]), replace=T,prob=c(pe,1-pe))
-	       resultado[t-1]=sum(paisag[,,t])/(c*l)
+           paisag[,,t][paisag[,,(t-1)]==0]<-sample(c(0,1),c*l-sum(paisag[,,(t-1)]), replace=T,prob=c(1-pi,pi))
-	      }
+           resultado[t-1]=sum(paisag[,,t])/(c*l)
+           }
+         F=1-(pe/i)
+         plot(1:tf,c(fi,resultado),type="l",xlab="Tempo",ylab="Fração de manchas ocupadas",
+ ylim=c(0,1),main=paste("Colonização Interna","\n c=",c," l=",l," fi=",fi," i=",i," pe=",pe),font.lab=2,lwd=2)
+ abline(h=F,col=2,lwd=2,lty=2)
+       return(paisag)
+}
-	F=1-(pe/i)
-	plot(1:tf,c(fi,resultado),type="l",xlab="Tempo",ylab="Fração de manchas ocupadas",
-	ylim=c(0,1),main=paste("Colonização Interna","\n c=",c," l=",l," fi=",fi," i=",i," pe=",pe),font.lab=2,lwd=2)
-	abline(h=F,col=2,lwd=2,lty=2)
-      return(paisag)
-	}
 </code>
-E agora você pode simular o modelo com os valores que escolher para os valores dos argumentos da função, como:
+E agora você pode simular o modelo com os valores que escolher para os argumentos da função, como:
 <code>
@@ Linha 55: / Linha 56: @@
- Para finalizar, uma última animaçãozinha:
+ Para finalizar, uma última animaçãozinha, antes salvo o resultado de uma simulação em um arquivo, por exemplo:
 <code>
-meta.inter2=function(tf,c,l,fi,i,pe){
-	paisag=array(0,dim=c(l,c,tf))
-	paisag[,,1]=sample(c(rep(0,round(c*l-fi*c*l)),rep(1,round(fi*c*l))))
-	resultado=numeric()
-	for(t in 2:tf){
-		 pi=i*sum(paisag[,,t-1])/(c*l)
-	       paisag[,,t][paisag[,,(t-1)]==1]<-sample(c(0,1),sum(paisag[,,1]),replace=T,prob=c(pe,1-pe))
-	       paisag[,,t][paisag[,,(t-1)]==0]<-sample(c(0,1),c*l-sum(paisag[,,1]),replace=T,prob=c(1-pi,pi))
-	       resultado[t-1]=sum(paisag[,,t])/(c*l)
-	      }
-	F=1-(pe/i)
+sim.int1 <- meta.inter(20,10,10,1, 0.4,0.2)
-      return(paisag)
+</code>
-	}
-anima2=function(tf,c,l,fi,i,pe){
+Agora passe a função abaixo para o programa
-	dados=meta.inter2(tf,c,l,fi,i,pe)
+<code>
+anima2=function(dados){
+	tf=dim(dados)[3]
 	for(i in 1:tf){
 	image(dados[,,i], main=("Ocupação de manchas"),col=c("white","red"),bty="n",xaxt='n',yaxt='n')
-	grid(c,l)
+	grid(dim(dados)[1],dim(dados)[2])
-	Sys.sleep(.5)
+	Sys.sleep(.2)
 	}
 	}
-anima2(tf=25,c=10,l=10,fi=.1,i=1,pe=0.1)
 </code>
-Abaixo o script das funções. A extensão está como pdf, mas na realidade é um arquivo de texto simples. Para abrir deve salvar em seu computador e remover a extensão //pdf//, deixando a extensão //r//. Depois disso pode abrir em qualquer editor de texto (Word, Notepad, ou, preferência com o próprio editor do R)
+Agora é só rodar a função acima com o resultado da simulação:
-  * {{mod1:mat_apoio:metapop.r.pdf}}
+<code>
+anima2(dados=sim.int1)
+</code>

mod1/mat_apoio/meta_inter.1281985015.txt.gz · Última modificação: 2024/01/11 15:21 (edição externa)

Mostrar página Revisões anteriores

[unknown link type]Voltar ao topo