日本の交通事故死亡者数の推移を性別・年齢別に棒グラフ化する

from pandas import *

df = read_csv('caraccident.csv',encoding='cp932')
df.head(2)

df.drop(['表章項目','死因年次推移分類','性別','/時間軸(年次)'],axis=1,inplace=True)
df.head(2)

df = df.convert_objects(convert_numeric=True)

from matplotlib.pyplot import *
from matplotlib.font_manager import FontProperties
from matplotlib import rcParams
style.use('ggplot')

fp = FontProperties(fname='/usr/share/fonts/opentype/ipaexfont-gothic/ipaexg.ttf', size=54)
rcParams['font.family'] = fp.get_name()
rcParams["font.size"] = "20"
fig, ax = subplots(figsize=(25,20))
df[:23][1:].set_index('年齢(5歳階級)').T.plot(kind='barh',ax=ax,stacked=True)
ax.legend(loc='lower right', fancybox=True, framealpha=0.5);
xticks(np.arange(0,2.6e4,1e4/5),
    ['{}万'.format(float(x/1e4)) if x > 0 else 0 for x in np.arange(0,2.6e4,1e4/5)]);

df1 = df[:5][1:].groupby(np.arange(len(df[:5][1:]))//4).sum()
df1.index = df[:5][1:].loc[1::4, '年齢(5歳階級)']
df2 = df[:9][5:].groupby(np.arange(len(df[:9][5:]))//4).sum()
df2.index = df[:9][5:].loc[1::4, '年齢(5歳階級)']
df3 = df[:13][9:].groupby(np.arange(len(df[:13][9:]))//4).sum()
df3.index = df[:13][9:].loc[1::4, '年齢(5歳階級)']
df4 = df[:22][13:].groupby(np.arange(len(df[:22][13:]))//9).sum()
df4.index = df[:22][13:].loc[1::9, '年齢(5歳階級)']
df5 = concat([df1,df2,df3,df4], axis=0)
df5.rename(index={'0～4歳':'0〜19歳','20～24歳':'20〜39歳','40～44歳':'40〜59歳','60～64歳':'60歳〜'},inplace=True)
df5.reset_index(inplace=True)

style.use('fivethirtyeight')

fp = FontProperties(fname='/usr/share/fonts/opentype/ipaexfont-gothic/ipaexg.ttf', size=54)
rcParams['font.family'] = fp.get_name()
rcParams["font.size"] = "18"
fig, ax = subplots(figsize=(18,10))
df5.set_index('年齢(5歳階級)').T.plot(kind='barh',ax=ax,stacked=True)
ax.legend(loc='lower right',fancybox=True,framealpha=0.5,prop={'size': 26})
xticks(np.arange(0,2.6e4,1e4/5),
    ['{}万'.format(float(x/1e4)) if x > 0 else 0 for x in np.arange(0,2.6e4,1e4/5)]);

df6 = df[['年齢(5歳階級)','2017年']][:44][24:]
df6.rename(columns={'2017年':'Male'}, inplace=True)
df6

df7 = df[['年齢(5歳階級)','2017年']][:67][47:]
df7.rename(columns={'2017年':'Female'}, inplace=True)
df7

df8 = merge(df6,df7,how='inner',on='年齢(5歳階級)')

df8['Total'] = df8['Male']+df8['Female']

df8['Female_Pct'] = df8['Female'] / df8['Total']
df8['Male_Pct'] = df8['Male'] / df8['Total']

rcParams['figure.figsize'] = 20,12
rcParams["font.size"] = "16"

for (i, row) in df8.iterrows():
    bar([row['年齢(5歳階級)']], [row['Male_Pct']],
        color=['#008AB8'], width=0.9, align='center', edgecolor='none')
    bar([row['年齢(5歳階級)']], [row['Female_Pct']], bottom=[row['Male_Pct']],
            color=['#CC6699'], width=0.9, align='center', edgecolor='none')
    
xlim(-0.6, 19.6)
ylim(0, 1)
xticks(range(0, 20), df8['年齢(5歳階級)'], fontsize=13.5)
grid(False)
yticks(np.arange(0, 1.01, 0.25),
       ['{}%'.format(int(x * 100)) for x in np.arange(0, 1.01, 0.25)])
xlabel('Age group')
ylabel('Percentage of age group')

plot([-0.425, 20.425], [0.5, 0.5], lw=2, color='black')
plot([-0.425, 20.425], [0.25, 0.25], lw=2, color='black')
plot([-0.425, 20.425], [0.75, 0.75], lw=2, color='black');