Income Distribution ECDF

import matplotlib.pyplot as plt
import numpy as np

fig, ax = plt.subplots(figsize=(10, 6), facecolor='#0a0a0f')
ax.set_facecolor('#0a0a0f')

np.random.seed(42)
income = np.random.lognormal(10.5, 0.8, 2000)

ax.ecdf(income, color='#22c55e', linewidth=2.5)

median = np.median(income)
ax.axhline(y=0.5, color='#ef4444', linestyle='--', alpha=0.7, label=f'Median: ${median:,.0f}')
ax.axvline(x=median, color='#ef4444', linestyle='--', alpha=0.7)

ax.set_xlabel('Annual Income ($)', fontsize=12, color='white')
ax.set_ylabel('Cumulative Probability', fontsize=12, color='white')
ax.set_title('Income Distribution ECDF', fontsize=16, fontweight='bold', color='white', pad=20)
ax.tick_params(colors='white')
ax.set_xlim(0, 200000)
ax.legend(facecolor='#1a1a2f', edgecolor='#333', labelcolor='white')
ax.grid(True, alpha=0.2, color='white')

for spine in ax.spines.values():
    spine.set_color('#333')

plt.tight_layout()