diff --git a/src/pop_models/powerlaw/joint_dist_plot.py b/src/pop_models/powerlaw/joint_dist_plot.py
index 1e2cf3a026613d2f022c18241af19016cd849c29..295dfbb9c5718b3648a13b23fcfcd72e34eda915 100644
--- a/src/pop_models/powerlaw/joint_dist_plot.py
+++ b/src/pop_models/powerlaw/joint_dist_plot.py
@@ -140,7 +140,7 @@ def _main(raw_args=None):
         log10_rate_smooth_samples = log10_rate_grid.ravel()
 
         kde = utils.gaussian_kde(
-            [log10_rates, alphas],
+            [alphas, log10_rates],
             bw_method="scott", weights=weights,
         )
         dP_dalpha_dlog10_rate_grid = kde(