knitr::opts_chunk$set(echo = TRUE)
library("xlsx")
library("readxl")

Se analiza la evolución del coronavirus desde el 27 de abril de 2020, cuando el Gobierno cambió el criterio de contabilización de hospitalizados y nuevos casos. Antes de esa fecha, el Gobierno contabilizaba los nuevos positivos tanto por anticuerpos cómo por PCR, y cómo hospitalizados sólo los positivos por PCR. El 26 de abril se incorporaron cerca de 24.000 hospitalizados a las estadísticas y se empezó a contar cómo positivos sólo los verificados por PCR.

En este estudio contaremos cómo casos tanto los confirmados por anticuerpos cómo por PCR, aprovechando que, aunque el Gobierno no computa los primeros, el Instituto de Salud Carlos III facilita ambos datos.

Selecciono mi archivo excel, confeccionado con datos del Instituto de Salud Carlos III, los datos que utilizaremos en el análisis, y lo guardo en la variable “coronadata”.

knitr::opts_chunk$set(echo = TRUE)
coronadata<-read_excel("C:/Users/JoséIgnacio/AppData/Roaming/Microsoft/Windows/Network Shortcuts/coronavirus 12-05-pcr+ag datos claros.xlsx")

Hacemos gráfico del ratio de Total hospitalizados/total detectados, siendo ambas magnitudes datos acumulados desde el inicio de la pandemia en España, donde se puede apreciar una clara tendencia descendente.

plot(coronadata$Fecha,coronadata$`hosp/casos acum.`,
      type='s',
      pch=19,
     col="blue",
      main = "Evolución Total hosp/total detectados",
     xlab = "Fecha")

Hacemos gráfico de los casos activos, dónde se puede observar también una tendencia descendente.

knitr::opts_chunk$set(echo = TRUE)
plot(coronadata$Fecha,coronadata$`casos act.`,
     type='s',
      pch=19,
    col="orange",
      main = "Evolución Casos Activos",
      xlab = "Fecha",
      ylab = "Casos Activos")

Hacemos gráfico del ratio casos nuevos hospitalizados/casos nuevos totales diarios.

knitr::opts_chunk$set(echo = TRUE)
plot(coronadata$Fecha,coronadata$`Hosp/casos (nuevos)`,
     type='b',
      pch=19,
    col="blue",
      main = "Evolución Hospitalizados/Casos(nuevos)",
      xlab = "Fecha",
      ylab = "Hospitalizados/Casos(nuevos)")
      abline(h=mean(coronadata$`Hosp/casos (nuevos)`), col="red")

Podemos ver que la media (línea horizontal) es ligeramente mayor del 30%. Se ha estimado que el 20% de los infectados requieren hospitalización. Atribuímos este mayor porcentaje a que se realizan test mayoritariamente a personas con síntomas, cómo efectivamente sabemos que ocurre.

También se sabe que aproximadamente el 25% de los ingresados requieren UCI. Veamos cual es el ratio en España en estas últimas semanas.

knitr::opts_chunk$set(echo = TRUE)
plot(coronadata$Fecha,coronadata$`UCI/Hosp (nuevos)`,
                 type='b',
                 pch=19,
                  col="orange",
                  main = "Evolución UCI/Hospitalizados(nuevos)",
                 xlab = "Fecha",
                ylab = "UCI/Hospitalizados(nuevos)")
 abline(h=mean(coronadata$`UCI/Hosp (nuevos)`, col="red"))

Podemos ver que la media (línea horizontal) está ligeramente encima del 10%, lo que indica que la gravedad de la enfermedad es ahora bastante menor.

Concluimos del estudio que el coronavirus está disminuyendo tanto en incidencia (casos activos), como en gravedad, por lo que es posible que llegue a desaparecer o se convierta en una enfermedad con muy poca mortalidad sin necesidad de vacuna.