Fuentes:

https://web.stanford.edu/~jhj1/teachingdocs/Jones-on-R0.pdf http://indico.ictp.it/event/7960/session/3/contribution/19/material/slides/0.pdf https://en.wikipedia.org/wiki/Compartmental_models_in_epidemiology

Hice un modelito SEIR , compartimente por simplicidad en 3 grupos etarios:

• [0] = 0-19
• [1] = 20-49
• [2] = 50+

(el modelo tiene mas modificaciones ver comments en el codigo)

En el sistema de ecuaciones pienso que:

• la poblacion 0 y 1 estan en contacto y no con los viejos
• la poblacion 0 no muere
• la poblacion 1 esta en contacto con todos
• la poblacion 2 esta en contacto solo con la 1

detalles:

• de los contagiados asintomaticos que pasan a recuperarse solos
• los contagiados, contagian 1/4 del rate que los sanos
• los contagiados viejos no se sanan, (si o si se infectan)

parametros:

t = entra la secuencia de tiempo en dias, en este caso no es usado, pero uno puede cambiar cosas en funcion del tiempo en que se llama a la funcion. Idea original: cambiar la porpagacion o las muertes segun las camas o aislamiento social en el tiempo

x = vector de poblaciones,ver notacion mia:

  S: Sanos. es decir poblacion original. a t-> inf esto tiene que ser cero
  E: Contagiados. poblacion que tiene el virus asintomaticamente
  I: Infectados. Gente que se le detecto el virus
  R: recuperados. Gente que sobrevino a la epidemia
  M: Muertos. gente que NO sobrevino a la epidemia

params = parametros generales del modelo:

  bet: 	beta en wikipedia, multiplicado por los contagiados es el rate de contagio, unidades 1/dias. Paso de S -> E
  a: 	1/tiempo de incuvacion, esto es propio de la sepa de virus (y del cuerpo de la persona).
      yo lo fije 1.0/5.2 dias (Whutan), como propio del virus. Paso de E -> I
  gam:    1/tiempo de reuperacion (con esto tmb calculamos el R0), esto si varia en el tiempo y grupo poblacional.
      En mi "modelo", permito que tanto los contagiados como los infectados se recuperen. Pasa de E -> R y I-> R
  mu:	1/Tiempo de muertos. Paso de I-> M
  N: 	Numero total de la poblacion. 3.5E6 es como la pobacion de la provincia de cordoba

** Ejemplo de modo de uso **

R> ddo2()