Fixed a bug in the PPD sampler

src/pop_models/ppd.py

@@ -31,7 +31,7 @@ def sample_ppd(
     # If ``post_weights`` is None, this is just 1/n_post, otherwise it will be
     # proportional to ``post_weights``, but properly normalized to 1.0.
     if post_weights is None:
-        weights = numpy.empty(n_samples, dtype=post_samples.dtype)
+        weights = numpy.empty(n_post, dtype=post_samples.dtype)
         weights[:] = 1.0 / n_post
     else:
         weights = post_weights / numpy.sum(post_weights)