[R-es] simulación1

[Carlos Ortega]

Hola José,
He incluido la línea como ejemplo, con fines didácticos para que veas cómo
se hace. Sin ningún otro objetivo que pudiera estar asociado a tu estudio.

Para añadir una línea/curva o puntos sueltos, tienes primero que conseguir
esos puntos en un data.frame y luego representarles como indico a través de
las funciones "points()", "lines()"...

Gracias,
Carlos Ortega
www.qualityexcellence.es

El dom, 7 mar 2021 a las 20:54, Jose Betancourt Bethencourt (<
betanster using gmail.com>) escribió:

> El anterior no se veia, esta linea no tiene sentido
>
>
>
>
> El 7/3/21, Carlos Ortega <cof using 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 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
> >
>
>
> --
> Dr. Jose A. Betancourt Bethencourt
> Universidad de Ciencias Medicas Carlos j. Finlay
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]