大阪市の超過死亡?

[2020-06-24] 2020年5月分のデータを追加した。

東京の超過死亡?に倣って大阪市を調べてみる。毎月の死者数を納めた osakadeaths.csv は,急いで作ったため,どなたかチェックしていただければありがたい(1箇所データの抜けがあったので修正)。

import pandas as pd
import matplotlib.pyplot as plt
import numpy as np

df = pd.read_csv('https://oku.edu.mie-u.ac.jp/~okumura/python/data/osakadeaths.csv')

def days(year, month):
    m = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31]
    return m[month - 1] + ((month == 2) and
                           ((year % 4 == 0) and (year % 100 != 0) or (year % 400 == 0)))

perday = np.array([r[2] / days(r[0], r[1]) for i, r in df.iterrows()])

for y in range(2015, 2021):
    plt.plot(df[df['year'] == y]['month'], perday[df['year'] == y], 'o-', label=y)

plt.legend()
plt.savefig('osakadeaths.svg', bbox_inches="tight")
Deaths in Osaka
plt.clf()
plt.plot(df[df['month'] == 1]['year'], perday[df['month'] == 1], 'o-', label='Jan')
plt.plot(df[df['month'] == 2]['year'], perday[df['month'] == 2], 'o-', label='Feb')
plt.plot(df[df['month'] == 3]['year'], perday[df['month'] == 3], 'o-', label='Mar')
plt.plot(df[df['month'] == 4]['year'], perday[df['month'] == 4], 'o-', label='Apr')
plt.plot(df[df['month'] == 5]['year'], perday[df['month'] == 5], 'o-', label='May')
plt.xticks(range(2016, 2021))
plt.legend()
plt.savefig('osakadeaths1.svg', bbox_inches="tight")
Deaths in Osaka

[TODO] 大阪府にも大阪市の過去の死亡数が載っているが微妙に数値が異なるらしい。要チェック。


Last modified: