Buenas! Hoy traigo algo complicado de entender, pero creo que muy útil. Lean con calma, vale la pena, y esperen al final para comentar, que es denso. Les voy a explicar cómo calcular aproximadamente qué día acabará el coronavirus... si todo va bien. Dentro hilo!
Como recordarán si hicieron mates "serias" en la universidad, los sistemas complejos se rigen por ecuaciones diferenciales. Se llaman así porque, en lugar de describir relaciones estáticas, describen relaciones cambiantes en base a su diferencia momento a momento.
Por ejemplo, si yo digo y=3x, eso NO es diferencial, es una verdad invariante. En cambio, si yo digo que el número de enfermos x crece según los que había antes, multiplicado por una constante (el contagio, R) por cada unidad de tiempo, me sale: R*x=dx/dt. O sea, cada enfermo x..
...multiplicado por su capacidad de contagiar R es igual al DIFERENCIAL (incremento) de enfermos (dx) por unidad de tiempo (dt). Lo sé, un lío. Háganme un favor, paren. Léanlo otra vez. Si lo leen con calma y se hacen un dibujo, se entiende, y es la base de lo que voy a explicar.
Las ecuaciones diferenciales complican todo: resolverlas no es tan fácil como aislar la X y ver su valor: hay que hacer un cálculo más avanzado (una integral), y ver a dónde nos lleva eso. Como son complejas, hay muchos "modelos" de ecuaciones diferenciales que están estudiadas.
Bien, si han entendido esto, el resto del hilo es fácil. Resulta que para epidemias está estudiado qué ecuación toca. Uno de los modelos de ecuación diferencial que se usa es la SIR (Suspected-Infected-Recovered). Son 3 curvas: una de sospechosos (o sea, sanos). Una de...
... infectados, y una de recuperados. La forma típica es así: 
Como pueden ver, la curva AZUL (sospechosos) empieza con que todos estamos sanos, y va bajando a medida que sube la VERDE (enfermos). La gente pasa de estar sana a infectada. Y, con un cierto retraso, aparece la curva ROJA (curados), como resultado de que la gente pasa de Azul...
... a verde y luego se cura. Bien, la curva que les he puesto es genérica. Evidentemente, cada epidemia tiene la suya. De ahí dije yo "es una exponencial" y me dijeron que estaba loco. Como pueden ver, la curva verde efectivamente es una exponencial, que luego frena y decae.
Bien. Por qué digo esto? Pues porque unos señores majísimos de la Universidad de Tecnología de Singapur han cogido todos los datos del Coronavirus y los han metido en ese modelo. Es francamente complicado, no me voy a enrollar, pero dejémoslo en que hay métodos numéricos que...
...ajustan curvas a conjuntos de datos "ruidosos", de forma que puedes coger los datos de España, y ver qué curva SIR es la que más se les parece. Evidentemente, lo que simularás es la curva de INFECTADOS (la verde). De esa se deducen las otras dos. Me han seguido de momento?
No se quejarán, si quieren ciencia y datos, se van a hartar conmigo. Como hay quien dice "Dani no sabes de qué hablas", seguiré hablando ;-) Venga, esta es la web, esto es una mina de oro. Si yo fuese @joethebrew y @NewsReputation, estaría babeando:
ddi.sutd.edu.sg
ddi.sutd.edu.sg
Paréntesis: a mis seguidores matemáticos (esto va por @PolarBearby), en la web misma les explican cómo han hecho el fitting (=el ajuste de datos). Lo digo por si hay errores, no soy TAN ninja para verlos. Parece que saben de lo que hablan, lo tienen aquí: mathworks.com/matlabcentral/…
Vale, volvamos al tema que nos ocupa. España. Esta es su curva, que como pueden ver ajusta bastante bien al modelo, esto es lo deseable.
Esto como se interpreta? Bien, esto lo que nos dice es que, hacia el día 7 de Mayo, el 97% de nuevos casos en España ya se habrá producido. Y que para el 22 de Mayo, el 99%. A ver, no abran el champán todavía, que no hemos acabado. Les cuento:
Primero: las curvas son tan buenas con los datos en los que se basan. Si los datos son incompletos, por ejemplo porque se a contado de más, o de menos, o porque un día se contó mal, todo eso afecta al modelo. Vamos, que esto es APROXIMADO. Pero hay dos motivos más:
Recuerden que el SIR son 3 curvas: esto nos dice que casi no va a enfermar gente NUEVA. De aquí a la victoria, faltan dos cosas: a) que los que enfermen estos días cubran su ciclo de enfermedad, y se curen b) que no se produzca un rebote. Les voy a enseñar un rebote "de libro":
Comparen Italia con Canadá. Se las pongo de lado y se ve mejor. Con Italia es fácil ver la tencencia descendente, en Canadá aún se producen estallidos y es difícil ver la tendencia.
Ven como Italia se ajusta mejor a la curva, pero Canadá tiene estallidos más o menos continuos? Eso es porque la epidemia aún no retrocede, y se producen rebrotes continuamente. De lo que hagamos los próximos días dependerá que seamos Italia, y pronto habremos vencido, o Canadá.
Matiz: evidentemente el coronavirus ha sido peor en Italia que en Canadá, no nos equivoquemos. Lo que la gráfica anterior explica es que, dentro de su TAMAÑO DE EPIDEMIA (mayor en Italia que en Canadá), Italia lo tiene controlado, y Canadá no tanto. Pero volvamos a España.
Esta gráfica es bonita. Hermosa. Esta gráfica huele a victoria si no la cagamos. Miren qué hermosura, qué campanita más mona:
En resumidas cuentas: esto sigue pintando BIEN. Aún hay infectados nuevos, pero pocos. Los que se infecten evidentemente seguirán su curso, con lo cual aún habrán muertes, tristemente, pero menos. Y a esto le quedan semanas antes de que la epidemia remita ya del todo.
Pero claro, eso es cierto bajo dos hipótesis: la primera, que estos señores de Singapur no la hayan cagado en su análisis. No lo parece, la verdad. La segunda, que NOSOTROS no la caguemos con nuestra parte. Porque recuerdan que he empezando diciendo "cada persona contagia a N"?
Pues la única manera de que no las contagie es... evidentemente, no contagiarlas. O sea no acercarse. No abrazarse. No salir fuera de horas. No saltarse las normas. No tocarlo todo y dejarlo lleno de virus. Joder, obedecer, que no es tan difícil señores.
Una reflexión final: si explico cosas de forma coloquial es para llegar a un público amplio, esa es la base de la divulgación científica. Que NADIE crea que me tomo esto a broma. 25.000 muertos es un drama. Es 25 veces los muertos en carretera en un año. Poca broma.
Ya llegará el momento de las críticas y las denuncias. Pero de momento centrémonos en vencer. Espero que les haya gustado el hilo, como ven, intento dar información de calidad que nos ayude a entender el mundo que nos toca vivir. Si les gustó, un RT hace milagros.
d
d
