INTRODUCCIÓN
La Organización Mundial de la Salud declara el 11 de marzo del 2020 la pandemia de COVID-19, instando a la imposición de medidas por parte de los gobiernos de cada país. Es de vital importancia contar con información sobre los casos y poder hacer un análisis estadístico que permita apoyar el estudio de la pandemia para poder aplicar las medidas necesarias con los recursos de que se dispone.(1)
El nuevo coronavirus Sars-CoV-2, causante de la enfermedad COVID-19, presenta una rápida propagación y elevada mortalidad, comparado con otros coronavirus.(2) Cuba adoptó medidas, con el fin de contrarrestar su propagación en la población. Sin embargo, como ha sucedido en otros países afectados por esta enfermedad, se produjo un crecimiento exponencial de casos confirmados.
Existen investigaciones científicas que antecedieron este estudio, realizadas al inicio de la pandemia en Cuba por autores como Medina Mendieta, Cortés Cortés, Cortés Iglesias,(3,4,5) donde exponen siete modelos de crecimiento poblacionales que han sido utilizados por diferentes autores en diferentes países para el pronóstico a corto plazo de picos en la pandemia provocada por la COVID-19.(6,7)
En el mundo se presenta una crisis de salud asociada al coronavirus SARS.CoV.2(3) y los gobiernos necesitan investigar sobre el comportamiento de la enfermedad y para ello han iniciado estudios de modelos matemáticos. De vital importancia es conocer, aproximadamente, la cantidad máxima de casos contaminados. Por ello se hace necesario contar con estimaciones que puedan servir para la toma de decisiones con los recursos disponibles.
Para pronosticar esta enfermedad se utilizan las modelaciones matemáticas. Se destacan: la modelación clásica con los modelos SIR (Susceptibles – Infestados – Recuperados), basada en ecuaciones diferenciales ordinarias, utilizada con éxito en algunas regiones.(7) En Cuba, varios autores han aplicado esta modelación a la epidemia provocada por la COVID-19.(3,4,5)También se han aplicado los modelos estadísticos de series cronológicas,(8) los modelos predictivos,(9) los modelos basados en inteligencia artificial y Machine Learning.(10)
Con el empleo de los métodos de la estadística se ofrece información relevante sobre el comportamiento de la pandemia y se identifican los indicadores más importantes que le permitan a los decisores de salud desarrollar las estrategias pertinentes en el momento oportuno.
Esta investigación permite demostrar el vínculo necesario e imprescindible que debe existir entre las ciencias y la relación intersectorial indispensable para el logro exitoso de los planes de desarrollo territoriales en los momentos matizados por la pandemia.
Los resultados de esta investigación forman parte de un proyecto nacional del grupo de investigadores del Departamento de Matemática de la Facultad de Economía de la Universidad de Cienfuegos avalado por el grupo nacional de trabajo de la Facultad de Matemática y Computación de la Universidad de la Habana.
Con el objetivo de realizar estudios estadísticos sobre los datos de casos confirmados en la provincia de Cienfuegos, se realizó este estudio.
MÉTODOS
Se realizó un estudio observacional descriptivo, transversal con un análisis estadístico matemático de los casos confirmados de COVID -19 en los 8 municipios de la provincia de Cienfuegos, desde marzo del 2020 hasta agosto del 2021.
Los estudios estadísticos descriptivos abarcan las variables: número de casos confirmados, acumulados, la edad, el sexo, las dosis recibidas de la vacuna cubana Abdala y las fechas en que se presentan los confirmados. Se obtienen ecuaciones de ajuste de Gompertz, Weibull y Loglogistic. Se calcularon los números de reproducción básica Ro y efectivo Rt.
El siguiente informe se adecuó al tiempo de estudio para la pertinencia de los modelos y se valoró en cada municipio de la provincia de Cienfuegos, el que mejor se comportaba. Se describen los modelos usados en la investigación:
Ecuación 1. Modelo logístico de cuatro-parámetros.(11,12)
Ecuación 2. Modelo logístico de cinco-parámetros.(11,13)
Ecuación 3. Modelo de medias para la curva de crecimiento de Gompertz de dosis-respuesta.(11,14,15)
Ecuación 4. Modelo Weibull de cuatro-parámetros.(11,16)
Ecuación 5. Los modelos de hormesis de Brain-Cousens.(17)
Ecuación 6. Modelo de dos fases de dosis-respuesta.(18)
Ecuación 7. Modelo de dosis-respuesta de cinco parámetros baroreflex.(19)
Número de reproducción básico (Ro) y efectivo (Rt)
Se trabajó en la estimación del número de reproducción básico (Ro), aplicando el método de máxima verosimilitud,(16) y el número reproductivo efectivo (Rt), con el método de Anne Cori,(6,20) se utilizan los paquetes \Ro" y \EpiEstim" del lenguaje de programación R.
Aplicación
(Ro) es una medida epidemiológica fundamental en tiempos de pandemia. Permite medir la tasa de contagio o de transmisión de la enfermedad en el territorio. Ro > 1 indica un incremento de casos y por ende una expansión de la pandemia. Ro < 1 es un indicio de reducción del número de casos. (Rt) mide la tasa de transmisión en cada momento (t). El conocimiento de estas medidas, permite a las autoridades de salud y del gobierno diseñar estrategias para reducir Ro y su seguimiento permite medir la efectividad de las medidas tomadas.
La información se obtuvo del registro estadístico de pacientes confirmados de COVID-19, proporcionada por el departamento de epidemiología de la Dirección Provincial de Salud de Cienfuegos.
Con los datos obtenidos y la relación que se presenta se estima la tasa de reproducción R de la epidemia, con información diaria disponible se obtiene de:(20,21)
donde:
t el día actual a días transcurridos desde el diagnóstico.
u y v los límites del periodo en que un caso es contagioso.
C(t) número de casos diagnosticados.
W(a) distribución Gamma es un ponderador de los casos ocurridos a días atrás.
El número 6 significa el rezago, 4 días para que aparezcan los síntomas y 2 para que se comuniquen.
La investigación realizada presenta limitaciones en su desarrollo como son: la dinámica de la búsqueda de la información diaria de la epidemia, la creación de la base de datos, la aplicación de estadística descriptiva, por un paquete de programa. Además, la obtención de los mejores modelos matemáticos de ajustes, el cálculo de los números básicos y reproductivos y la variedad de paquetes de programas.
La investigación fue aprobada por el consejo científico de la Universidad de Cienfuegos.
RESULTADOS
La distribución de los pacientes confirmados con COVID-19 muestra que más del 50 % pertenecía al municipio de Cienfuegos (55,43 %). El promedio de contactos por confirmados osciló entre 8 y 9 pacientes. La edad entre 20 y 59 años fue la más representada, mientras que las edades extremas de la vida (menor de un año y más de 80) mostraron el menor porciento de confirmados. La vacunación preventiva mostró los primeros resultados, aún con coberturas bajas (en segunda y tercera dosis), teniendo en cuenta el período estudiado. (Tabla 1).
Ecuaciones de pronóstico por municipios en el período analizado.
Municipio Abreus:
Una vez corridos los modelos estadísticos en estudio, en el período estudiado, se observó que el que mejor se adecuó a los datos de los casos confirmados reportados acumulados en el municipio de Abreus fue el modelo de Weibull de cuatro parámetros. La curva de contagiados confirmados puede verse creciente en el período analizado y se pronosticó que el pico máximo de contagios fuera de 524.54 y se alcanzara el 17 de agosto 2021. El pico real se obtuvo el día pronosticado. (Gráfico 1).
Municipio Aguada de Pasajeros:
El modelo estadístico que mejor se adecuó a los datos, en el período estudiado, de los casos confirmados reportados acumulados en el municipio de Aguada de Pasajeros fue el modelo de Gompertz de cuatro parámetros. La curva de contagiados confirmados puede verse creciente en el período analizado y se pronosticó que el pico máximo de contagios fuera de 519.89 y se alcanzaría el 12 de agosto 2021. El pico real se obtuvo unos días posteriores. (Gráfico 2).
Municipio Cienfuegos:
El modelo estadístico que mejor se adecuó a los datos, en el período estudiado, de los casos confirmados reportados acumulados en el municipio de Cienfuegos fue el modelo de Weibull de cuatro parámetros. La curva de contagiados confirmados puede verse creciente en el período analizado y se pronosticó que el pico máximo de contagios fuera de 522.82 y se alcanzaría el 15 de agosto 2021. El pico real fue obtenido un día después del pronosticado. (Gráfico 3).
Municipio Cumanayagua:
El modelo que mejor se adecuó a los datos, en el período estudiado, de los casos confirmados reportados acumulados en el municipio de Cumanayagua fue el modelo de Weibull de cuatro parámetros. La curva de contagiados confirmados puede verse creciente en el período analizado y se pronosticó que el pico máximo de contagios fuera de 524.79 y se alcanzaría el 17 de agosto 2021. El pico real se obtuvo el mismo día pronosticado. (Gráfico 4).
Municipio Cruces:
El modelo que mejor se adecuó a los datos, en el período estudiado, de los casos confirmados reportados acumulados en el municipio de Cruces fue el modelo logístico de cinco parámetros. La curva de contagiados confirmados puede verse creciente en el período analizado y se pronosticó que el pico máximo de contagios fuera de 523.11 y se alcanzara el 16 de agosto 2021. El pico real se obtuvo unos días después del pronosticado. (Gráfico 5).
Municipio Lajas:
El modelo que mejor se adecuó a los datos de los casos, en el período estudiado, de los casos confirmados reportados acumulados en el municipio de Lajas fue el modelo logístico de cinco parámetros. La curva de contagiados confirmados puede verse creciente en el período analizado y se pronosticó que el pico máximo de contagios fuera de 524.42 y se alcanzara el 17 de agosto 2021. El pico real se obtuvo unos días posteriores. (Gráfico 6).
Municipio Palmira:
El modelo que mejor se adecuó a los datos, en el período estudiado, de los casos confirmados reportados acumulados en el municipio de Palmira fue el modelo Weibull de cuatro parámetros. (Gráfico 7).
Municipio Rodas:
El modelo que mejor se adecuó a los datos, en el período estudiado, de los casos confirmados reportados acumulados en el municipio de Rodas fue el modelo logístico de cinco parámetros. La curva de contagiados confirmados puede verse creciente en el período analizado y se pronosticó que el pico máximo de contagios fuera de 513.26 y se alcanzara el 6 de agosto 2021. El pico real se obtuvo unos días anteriores al pronosticado. (Gráfico 8).
Con el cálculo del número básico Ro se pasó a obtener el número reproductivo Rt por municipios, este número representa la tasa de trasmisión en función de su comparación con la unidad. Un Rt < 1 implica reducción del número de casos, representados en tonos de verde y un Rt > 1 representa un incremento de casos, representados en tonos de rojo. En la figura 1 se observa la distribución de la trasmisión en cada uno de los municipios de la provincia. (Fig. 1).
DISCUSIÓN
En el presente estudio se analizaron siete modelos predictivos, Gompertz, Weibull y Loglogistic, de diferentes parámetros, para la base de datos de los confirmados en el período de inicio de marzo del 2020 hasta julio del 2021, en los municipios de la provincia de Cienfuegos, Cuba. Se calculó la ecuación del mejor modelo predictivo para cada municipio de la provincia de Cienfuegos, según las pruebas de significación realizadas.
De los modelos predictivos encontrados para cada municipio se calculó el valor y la fecha del pico de la pandemia. La fecha pico pronosticada en los municipios de Abreus, Cumanayagua, Palmira y Cienfuegos coincidieron totalmente con la fecha real; en el resto de los municipios el pronóstico tuvo pocos días de diferencia con la realidad.
Los modelos predictivos poblacionales logístico, de Weibull y de Gompertz utilizados para la obtención de ecuaciones pronósticos en los municipios de la provincia de Cienfuegos de casos positivos de la COVID-19, han sido validados estadísticamente. Investigaciones similares, realizadas y publicadas por los autores han aplicado estas técnicas para obtener pronósticos sobre la epidemia provocada por la COVID-19 en otros países del mundo, con resultados positivos.(3,4,5)
Se calcularon los números básicos Ro y reproductivo Rt por municipios, para medir la tasa de reproducción y se obtuvo el mapa de la distribución de la trasmisión por colores de la provincia de Cienfuegos.
La investigación contribuyó a la toma de decisiones en la provincia de Cienfuegos, para proyectar, controlar y evaluar el comportamiento de la COVID-19 en sus municipios.
Se obtuvieron las ecuaciones predictivas de cada municipio de Cienfuegos en el período estudiado utilizando los modelos de Weibull, Gompertz y Loglogistic con varios parámetros, con resultados estadísticos significativos.
Para cada modelo municipal se obtuvo el valor del pico de la pandemia y la fecha en que se obtuvieron, con resultados adecuados.
El mapa de distribución de la trasmisión mediante el cálculo del número reproductivo Rt constituye una valiosa herramienta gráfica para la toma de decisiones en la provincia de Cienfuegos.
Conflicto de intereses
Los autores de este trabajo no presentan conflicto de intereses.
Contribuciones de los autores
Conceptualización: Manuel Cortés Cortés.
Metodología: Manuel Cortés Iglesias.
Investigación: Felipe Medina Mendieta.
Análisis formal: Felipe Medina Mendieta, Miguel Santana Justiz.
Supervisión: Ridelio Miranda Pérez.
Redacción: Manuel Cortés Cortés,
Financiación
Universidad Carlos Rafael Rodríguez. Cienfuegos. Cuba.