Fit a random-effects meta-regression with one moderator. Paule-Mandel τ² + Knapp-Hartung t(k-2) CIs (with max(1, …) floor per advanced-stats.md). Bubble plot with fitted line.
Studies
Bubble plot
τ² via Paule-Mandel bisection; CIs use Knapp-Hartung t(k−2) with q = max(1, RSS/(k−2)) floor. Export the R script to replicate with metafor::rma(yi, vi, mods=~x, method="PM", test="knha").