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