[R-es] R pred

[Javier Marcuzzi]

Estimado José

Ese último código me da esta gráfica

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: Captura de pantalla 2023-05-24 a la(s) 08.42.49.png
Type: image/png
Size: 251278 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20230524/4ea738b0/attachment-0001.png>

------------ próxima parte ------------


Javier Rubén Marcuzzi

> El 24 may. 2023, a las 08:39, Jose Betancourt Bethencourt <betanster using 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 using 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 using 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