El anterior no se veia, esta linea no tiene sentido
El 7/3/21, Carlos Ortega <cof en qualityexcellence.es> escribió:
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
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