レインクラウドプロットで分布の全体像を一枚に

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53

import numpy as np
import matplotlib.pyplot as plt

rng = np.random.default_rng(42)
groups = ["Group A", "Group B", "Group C"]
data = [
    rng.normal(loc=50, scale=8, size=120),
    rng.normal(loc=62, scale=12, size=100),
    rng.exponential(scale=10, size=110) + 40,
]

fig, ax = plt.subplots(figsize=(8, 5))
positions = [1, 2, 3]

# --- half-violin (cloud) ---
vp = ax.violinplot(data, positions=positions, showextrema=False, widths=0.6)
for body in vp["bodies"]:
    m = np.mean(body.get_paths()[0].vertices[:, 0])
    body.get_paths()[0].vertices[:, 0] = np.clip(
        body.get_paths()[0].vertices[:, 0], -np.inf, m
    )
    body.set_facecolor("#2563eb")
    body.set_alpha(0.35)

# --- strip (rain) ---
for pos, d in zip(positions, data):
    jitter = rng.uniform(-0.08, 0.08, size=len(d))
    ax.scatter(
        pos + 0.15 + jitter, d, s=8, alpha=0.45, color="#1e40af", edgecolors="none"
    )

# --- boxplot (puddle) ---
bp = ax.boxplot(
    data,
    positions=positions,
    widths=0.12,
    patch_artist=True,
    showfliers=False,
    zorder=3,
)
for box in bp["boxes"]:
    box.set(facecolor="white", edgecolor="#1e40af", linewidth=1.2)
for med in bp["medians"]:
    med.set(color="#ef4444", linewidth=2)

ax.set_xticks(positions)
ax.set_xticklabels(groups)
ax.set_ylabel("スコア")
ax.set_title("レインクラウドプロット — グループ別スコア分布")
ax.grid(axis="y", alpha=0.2)

fig.tight_layout()
plt.show()