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 (<
betanster en 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 en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]