Heart Rate by Exercise Intensity

import numpy as np
import matplotlib.pyplot as plt

np.random.seed(305)
BG_COLOR = '#0a0a0f'
TEXT_COLOR = 'white'
COLORS = ['#6CF527', '#F5B027', '#F5276C', '#9C2007']

def simple_beeswarm(y, nbins=None, width=1.):
    y = np.asarray(y)
    if nbins is None:
        nbins = np.ceil(len(y) / 6).astype(int)
    nn, ybins = np.histogram(y, bins=nbins)
    nmax = nn.max()
    x = np.zeros(len(y))
    ibs = []
    for ymin, ymax in zip(ybins[:-1], ybins[1:]):
        i = np.nonzero((y > ymin) * (y <= ymax))[0]
        ibs.append(i)
    dx = width / (nmax // 2) if nmax > 1 else width
    for i in ibs:
        yy = y[i]
        if len(i) > 1:
            j = len(i) % 2
            i = i[np.argsort(yy)]
            a = i[j::2]
            b = i[j+1::2]
            x[a] = (0.5 + j / 3 + np.arange(len(b))) * dx
            x[b] = (0.5 + j / 3 + np.arange(len(b))) * -dx
    return x

intensities = ['Rest', 'Light', 'Moderate', 'High']
data = {
    'Rest': np.random.normal(65, 8, 40),
    'Light': np.random.normal(95, 12, 45),
    'Moderate': np.random.normal(135, 15, 50),
    'High': np.random.normal(170, 12, 42)
}

fig, ax = plt.subplots(figsize=(10, 6), facecolor=BG_COLOR)
ax.set_facecolor(BG_COLOR)

boxplot_data = []
for i, (intensity, values) in enumerate(data.items()):
    x = simple_beeswarm(values, width=0.3)
    ax.scatter(x + i + 1, values, c=COLORS[i], alpha=0.7, s=45, edgecolors='white', linewidth=0.5)
    boxplot_data.append(values)

bp = ax.boxplot(boxplot_data, positions=range(1, len(intensities)+1), widths=0.5, patch_artist=True)
for patch in bp['boxes']:
    patch.set_facecolor('none')
    patch.set_edgecolor('#666')
for element in ['whiskers', 'caps', 'medians']:
    for item in bp[element]:
        item.set_color('#666')

ax.set_xticks(range(1, len(intensities)+1))
ax.set_xticklabels(intensities, color=TEXT_COLOR)
ax.set_xlabel('Exercise Intensity', fontsize=12, color=TEXT_COLOR, fontweight='500')
ax.set_ylabel('Heart Rate (BPM)', fontsize=12, color=TEXT_COLOR, fontweight='500')
ax.set_title('Heart Rate by Exercise Level', fontsize=14, color=TEXT_COLOR, fontweight='bold', pad=15)

ax.tick_params(colors='#888', labelsize=10)
for spine in ax.spines.values():
    spine.set_color('#333')

plt.tight_layout()
plt.show()