Modelos compartimentales en epidemiología

Los modelos compartimentales constituyen una técnica utilizada para simplificar la modelización matemática de las enfermedades infecciosas. La población se divide en compartimentos, asumiendo que cada individuo en un mismo compartimento tiene las mismas características.

Su origen se encuentra a principios del siglo XX, siendo un trabajo importante la primera investigación de Kermack y McKendrick en 1927. ​ Los modelos suelen investigarse mediante ecuaciones diferenciales ordinarias (las cuales son determinísticas), pero también se pueden observar en un marco estocástico, el cual es más realista pero también más complicado de analizar.

Los modelos compartimentales se pueden utilizar para predecir las características de propagación de una enfermedad, por ejemplo, la prevalencia (número total de infectados) o la duración de una epidemia.

Además, el modelo permite comprender cómo diferentes situaciones pueden afectar el resultado de la epidemia, por ejemplo, cuál es la técnica más eficiente para emitir un número limitado de vacunas en una determinada población.

El modelo SIR

El modelo SIR es uno de los modelos compartimentales más simples y muchos modelos derivan de este modelo básico. El modelo consta de tres compartimentos: S para el número de susceptibles, I para el número de infectados y R para el número de recuperados (o inmunes). Este modelo es razonablemente predictivo para enfermedades infecciosas que se transmiten de humano a humano y donde la recuperación confiere resistencia duradera, como el sarampión, las paperas y la rubéola.

Estas variables (S, I y R) representan el número de personas en cada compartimiento en un momento concreto. Para representar que el número de individuos susceptibles, infectados y recuperados puede variar con el tiempo (incluso si el tamaño de la población total permanece constante), los números precisos se convierten en una función de t (tiempo): S(t), I(t) y R(t). Para una enfermedad específica en una población específica, estas funciones pueden calcularse para predecir posibles brotes y controlarlos.

El modelo SIR es dinámico en tres sentidos

En cuanto a la función variable de la t, el modelo es dinámico en tanto que los números en cada compartimento pueden fluctuar con el tiempo. La importancia de este aspecto dinámico es más evidente en una enfermedad endémica con un corto período infeccioso, como el sarampión en el Reino Unido antes de la introducción de la vacuna en 1968. Estas enfermedades tienden a ocurrir en ciclos de brotes debido a la variación en el número de casos susceptibles (S(t)) a lo largo del tiempo. Durante una epidemia, el número de individuos susceptibles decrece rápidamente a medida que se infectan y, en consecuencia, entran en los compartimentos de infectados y recuperados. La enfermedad no puede brotar de nuevo hasta que el número de susceptibles haya vuelto a acumularse, por ejemplo, como resultado de la descendencia, que nace en el compartimento susceptible.

Cada miembro de la población suele pasar de susceptible a infectado, y a recuperado. Esto se puede mostrar como un diagrama de flujo en el que los cuadros representan los diferentes compartimentos y las flechas la transición entre compartimentos, por ejemplo:

Tasas de transición

Para que el modelo sea completamente específico, las flechas se tienen que etiquetar con las tasas de transición entre los compartimentos. Entre S y I, la tasa de transición es βI, donde β es el número promedio de contactos por persona por tiempo, multiplicado por la probabilidad de transmisión de la enfermedad en un contacto entre un sujeto susceptible y uno infectado. Entre I y R, la tasa de transición es γ (simplemente es la tasa de recuperación o muerte, es decir, el número de recuperados o muertos durante un día, dividido entre el número total de infectados de ese mismo día, suponiendo que el «día» es la unidad de tiempo). Si la duración de la infección es indicada como D, entonces γ = 1/D, ya que un individuo experimenta una recuperación en D unidades de tiempo. Se asume que la permanencia de cada sujeto en los estados epidémicos es una variable aleatoria con distribución exponencial. Se pueden utilizar distribuciones más complejas y realistas (como la distribución de Erlang) con algunas modificaciones.

Tratamiento determinístico biomatemático del modelo SIR

El modelo SIR sin dinámica vital.

La dinámica de una epidemia, por ejemplo la gripe, a menudo es mucho más rápida que la dinámica del nacimiento y la muerte, por lo tanto, el nacimiento y la muerte suelen omitirse en modelos compartimentales simples. El sistema SIR sin la denominada dinámica vital (nacimiento y muerte, a veces denominada demografía), descrito anteriormente, puede expresarse mediante el siguiente conjunto de ecuaciones diferenciales ordinarias:​

${\displaystyle {\begin{aligned}&{\frac {dS}{dt}}=-{\frac {\beta IS}{N}},\\[6pt]&{\frac {dI}{dt}}={\frac {\beta IS}{N}}-\gamma I,\\[6pt]&{\frac {dR}{dt}}=\gamma I,\end{aligned}}}$

donde ${\displaystyle S}$ es el conjunto de la población susceptible, ${\displaystyle I}$ es el conjunto de infectados, ${\displaystyle R}$ el de la población sin la enfermedad (por muerte o recuperación) y ${\displaystyle N}$ la suma de los tres.

Este modelo fue propuesto por primera vez por O. Kermack y Anderson Gray Mckendrick, como un caso especial de lo que ahora llamamos la teoría de Kermack y McKendrick, y dio seguimiento al trabajo de McKendrick con Ronald Ross. Este sistema es no lineal, sin embargo, es posible derivar su solución analítica de forma cerrada. ​ Otras herramientas numéricas incluyen métodos de Monte Carlo, como el algoritmo Gillespie. Primero, nótese que:

${\displaystyle {\frac {dS}{dt}}+{\frac {dI}{dt}}+{\frac {dR}{dt}}=0,}$

a lo que sigue:

${\displaystyle S(t)+I(t)+R(t)={\text{constant}}=N,}$

expresando en términos matemáticos la constante de población ${\displaystyle N}$. Nótese que la relación anterior implica que tan solo se tiene que estudiar dos de las tres ecuaciones.

En segundo lugar, observamos que la dinámica de la clase infecciosa depende de la siguiente relación:

${\displaystyle R_{0}={\frac {\beta }{\gamma }},}$

el llamado número básico de reproducción (también llamada ratio básica reproductiva). Esta ratio está derivada como el número de nuevas infecciones esperadas (estas nuevas infecciones se llaman secundarias) de una sola infección en una población donde todos los sujetos son susceptibles. ​​ Esta idea se puede comprobar mejor si se establece el tiempo típico entre contactos como ${\displaystyle T_{c}=\beta ^{-1}}$ y el tiempo típico hasta que no se tiene la enfermedad como ${\displaystyle T_{r}=\gamma ^{-1}}$. Desde aquí sigue que, en general, el número de contactos con una persona infecciosa antes de que ya no tenga la enfermedad es: ${\displaystyle T_{r}/T_{c}.}$

Dividiendo la primera ecuación diferencial entre un tercio, separando las variables e integrando los datos obtenemos

${\displaystyle S(t)=S(0)e^{-R_{0}(R(t)-R(0))/N},}$

donde ${\displaystyle S(0)}$ y ${\displaystyle R(0)}$ son los números iniciales de, respectivamente, susceptibles y sin enfermedad. Si escribimos ${\displaystyle s_{0}=S(0)/N}$ para la proporción inicial de susceptibles y ${\displaystyle s_{\infty }=S(\infty )/N}$ and ${\displaystyle r_{\infty }=R(\infty )/N}$ para los susceptibles y sin enfermedad, respectivamente, en el límite ${\displaystyle t\to \infty ,}$ se obtiene

${\displaystyle s_{\infty }=1-r_{\infty }=s_{0}e^{-R_{0}(r_{\infty }-r_{0})}}$

(nótese que el compartimento de infecciosos se vacía en este límite).

Esta ecuación trascendente tiene una solución en términos de la función W de Lambert,​ que es:

${\displaystyle s_{\infty }=1-r_{\infty }=-R_{0}^{-1}\,W(-s_{0}R_{0}e^{-R_{0}(1-r_{0})}).}$

Esto muestra que al final de una epidemia, a no ser que ${\displaystyle s_{0}=0}$, no todos los individuos de una población no tienen la enfermedad, por lo tanto algunos pueden ser susceptibles. Esto quiere decir que el final de una epidemia está causada por el descenso en el número de infecciosos y no por la falta de sujetos susceptibles.

El rol del número básico de reproducción es extremadamente importante. De hecho, si reescribimos la ecuación para un individuo infeccioso de la siguiente manera:

${\displaystyle {\frac {dI}{dt}}=\left(R_{0}{\frac {S}{N}}-1\right)\gamma I,}$

vemos que si:

${\displaystyle R_{0}>{\frac {N}{S(0)}},}$

entonces:

${\displaystyle {\frac {dI}{dt}}(0)>0,}$

habrá un brote epidémico con un número más alto de casos (lo que puede alcanzar a una gran parte de la población). Por el contrario:

${\displaystyle R_{0}<{\frac {N}{S(0)}},}$

entonces:

${\displaystyle {\frac {dI}{dt}}(0)<0,}$

independientemente del tamaño inicial de la población susceptible, la enfermedad nunca puede causar un brote epidémico. Como consecuencia, está claro que el [[número básico reproductivo es muy importante.

La fuerza de la infección

Hay que tener en cuenta que en el modelo anterior la función:

${\displaystyle F=\beta I,}$

modela la ratio de transición del compartimento susceptible al infeccioso, así que se llama fuerza de la infección. Sin embargo, para clases de enfermedades más comunicables, es más realista considerar la fuerza de la infección que no depende del número absoluto de infecciosos, pero de la fracción (con respecto a la constante total de población ${\displaystyle N}$):

${\displaystyle F=\beta {\frac {I}{N}}.}$

Capasso y, después, otros autores han propuesto fuerzas no-lineares para modelar de manera más realista el proceso de contagio.

Soluciones analíticas exactas al modelo SIR

En el 2014, Harko y coautores derivaron una solución analítica exacta al modelo SIR. ​ En el caso sin dinámica vital, para ${\displaystyle {\mathcal {S}}(u)=S(t)}$, etc. , corresponde a la siguiente parametización:

${\displaystyle {\mathcal {S}}(u)=S(0)u}$

${\displaystyle {\mathcal {I}}(u)=N-{\mathcal {R}}(u)-{\mathcal {S}}(u)}$

${\displaystyle {\mathcal {R}}(u)=R(0)-\rho \ln(u)}$

para

${\displaystyle t={\frac {N}{\beta }}\int _{u}^{1}{\frac {du^{*}}{u^{*}{\mathcal {I}}(u^{*})}},\quad \rho ={\frac {\gamma N}{\beta }},}$

con condiciones iniciales

${\displaystyle ({\mathcal {S}}(1),{\mathcal {I}}(1),{\mathcal {R}}(1))=(S(0),N-R(0)-S(0),R(0)),\quad u_{T}<u<1,}$

donde ${\displaystyle u_{T}}$ satisface ${\displaystyle {\mathcal {I}}(u_{T})=0}$. Con la ecuación trascendental ${\displaystyle R_{\infty }}$ anterior, sigue que ${\displaystyle u_{T}=e^{-(R_{\infty }-R(0))/\rho }(=S_{\infty }/S(0)}$, si ${\displaystyle S(0)\neq 0)}$ y${\displaystyle I_{\infty }=0}$.

Se encontró una solución analítica equivalente a los rendimientos de Miller:​​

${\displaystyle {\begin{aligned}S(t)&=S(0)e^{-\xi (t)}\\[8pt]I(t)&=N-S(t)-R(t)\\[8pt]R(t)&=R(0)+\rho \xi (t)\\[8pt]\xi (t)&={\frac {\beta }{N}}\int _{0}^{t}I(t^{*})\,dt^{*}\end{aligned}}}$

Aquí ${\displaystyle \xi (t)}$ se puede interpretar como el número esperado de transmisiones que un individuo ha recibido por tiempo ${\displaystyle t}$. Estas dos soluciones están relacionadas por ${\displaystyle e^{-\xi (t)}=u}$

Efectivamente, el mismo resultado puede encontrarse en el trabajo original de Kermack y Mckendrick. ​

Estas soluciones se entienden mejor si vemos que todos los términos en la parte de la derecha de las ecuaciones diferenciales son proporcionales a ${\displaystyle I}$. Las ecuaciones pueden ser divididas por ${\displaystyle I}$, y el tiempo reescalado para que el operador en la izquierda se convierta simplemente en ${\displaystyle d/d\tau }$, donde ${\displaystyle d\tau =Idt}$, por ejemplo ${\displaystyle \tau =\int Idt}$. Las ecuaciones diferenciales son ahora lineares y la tercera, con forma ${\displaystyle dR/d\tau =}$ const. , muestra que ${\displaystyle \tau }$ y ${\displaystyle R}$ (y arriba ${\displaystyle \xi }$) se crean simplemente lineares.

El modelo SIR con dinámica vital y población constante

Considere una población caracterizada por una tasa de mortalidad ${\displaystyle \mu }$ y una ratio de nacimiento ${\displaystyle \Lambda }$, y donde enfermedades comunicables se propagan. El modelo de transmisión es:

${\displaystyle {\begin{aligned}{\frac {dS}{dt}}&=\Lambda -\mu S-{\frac {\beta IS}{N}}\\[8pt]{\frac {dI}{dt}}&={\frac {\beta IS}{N}}-\gamma I-\mu I\\[8pt]{\frac {dR}{dt}}&=\gamma I-\mu R\end{aligned}}}$

para el cual el equilibrio sin enfermedad (DFE, en inglés) es:

${\displaystyle \left(S(t),I(t),R(t)\right)=\left({\frac {\Lambda }{\mu }},0,0\right).}$

En este caso, podemos derivar el número de reproducción básico:

${\displaystyle R_{0}={\frac {\beta \Lambda }{\mu (\mu +\gamma )}},}$

que tiene propiedades de entrada. De hecho, independientemente de los valores biológicos iniciales, podemos ver que:

${\displaystyle R_{0}\leq 1\Rightarrow \lim _{t\to \infty }(S(t),I(t),R(t))={\textrm {DFE}}=\left({\frac {\Lambda }{\mu }},0,0\right)}$

${\displaystyle R_{0}>1,I(0)>0\Rightarrow \lim _{t\to \infty }(S(t),I(t),R(t))={\textrm {EE}}=\left({\frac {\gamma +\mu }{\beta }},{\frac {\mu }{\beta }}\left(R_{0}-1\right),{\frac {\gamma }{\beta }}\left(R_{0}-1\right)\right).}$

El punto EE, es el Equilibrio Endémico (la enfermedad no está erradicada y se mantiene en la población). Con argumentos heurísticos, uno puede ver que ${\displaystyle R_{0}}$ se puede leer como el promedio de infecciones causadas por un solo sujeto infeccioso en una población susceptible. La relación biológica anterior implica que si el número es igual o menor que 1 la enfermedad se extingue; por el contrario, si es mayor, la enfermedad se quedara permanentemente endémica en la población.

Modelo SIS

Algunas infecciones, por ejemplo las de la gripe común y la influenza, no confieren inmunidad duradera. Estas infecciones no dan inmunización al recuperarse de la infección, y los individuos vuelven a ser susceptibles.

Tenemos el modelo:

${\displaystyle {\begin{aligned}{\frac {dS}{dt}}&=-{\frac {\beta SI}{N}}+\gamma I\\[6pt]{\frac {dI}{dt}}&={\frac {\beta SI}{N}}-\gamma I\end{aligned}}}$

Nótese que denominar con N al total de la población:

${\displaystyle {\frac {dS}{dt}}+{\frac {dI}{dt}}=0\Rightarrow S(t)+I(t)=N}$.

A lo que sigue:

${\displaystyle {\frac {dI}{dt}}=(\beta -\gamma )I-{\frac {\beta }{N}}I^{2}}$,

Las dinámicas infecciosas se rigen por una logística funcional que implica ${\displaystyle \forall I(0)>0}$:

${\displaystyle {\begin{aligned}&{\frac {\beta }{\gamma }}\leq 1\Rightarrow \lim _{t\to +\infty }I(t)=0,\\[6pt]&{\frac {\beta }{\gamma }}>1\Rightarrow \lim _{t\to +\infty }I(t)=\left(1-{\frac {\gamma }{\beta }}\right)N.\end{aligned}}}$

Es posible encontrar una solución analítica a este modelo (transformando las variables: ${\displaystyle I=y^{-1}}$ y sustituyéndolas por ecuaciones de campo),​ para que la ratio de reproducción básica sea mayor que la unidad. La solución se presenta como:

${\displaystyle I(t)={\frac {I_{\infty }}{1+Ve^{-\chi t}}}}$.

where ${\displaystyle I_{\infty }=(1-\gamma /\beta )N}$ is the endemic infectious population, ${\displaystyle \chi =\beta -\gamma }$, and ${\displaystyle V=I_{\infty }/I_{0}-1}$. As the system is assumed to be closed, the susceptible population is then ${\displaystyle S(t)=N-I(t)}$.

Como caso especial, se obtiene la función logística usual asumiendo que ${\displaystyle \gamma =0}$. Esto también se puede considerar en el modelo SIR con ${\displaystyle R=0}$, en este caso no se quitará nada. Ese es el modelo SI. ​ La ecuación diferencial usando ${\displaystyle S=N-I}$ se reduce a:

${\displaystyle {\frac {dI}{dt}}\propto I\cdot (N-I).}$

A la larga, en este modelo, todos los individuos se infectan.

Explicaciones sobre el modelo SIR básico

El modelo MSIR

En muchas infecciones, incluyendo el sarampión, los bebés no nacen en el compartimento susceptible, sino que son inmunes a la enfermedad durante los primeros meses de vida debido a la protección de los anticuerpos maternos (transmitidos a través de la placenta y adicionalmente a través del calostro). Esto se denomina inmunidad pasiva. Este detalle adicional puede demostrarse con una clase M (para la inmunidad maternal) al inicio del modelo

Para indicar esto matemáticamente, se añade un compartimento adicional, M(t), que tiene como resultado las siguientes ecuaciones diferenciales:

Para indicar esto matemáticamente, se añade un compatimento adicional M(t). Esto resulta en las siguientes ecuaciones:

${\displaystyle {\begin{aligned}{\frac {dM}{dT}}&=\Lambda -\delta M-\mu M\\[8pt]{\frac {dS}{dT}}&=\delta M-{\frac {\beta SI}{N}}-\mu S\\[8pt]{\frac {dI}{dT}}&={\frac {\beta SI}{N}}-\gamma I-\mu I\\[8pt]{\frac {dR}{dT}}&=\gamma I-\mu R\end{aligned}}}$

Estado de portador

Algunas personas que han tenido una enfermedad infecciosa como la tuberculosis nunca se recuperan completamente y siguen llevando la infección, pese a que no sufren la enfermedad. Pueden volver a entrar al compartimento de infectados y sufrir síntomas (como en la tuberculosis) o pueden seguir infectando a otros en su estado de portador, a pesar de no padecer síntomas. El ejemplo más famoso de esto es, probablemente, Mary Mallon, que infectó a 22 personas con fiebre tifoidea. El compartimento de portadores está marcado con C.

Modelo SEIR

Para muchas infecciones importantes existe un período de incubación significativo durante el cual los individuos han sido infectados, pero aún no son infecciosos por sí mismos. Durante este periodo el individuo se encuentra en el compartimento E (de exposición).

Asumiendo que el período de incubación es una variable aleatoria con distribución exponencial con parámetro${\displaystyle a}$ (el periodo medio de incubación es ${\displaystyle a^{-1}}$) y asumiendo la presencia de dinámicas vitales en la ratio de nacimiento ${\displaystyle \Lambda }$ iguales a la ratio de muertes ${\displaystyle \mu }$, obtenemos este modelo:

${\displaystyle {\begin{aligned}{\frac {dS}{dt}}&=\mu N-\mu S-\beta {\frac {I}{N}}S\\[8pt]{\frac {dE}{dt}}&=\beta {\frac {I}{N}}S-(\mu +a)E\\[8pt]{\frac {dI}{dt}}&=aE-(\gamma +\mu )I\\[8pt]{\frac {dR}{dt}}&=\gamma I-\mu R.\end{aligned}}}$

Tenemos ${\displaystyle S+E+I+R=N,}$ pero tan solo es constante porque la asunción (degenerativa) de que las muertes y nacimientos son iguales; en general ${\displaystyle N}$ es variable.

Para este modelo, el número de reproducción básico es:

${\displaystyle R_{0}={\frac {a}{\mu +a}}{\frac {\beta }{\mu +\gamma }}.}$

Similar al modelo SIR, tenemos el equilibrio sin enfermedad (DFE) (N,0,0,0) y el Equilibrio Endémico (EE); se puede ver que, de manera independiente a las condiciones biológicas iniciales:

${\displaystyle \left(S(0),E(0),I(0),R(0)\right)\in \left\{(S,E,I,R)\in [0,N]^{4}:S\geq 0,E\geq 0,I\geq 0,R\geq 0,S+E+I+R=N\right\}}$

esto mantiene:

${\displaystyle R_{0}\leq 1\Rightarrow \lim _{t\to +\infty }\left(S(t),E(t),I(t),R(t)\right)=DFE=(N,0,0,0),}$

${\displaystyle R_{0}>1,I(0)>0\Rightarrow \lim _{t\to +\infty }\left(S(t),E(t),I(t),R(t)\right)=EE.}$

En el caso de ratios de contacto que varían periódicamente ${\displaystyle \beta (t)}$, la condición para la atracción global del DFE es siguiendo un sistema linear con coeficientes periódicos:

${\displaystyle {\begin{aligned}{\frac {dE_{1}}{dt}}&=\beta (t)I_{1}-(\gamma +a)E_{1}\\[8pt]{\frac {dI_{1}}{dt}}&=aE_{1}-(\gamma +\mu )I_{1}\end{aligned}}}$

es estable (tiene los valores de Floquet dentro de la unidad círculo en el plano complejo)

Modelo SEIS

El modelo SEIS contempla el período de exposición, o periodo latente de la enfermedad, dando lugar a un compartimento adicional, E(t).

${\displaystyle {\color {blue}{{\mathcal {S}}\to {\mathcal {E}}\to {\mathcal {I}}\to {\mathcal {S}}}}}$

En este modelo, una infección no deja inmune a los individuos que se han recuperado; estos vuelven a ser suscptibles, volviendo al compartimentoS(t). Le ecuación diferencial siguiente describe el modelo:

${\displaystyle {\begin{aligned}{\frac {dS}{dT}}&=\Lambda -{\frac {\beta SI}{N}}-\mu S+\gamma I\\[6pt]{\frac {dE}{dT}}&={\frac {\beta SI}{N}}-(\epsilon +\mu )E\\[6pt]{\frac {dI}{dT}}&=\varepsilon E-(\gamma +\mu )I\end{aligned}}}$

Modelo MSEIR

En el caso de una enfermedad con los factores de inmunidad pasiva y un período de latencia, existe el modelo MSEIR.

${\displaystyle \color {blue}{{\mathcal {M}}\to {\mathcal {S}}\to {\mathcal {E}}\to {\mathcal {I}}\to {\mathcal {R}}}}$

${\displaystyle {\begin{aligned}{\frac {dM}{dT}}&=\Lambda -\delta M-\mu M\\[6pt]{\frac {dS}{dT}}&=\delta M-{\frac {\beta SI}{N}}-\mu S\\[6pt]{\frac {dE}{dT}}&={\frac {\beta SI}{N}}-(\varepsilon +\mu )E\\[6pt]{\frac {dI}{dT}}&=\varepsilon E-(\gamma +\mu )I\\[6pt]{\frac {dR}{dT}}&=\gamma I-\mu R\end{aligned}}}$

Modelo MSEIRS

Un modelo MSEIRS es similar al modelo MSEIR, pero la inmunidad en la clase R sería temporal, por lo que los individuos recuperarían su susceptibilidad cuando terminara la inmunidad temporal.

${\displaystyle {\color {blue}{{\mathcal {M}}\to {\mathcal {S}}\to {\mathcal {E}}\to {\mathcal {I}}\to {\mathcal {R}}\to {\mathcal {S}}}}}$

Tasas de contacto variables y epidemias plurianuales o caóticas.

Es bien sabido que la probabilidad de contraer una enfermedad no es constante en el tiempo. Algunas enfermedades son estacionales, tales como los virus del resfriado común, que son más frecuentes durante el invierno. Con las enfermedades de la infancia, tales como el sarampión, las paperas y la rubéola, existe una fuerte correlación con el calendario escolar, por lo que durante las vacaciones escolares la probabilidad de contraer dicha enfermedad disminuye notablemente. Como consecuencia, para muchas clases de enfermedades se debe considerar una fuerza de infección con tasa de contacto periódica y variable ('estacional').

Como consecuencia, para muchas clases de enfermedades, uno debe considerar la fuerza de infección con una ratio de contacto variable periódicamente ('estacional')

${\displaystyle F=\beta (t){\frac {I}{N}},\quad \beta (t+T)=\beta (t)}$

con un periodo T igual a un año.

Por lo tanto, nuestro modelo se convierte en

${\displaystyle {\begin{aligned}{\frac {dS}{dt}}&=\mu N-\mu S-\beta (t){\frac {I}{N}}S\\[8pt]{\frac {dI}{dt}}&=\beta (t){\frac {I}{N}}S-(\gamma +\mu )I\end{aligned}}}$

(las dinámicas de los recuperados siguen la forma ${\displaystyle R=N-S-I}$). Una serie de ecuaciones diferenciales no-lineares con parámetros variables periódicos. Es sabido que esta clase de sistema dinámico puede sufrir interesantes y complejos fenómenos de resonancia paramétrica no-linear. Es fácil de ver si:

${\displaystyle {\frac {1}{T}}\int _{0}^{T}{\frac {\beta (t)}{\mu +\gamma }}\,dt<1\Rightarrow \lim _{t\to +\infty }(S(t),I(t))=DFE=(N,0),}$

donde si la integral es mayor que uno, la enfermedad no morirá y puede haber estas resonancias. Por ejemplo, considerando la ratio de contacto variable periódica como la entrada del sistema, la salida es una función periódica cuyo periodo es un múltiple del periodo de entrada.

Esto permitió ayudar a explicar los brotes epidémicos polianuales (típicamente bienales) de algunas enfermedades infecciosas como interacción entre el período de las oscilaciones de la tasa de contacto y el seudo-período de las fluctuaciones que se encuentran cerca de la región endémica. Notablemente, el comportamiento también puede ser cuasi-periódico o incluso caótico.

Modelización de la vacunación

El modelo SIR se puede modificar para modelar la vacunación. Normalmente estos introducen un compartimento adicional al modelo SIR, ${\displaystyle V}$, para individuos vacunados. ^{[cita requerida]} A continuación, algunos ejemplos.

Vacunando recién nacidos

En presencia de enfermedades transmisibles, una de las principales tareas es la erradicación mediante medidas de prevención y, si es posible, mediante el establecimiento de un programa de vacunación masiva. Considérese una enfermedad por la cual el recién nacido está vacunado (con una vacuna que administra inmunidad de por vida) a un ritmo ${\displaystyle P\in (0,1)}$:

${\displaystyle {\begin{aligned}{\frac {dS}{dt}}&=\mu N(1-P)-\mu S-\beta {\frac {I}{N}}S\\[8pt]{\frac {dI}{dt}}&=\beta {\frac {I}{N}}S-(\mu +\gamma )I\\[8pt]{\frac {dV}{dt}}&=\mu NP-\mu V\end{aligned}}}$

donde ${\displaystyle V}$ es la clase de sujetos vacunados. Esto muestra inmediatamente que::

${\displaystyle \lim _{t\to +\infty }V(t)=NP,}$

entonces tenemos que tratar con un comportamiento a mayor plazo de ${\displaystyle S}$ e${\displaystyle I}$, por lo que:

${\displaystyle R_{0}(1-P)\leq 1\Rightarrow \lim _{t\to +\infty }\left(S(t),I(t)\right)=DFE=\left(N\left(1-P\right),0\right)}$

${\displaystyle R_{0}(1-P)>1,\quad I(0)>0\Rightarrow \lim _{t\to +\infty }\left(S(t),I(t)\right)=EE=\left({\frac {N}{R_{0}(1-P)}},N\left(R_{0}(1-P)-1\right)\right).}$

En otras palabras, si

${\displaystyle P<P^{*}=1-{\frac {1}{R_{0}}}}$

el programa de vacunación no es capaz de erradicar la enfermedad, al contrario, permanecerá endémica, aunque a niveles inferiores que si no hubiera vacunas. Esto quiere decir que el modelo matemático sugiere que para una enfermedad cuyo número de reproducción básico puede llegar a ser 18, se deberían vacunar al menos al 94. 4% de los recién nacidos para erradicar la enfermedad.

Vacunación e información

Las sociedades modernas están enfrentándose al desafío de la exención “racional", es decir, la decisión de la familia de no vacunar a los niños como consecuencia de una comparación "racional" entre el riesgo percibido por la infección y los posibles daños por la vacuna. Para valorar si este comportamiento es realmente racional, es decir, si esto puede igualmente llevar a la erradicación de la enfermedad, simplemente hay que asumir que la tasa de vacunación es una función creciente del número de sujetos infecciosos:

${\displaystyle P=P(I),\quad P'(I)>0.}$

En este caso la condición de erradicación se convierte en:

${\displaystyle P(0)\geq P^{*},}$

La línea de la ratio de vacunación debería ser más alta que el umbral de "vacunación obligatoria" que, en caso de exención, no puede mantenerse. Entonces, la exención "racional" puede ser miope puesto que está basada en la baja incidencia debido a la alta cobertura de la vacuna, en lugar de tener en cuenta un futuro rebrote de la infección debido a un descenso en la cobertura.

Vacunación de los no recién nacidos

En el caso de que además haya vacunas a no recién nacidos a una tasa ρ, la ecuación para el sujeto susceptible y vacunado se debe modificar de la siguiente manera :

${\displaystyle {\begin{aligned}{\frac {dS}{dt}}&=\mu N(1-P)-\mu S-\rho S-\beta {\frac {I}{N}}S\\[8pt]{\frac {dV}{dt}}&=\mu NP+\rho S-\mu V\end{aligned}}}$

que lleva a la siguiente condición de erradicación:

${\displaystyle P\geq 1-\left(1+{\frac {\rho }{\mu }}\right){\frac {1}{R_{0}}}}$

Estrategia de vacunación frecuente

Esta estrategia vacuna repetidamente a un grupo de edad definido (como niños pequeños o ancianos) en una población susceptible a lo largo del tiempo. Con esta estrategia, se retira inmediatamente al grupo de individuos susceptibles, lo que permite la eliminación de una enfermedad infecciosa, (como el sarampión), de toda la población. Cada unidad de tiempo T, una fracción constante p de sujetos susceptibles se vacunan en un tiempo relativamente corto (respecto a la dinámica de la enfermedad). Esto lleva a las siguientes ecuaciones diferenciales impulsivas para los sujetos susceptibles y vacunados:

${\displaystyle {\begin{aligned}{\frac {dS}{dt}}&=\mu N-\mu S-\beta {\frac {I}{N}}S,\quad S(nT^{+})=(1-p)S(nT^{-}),&&n=0,1,2,\ldots \\[8pt]{\frac {dV}{dt}}&=-\mu V,\quad V(nT^{+})=V(nT^{-})+pS(nT^{-}),&&n=0,1,2,\ldots \end{aligned}}}$

Es fácil de ver que si proponemos I = 0 se obtiene que la dinámica de los sujetos susceptibles se da por:

${\displaystyle S^{*}(t)=1-{\frac {p}{1-(1-p)E^{-\mu T}}}E^{-\mu MOD(t,T)}}$

y la condición de erradicación es:

${\displaystyle R_{0}\int _{0}^{T}S^{*}(t)\,dt<1}$

La influencia de la edad: modelos estructurados por edad

La edad influye profundamente en la propagación de la enfermedad en una población, especialmente en la tasa de contacto. Esta tasa resume la efectividad de los contactos entre sujetos susceptibles e infecciosos. Teniendo en cuenta las edades de las clases epidémicas ${\displaystyle s(t,a),i(t,a),r(t,a)}$ (to limit ourselves to the susceptible-infectious-removed scheme) such that:

${\displaystyle S(t)=\int _{0}^{a_{M}}s(t,a)\,da}$

${\displaystyle I(t)=\int _{0}^{a_{M}}i(t,a)\,da}$

${\displaystyle R(t)=\int _{0}^{a_{M}}r(t,a)\,da}$

(donde${\displaystyle a_{M}\leq +\infty }$ es la edad máxima admisible) y las dinámicas no se describen, como uno pensaría, con una "simple" ecuación diferencial parcial, sino por ecuaciones integro-diferenciales:

${\displaystyle \partial _{t}s(t,a)+\partial _{a}s(t,a)=-\mu (a)s(a,t)-s(a,t)\int _{0}^{a_{M}}k(a,a_{1};t)i(a_{1},t)\,da_{1}}$

${\displaystyle \partial _{t}i(t,a)+\partial _{a}i(t,a)=s(a,t)\int _{0}^{a_{M}}{k(a,a_{1};t)i(a_{1},t)da_{1}}-\mu (a)i(a,t)-\gamma (a)i(a,t)}$

${\displaystyle \partial _{t}r(t,a)+\partial _{a}r(t,a)=-\mu (a)r(a,t)+\gamma (a)i(a,t)}$

donde:

${\displaystyle F(a,t,i(\cdot ,\cdot ))=\int _{0}^{a_{M}}k(a,a_{1};t)i(a_{1},t)\,da_{1}}$

es la fuerza de infección que, por supuesto, dependerá del contactp kernel ${\displaystyle k(a,a_{1};t)}$ en la interacción entre edades.

La complejidad se agrega por las condiciones iniciales de los recién nacidos (es decir, por a=0), que son fáciles de infectar y eliminar:

${\displaystyle i(t,0)=r(t,0)=0}$

pero que son no-locales para la densidad de los recién nacidos susceptibles:

${\displaystyle s(t,0)=\int _{0}^{a_{M}}\left(\varphi _{s}(a)s(a,t)+\varphi _{i}(a)i(a,t)+\varphi _{r}(a)r(a,t)\right)\,da}$

donde ${\displaystyle \varphi _{j}(a),j=s,i,r}$ son las fertilidades de los adultos.

Además, definiendo ahora la densidad de la población total ${\displaystyle n(t,a)=s(t,a)+i(t,a)+r(t,a)}$ se obtiene:

${\displaystyle \partial _{t}n(t,a)+\partial _{a}n(t,a)=-\mu (a)n(a,t)}$

En el caso más simple de fertilidades iguales en las tres clases epidémicas, se obtiene que para tener equilibrio demográfico la siguiente necesaria y suficiente condición vinculando la fertilidad ${\displaystyle \varphi (.)}$ con la mortalidad ${\displaystyle \mu (a)}$ se tiene que mantener:

${\displaystyle 1=\int _{0}^{a_{M}}\varphi (a)\exp \left(-\int _{0}^{a}{\mu (q)dq}\right)\,da}$

y el equilibrio demográfico es:

${\displaystyle n^{*}(a)=C\exp \left(-\int _{0}^{a}\mu (q)\,dq\right),}$

automáticamente asegurando la existencia de una solución libre de enfermedad:

${\displaystyle DFS(a)=(n^{*}(a),0,0).}$

Un número de reproducción básico se puede calcular como el radio espectral de un operador funcional apropiado.

Otras consideraciones dentro de los modelos compartimentales epidémicos

Transmisión vertical

En el caso de algunas enfermedades como el SIDA y la Hepatitis B, es posible que la descendencia de los padres infectados nazca con la infección. La transmisión de la enfermedad por medio de la madre se denomina transmisión vertical. La afluencia de miembros adicionales en la categoría infectada puede considerarse dentro del modelo al incluir una fracción de los miembros recién nacidos en el compartimento infectado.

Transmisión vectorial

Las enfermedades transmitidas indirectamente entre seres humanos, - por ejemplo, la malaria que se transmite por medio de mosquitos - se transmiten a través de un vector. En estos casos, las transferencias de infección de humanos al insecto y un modelo epidémico deben incluir ambas especies, generalmente requiriendo más compartimentos que un modelo de transmisión directa.

Otros

Otros casos que necesitan considerarse al modelar una epidemia incluyen cosas como las siguientes: Mezclas no homogéneas Infectividad variable Distribuciones que son espacialmente no uniformes Enfermedades causadas por macroparásitos

Modelos epidémicos determinísticos contra estocásticos

Es importante destacar que los modelos determinísticos que aquí se presentan son válidos únicamente en poblaciones suficientemente grandes y en consecuencia deben ser utilizados con precaución. Para ser más precisos, estos modelos sólo son válidos en el límite termodinámico, donde la población es efectivamente infinita. En modelos estocásticos, el equilibrio endémico a largo plazo obtenido previamente, no se mantiene, ya que existe una probabilidad finita de que el número de individuos infectados disminuya por debajo de uno en un sistema. En un sistema real, el patógeno puede no propagarse, ya que ningún huésped estaría infectado. Sin embargo, en modelos de campo determinísticos, el número de infectados puede asumir valores reales, concretamente, valores no enteros de huéspedes infectados, y el patógeno aún puede persistir en el sistema con un número de huéspedes que es menor a uno, pero más de cero.

Véase también

• Modelización matemática en epidemiología
• Problema de unidad de área modificable
• Matriz de última generación
• Asesoramiento de riesgo
• Ratio de ataque

Lectura adicional

• May, Robert M.; Anderson, Roy M. (1991). Infectious diseases of humans: dynamics and control. Oxford: Oxford University Press. ISBN 0-19-854040-X. 
• Capasso, V. (1993). Mathematical Structure of Epidemic Systems. Berlin: Springer. ISBN 3-540-56526-4. 
• Vynnycky, E.; White, R. G., eds. (2010). An Introduction to Infectious Disease Modelling. Oxford: Oxford University Press. ISBN 978-0-19-856576-5. 
• N. Bacaër, J. Ripoll, R. Bravo de la Parra, X. Bardina, S. Cuadrado: Matemáticas y epidemias (2021) Paris, ISBN 979-10-343-8464-8. Pdf

Enlaces externos

Modelo SIR: Experimentos en línea con JSXGraph

• SIR model: Online experiments with JSXGraph
• «Simulating an epidemic». 3Blue1Brown. 27 de marzo de 2020 – via YouTube.