Beeswarm Plot

Daily Commute Duration by Transport

Travel time distribution across commute methods

Output
Daily Commute Duration by Transport
Python
import numpy as np
import matplotlib.pyplot as plt

np.random.seed(312)
BG_COLOR = '#ffffff'
TEXT_COLOR = '#1f2937'
COLORS = ['#6CF527', '#27D3F5', '#F54927', '#5314E6']

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

modes = ['Bike', 'Transit', 'Car', 'Walk']
data = {
    'Bike': np.random.normal(22, 6, 45),
    'Transit': np.random.normal(38, 12, 55),
    'Car': np.random.normal(32, 15, 60),
    'Walk': np.random.normal(18, 5, 35)
}

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

boxplot_data = []
for i, (mode, values) in enumerate(data.items()):
    values = np.maximum(values, 5)
    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(modes)+1), widths=0.5, patch_artist=True)
for patch in bp['boxes']:
    patch.set_facecolor('none')
    patch.set_edgecolor('#9ca3af')
for element in ['whiskers', 'caps', 'medians']:
    for item in bp[element]:
        item.set_color('#9ca3af')

ax.set_xticks(range(1, len(modes)+1))
ax.set_xticklabels(modes, color=TEXT_COLOR)
ax.set_xlabel('Transportation Mode', fontsize=12, color=TEXT_COLOR, fontweight='500')
ax.set_ylabel('Time (min)', fontsize=12, color=TEXT_COLOR, fontweight='500')
ax.set_title('Daily Commute by Transport', fontsize=14, color=TEXT_COLOR, fontweight='bold', pad=15)

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

plt.tight_layout()
plt.show()
Library

Matplotlib

Category

Statistical

Did this help you?

Support PyLucid to keep it free & growing

Support