Income Distribution

import matplotlib.pyplot as plt
import numpy as np

COLORS = {'bars': '#F59E0B', 'median': '#EF4444', 'mean': '#6366F1', 'background': '#FFFFFF', 'text': '#1E293B', 'text_muted': '#64748B', 'grid': '#F1F5F9'}

np.random.seed(42)
income = np.random.lognormal(10.5, 0.8, 2000)
income = income[income < 300000]

fig, ax = plt.subplots(figsize=(10, 6), dpi=100)
ax.set_facecolor(COLORS['background'])
fig.patch.set_facecolor(COLORS['background'])

ax.hist(income/1000, bins=40, color=COLORS['bars'], edgecolor='white', linewidth=0.5, alpha=0.8)

median, mean = np.median(income)/1000, np.mean(income)/1000
ax.axvline(median, color=COLORS['median'], linewidth=2, linestyle='--', label=f'Median: ${median:.0f}K')
ax.axvline(mean, color=COLORS['mean'], linewidth=2, linestyle='--', label=f'Mean: ${mean:.0f}K')

ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.spines['left'].set_color(COLORS['grid'])
ax.spines['bottom'].set_color(COLORS['grid'])
ax.yaxis.grid(True, color=COLORS['grid'], linewidth=1, zorder=0)
ax.set_axisbelow(True)
ax.tick_params(axis='both', colors=COLORS['text_muted'], labelsize=9, length=0, pad=8)
ax.set_xlabel('Income ($K)', fontsize=10, color=COLORS['text'], labelpad=10)
ax.set_ylabel('Count', fontsize=10, color=COLORS['text'], labelpad=10)
ax.legend(loc='upper center', bbox_to_anchor=(0.5, -0.12), ncol=2, frameon=False, fontsize=10, labelcolor=COLORS['text_muted'])

plt.tight_layout()
plt.show()