| Ambos lados da revisão anteriorRevisão anteriorPróxima revisão | Revisão anterior |
| mod1:mat_apoio:meta_resgate [2010/08/17 23:38] – adalardo | mod1:mat_apoio:meta_resgate [2024/01/11 15:21] (atual) – edição externa 127.0.0.1 |
|---|
| "\n c=",c," l=",l," fi=",fi," pi=",pi," e=",e), | "\n c=",c," l=",l," fi=",fi," pi=",pi," e=",e), |
| font.lab=2,lwd=2) | font.lab=2,lwd=2) |
| abline(h=F,col=2,lwd=2,lty=2) | abline(h=F,col=2,lwd=2,lty=2) # equilibrio F |
| |
| | points(1:tf,c(e*(1-fi),res),type='l',lwd=2,col="blue") # pe observado |
| points(1:tf,c(e*(1-fi),res),type='l',lwd=2,col=4) | abline(h=e-pi,col="green",lwd=2,lty=2) # pe equilibrio |
| abline(h=e-pi,col=3,lwd=2,lty=2) | legend("topright", legend=c("proporção ocupada", "equilíbrio F", "prob. extinção (pe)", "equilíbrio pe"), lty=c(1,2,1,2), col=c("black","red","blue", "green"), bty="n") |
| abline(h=0,lty=2) | |
| | |
| return(paisag) | return(paisag) |
| pe=e*(1-sum(paisag[,,t-1])/(c*l)) | pe=e*(1-sum(paisag[,,t-1])/(c*l)) |
| pi=i*sum(paisag[,,t-1])/(c*l) | 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)]==1]<-sample(c(0,1),sum(paisag[,,t-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)) | 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) | resultado[t-1]=sum(paisag[,,t])/(c*l) |
| rese[t-1]=pe | rese[t-1]=pe |
| resi[t-1]=pi | resi[t-1]=pi |
| } | } |
| plot(1:tf,c(fi,resultado),type="l",xlab="Tempo",ylab="Fração de manchas ocupadas", | plot(1:tf,c(fi,resultado),type="l",xlab="Tempo",ylab="Proporção/Probabilidade", |
| ylim=c(0,1),main=paste("Colonização Interna","\n c=",c," l=",l," fi=",fi," i=",i," e=",e),font.lab=2,lwd=2) | ylim=c(0,1),main=paste("Colonização Interna","\n c=",c," l=",l," fi=",fi," i=",i," e=",e),font.lab=2,lwd=2) |
| abline(h=0,lty=2) | abline(h=0,lty=2) |
| | |
| points(1:tf,c(i*fi,resi),type='l',lwd=2,col=6,lty=3) | points(1:tf,c(i*fi,resi),type='l',lwd=2,col=6,lty=3) |
| | legend("topright", legend=c("manchas ocupadas", "prob.colonização", "prob.extinção"), lty=c(1,3,3), col=c(1,6,4), bty="n") |
| | |
| | |
| return(paisag) | return(paisag) |
