José
Envié antes el correo
El código anda, yo tendrÃa que estudiarlo, pero, posiblemente no funcione
predict porque no está pensado para eso, pero podrÃa calcularlo, sin saber
como, posiblemente desde tiempo, su gráfica da 500 en tiempo. Pero no
trabajo sobre esa librerÃa, opinar de mi parte serÃa absurdo.
Javier Marcuzzi
El 24 may. 2023, a las 08:39, Jose Betancourt Bethencourt
<betanster en gmail.com> escribió:
Quisiéramos agregar predicción para siete dÃas a este modelo :
library(EpiModel)
param <- param.dcm(inf.prob = 0.2, act.rate = 5,
rec.rate = 1/3, a.rate = 1/90, ds.rate = 1/100,
di.rate = 1/35, dr.rate = 1/100)
init <- init.dcm(s.num = 500, i.num = 1, r.num = 0)
control <- control.dcm(type = "SIR", nsteps = 500)
mod2 <- dcm(param, init, control)
mod2
plot(mod2)
El 24/5/23, Javier Marcuzzi <javier.ruben.marcuzzi en gmail.com> escribió:
Estimado José Betancourt
Corrà si código en Mac, el mensaje que me da a mà es el siguiente:
Error in UseMethod("predict") :
no applicable method for 'predict' applied to an object of class "icm"
No sabrÃa más que decir.
Javier Rubén Marcuzzi
El 24 may. 2023, a las 08:23, Jose Betancourt Bethencourt
<betanster en gmail.com> escribió:
library(EpiModel)
# Definir los parámetros iniciales
param <- param.icm(inf.prob = 0.2, act.rate = 0.25, rec.rate = 1/50)
# Definir las condiciones iniciales
init <- init.icm(s.num = 500, i.num = 1, r.num = 0)
# Definir las opciones de control
control <- control.icm(type = "SIR", nsteps = 500, nsims = 10)
# Crear el modelo original
mod1 <- icm(param, init, control)
# Ajustar los parámetros para reflejar la situación actual
param2 <- param.icm(inf.prob = 0.3, act.rate = 0.3, rec.rate = 1/40)
# Crear un nuevo modelo conlos nuevos parámetros
mod2 <- icm(param2, init, control)
# Generar predicciones para los próximos 7 dÃas
pred <- predict(mod2, times = 501:508, nsims = 100)
# Graficar las predicciones
plot(pred)
--
Dr. Jose A. Betancourt Bethencourt
Universidad de Ciencias Medicas Carlos j. Finlay