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 (< betans...@gmail.com>) escribió: > El anterior no se veia, esta linea no tiene sentido > > > > > El 7/3/21, Carlos Ortega <c...@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 (< > > betans...@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@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]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es