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
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 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
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