[R-es] simulación1

[Jose Betancourt Bethencourt]

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 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

------------ 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>