Modelos de crecimiento poblacional: Malthus, Verhulst, Gompertz, Von Bertalanffy

[jose romero]

Saludos a todos:

Estoy re-estrenando este grupo Google que como pueden ver, data del último trimestre del 2008 cuando dí este curso en la UNEFA.  Si revisan los temas de debate verán como eran nuestras "parrilladas" en aquel entonces.  Una vez inscritos, ustedes pueden abrir nuevos temas bien sea desde la interfáz web del grupo o enviando el tema por correo electrónico (a la dirección simulacion...@googlegroups.com).  Traten de no abrir muchos temas nuevos- busquen primero entre los temas si hay alguno ya creado en el cual puedan plantear sus preguntas y comentarios.  Por ejemplo, en este tema podemos debatir todo lo relacionado a los modelos de crecimiento poblacional (dentro del tema de la dinámica de sistemas).

Para comenzar el debate, le anexo dos archivos: uno con la data del censo como un archivo csv y el otro es el script en R para simular el crecimiento poblacional según el modelo de Malthus usando simecol.  Observen cómo se utiliza la función fitOdeModel para ajustar los parametros del modelo, retroalimentarlo, volver a simular y comparar los resultados entre un conjunto y otro de valores de parámetros.

Quiero que investiguen sobre los otros modelos de crecimiento poblacional (Varhulst, Gompertz, Von Bertalanffy y otros) para ajustarlos a la data del censo venezolano utilizando la función fitOdeModel.  Practiquen con esto y planteen sus dudas aqui.

[Romero, José Loreto]

Unas preguntas para debatir:

1) ¿Cuál de los modelos de crecimiento poblacional se ajustan mejor al caso venezolano? ¿Qué interpretación tiene esto?
2) ¿Cuál métrica de bondad de ajuste utilizamos para evaluar los distintos modelos?
3) Hay años para los cuales algunos modelos sobre-estiman (o sub-estiman?) la población y hay otras décadas para los cuales la población venezolana según el censo coincide más o menos con la proyección según el modelo.  ¿Porqué el ajuste es cási siempre mejor en cierto rango de fechas y en otros no?
4) ¿Podríamos plantear un modelo hibrido que utilice un modelo determinado para ajustar la población en cierto rango de fechas y otro modelo para el otro rango de fechas?

[Romero, José Loreto]

Hola:

Esto parece más bien un monologo... Quisiera saber si tienen problemas para intervenir en este foro, porque su silencio es preocupante.  Recuerden que habia mandado como tarea hacer el código R para simular el crecimiento demográfico venezolano según los otros modelos de crecimiento poblacional.  En cada caso, tienen que ajustar los parametros del modelo usando fitOdeModel.  Recuerden que esta asignación es para ser entregada este miercoles.

Aprovechen para plantear sus dudas por aqui.

[seccio...@gmail.com]

Buenas noches profesor. Soy Sandoval y estoy al tanto de la materia. Como hacemos aquellos que aun no tiene R en sus computadores y no han podido realizar nada de codigos. :S

Enviado desde mi BlackBerry de Movistar

---

From: Romero, José Loreto jlaur...@yahoo.com

Sender: simulacion...@googlegroups.com

Date: Mon, 21 May 2012 16:58:07 -0700 (PDT)

To: <simulacion...@googlegroups.com>

ReplyTo: simulacion...@googlegroups.com

Cc: jose romero<jlaur...@yahoo.com>

Subject: [Simulación y Modelos:73] Re: Modelos de crecimient o poblacional: Malthus, Verhulst, Gompertz, Von Bertalanffy

Hola:

Esto parece más bien un monologo... Quisiera saber si tienen problemas para intervenir en este foro, porque su silencio es preocupante.  Recuerden que habia mandado como tarea hacer el código R para simular el crecimiento demográfico venezolano según los otros modelos de crecimiento poblacional.  En cada caso, tienen que ajustar los parametros del modelo usando fitOdeModel.  Recuerden que esta asignación es para ser entregada este miercoles.

Aprovechen para plantear sus dudas por aqui.

--
Has recibido este mensaje porque estás suscrito al grupo "Simulación y Modelos- Unefa San Tomé" de Grupos de Google.
Para ver este debate en la Web, visita https://groups.google.com/d/msg/simulacion-y-modelos/-/91znvmZmQL8J.
Para publicar una entrada en este grupo, envía un correo electrónico a simulacion...@googlegroups.com.
Para anular tu suscripción a este grupo, envía un correo electrónico a simulacion-y-mod...@googlegroups.com
Para tener acceso a más opciones, visita el grupo en http://groups.google.com/group/simulacion-y-modelos?hl=es.

[Romero, José Loreto]

Hola Jesus:

Al principio del curso aclaré que ibamos a trabajar con lenguaje R y que debian de instalarlo o copiar el live-cd que les proporcione. 

[karmen...@gmail.com]

Buenas
profesor me puede eliminar del grupo

Este mensaje ha sido enviado gracias al servicio BlackBerry de Movilnet

---

From: Romero, José Loreto jlaur...@yahoo.com

Sender: simulacion...@googlegroups.com

Date: Wed, 23 May 2012 05:00:37 -0700 (PDT)

To: <simulacion...@googlegroups.com>

ReplyTo: simulacion...@googlegroups.com

Cc: <seccio...@gmail.com>

Subject: Re: [Simulación y Modelos:75] Re: Modelos de crecim ient o poblacional: Malthus, Verhulst, Gompertz, Von Bertal anffy

Hola Jesus:

Al principio del curso aclaré que ibamos a trabajar con lenguaje R y que debian de instalarlo o copiar el live-cd que les proporcione. 

--
Has recibido este mensaje porque estás suscrito al grupo "Simulación y Modelos- Unefa San Tomé" de Grupos de Google.

Para ver este debate en la Web, visita https://groups.google.com/d/msg/simulacion-y-modelos/-/hP0YG0piEWIJ.

[freddy]

tengan un bue día todos los presentes. profesor existimos un grupo de
personas que tenemos una duda y espero no las aclare y saque de
nuestra ignorancia

la sigueiente formula es la que tenemos que trabajar en el modelo de
Gompertz, o estoy cometiendo un gran error.?
X(t)=alfa v*(1-(X(t)/K)^1/v)*X(t)

[Rafael Rosas]

Con respecto a la Ecuacion de Van bertanlanffy
dpoblacion= poblacion_maxima (1-euler^-alfa-(tiempo - tiempo_inicial))

Como Introduzco el vector tiempo en esta parte ==> (tiempo_actual -
tiempo_inicial))??????????????

Aca esta el codigo en cuestion... ¿Como declaro los tiempo?

library(simecol)
censo <- read.table("censos.csv",header=TRUE)
censo
bertalanffy <-new("odeModel",
main=function(time,init,parms)
{
x<- init
alfa<- parms["alfa"]
poblacion_maxima<- parms["pob"]
dpoblacion= poblacion_maxima (1-euler^-alfa-
(tiempo_actual - tiempo_inicial))
list(dp)
},
parms=c(alfa=0.022696924,pob=100000000),
init=c(pob=2221572),
times=seq(from=1891,to=2011,by=0.5),
solver="rk4"
)
salida <- out(sim(gompertz))

[Rafael Rosas]

para el Modelo de SIR
Necesito el DATA FRAME de infectados por dias
Quen lo Tiene?

[freddy]

rafael hay un hilo de conversación en el que esta eso que pides.
revisa

[ing.pet...@gmail.com]

Disculpen pero estan equivocados de correo! Dejen de comunicarse por aqui.. Gracias!

Este mensaje ha sido enviado gracias al servicio BlackBerry de Movilnet

-----Original Message-----
From: freddy <alejan...@hotmail.com>
Sender: simulacion...@googlegroups.com
Date: Sat, 26 May 2012 16:05:31
To: Simulación y Modelos- Unefa San Tomé<simulacion...@googlegroups.com>
Reply-To: simulacion...@googlegroups.com
Subject: [Simulación y Modelos:80] Re: Modelos de crecim ien
t o poblacional: Malthus, Verhulst, Gompertz, Von Bertal anf
fy

--
Has recibido este mensaje porque estás suscrito al grupo "Simulación y Modelos- Unefa San Tomé" de Grupos de Google.

[freddy]

según:
modelo de Gompertz
trabajando vajo la creencia de que la ecuacion diferencial es
X(t)=alfa v*(1-(X(t)/K)^1/v)*X(t)




#modelo de "Gompertz"
library(deSolve)
library(lattice)

library(simecol)
censo <- read.table("censos.csv",header=TRUE)
censo

Gompertz <-new("odeModel",
                        main=function(time,init,parms)
                        {
              X<- init
              alfa<- parms["alfa"]
              v <- parms
              dv= alfa*v["pob"]*(1-(X/v["k"])^1/v["pob"])*X
                                list(dv)
                        },
                        parms=c(alfa=0.022696924,k=1000000),
                        init=c(pob=2221572),
                        times=seq(from=1891,to=2011,by=0.3),
                        solver="lsoda"
                 )




resultado <- fitOdeModel(Gompertz,obstime=Gompertz$a,
yobs=Gompertz[2],fn=ssqOdeModel,
method = "PORT", lower = 0 )
resultado
parms(Gompertz) <- resultado$par
salida <- out(sim(censo1))












by:




freddy
migdalis zacarias
nordis guzman
yolimar cabrera
leidis salazar








continuara.? xD

[Romero, José Loreto]

Veo mucha confusión aqui... Algunos trozos de código colgados ni siquiera correrían porque faltan los respectivos operadores aritméticos en las fórmulas, identificadores como "euler" no estan definidos, etc. También pareciera que a veces trabajan con la solución analítica del modelo porque todo eso del tiempo actual menos el tiempo inicial en una ecuación diferencial no me suena.

En fin, un desastre.

Para comenzar a desenredar esto, podrían suministrarme los documentos o la dirección donde investigaron las fórmulas de las ecuaciones diferenciales que están usando? Yo sé que en clase, les suministré las ecuaciones diferenciales del modelo de Gompertz...

[leidisss]

modelo de Von Bertalanffy

si existe algún error por favor comuníquelo. que no estamos seguras

library(simecol)
censo <- read.table("censos.csv",header=TRUE)
censo

von <-new("odeModel",
main=function(time,init,parms)
{
x<- init
r<- parms["r"]
k<- parms["k"]
dp= r*(k-x)
list(dp)
},
parms=c(r=0.022696924,k=100000000),

init=c(pob=2221572),
times=seq(from=1891,to=2011,by=0.5),
solver="rk4"
)

salida <- out(sim(von))
plot(salida)
resultado <- fitOdeModel( von,
obstime=censo$a,
yobs=censo[2],

fn=ssqOdeModel,
method = "PORT",

lower =c(r=0,k=2221572),
scale=c(1/0.022696924,1/100000000)
)
resultado

plot(resultado)



by
leidis salazar
nordis guzman

sea dulce con sus criticas :yaoming:

[Migdalys Zacarias]

hola profe.. creo que pudimos correjir el ejercicio de Gompertz y
creemos que esta es la respuesta, pero no estamos seguro. por fa
esperamos respuestas

library(simecol)
censo <- read.table("censos.csv",header=TRUE)
censo

gompertz <-new("odeModel",
main=function(time,init,parms)
{
x<- init
alfa<- parms["alfa"]
k<- parms["k"]
dp<-alfa*log(k/x["pob"])*x["pob"]
list(dp)
},
parms=c(alfa=0.022696924,k=100000000),

init=c(pob=2221572),
times=seq(from=1891,to=2011,by=0.5),
solver="rk4"
)

salida <- out(sim(gompertz))
plot(salida)
resultado <- fitOdeModel( gompertz,

obstime=censo$a,
yobs=censo[2],
fn=ssqOdeModel,
method = "PORT",

lower =c(alfa=0,k=2221572),
scale=c(1/0.022696924,1/100000000)
)
resultado

plot(resultado)


by:
yolimar cabrera
freddy jimenez


porfa. este ultimo esta bien.?? tiene errores. es que en la primera
vercion teníamos dudas con respecto a la ecuacion diferencial. pero ya
creo que salimos de esa duda.. esta mas o menos bien o mal.??? emos
echo nuestro mejor esfuerzo :(