Examiner la plage d'erreur dans le nombre de décès dus à la pneumonie

Si les décès de Covid19 étaient comptabilisés pour une pneumonie, nous aimerions savoir à quel point cela ne peut plus être considéré comme une erreur. Puisqu'il existe clairement des fluctuations et des tendances saisonnières, nous examinerons la variation des résidus après avoir retiré chacun d'eux.

Analyse des séries chronologiques

Les données proviennent des statistiques de l'état civil. Nous confirmerons l'évolution mensuelle du nombre de décès de 2,015 à 2,018. L'axe vertical représente le nombre de décès et l'axe horizontal le mois.

4nen.png

Cycle de 12 mois sans regarder l'autocorrélation partielle. Chaque année, il y a une tendance générale à la hausse.

Cette fois, le nombre de décès dus à la pneumonie est considéré comme 9% du total (à partir du dernier document du poste), et il est décomposé en tendances, fluctuations saisonnières et résidus à l'aide d'un modèle additif. Utilisez season_decompose de StatsModels dans un modèle additif (il doit répéter la moyenne mobile pour calculer la variation saisonnière)

Y[t] = T[t] + S[t] + e[t]

T est la tendance, S est la fluctuation saisonnière, e est le résidu


import statsmodels.api as sm

#
#les décès sont des données qui ne contiennent que le nombre de décès dans une colonne
# 

dec = sm.tsa.seasonal_decompose(death,period=12)
trend = dec.trend
seasonal = dec.seasonal
residual = dec.resid

Graphique de données

tendance

trend.png

Fluctuation saisonnière

seasonal.png

Erreur résiduelle

residual2.png

Conclusion

La valeur résiduelle minimale est -409,1 et la valeur maximale est 282,8. L'écart type est de 136,5. La moyenne est de -4,59. Lorsque le nombre augmente de 1 000 ou plus, on a l'impression qu'il contient probablement beaucoup de Covid19.

Si le taux de létalité est de 2% et le nombre de personnes infectées est d'environ 50 000, la ligne de 1 000 décès.

Les patients Covid19 étant peu répartis dans tout le pays, il est difficile de voir l'augmentation des décès par pneumonie dans les hôpitaux (le nombre d'hôpitaux sans hôpitaux psychiatriques devrait être supérieur à 7 000, environ un de plus. Pas même).

Étant donné que les données par cause de décès à l'échelle nationale sont retardées de 5 mois, Afin de saisir la situation actuelle, il est nécessaire de chercher une place rapide dans le bulletin du gouvernement local. Cela peut ne pas exister.

À propos de la situation actuelle

À l'heure actuelle, deux semaines se sont écoulées depuis que la réunion d'experts a désigné les semaines restantes d'une à deux semaines au Japon. Une déclaration d'urgence a été émise à New York et les dommages dans la zone de l'UE deviennent graves.

Au Japon, aucun dommage n'a été signalé autant que dans l'UE et en Corée du Sud, Les caractéristiques du principal virus dominant peuvent être différentes, le système médical peut être différent, etc. Il existe plusieurs facteurs possibles et il n'y a pas suffisamment de données pour établir un lien avec un scénario spécifique.

Si l'environnement est équivalent à celui du virus de Wuhan et similaire à celui de l'Europe et des États-Unis, Le nombre de décès est inférieur aux prévisions, nous avons donc examiné si les données pouvaient détecter si les décès de Covid 19 faisaient partie des décès par pneumonie.

(Je fixe simplement le virus et l'environnement en tant que variables, je ne juge pas cela, et la vérification d'autres scénarios est une autre histoire. J'ai beaucoup confirmé cette fois)

Les données

Données de mortalité mensuelles (taux annuel de mortalité et de mortalité par mois de décès (pour 1000 habitants)): https://www.e-stat.go.jp/stat-search/files?page=1&layout=datalist&toukei=00450011&tstat=000001028897&cycle=7&year=20180&month=0&tclass1=000001053058&tclass2=000001053061&tclass3=000001053065&result_back=1

Il existe un rapport préliminaire sur le nombre total de décès par mois, et le dernier date d'il y a 2 mois (Bulletin des statistiques de l'état civil du 19 novembre): https://www.mhlw.go.jp/toukei/saikin/hw/jinkou/geppo/s2019/11.html

Nombre mensuel approximatif par cause de décès (il y a 5 mois au plus tard) https://www.e-stat.go.jp/stat-search/files?page=1&layout=datalist&toukei=00450011&tstat=000001028897&cycle=1&year=20190&month=23070909&tclass1=000001053058&tclass2=000001053060

Pneumonie (Décès par pneumonie au Japon du point de vue de l'état civil): http://www.tokyo-eiken.go.jp/sage/sage2018/

Recommended Posts

Examiner la plage d'erreur dans le nombre de décès dus à la pneumonie
L'histoire d'une erreur dans PyOCR
Sortie du nombre de cœurs de processeur en Python
Trouvez le nombre de jours dans un mois
Utilisé depuis l'introduction de Node.js dans l'environnement WSL
Comment obtenir le nombre de chiffres en Python
Comptez le nombre de paramètres dans le modèle d'apprentissage en profondeur
Calcul du nombre minimum de voix requis à partir du taux de vote
Obtenir la taille (nombre d'éléments) de Union Find en Python
Examinez les paramètres de RandomForestClassifier dans le didacticiel Kaggle / Titanic
Visualisez le nombre de plaintes des compagnies d'assurance-vie
Examinons le temps de convergence de la tendance mondiale du nombre de reproduction effectif du nouveau virus corona
[Python] Représentation du nombre de plaintes des compagnies d'assurance-vie dans un graphique à barres
Comment trouver le nombre optimal de clusters pour les k-moyennes
Maya | Découvrez le nombre de polygones dans l'objet sélectionné
Obtenez le nombre de chiffres
Comptez le nombre de caractères dans le texte dans le presse-papiers sur Mac
Obtenez le nombre d'éléments spécifiques dans la liste python
Python --Trouvez le nombre de groupes dans l'expression regex
VisibleDeprecation Avertissement: l'utilisation d'un nombre non entier au lieu d'un entier entraînera une erreur dans le futur
[Homologie] Comptez le nombre de trous dans les données avec Python
Obtenez le nombre d'occurrences pour chaque élément de la liste
Calculez le nombre de changements
J'ai essayé de trouver la tendance du nombre de navires dans la baie de Tokyo à partir d'images satellites.
[Python] Réduisons le nombre d'éléments dans le résultat dans le fonctionnement de l'ensemble
Erreur dans l'ordre lors de l'importation de Numpy dans l'environnement Windows Anaconda
Obtenez le nombre de lecteurs d'articles sur Mendeley en Python
Obtenez le nombre de vues de Qiita
L'histoire de la participation à AtCoder
Existence du point de vue de Python
Calcul du nombre d'associations de Klamer
L'histoire du "trou" dans le fichier
Obtenez le nombre d'abonnés Youtube
Examiner l'erreur de coupure de la formule de Lie-Trotter
Examiner la classe d'un objet avec python
[Question] Dans la régression de forêt aléatoire de sk-learn, une erreur se produit lorsque le nombre de parallèles est défini sur -1.
[Version terminée] Essayez de connaître le nombre d'habitants de la ville à partir de la liste d'adresses avec Python
Une histoire sur la création d'un programme qui augmentera le nombre d'abonnés Instagram de 0 à 700 en une semaine
"Le gars qui prédit le nombre de vues à partir du titre de la vidéo de Jar Jal"
Graphique de l'historique du nombre de couches de deep learning et du changement de précision
Générez une liste contenant le nombre de jours du mois en cours.
Vérifions la chaîne d'octets en mémoire du nombre flottant flottant en Python
[Python] Calculez le nombre de chiffres requis lors de la saisie de 0 [Note]
Représentez graphiquement l'évolution du nombre d'apparitions de mots clés par mois à l'aide de pandas
Tendances des langages de programmation du point de vue de GitHub (mises à jour semestrielles)
J'ai essayé de prédire le nombre de personnes infectées par le virus corona en tenant compte de l'effet de s'abstenir de sortir
[Compétition Pro] Résolvez le nombre de M cartes extraites de N cartes en utilisant un itinéraire [Explication sur la figure]
[Comprendre en 3 minutes] Le début de Linux
Vérifiez le comportement du destroyer en Python
Compter / vérifier le nombre d'appels de méthode.
Implémenter une partie du processus en C ++
Notes d'apprentissage depuis le début de Python 1
Omettre la nomenclature depuis le début de la chaîne
Quelle est la cause de l'erreur suivante?
Résoudre l'erreur d'adresse déjà utilisée
Le résultat de l'installation de python sur Anaconda