Javier Milei - Pandenomics

En el modelo SIR, como en todo modelo matemático de algún fenómeno natural, se requiere incluir en él condiciones que representen, a grandes rasgos, los mecanismos involucrados. En el caso del desarrollado por Kermack y McKendrick (1927, 1933) se divide en tres clases:

(i) el número de los individuos susceptibles S ( t ), es decir aquellos individuos que pueden contagiarse;

(ii) los infecciosos I(t) , aquellos individuos que al estar enfermos pueden transmitir la enfermedad;

(iii) los removidos R(t) , aquellos individuos que luego de haberse enfermado se alivian y quedan en un estado de inmunidad o bien están aislados en algún sitio (como un hospital) o fallecieron.

De este modo, el escenario que se le presenta a un individuo descrito por el modelo SIR estándar será:

S → I → R

El modelo presenta los siguientes postulados:

a) En la epidemia, una sola infección es la responsable de ocasionar un proceso infeccioso en el huésped.

b) El desenlace de la enfermedad es la muerte o la inmunidad completa.

c) La tasa de contagio c ( t ) es proporcional al número de enfermos: c ( t ) = a.I ( t ), siendo a > 0 constante.

d) Número de individuos susceptibles que se contagian por unidad de tiempo está dado por c ( t ). S ( t ), por lo tanto, la cantidad de personas que se contagian por unidad de tiempo es:

a.I ( t ) .S ( t ) con a >0

En otras palabras el número de susceptibles S ( t ) al trascurrir un tiempo Δ t será:

S ( t +∆ t )= S ( t )- a . I ( t ). S ( t ).∆ t (1)

e) Todos los individuos sanos son susceptibles.

La población es cerrada para cualquier tiempo S + I + R = N . Es decir, no se toman en cuenta nacimientos ni migraciones.

El periodo de incubación que contempla el modelo más general de Kermack y McKendrick es muy corto. El individuo susceptible se enferma y puede contagiar a otro.

La velocidad de decremento de infectados es proporcional a su número. En términos matemáticos:

b.I ( t ) con b >0

donde b es la razón de removidos por unidad de tiempo y es constante. Esto es, el número de infecciosos I ( t ) al trascurrir un tiempo Δ t será:

I ( t +∆ t )= I ( t )+ a.I ( t ). S ( t ).∆ t - b.I ( t ).∆ t (2)

y el de removidos R ( t ) al transcurrir un tiempo Δ t será:

R ( t +∆ t )= R ( t )+ b.I ( t ).∆ t (3)

El modelo supone además que los individuos de todas las clases están mezclados y entran en contacto (excepto los que están aislados en hospitales y los muertos). Al tomar estas condiciones y al usar variables normalizadas en función del número total de la población, N :

s ≡ S ⁄ N , i ≡ I ⁄ N , r ≡ R ⁄ N

Por lo que entonces, el modelo estará descripto por el conjunto de las siguientes ecuaciones diferenciales, que resultan del sistema dinámico dado por las tres relaciones dinámicas indicadas por (1), (2) y (3). Así, la ecuación (4) indicará la velocidad a la que disminuye el número de susceptibles, debido a que se están infectando, de acuerdo con el postulado (d):

ds ⁄ d τ = - s ( τ ). i ( τ ) (4)

A su vez, la ecuación (5) nos presenta la velocidad a la que se incrementan los infecciosos de acuerdo con los postulados (d) y (h):

di ⁄ d τ = s ( τ ). i ( τ )- ρ . i ( τ ) (5)

En último lugar, la velocidad de incremento de los removidos está contemplada en la ecuación (6) de acuerdo con el postulado (h):

dr ⁄ d τ = ρ . i ( τ ) (6)

donde la constante:

ρ = b ⁄ a.N

es la razón de los removidos a los infectados, mientras que la variable temporal:

τ ≡ N.a.t

arroja el tiempo en unidades del periodo de infección.

La formulación matemática del modelo epidemiológico se completa fijando las condiciones iniciales del sistema, por ejemplo al tiempo τ = 0, que marca el inicio de la epidemia:

s ( τ =0) = s 0> 0

i ( τ =0) = i 0> 0

r ( τ =0) = r 0= 0

Éste es un modelo muy rudimentario, al cual se le han hecho algunas modificaciones a lo largo del tiempo, adecuándolo a situaciones específicas. Algunos modelos han relajado las condiciones señaladas arriba con los postulados (e) y (h) al introducir dinámicas más complicadas. Por ej. el periodo latente del HIV es muy variable. Hay enfermos que manifiestan sus primeros síntomas pocos años después de haber sido infectados, algunos nunca. En una epidemia muy larga como el SIDA los hijos de padres infectados como no infectados alteran el número total de pobladores. Esto último no está contemplado en el modelo original ya que lleva implícito conservar la población. Simplemente la suma de las tres ecuaciones diferenciales en (4), (5) y (6) es igual a cero, así:

s ( τ ) + i ( τ ) + r ( τ ) = 1

donde la unidad representa el total de la población, esto es, la suma de susceptibles, infectados y removidos es constante. Otras variaciones al modelo SIR contemplan casos de epidemias donde la población de susceptibles es heterogénea, ya que pueden integrarla también niños y ancianos. Mollison (1995) detalla más sobre las modificaciones al modelo SIR. Sin embargo, aun cuando el modelo estándar es muy elemental, permite hacer algunos señalamientos relevantes sobre las epidemias.

A su vez, una de las características de las ecuaciones diferenciales que representan al modelo SIR, si bien no tienen una solución analítica, es posible encontrar soluciones numéricas. En el gráfico que se presenta a continuación nos muestra la solución numérica sobre cómo evolucionan las tres variables del modelo SIR para un caso de epidemia, donde se observa que el número de infectados aumenta hasta un máximo y después decrece.

EVOLUCIÓN TÍPICA DE LAS VARIABLES SIR EN UNA EPIDEMIA

Sin embargo, sin necesidad de recurrir a métodos numéricos, inferimos algunas de las características cualitativas de las soluciones. Para comenzar, el modelo sólo tiene sentido si las funciones s(τ) e i(τ) se mantienen no negativas porque representan un número de personas de alguna de las dos clases. En otras palabras, si alguna de las dos funciones alcanza el valor 0, el sistema termina, es decir, con cero susceptibles o cero infecciosos.

Subrayamos además que el número de susceptibles disminuye todo el tiempo en el modelo, esto es, ds / d τ < 0 . Asimismo, observamos que para que exista una epidemia debe cumplirse que el número de infecciosos i( τ ) al inicio aumente con el tiempo. Queremos remarcar que por epidemia se entiende que i ( τ ) > i 0para algún tiempo τ . En términos matemáticos, hay epidemia si di / d ô > 0 y esto ocurre si y sólo si s > ñ (véase ecuación 5). Sin embargo, acabamos de mencionar que s disminuye todo el tiempo, por lo que en algún momento la desigualdad se puede invertir y causaría que ahora di / d τ < 0: el número de infecciosos disminuirá con el tiempo. En pocas palabras, eventualmente se acaba la epidemia. La discusión anterior significa que existe un umbral, dentro del modelo, que determina si se desarrolla o no epidemia.