Raw Message

Message-ID: <CANuuowMAb6o5FFKRVm7_u6-s2YwoNeHsJXGVsZHPw83v6L2d_A@mail.gmail.com>
Date: 2023-05-24T11:39:31Z
From: Jose Betancourt B.
Subject: [R-es] R pred
In-Reply-To: <60F46D78-7B3B-411B-A5A5-72A0627ACB6F@gmail.com>

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