Introducción
La enfermedad por coronavirus del 2019 conocida como COVID-19 es una enfermedad infecciosa causada por el virus coronavirus de tipo 2 causantes del síndrome respiratorio agudo severo(SARS-CoV-2) que se caracteriza por causar diversas complicaciones en la salud de los humanos. La partícula viral está compuesta por una envoltura lipídica, la nucleocápside y las proteínas de membrana (M), envoltura (E) y la espiga (S), siendo esta última estructura la que otorga mayor variabilidad genética al virus.1
La aparición del virus en la población humana pudo darse probablemente por dos escenarios. El primero corresponde a un evento de selección natural en un hospedero animal antes de la transmisión zoonótica y el segundo corresponde a una selección natural en humanos luego de la transferencia zoonótica.2
La Organización Mundial de la Salud declaró el 11 de marzo de 2020 al COVID-19 como enfermedad epidémica y a nivel mundial hasta el 08 de febrero del 2021 se han confirmado tres variantes de preocupación de SARS-CoV-2 produciendo un total de 105.658.476 de casos acumulados llegando a producir una cantidad de 2.309.370 muertes.3
El primer caso de SARS-CoV-2 en Paraguay se reportó el 07 de marzo del 2020 y, luego de dos meses de aislamiento, las autoridades de salud pública aplicaron una estrategia de cuarentena inteligente, consistente en 4 fases que tenían varios niveles de restricciones para la realización de actividades. La última fase de la cuarentena inteligente se decretó el 18 de julio del 2020,4donde si bien se levantó el aislamiento en casi todo el país, hubo excepción en la Capital y los departamentos de Alto Paraná y Central debido principalmente a la elevada cantidad de casos reportados en estas localidades. En Paraguay se reportaron hasta el 01 de febrero del 2021 una cantidad de 2734 fallecidos y 21820 casos activos.
Los modelos epidemiológicos como el de susceptibles, infectados y recuperados (SIR) de Kermack y McKendrickse aplicaron a nivel mundial para predecir futuros escenarios de la propagación del virus,5 para recrear dicho modelo se pueden usar diversas herramientas computacionales entre ellas Vensim.
El software Vensim permite realizar diagramas de los componentes del sistema a estudiar, las relaciones entre las variables e introducir ecuaciones o valores numéricos a las distintas variables ingresadas, este software ha sido utilizado en diferentes estudios de epidemiología,6 ayudando a la toma de decisiones, es por ello que se establece como objetivo utilizar una modificación del modelo básico de epidemiología(SIR) y simular la dinámica del COVID en febrero del 2021 en Paraguay mediante el software Vensim, evaluar la cantidad de contagios activos, muertes y recuperados a nivel país, y mostrar si lo obtenido por simulación tiene correspondencia con los datos reales presentados en los informes del Ministerio de Salud del Paraguay. 7
Métodos
Se realizó una investigación observacional de tipo correlacional, en todos los pacientes que se enfermaron de Covid-19 durante el mes de febrero del 2021, en Paraguay. Primeramente se realizaron las simulaciones con el software Vensim y luego se compararon con los datos publicados por el Ministerio de Salud de Paraguay.
Modelo matemático.
El modelo utilizado corresponde a una variación del modelo SIR.
Donde β es la tasa de transmisión, γ es la inversa del periodo infeccioso, N es el tamaño de población, I son los infectados, S son los susceptibles, y Tm es la tasa de muerte. Las tres primeras ecuaciones corresponden al modelo SIR clásico y la cuarta ecuación es para obtener la cantidad de individuos muertos.
Se tienen además en cuenta las siguientes ecuaciones Re = R0x (S/N) donde Re es el número de reproducción efectiva y R0el número reproductivo básico cuya ecuación es R0= β/γ.
Herramienta computacional.
La simulación se realizó a través del software Vensim que pudo ser descargado de manera gratuita a través del sitio web del programa.
Las ecuaciones utilizadas con Vensim son las siguientes.
Individuos sanos= INTEG(-Nuevos contagiados). Valor Inicial= Población total.
Tasa de contagio= R0* γ
Nuevos contagiados= Individuos contagiados*Tasa de contagio
Individuos contagiados= INTEG(Nuevos contagiados-Nuevos muertos-Nuevos recuperados). Valor Inicial=21820
Nuevos recuperados= Individuos contagiados*Tasa de recuperados
Individuos recuperados= INTEG(Nuevos recuperados). Valor Inicial=109227
Tasa de recuperados= γ
Individuos muertos= INTEG(Nuevos muertos). Valor Inicial=2734
Nuevos muertos= Tasa de muertes*Individuos contagiados
Tasa de muertes= Muertos/Población total.
Resultados
Se realizó un diagrama causal del tipo Forrester, teniendo en cuenta variables que permitan la interacción de contagiados, muertos y recuperados como se indica en la Figura 1. Esto permitió observar de manera gráfica la relación entre los flujos y las variables.
La dinámica de sistemas simulada a través de Vensim permitió obtener valores para los distintos tiempos de simulación. Luego de 27 días se obtiene que los individuos contagiados considerados como activos son de 23616 casos como lo observado en la Figura 2.
Al realizar la simulación para la cantidad de decesos se obtiene un total de 2965 muertes como se visualiza en la Figura 3.
Así mismo se pudo realizar una estimación de la cantidad de recuperados, obteniendo una cantidad de 170444 recuperados, como se visualiza en la Figura 4.
Para evaluar la precisión de la simulación de manera visual, se graficaron los desvíos como el valor absoluto de la diferencia entre los valores simulados y valores reales de cada categoría, contra el transcurso del tiempo como se observa en la Figura 5. Debido a la existencia de rangos de órdenes de magnitud distintos, estos desvíos fueron normalizados, dividiendo con cada media respectiva.
Discusión
La complejidad de los distintos modelos epidemiológicos puede variar y de esta manera ofrecer resultados más cercanos a la realidad conforme a las variables que se lleguen a utilizar. En este caso no se utilizaron algunas variables que se tienen en cuenta en otros modelos como las variantes del virus SARS-CoV-2 o el número de personas hospitalizadas,8 sin embargo, la simulación realizada posee aproximaciones importantes respecto a la cantidad de casos activos reportados y los obtenidos a través del modelo. Se obtuvo que solo existe una diferencia de 758 casos entre los casos activos obtenidos por la simulación y los casos reportados del mes acertando en un 96,79%. Esta cantidad de casos activos también se ve reflejada con el número de reproducción efectiva reportado por el ministerio de salud del Paraguay que era de 1,03 para el 01 de febrero y fue el dato inicial con el que se inició el modelado. El número de reproducción efectiva indica el tipo de propagación que existe, donde valores superiores a 1 indican que la enfermedad seguirá propagándose,9 en este caso el número fue ligeramente superior a 1 por lo cual hubo un aumento en cuanto a los casos reportados inicialmente.
Los datos de tasa de transmisión utilizados para el modelo fueron obtenidos a partir del número reproductivo básico, que es igual a la razón entre la tasa de contagio y la tasa de recuperación. La tasa de recuperación establecida fue de 0,1 que corresponde a la inversa del periodo infeccioso, dicho valor fue escogido debido a que las personas confirmadas con COVID-19 en Paraguay deben presentar aislamiento por 10 días a partir del primer día en el que se tuvieron los síntomas. Esto coincide con otros estudios que refieren que los pacientes con síntomas leves pueden llevar 10 días de recuperación y que los que poseen una enfermedad crítica o inmunodeficiencia podrían tardar 15 días.10
Para los recuperados existe una mayor brecha entre lo simulado que corresponde a 170.444 y lo que se reportó 133.435 teniendo de esta manera una diferencia de 37.009 casos. Esto se debería a que cada persona en particular tiene una evolución distinta durante su recuperación. Existen individuos que por presentar enfermedades de base son consideradas de riesgo y son más propensos a tener una recuperación tardía incluso superior a 15 días. En Paraguay existe una prevalencia de 46,8% personas hipertensas.11Otras enfermedades de consideración que no se tienen en cuenta durante el modelado son la cantidad de personas con diabetes y cardiopatías donde se observa que estos pacientes necesitan de mayores cuidados.12
El modelo sencillo de SIR está representado por los individuos susceptibles, infectados y recuperados y no tiene en cuenta la cantidad de muertes, 13 sin embargo, para realizar la simulación a través de Vensim se integró este flujo de manera a obtener la cantidad estimada de muertes durante el periodo de simulación. Las complicaciones inherentes de padecer COVID-19 con cuadros de base podría aumentar la tasa de fallecidos, la diferencia entre lo esperado y lo reportado fue de 216. Según los informes desde el 01 al 28 de febrero de 2021 se produjeron 447 muertes, de esta manera la simulación acertó en un 51,67% sobre la cifra total de decesos ya que la cantidad final de la simulación fue de 2965 muertos. Algunos aspectos a tener en cuenta y que influyen en el número de pacientes muertos es que no todos poseen el fácil acceso a medicamentos por su costo y por la disponibilidad en el mercado y los hospitales, así mismo existen ocasiones en las que no hay un rápido ingreso a unidades de terapia intensiva debido a que los mismos se encuentran llenos o viven lejos de los hospitales con dicho servicio. Cabe destacar que también existen muchas personas contagiadas que no acuden a los centros de salud y consumen medicamentos sin prescripción médica lo cual desemboca en la complicación del cuadro clínico.
El modelo propuesto y simulado a través de Vensim podría ayudar a la correcta toma de decisiones en cuanto a la estimación de individuos contagiados. Esto se visualizó mediante los valores obtenidos de los desvíos normalizados de las distintas categorías utilizadas en el modelo, sin embargo las aproximaciones respecto al escenario de personas recuperadas y la cantidad de decesos deberían ser mejoradas incluyendo variables como capacidades de los hospitales para pacientes con COVID-19, 14 se debe tener en cuenta las afecciones de base que presenten las personas contagiadas y también a aquellos pacientes que presentan la enfermedad y no quieren recibir tratamientos médicos.15
Conclusiones
La simulación de la dinámica de sistemas es compleja y depende en gran medida de los datos de partida, en tal sentido si bien la simulación tuvo buena aproximación respecto a los casos activos hubo menor precisión en los valores de recuperados y muertes lo cual fue observado por la mayor diferencia entre lo esperado y lo reportado. Esta diferencia plantea la necesidad de realizar ajustes de los parámetros del modelo para futuras simulaciones. El tiempo de simulación se aplicó para un intervalo de 27 días sin embargo este tipo de simulaciones se puede realizar en otros periodos más amplios, no obstante, al tener un mayor intervalo de tiempo los resultados obtenidos podrían tener mayores variaciones debido a la aparición de nuevas variables que no se tienen en cuenta al principio del modelado como las características individuales relacionados a la salud de las personas, hábitos sociales y la toma de decisiones a nivel país.