Pharmacokinetic Half-Life

import matplotlib.pyplot as plt
import numpy as np

np.random.seed(42)
drugs = ['Antibiotics', 'NSAIDs', 'Statins', 'SSRIs', 'Beta-blockers']
data = [
    np.random.lognormal(1, 0.5, 120),
    np.random.lognormal(1.5, 0.6, 120),
    np.random.lognormal(2.5, 0.4, 120),
    np.random.lognormal(3, 0.5, 120),
    np.random.lognormal(2, 0.7, 120)
]

fig, ax = plt.subplots(figsize=(10, 6), dpi=100)
ax.set_facecolor('#0d1117')
fig.patch.set_facecolor('#0d1117')

colors = ['#F5276C', '#F5B027', '#27D3F5', '#4927F5', '#6CF527']
bp = ax.boxplot(data, widths=0.5, patch_artist=True, showfliers=False, notch=True,
                medianprops=dict(color='#F5D327', linewidth=2.5))

for patch, color in zip(bp['boxes'], colors):
    patch.set_facecolor(color)
    patch.set_alpha(0.8)
    patch.set_edgecolor('white')
for i, color in enumerate(colors):
    bp['whiskers'][i*2].set_color(color)
    bp['whiskers'][i*2+1].set_color(color)
    bp['caps'][i*2].set_color(color)
    bp['caps'][i*2+1].set_color(color)

ax.set_xticklabels(drugs, rotation=15, ha='right')
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.spines['left'].set_color('#333333')
ax.spines['bottom'].set_color('#333333')
ax.yaxis.grid(True, color='#1a1a3f', linewidth=0.5, zorder=0)
ax.set_axisbelow(True)
ax.tick_params(axis='both', colors='#888888', labelsize=9, length=0, pad=8)
ax.set_ylabel('Half-Life (hours)', fontsize=11, color='white', fontweight='500')
ax.set_title('Drug Pharmacokinetic Half-Life', fontsize=14, color='white', fontweight='bold', pad=15)

plt.tight_layout()
plt.show()