Foro de debate

Ayuda con algoritmo en R

Ayuda con algoritmo en R

de Mauricio Bernal -
Número de respuestas: 3
 si no es mucha molestia me podrían ayudar con este 
Hacer un programa que simule el crecimiento exponencial de una población con 
tasa finita de crecimiento  λ y población inicial N0. Hacer un gráfico con el 
tamaño poblacional hasta un tiempo t = 100.

gracias

En respuesta a Mauricio Bernal

Re: Ayuda con algoritmo en R

de Juan Miguel Requena Mullor -
Hola Mauricio, las cosas que tiene esta vida, justo este fin de semana estuve trabajando cosas parecidas a las que necesitas, así que aquí tienes una opción con gráfico incluido, espero te sirva. (debes tener instalado el paquete "abind").

No=10
λ=1
t=100
library("abind")
par(mfrow=c(1,1))

j<-c(1:100)
n1<-NULL
a<-function(x){
  b<-(1*(x/100))*(1-(x/100))
  c<-b*100
}
for (i in j){
  if (1==i){ n1<-as.vector(a(10))       
  }else{m<-a(n1[i-(i-1)])
  n1<-abind(m,n1,along=)
  }
}
#n1
par(bg="white",bty="l")
plot(j,n1,main="Crecimiento poblacional",type="l",lwd=2.5,lty=2,
xlab="Generaciones",ylab="Tamaño Poblacional",pch=21,bg="grey")
En respuesta a Mauricio Bernal

Re: Ayuda con algoritmo en R

de Juan Miguel Requena Mullor -
Hola, el script que te pasé ayer tiene algunos fallos, sustitúyelo por este:

library("abind")
par(mfrow=c(1,1))

j<-c(1:100)
n1<-NULL
a<-function(x){
  b<-(1*(x/100))*(1-(x/100))
  c<-b*100
}
for (i in j){
  if (1==i){ n1<-as.vector(a(10))       
  }else{m<-a(n1[i-(i-1)])
        n1<-abind(m,n1,along=)
  }
}
n2<-n1[100:1]

par(bg="white",bty="l")
plot(j,n2,main="Crecimiento poblacional",type="l",lwd=2.5,lty=2,
     xlab="Generaciones",ylab="Tamaño Poblacional",pch=21,bg="grey")

#se me olvidó decirte que el modelo tiene capacidad de carga y es de 100 individuos.