[R-es] simulación1

[Carlos Ortega]

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 using 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 using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]