(1) Ich habe eine Trendanalyse der neuen Korona mit Pythons Statistikmodell versucht (2) Die Analyse der Daten des Ministeriums für Gesundheit, Arbeit und Soziales ergab, dass der erste Höchststand im April und dieser Höchststand im August lag. (3) Das bedeutet, dass der nächste Höhepunkt im Dezember sein wird.
(1) Informationen zu neuen mit Korona infizierten Personen erhalten Sie auf der Homepage des Ministeriums für Gesundheit, Arbeit und Soziales (2) StatsModel of Python zerlegt sich in Trends, saisonale Faktoren und Residuen
(Special Thanks to) Ich bezog mich auf Beginnen wir mit der Datenanalyse mit Momoki. Vielen Dank.
Laden Sie die Anzahl der positiven Personen von der [Website des Ministeriums für Gesundheit, Arbeit und Soziales] herunter (https://www.mhlw.go.jp/stf/covid-19/open-data.html). Ich war sehr beeindruckt, wie einfach es war, CSV-Daten herunterzuladen. Ministerium für Gesundheit, Arbeit und Soziales, unglaublich!
Ausführliche Analysemethoden finden Sie oben unter "Erste Schritte mit Momoki und Datenanalyse".
Zunächst vorbereitende Arbeiten.
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import statsmodels.api as sm
% matplotlib inline
Dann werden die heruntergeladenen Daten gelesen. Enthält Daten vom 16. Januar.
df=pd.read_csv('pcr_positive_daily.csv')
df.head()
Die neuesten Daten stammen übrigens von gestern. Grundsätzlich können Sie die Daten des Vortages abrufen, diese werden jedoch abends aktualisiert. Erai wird sogar samstags und sonntags aktualisiert!
df.tail()
Lassen Sie uns nun einen Blick auf die Veränderungen der Anzahl infizierter Personen bis zur Gegenwart seit der zweiten Januarhälfte werfen, als die Daten veröffentlicht wurden.
%matplotlib inline
df.plot()
Übrigens, wenn Sie Matplotlib so wie es ist verwenden, wird es sogenannter "Tofu" (lacht) Die folgende Seite war für die Lösung sehr hilfreich. Vielen Dank. So übersetzen Sie Google Collaboration Graph (matplotlib) ins Japanische!
Nun, endlich das Hauptthema. Zerlegen Sie die Daten des Ministeriums für Gesundheit, Arbeit und Soziales mit Pythons StatsModel.
numbers = pd.Series(df['Anzahl positiver PCR-Tests(Ein Tag)'], dtype='float')
numbers.index = pd.to_datetime(df['Datum'])
res = sm.tsa.seasonal_decompose(numbers)
original = numbers #Originale Daten
trend = res.trend #Trenddaten
seasonal = res.seasonal #Saisonale Daten
residual = res.resid #Restdaten
plt.figure(figsize=(8, 8)) #Erstellung eines Diagrammzeichnungsrahmens, Größenangabe
#Originaldatenplot
plt.subplot(411) #Grafik 4 Zeilen 1 Spalte 1. Position (oben)
plt.plot(original)
plt.ylabel('Original')
#Darstellung der Trenddaten
plt.subplot(412) #Zweite Position in 4 Zeilen und 1 Spalte des Diagramms
plt.plot(trend)
plt.ylabel('Trend')
#Darstellung der saisonalen Daten
plt.subplot(413) #3. Position in 4 Zeilen und 1 Spalte des Diagramms
plt.plot(seasonal)
plt.ylabel('Seasonality')
#Darstellung der Restdaten
plt.subplot(414) #4. Position in 4 Zeilen und 1 Spalte Grafik (unten)
plt.plot(residual)
plt.ylabel('Residuals')
plt.tight_layout() #Automatische Anpassung des Diagrammabstands
Die Ergebnisse sind wie folgt. Ich habe es in den Code geschrieben, aber in der Reihenfolge von oben, ① Originaldaten ② Trenddaten ③ Saisonale Daten ④ Restdaten Es wird sein.
Bitte beachten Sie die zweiten Trenddaten. Es ist zu sehen, dass der erste Gipfel Anfang April und dieser Gipfel Anfang August ist. Ist aus dieser Sicht der nächste Höhepunkt Anfang Dezember?
Ich bete nur, dass die neue Corona vorzeitig endet. Wir erkennen jedoch, dass die Realität hart ist.
Der Krankheitsname der neuen Korona lautet COVID-19, der Virusname jedoch SARS-CoV-2. Es scheint, dass dieser Virenname SARS ähnlich ist.
Das SARS wurde im Jahr 2002 populär, also war es vor fast 20 Jahren. Es scheint jedoch, dass der SARS-Impfstoff noch nicht hergestellt wurde.
Ich möchte ruhig tun, was ich kann, weil es eine schwierige Zeit ist.
Ich fühle mich einsam, dass die Trinkparty aufgrund des Einflusses der neuen Corona verschwunden ist, Es war auch eine gute Gelegenheit, rationale Arbeitsstile wie die Förderung der Telearbeit zu fördern.
Lass uns in solchen Zeiten ruhig leben: entspannt:
Last but not least möchte ich mich bei allen an der Website beteiligten Personen für ihre Referenz bedanken.
Recommended Posts