Fotografía de Giacomo Carra, tomada de https://unsplash.com/photos/gf6UDwpl0ac
1. Representación gráfica y pronósticos
Se usó una calculadora en línea [1] para la solución numérica y la representación gráfica del modelo epidemiológico SEIR para México y Michoacán: Susceptible (Susceptibles) → Exposed (Expuestos) → Infected (Infectados) → Removed (Removidos). El modelo también representa Hospitalized (Hospitalizados) y Fatalities (Fallecimientos). Se consideraron solamente pacientes hospitalizados de terapia intensiva (UCI y/o intubación).
Los modelos y la metodología están explicados en las secciones siguientes. En la Sección 4 se aclara por qué se multiplican los números de los gráficos con un factor 0.04, en el caso de México y con un factor 0.05, en el caso de Michoacán para obtener estimaciones realistas.
Conclusiones para México:
Día 0: 25 de marzo |
Duración: mayor que 200 días |
Pico infectados: 13-15 de julio |
Pico fallecidos: 18-20 de julio |
Infectados: 3 300 000 |
Infectados reportados: 660 000 |
Fallecidos: 110 000 |
Fallecidos reportados: 98 000 |
Pico fallecidos: 1 400 |
Pico fallecidos reportado: 1 200 |
Máximo de pacientes de UCI: 22 000 |
|
Conclusiones para Michoacán:
Día 0: 29 de marzo |
Duración: mayor que 200 días |
Pico infectados: 28-30 de julio |
Pico fallecidos: 2-4 de agosto |
Infectados: 137 000 |
Infectados reportados: 27 000 |
Fallecidos: 4 500 |
Fallecidos reportados: 4 000 |
Pico fallecidos: 60 1000 |
Pico fallecidos reportado: 50 |
Máximo de pacientes de UCI: 1000 |
|
2. El modelo matemático
Los cálculos están basados en el modelo epidemiológico conocido como SEIR con compartimentos. Se considera una división de una población de N personas en los conjuntos disjuntos Susceptibles, Expuestos, Infectados, Recuperados, Hospitalizados de terapia intensiva y Fallecidos con la siguiente dinámica indicada por las flechas:
Esta dinámica se puede expresar en ecuaciones de diferencias contando los días n desde un día inicial:
1.
El factor r0 es conocido como factor reproductivo básico y representa el número de nuevos contagios causados por I(n) infectados en el periodo de días infecciosos Tinf . Si cada persona infectada contagia al inicio de la epidemia a r0 personas, entonces en un día se infectan en promedio personas. Cuando el número de personas susceptibles se reduce, el virus encuentra con probabilidad de a una persona susceptible, por lo tanto .
2.
Las personas que se infectan según la Ecuación 1, se encuentran primero en un periodo de incubación y son consideradas \expuestos". Después de Tinc se vuelven \infectados", entonces en un día personas dejan de ser expuestos.
3.
Los expuestos de la Ecuación 2 se vuelven infectados y dejan de ser infectados después de Tinf .
4.
Se considera que la p parte de los infectados son casos graves y entran por Th días en un hospital.
5.
La 1—p parte de los infectados y la q parte de los hospitalizados se recuperan.
6.
La 1—p parte de los hospitalizados mueren.
Aunque se podrían resolver numéricamente las ecuaciones de diferencia, los matemáticos prefieren convertirlas en ecuaciones diferenciales obteniendo el siguiente sistema dinámico:
donde el parámetro continuo t denota el tiempo a partir de un \día 0". En este modelo, la unidad de tiempo es día.
Para resolver este sistema de ecuaciones diferenciales, se tienen que conocer los valores de los parámetros y las condiciones iniciales, es decir, los valores en un día específico. En la próxima sección se explica cómo se determinan estos valores.
3. Metodología
Días Infecciosos: Tinf = 7.2
Para determinar este valor, se calculó el promedio de los días entre presentar síntomas y el ingreso como paciente, considerando que un paciente no contagia más personas por estar en un hospital o en cuarentena.
Cabe mencionar que 80% de las personas infectadas presentan pocos síntomas y no son detectadas por pruebas. Por otro lado, un artículo en Science [3] afirma que el Factor de Dispersión es muy bajo, concluyendo que 10% de los infectados causan 80% de los contagios. Si se aplican consistentemente pruebas a las personas con síntomas, se puede retirar a las personas de alto contagio (super spreaders) controlando así el valor de Tinf y la propagación del virus. De este modo, el valor Tinf refleja directamente las intervenciones gubernamentales.
Se consideraron 177 pacientes positivos de Morelia y se calculó un promedio de 5.2 días. A este valor, se sumaron 2 días contagiosos antes de presentar síntomas obteniendo el valor Tinf = 7.2. El mismo valor se aplica para México y Michoacán.
Factor reproductivo básico: r0 > 1.6
Al inicio de la epidemia se tiene y el factor reproductivo básico representa las nuevas infecciones causadas por I(n) infectados después de Tinf días, entonces , por lo tanto .
Suponiendo que el factor reproductivo diario es constante, se obtiene
El problema es que el número de infectados confirmados está sujeto a la disponibilidad de pruebas y las capacidades de los laboratorios. Por otro lado, la mayoría de las personas gravemente afectadas van a acudir a un hospital, entonces el número de defunciones refleja mejor la dinámica de la epidemia. Considerando que un cierto ratio de los infectados muere, se obtiene , donde denota el número total de fallecidos hasta el día n. Para mitigar las fluctuaciones diarias, se consideró el
promedio semanal , entonces .
Esta fórmula para calcular el factor reproductivo es más precisa que la fórmula a base de los infectados confirmados, por los problemas mencionados. Las cifras de fallecidos están disponibles en los datos abiertos de la Secretaría de Salud [2]. Para determinar rd para México, se tomó el promedio de los factores reproductivos diarios del periodo 1 de abril hasta 6 de mayo, y para Michoacán se consideró el periodo del 1 de mayo hasta el 26 de mayo. En ambos casos, el resultado fue rd = 1.067. Tomando en cuenta que las personas que no son detectadas por la pruebas son infecciosas por un periodo mayor que 7.2 días, se obtiene .
Tiempo de incubación: Tinc = 3.2
Se consideró un promedio de 5.2 días hasta presentar síntomas y que una persona es infecciosa 2 días antes de presentar síntomas [4].
Tiempo en terapia intensiva: Tinc = 14
Se usó el valor mencionado en El Universal [5].
Ratio de sobrevivientes de terapia intensiva: q = 0.4
Hasta el día 17 de junio se registraron 2961 pacientes sobrevivientes intubados y/o de Unidades de Cuidado Intensivo (UCI). Además, del los 19 080 fallecidos confirmados, solamente 3 965 ingresaron a una UCI y/o fueron intubados [2]. Entonces el ratio de sobrevivientes de terapia intensiva es
Ratio de infectados que requieren terapia intensiva: p = 0.04
Un estudio publicado en El Lancet [6] estimó que la mortalidad es del 3 %. Considerando que todos los fallecidos son casos graves y el ratio de sobrevivientes de terapia intensiva es 0.4, se estima que 3% + 0.4 x 3% ≈ 4% de los infectados requieren una terapia intensiva (UCI y/o intubación).
Días entre prueba positiva y defunción: T∆ = 5
Usando los datos abiertos [2], se calculó el promedio de días entre el ingreso como paciente y la defunción para todos los fallecidos en Michoacán hasta el 28 de mayo.
Condiciones iniciales: I(25 de marzo) = 6200 para México, I(29 de marzo) = 125 para Michoacán
Como se mencionó al final de la Sección 2, se tiene que determinar en un día especifico para poder resolver las ecuaciones diferenciales. Desafortunadamente se desconoce el número verdadero de infectados porque las pruebas solamente detectan los casos graves que acuden a un hospital. De modo similar al cálculo del factor reproductivo, es más confiable basar el análisis en el número de defunciones.
El 19 de junio hubo 22 874 fallecidos en México, de los cuales 2 425 no entraron en un hospital. Por eso se estima que aproximadamente el 10% de los fallecidos por Covid-19 no acuden a un hospital y por lo tanto nunca son registrados. Además hubo 1 567 defunciones con el resultado pendiente por más que 5 días, entonces aproximadamente el 7% de los fallecidos que acudieron a un hospital se quedan con el resultado pendiente. Por otro lado, se tienen 191 410 pacientes positivos y 251 355 pacientes negativos registrados, entonces se concluye que menos que la mitad de los fallecidos con resultado pendiente murieron por Covid-19. Como una estimación aproximada, se considera que se debe sumar una 1/8 parte de los fallecidos por Covid-19 que no son detectados por las pruebas. Anteriormente se observó un promedio de T∆ = 5 días entre el ingreso como paciente y la defunción. Suponiendo una mortalidad de 3% [6], los 1.125 x 22 874 fallecidos hasta el 19 de junio corresponden a casos positivos hasta el 14 de junio. Por otro lado, hubo 169 139 casos confirmados hasta el 14 de junio, entonces las pruebas detectan solamente de los infectados. El 80% de los infectados sin aplicar una prueba corresponde principalmente a los casos leves.
Para determinar las condiciones iniciales, se considera que el virus fue \importado" a México y Michoacán por personas infecciosas. Considerando un día temprano de la epidemia, se puede suponer que aún no hubo recuperados, hospitalizados ni fallecidos, entonces . Por simplicidad, se supone que todas las personas que importaron el virus se encontraron en el periodo infeccioso, es decir, . Debido a que nadie es inmune a este virus, el número de susceptibles es aproximadamente la población total de México y Michoacán.
Para Michoacán se consideró el día = 29 de marzo con 25 casos confirmados, lo que corresponde a . Para México se consideró = 25 de marzo con 1246 casos confirmados lo que corresponde a .
Intervenciones gubernamentales: 22 de abril
El día 21 de abril fue declarada la fase 3 en todo México. Se considera que la aplicación de medidas de aislamiento causó después del 21 de abril un decremento gradual del factor reproductivo al valor calculado. Sin intervención gubernamental, el factor reproductivo puede ser mayor que 3.28 [7].
4. Validez del modelo
Una gran parte de la población no participa en la propagación del virus por evitar exitosamente un contagio o por vivir en lugares remotos. Se observa que el modelo es invariante bajo un cambio de escala, es decir, si es una solución de las ecuaciones diferenciales , entonces
por lo tanto es una solución con la condiciones iniciales . La parte de la población total representa a las personas que se quedan en cuarentena estricta (personas que evitan contagios o viven en un lugar remoto), mientras la representa la parte de la población que participa activamente en la vida pública.
Para estimar el valor de, basta comparar un valor confiable con la predicción del modelo. Hasta el día 14 de junio hubo 20 621 fallecidos confirmados en México pero el modelo pronosticó 522 329 fallecidos, entonces . Análogamente, hasta el 19 de junio hubo 373 fallecidos confirmados en Michoacán y el modelo pronosticó 8 087 fallecidos, entonces .
Agradecimiento
Se extiende un agradecimiento especial al Dr. Alfredo Raya, Director del Instituto de Física y Matemáticas por su paciencia y valiosas discusiones sobre el tema.
Referencias
[1] https://gabgoh.github.io/COVID/index.html
[2] https://www.gob.mx/salud/documentos/datos-abiertos-152127
[3] https://www.sciencemag.org/news/2020/05/why-do-some-covid-19-patients-infect-many-others-
whereas-most-don-t-spread-virus-all
[4] https://www.acc.org/latest-in-cardiology/journal-scans/2020/05/11/15/18/the-incubation-period-
of-coronavirus-disease
[5] https://www.eluniversal.com.mx/ciencia-y-salud/salud/coronavirus-cuanto-tiempo-lleva-
recuperarse-de-la-covid-19
[6] https://www.thelancet.com/journals/lancet/article/PIIS0140-6736(20)30185-9/fulltext#tbl1
[7] https://academic.oup.com/jtm/article/27/2/taaa021/5735319
ELMAR WAGNER
Instituto de Física y Matemáticas
Universidad Michoacana de San Nicolás de Hidalgo
Report