random_intercepts.png

def fig_random_intercepts(mixed) -> None:
    """Distribution of the estimated per-patient random intercepts (BLUPs)."""
    re_values = np.array([float(v.iloc[0]) for v in mixed.random_effects.values()])
    group_var = float(mixed.cov_re.iloc[0, 0])

    fig, ax = plt.subplots(figsize=(7, 4.2))
    ax.hist(re_values, bins=30, color=PRIMARY, alpha=0.85, edgecolor="white")
    ax.axvline(0, color=EVENT, lw=1.6, ls="--", label="Populationsmittel (0)")
    ax.set_xlabel("Individuelle Abweichung vom Populationsmittel (mmHg)")
    ax.set_ylabel("Anzahl Patient:innen")
    ax.set_title("Verteilung der geschätzten Random Intercepts\n"
                 f"(Between-Patient-Varianz = {group_var:.2f}, "
                 f"SD = {np.sqrt(group_var):.2f} mmHg)")
    ax.legend(loc="upper right")
    save(fig, ASSETS / "random_intercepts.png")