|
|
@@ -373,7 +373,6 @@ def get_parametric_boots_pars(lm, x, y, m, degree = 1, seed = 1977):
|
|
|
p = logit_poly_model(x, pars)
|
|
|
|
|
|
rng = np.random.default_rng(seed)
|
|
|
-
|
|
|
n = len(x)
|
|
|
|
|
|
# generate parameters
|
|
|
@@ -381,7 +380,7 @@ def get_parametric_boots_pars(lm, x, y, m, degree = 1, seed = 1977):
|
|
|
while True:
|
|
|
|
|
|
# Generate new binary outcomes from Bernoulli(p_i)
|
|
|
- y_sim = np.random.binomial(n = 1, p = p)
|
|
|
+ y_sim = rng.binomial(n = 1, p = p)
|
|
|
|
|
|
if np.sum(y_sim) in [0, n]: continue
|
|
|
|
