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]]
Saludos la linea que genera no tiene sentido para mi, deber ia ser una linea similar a la linea ascendente de los infectados con crecimiento exponencial 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 "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 >-- Dr. Jose A. Betancourt Bethencourt Universidad de Ciencias Medicas Carlos j. Finlay ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: plot_zoom_png Type: application/octet-stream Size: 8151 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20210307/5129cd66/attachment-0001.obj>
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 "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 >-- Dr. Jose A. Betancourt Bethencourt Universidad de Ciencias Medicas Carlos j. Finlay ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: plot_zoom_png.png Type: image/png Size: 8151 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20210307/fd61454c/attachment.png>