[PYTHON] Lassen Sie uns die medizinische Kollapshypothese des neuen Koronavirus testen

Einführung

In Bezug auf die Bedrohung durch die neue Coronavirus-Infektion (COVID-19) wurde auch in Japan eine Notfallerklärung abgegeben, und die Gefahr eines medizinischen Zusammenbruchs ist gefordert. Insbesondere als Reaktion auf den raschen Anstieg der Zahl infizierter Menschen in städtischen Gebieten haben Maßnahmen wie die Aufnahme milder und asymptomatischer Menschen in Unterkünfte begonnen. Darüber hinaus wurde, wie im Fernsehen berichtet, darauf hingewiesen, dass es in China, Italien, den Vereinigten Staaten usw. zu einem medizinischen Zusammenbruch gekommen ist, und es gibt Fälle, in denen aufgrund einer Lähmung der Krankenhausfunktionen keine ausreichende Behandlung erhalten werden kann. In diesem Artikel möchte ich den medizinischen Kollaps anhand der von der Johns Hopkins University veröffentlichten Daten ** quantitativ verifizieren ** und gleichzeitig konkret definieren, was ein medizinischer Kollaps ist.

Was ist ein medizinischer Zusammenbruch?

Laut Wikipedia Es scheint, dass es bedeutet, dass "ein stabiles und kontinuierliches medizinisches Versorgungssystem aufgrund einer Abnahme der Arbeitsmoral der Ärzte, einer Zunahme der verteidigungsmedizinischen Versorgung, einer Verschlechterung des Krankenhausmanagements usw. nicht eingerichtet wird", sondern in Bezug auf COVID-19 Es gibt einen Mangel an medizinischem Personal und medizinischer Ausrüstung, und es ist schwierig geworden, schwerkranke Menschen zu behandeln. " Mit anderen Worten, es scheint, dass sich die Bedeutung je nach Kontext, in dem sie verwendet wird, geringfügig ändert. Im früheren Sinne [diese Grafik](https://ja.wikipedia.org/wiki/%E3%83%95%E3%82%A1%E3%82%A4%E3%83%AB:Oecd- Mit Blick auf doctorconsult2013.svg) scheinen sich Japan und Südkorea als Länder mit einer großen Anzahl von Konsultationen im Vergleich zur Anzahl der Ärzte pro Bevölkerung hervorzuheben. Dieses Mal möchte ich die Definition des medizinischen Kollapses * im Zusammenhang mit * COVID-19 wie folgt definieren.

Quantifizieren Sie den medizinischen Kollaps

Von der Johns Hopkins University veröffentlichte Daten aggregiert tägliche kumulative Infektionen, Genesung und Todesfälle aus der ganzen Welt. Ich bin. Basierend auf diesen Daten werden die folgenden zwei Indikatoren berechnet.

  1. Mortalität rD: Anzahl der Todesfälle an diesem Tag, die sich aus der Anzahl der Personen ergeben, die am Vortag behandelt werden müssen (C)
  2. Wiederherstellungsrate rR: Anzahl der Wiederhersteller an diesem Tag, generiert aus der Anzahl der Personen, die am Vortag behandelt werden müssen (C)

Tatsächlich gibt es eine Zeitverzögerung vom Krankenhausaufenthalt bis zum Tod oder vom Krankenhausaufenthalt bis zur Genesung. Daher muss dies bei der Berechnung berücksichtigt werden. Da jedoch keine Daten zur für den Tod oder die Genesung pro Patient erforderlichen Zeit enthalten sind, ist dies einfach. Ich entschied mich zu rechnen. Obwohl es einfach ist, scheint die Anzahl der Personen, die eine Behandlung benötigen (C), Patienten mit unterschiedlichen Tagen nach dem Krankenhausaufenthalt zu umfassen, so dass dies in gewissem Maße als Index angesehen wird. Die Anzahl der behandlungsbedürftigen Personen wird übrigens nach folgender Formel berechnet.

Versuchen Sie mit Python zu rechnen

Dieser Code ist auf GitHub verfügbar. Es wird im Jupyter Notebook-Format gespeichert. (Dateiname: 03_R0_estimation-WLD-RDR-01c.ipynb)

Code

Ich werde nicht den gesamten Code im Artikel behandeln, da er lang sein wird, aber ich werde die wichtigsten Punkte in der Berechnung erläutern.

Die Daten der Johns Hopkins University werden nach Land und Bundesstaat aggregiert. Da jedoch die nach Bundesstaat aufgelisteten Daten und die nach Land insgesamt aufgeführten Daten gemischt sind, werden sie nach Ländern neu aggregiert. ..

def getSumOfCountry01(df, country):
    df1 = df[df["Country/Region"] == country]
    # summarise for each country/region
    if df1['Province/State'].isnull().any():
        df2 = df1.fillna({'Province/State':'SUM'})
        df3 = df2[df2['Province/State'] == 'SUM']
    else:
        df3 = pd.DataFrame(df1.sum(axis=0)).T
        df3['Province/State'] = 'SUM'
        df3['Country/Region'] = df1.iloc[0,1]    
    return df3

Dies ist der Teil, der die obige Sterblichkeitsrate rD und die Wiederfindungsrate rR berechnet. Der gleitende Durchschnitt wird unter Berücksichtigung der Zeitverzögerung ermittelt.

def calcCRD(df, keys):
    #
    nth = keys['nth']
    tf = keys['tf']
    #
    nrow  = len(df)
    getV  = lambda s, tag: df.loc[s, tag] if s < nrow else np.NaN
    # tf [days]Gleitender Durchschnitt bis zum Ende
    getV2 = lambda s, tag: np.nanmean([getV(u, tag) for u in range(s,s + tf)])
    df.loc[0, 'C'] = 0.
    for t in range(1, nrow):
        df.loc[t, 'C'] = getV(t, 'PS') - getV(t, 'RS') - getV(t, 'DS')
        if df.loc[t-1, 'C'] > nth: #Für Maßnahmen gegen kleinen Nenner
            df.loc[t, 'rR'] = getV2(t, 'R') / df.loc[t-1, 'C']
            df.loc[t, 'rD'] = getV2(t, 'D') / df.loc[t-1, 'C']
        else:
            df.loc[t, 'rR'] = np.NaN
            df.loc[t, 'rD'] = np.NaN
    return df

Mehrere Dateien werden zu einem Datenrahmen zusammengefasst.

def CRDinWorldArea01(area):
    keys = { 'nth':100, 'tf': 3 }
    df1 = makeCalcFrame(100) #
    df2 = readCsvOfWorldArea_Confirmed(area)
    df3 = readCsvOfWorldArea_Recovered(area)
    df4 = readCsvOfWorldArea_Deaths(   area)
    df = df1
    df = mergeCalcFrame(df, df2)
    df = mergeCalcFrame(df, df3)
    df = mergeCalcFrame(df, df4)
    df = calcCRD(df, keys)
    return df

Ergebnis der Analyse

Schauen wir uns nun die Berechnungsergebnisse an. Das Berechnungsziel sind Länder, in denen die Zahl der Infizierten zum 15. April 5.000 übersteigt.

Gesamttrend

Zeichnen Sie zunächst die durchschnittliche Wiederfindungsrate rR auf der horizontalen Achse und die durchschnittliche Mortalität rD auf der vertikalen Achse. CRD_COVID-19_Map_WLD.png Die Sterblichkeitsraten variieren erheblich von Land zu Land. Zu den Ländern mit besonders hohen Sterblichkeitsraten zählen EU-Länder wie Spanien, das Vereinigte Königreich, Italien, Frankreich sowie Belgien und der Iran. Umgekehrt gehören zu den Ländern mit niedrigen Sterblichkeitsraten Australien, Israel, Norwegen, Chile, Südkorea und Deutschland. Die Wiederherstellungsrate ist ebenfalls sehr unterschiedlich, aber es ist charakteristisch, dass der Iran und Peru hoch sind. China scheint eine hohe Erholungsrate zu haben, da es seit der Pandemie lange her ist.

Vergleich innerhalb von G7

Lassen Sie uns einen Vergleich innerhalb von G7 machen. Erstens ist die Sterblichkeitsrate. Das Land links hat eine relativ niedrige Sterblichkeitsrate und das Land rechts hat eine relativ hohe Sterblichkeitsrate. CRD_COVID-19_DTR_G7.png In Ländern mit relativ niedrigen Sterblichkeitsraten (links) sinken die Sterblichkeitsraten tendenziell, wenn die Zahl der behandlungsbedürftigen Personen zunimmt. Insbesondere die Zahl der Infizierten in den Vereinigten Staaten ist die höchste der Welt, wird jedoch auf etwa 0,5% unterdrückt. In Ländern mit relativ hohen Sterblichkeitsraten (rechts) steigt die Sterblichkeitsrate tendenziell mit zunehmender Anzahl der behandlungsbedürftigen Personen. Wenn die Anzahl der behandlungsbedürftigen Personen 1000 oder mehr erreicht, steigt die Sterblichkeitsrate schnell von etwa 1% auf maximal etwa 3%.

Schauen wir uns die Wiederherstellungsrate an. CRD_COVID-19_RTR_G7.png In Ländern mit relativ niedrigen Sterblichkeitsraten (links) mit Ausnahme der Vereinigten Staaten scheinen die Genesungsraten zu steigen, wenn die Zahl der behandlungsbedürftigen Personen relativ gering ist. In Kanada und Deutschland steigt die Genesungsrate sprunghaft an, wenn die Zahl der behandlungsbedürftigen Personen ein bestimmtes Niveau erreicht. In Ländern mit relativ hohen Sterblichkeitsraten (rechts) ** sinkt die Genesungsrate tendenziell, wenn die Zahl der behandlungsbedürftigen Personen steigt **.

Vergleich von Ländern mit hohen und niedrigen Sterblichkeitsraten weltweit

Schauen wir uns einen Vergleich zwischen Ländern mit hohen und niedrigen Sterblichkeitsraten weltweit an. Erstens ist die Sterblichkeitsrate. Das Land links hat eine relativ niedrige Sterblichkeitsrate und das Land rechts hat eine relativ hohe Sterblichkeitsrate. CRD_COVID-19_DTR_GoodBad.png In Ländern mit relativ geringer Sterblichkeit (links) sind flache Formen von 0,5% oder weniger zu sehen. In Ländern mit relativ hohen Sterblichkeitsraten (rechts) steigen ausnahmslos ** die Sterblichkeitsraten mit zunehmender Anzahl der behandlungsbedürftigen Personen **.

Schauen wir uns die Wiederherstellungsrate an. CRD_COVID-19_RTR_GoodBad.png In Ländern mit relativ niedrigen Sterblichkeitsraten (links) sind die Genesungsraten bei einer bestimmten Anzahl von Personen, die eine Behandlung benötigen, in die Höhe geschossen. In Ländern mit relativ hohen Sterblichkeitsraten (rechts) sinken die Genesungsraten tendenziell, wenn die Zahl der behandlungsbedürftigen Personen zunimmt.

BCG japanische Stamm geimpfte Länder und Nachbarländer

Ein Artikel darüber, ob eine routinemäßige BCG-Impfung COVID-19 beeinflussen kann wurde in Deutschland, den Niederlanden, Australien usw. veröffentlicht. Die klinische Studie hat begonnen. Schauen wir uns also die Auswirkungen von BCG an, insbesondere japanischer Aktien. Erstens ist die Sterblichkeitsrate. Das linke ist das Land, das mit japanischen BCG-Stämmen geimpft wurde, und das rechte ist das Nachbarland, in das japanische Stämme ohne BCG geimpft wurden. CRD_COVID-19_DTR_BCG.png Die Sterblichkeitsrate ist im Iran hoch, aber die Sterblichkeitsrate ist in anderen Ländern, die nicht mit BCG-japanischen Stämmen geimpft wurden, relativ unterdrückt.

Schauen wir uns die Wiederherstellungsrate an. CRD_COVID-19_RTR_BCG.png In Ländern, die mit BCG-Stämmen geimpft wurden, beträgt die Anzahl der behandlungsbedürftigen Personen weniger als 1000 und die Wiederherstellungsrate ist hoch. In Ländern, die nicht mit BCG-japanischen Stämmen geimpft wurden, steigt die Wiederherstellungsrate jedoch nur an, wenn die Anzahl der Personen, die eine Behandlung benötigen, auf 1000 oder mehr steigt. .. Mit anderen Worten, ** BCG-Inokulation mit japanischen Stämmen hat keinen Effekt auf die Senkung der Sterblichkeitsrate, kann jedoch den Effekt haben, dass die Wiederfindungsrate in einem frühen Stadium erhöht wird **.

Erwägung

Aus dem Obigen wurden die folgenden Trends aus den Daten als Verifizierungsergebnisse bezüglich des medizinischen Kollapses gefunden.

Außerdem ...

Referenzlink

Ich habe auf die folgende Seite verwiesen.

  1. Medizinischer Kollaps
  2. Daten veröffentlicht von der Johns Hopkins University
  3. Association of BCG vaccination policy with prevalence and mortality of COVID-19
  4. Correlation between universal BCG vaccination policy and reduced morbidity and mortality for COVID-19: an epidemiological study
  5. JSatoNotes

Recommended Posts

Lassen Sie uns die medizinische Kollapshypothese des neuen Koronavirus testen
Zeichnen Sie die Ausbreitung des neuen Koronavirus
Schätzen Sie die maximale Infektiosität des neuen Coronavirus ab
Berechnen wir den Übergang der Grundreproduktionszahl des neuen Koronavirus nach Präfektur
GUI-Simulation des neuen Koronavirus (SEIR-Modell)
Lassen Sie uns die Konvergenzzeit anhand des globalen Trends der effektiven Reproduktionszahl des neuen Koronavirus untersuchen
Lassen Sie uns eine Rangfolge der Anzahl der wirksamen Reproduktionen des neuen Koronavirus nach Präfektur erstellen
Stellen wir uns die Anzahl der mit Matplotlib mit dem Coronavirus infizierten Personen vor
Quantifizieren Sie den Grad der Selbstbeherrschung, der erforderlich ist, um das neue Koronavirus einzudämmen
Simulieren wir den Effekt der Einführung einer Kontaktverfolgungs-App als Gegenmaßnahme gegen den neuen Koronavirus
Bedeutung der Eingrenzung des Testziels des PCR-Tests auf neues Koronavirus nach dem Bayes-Theorem
Messen Sie die Testabdeckung von Push-Python-Code auf GitHub.
Testen Sie die Eignung der Verteilung
Hypothese, warum das neue Koronavirus in städtischen Gebieten wie Tokio nicht so häufig vorkommt
Werfen wir einen Blick auf die Infektionstendenz des neuen Coronavirus COVID-19 in jedem Land und den Status der medizinischen Reaktion (zusätzliche Informationen).
Analysieren Sie den altersspezifischen Schweregrad des Coronavirus
Ich habe versucht, PDF-Daten der medizinischen Online-Versorgung zu verwenden, die auf der Ausbreitung einer neuen Coronavirus-Infektion basieren
Lassen Sie uns den Gewinner des Bingo bestimmen
Die Geschichte des Studenten, der die neue Website für Gegenmaßnahmen gegen Koronaviren entwickelt hat (Ishikawa-Version)
Ich habe versucht, das Verhalten des neuen Koronavirus mit dem SEIR-Modell vorherzusagen.
Folding @ Home unter Linux Mint soll zur Analyse des neuen Corona-Virus beitragen
Untersuchen wir den Mechanismus von Kaijis Chinchirorin
[AWS] Lassen Sie uns einen Komponententest der Lambda-Funktion in der lokalen Umgebung durchführen
Die epidemische Vorhersage des neuen Koronavirus wurde mit explosiver Geschwindigkeit im Internet veröffentlicht
Ich habe versucht, die Literatur des neuen Corona-Virus mit Python automatisch an LINE zu senden
[Python] Testen Sie die Mondmatagi des relativen Deltas
Ich habe einige der neuen Funktionen von Python 3.8 touched angesprochen
[Django] Lassen Sie uns versuchen, den Teil von Django zu klären, der im Test irgendwie durch war
Überprüfen Sie die Wirkung des Urlaubs als Gegenmaßnahme gegen das neue Koronavirus mit dem SEIR-Modell
Die Theorie, dass der Schlüssel zur Kontrolle der Infektion des neuen Coronavirus die Hyperdispersion der Anfälligkeit ist.
Ich habe versucht, die Eigenschaften der neuen Informationen über mit dem Corona-Virus infizierte Personen mit Wordcloud zu visualisieren