[PYTHON] 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.

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

Fluctuation saisonnière

Erreur résiduelle

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/