Skip to content

[R-es] simulación1

5 messages · Jose Betancourt B., Carlos Ortega

Original
Jose Betancourt B.
# 
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)
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 en gmail.com>) escribió:

  
    

      
      
    

  


  


      

    

    

      
      

        


  

        

      Jose Betancourt B.
    

    

      
      #
    

  


  

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

            

      
      
      
    

        

      
      
    

            

      
      
      
    

        

  
    

      
      
    

  


  


      

    

    

      
      

        


  

        

      Jose Betancourt B.
    

    

      
      #
    

  


  

      
El anterior no se veia, esta linea no tiene sentido




El 7/3/21, Carlos Ortega <cof en qualityexcellence.es> escribió:

            

      
      
      
    

        

      
      
    

            

      
      
      
    

        

  
    

      
      
    

  


  


      

    

    

      
      

        


  

        

      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 en gmail.com>) escribió: