Modelo epidemiológico Covid-19



Recientemente, ha habido modelos bastante diferentes para el desarrollo de la epidemia, incluso en Habré. Este tema tampoco me ha pasado por alto. Apenas habría escrito aquí, pero dado lo que logré encontrar, la importancia de las dependencias descubiertas y su impacto en nuestras vidas, no puedo dejar de compartir el hallazgo. Habrá muchas fórmulas, gráficos y poco texto. Información básica y gráficos para Alemania, donde vivo.

Entonces, el modelo epidemiológico en la primera aproximación es descrito por la fórmula de crecimiento de los infectados.

Minuto de atención ovni


La pandemia COVID-19, una infección respiratoria aguda potencialmente grave causada por el coronavirus SARS-CoV-2 (2019-nCoV), se ha anunciado oficialmente en el mundo. Hay mucha información sobre Habré sobre este tema; recuerde siempre que puede ser confiable / útil, y viceversa.

Le instamos a que sea crítico con cualquier información publicada.


Fuentes oficiales

Si no vive en Rusia, consulte sitios similares en su país.

Lávese las manos, cuide a sus seres queridos, quédese en casa siempre que sea posible y trabaje de forma remota.

Leer publicaciones sobre: coronavirus | trabajo remoto

Nt=N02t/Td


Dónde Td- el tiempo de duplicación de infectados en nuestro caso en días,t número de díasN0 número de infectados en un momento determinado yNt - el número de casos hastat días Si dividimos ambas partes de la fórmula por la población total de la región, obtenemos la misma fórmula, pero en partes de la poblaciónP .

Pt=P02t/Td(1)


El problema con esta fórmula es que la fórmula no tiene en cuenta la población limitada y Ppronto será más de 1. Esto no sucede en la vida real.

Existe un factor epidemiológico que determina hasta qué nivel puede aumentar el número de casos. Se calcula en función delnúmero base de reproducciones. . Este número nos muestra cuántas personas infecta aproximadamente una persona infectada, es constante y específico para cada región en particular, dependiendo de la densidad de población y otras características de la región. Se puede determinar solo al comienzo de la epidemia, cuando no hay factores limitantes. La fórmula en sí se ve así:R0

Psat=11/R0


También hay un número efectivo de reproducciones. , que también nos permite saber cuántas personas infecta el paciente. A diferencia del número base, el efectivo está cambiando constantemente. Puede determinar este valor utilizando la fórmula anterior y conocer el número de infectados en un determinado momento:Rt

Rt=R0(1Pt)(2)


Si tomamos el modelo simplificado SEIR [1] de la epidemia, podemos encontrar factores adicionales que describen las características de la epidemia, como la tasa de crecimientoGr o tiempo de infectividad del pacienteD . Las siguientes fórmulas muestran la relación entre las cantidades.

Gr=Rt1D


Td=ln(2)Gr


Usando las fórmulas anteriores, podemos derivar la siguiente dependencia

Td=Dln(2)Rt1


y sustituyéndolo en (1) obtenemos:

Pt=P02t(Rt1)Dln(2)


o después de la simplificación

Pt=P0et(Rt1)D


Si necesitamos determinar el valor al día siguiente, entonces t=1

P1=P0eRt1D


El número efectivo de reproducciones para un tiempo específico currentpuede calcularse a partir de la fórmula (2) y luego conocer solo el número básico de reproduccionesR0 y el número de infectados en este momentoPcurrent podemos calcular fácilmente el porcentaje de personas infectadas al día siguiente.

Pnext=PcurrenteR0(1Pcurrent)1D


Solo hay un parámetro en esta fórmula , que puede calcularse a partir del tiempo de duplicación al comienzo de la epidemia. Tomando por ejemploR0Pcurrent=0,0001y tomando n pasos, obtendremos un estado epidemiológico en n días. El tiempo, la forma de la curva, el valor de saturación, el número de casos en un momento determinado y otros parámetros saltan de la fórmula "como un demonio de una caja de rapé".

¿Qué pasa con la cuarentena y otros factores que afectan el curso de la epidemia?

Cada una de las medidas tomadas corrige el número base de reproducciones por un determinado factor (factor)μ[0;1] de la siguiente manera:

Pnext=PcurrenteμR0(1Pcurrent)1D


Para simplificar, incluso puede definir un número de reproducción base "limitante":

R0=μR0


Además, en algunos momentos, simplemente puede reemplazar un número de reproducción con otro, ajustando así la propagación de la epidemia. La epidemia continúa propagándose con nuevas condiciones durante un cierto período de tiempo, hasta el momento de un nuevo cambio. Los puntos de cambio o los puntos de intervención están determinados por factores externos, como la cuarentena, el cierre de escuelas o la necesidad de usar máscaras. El tiempo y el alcance de la exposición a estos factores, como regla, no se pueden conocer de antemano. Sin embargo, si se conoce el valor del cambio del número de correcciónμ, que determina la efectividad de la cuarentena, puede averiguar cómo afectará su cancelación en un momento determinado en el futuro. Esto proporciona buenas oportunidades predictivas para este modelo. Esta es también una forma de probar su validez en la práctica.

Como el tiempo de infectividad del paciente.Dpara Covid-19, se tomó un valor de 10 [2] .

No hay otros parámetros en el modelo, así como grados adicionales de libertad.

¿Qué pasa con la verificación?

Gráficos basados ​​en datos de Alemania.

Solo se indicaron 3 puntos de intervención en la siguiente tabla: lo



que condujo a los siguientes resultados.

imagen



Los puntos de intervención y comparación de los datos del modelo con los valores reales extraídos de los datos públicos son visibles en el gráfico de cambios en el número efectivo de reproducciones.



La coincidencia de los datos y la calidad del modelo se pueden verificar en los cuadros de regresión:



el modelo y los cálculos para Alemania se publican en GitHub . No solo existen estos datos, sino también estudios sobre muertes.

Actualizar:
Verificaciones adicionales realizadas. Cada país tiene su propio factor de corrección.

Rusia:




Italia:



EE . UU .:



Actualización 2: se
agregó una versión de "simple" en GitHub que elimina todos los innecesarios, puede insertar valores de otros países, cambiar puntos de intervención y factores de corrección. Existe una alta probabilidad de que este mismo factor de corrección sea la proporción de personas diagnosticadas a infectadas. Pero necesita ser revisado. Un mayor desarrollo y finalización de la epidemia confirmará o refutará esta hipótesis.

En los gráficos, los valores de los casos detectados están en% y estos valores no se corrigen para el valor de Ratio (Factor de corrección). Dividiendo por este número obtenemos el porcentaje real de casos detectados y el pronóstico de la infección. En la versión simple, esta corrección se ha realizado.

Referencias
[1] JM Heffernan et al. Perspectivas sobre la relación reproductiva básica. doi.org/10.1098/rsif.2005.0042 JR Soc. Interface 2005 2, 281–293 (2005)

[2] Xi He, Eric HY Lau y cols. Dinámica temporal en la eliminación viral y la transmisibilidad de COVID-19. www.nature.com/articles/s41591-020-0869-5 Naturaleza (2020)

All Articles