Hola, Fíjate que usas "plot()", el sitema gráfico por defecto de "R". Sobre el plot que tienes, puedes usar las funciones "points()", o "lines()" para añadir elementos a ese gráfico.
En este ejemplo añado una línea horizontal a la altura de y = 100. #-------------------------------------- library(EpiModel) # SIR Model with Vital Dynamics (One-Group) param <- param.dcm(inf.prob = 0.2, act.rate = 10, 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 = 15) mod2 <- dcm(param, init, control) mod2 plot(mod2) *lines( 1:10, rep(100,10))* #------------------------------------- Gracias, Carlos Ortega www.qualityexcellence.es El dom, 7 mar 2021 a las 12:36, Jose Betancourt Bethencourt (< betans...@gmail.com>) escribió: > Estimados > Quisieramos saber como agregar al modelo que se presenta debajo una > linea de casos reales a la simulación, la cual ya tiene ya tiene una > linea de susceptibles, una de infecciosos y otra de rcuperados > > > saludos > José > > library(EpiModel) > SIR Model with Vital Dynamics (One-Group) > param <- param.dcm(inf.prob = 0.2, act.rate = 10, > 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 = 15) > mod2 <- dcm(param, init, control) > mod2 > plot(mod2) > > -- > 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 > -- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es