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
> <betans...@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.marcu...@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
>>> <betans...@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
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
