Violin with Box Plot

import matplotlib.pyplot as plt
import numpy as np

# Data
np.random.seed(42)
data = [np.random.exponential(2, 300) + 2,
        np.random.normal(6, 1.5, 300),
        np.random.gamma(3, 1, 300) + 1,
        np.random.normal(5, 2, 300)]
labels = ['Exponential', 'Normal', 'Gamma', 'Bimodal']

# Colors
colors = ['#10B981', '#3B82F6', '#F59E0B', '#EF4444']

# Create figure
fig, ax = plt.subplots(figsize=(10, 6), facecolor='white')

positions = [1, 2, 3, 4]

# Violin plot
vp = ax.violinplot(data, positions=positions, widths=0.7,
                   showmeans=False, showmedians=False, showextrema=False)

for i, body in enumerate(vp['bodies']):
    body.set_facecolor(colors[i])
    body.set_edgecolor('white')
    body.set_alpha(0.6)

# Box plot overlay
bp = ax.boxplot(data, positions=positions, widths=0.15,
                patch_artist=True, showfliers=False)

for i, (box, median) in enumerate(zip(bp['boxes'], bp['medians'])):
    box.set_facecolor('white')
    box.set_edgecolor(colors[i])
    box.set_linewidth(2)
    median.set_color(colors[i])
    median.set_linewidth(2)

for element in ['whiskers', 'caps']:
    for i, item in enumerate(bp[element]):
        item.set_color(colors[i // 2])
        item.set_linewidth(1.5)

# Customize axes
ax.set_xticks(positions)
ax.set_xticklabels(labels, fontsize=11, fontweight='500')
ax.set_ylabel('Value', fontsize=12, fontweight='500', color='#374151')
ax.set_xlim(0.4, 4.6)

# Clean styling
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.spines['left'].set_color('#E5E7EB')
ax.spines['bottom'].set_color('#E5E7EB')
ax.tick_params(colors='#6B7280', labelsize=10)
ax.yaxis.grid(True, linestyle='--', alpha=0.3, color='#9CA3AF')

plt.tight_layout()
plt.show()