[R-es] R pred
Jose Betancourt Bethencourt
bet@n@ter @end|ng |rom gm@||@com
Mie Mayo 24 13:23:33 CEST 2023
Estimados
Quisiera saber porque el modelo siguiente no me hace las predicciones
--
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)
saludos
José
Más información sobre la lista de distribución R-help-es