95% Bootstrap метод для медианы

Question

Я рассчитал мой Bootstrap метод для распределения δn∗ = Xn − xn для среднего μ двух наборов данных, которые у меня есть.

beer = [27, 19, 20, 20, 23, 17, 21, 24, 31, 26, 28, 20,27, 19, 25, 31, 24, 28, 24, 29, 21, 21, 18, 27, 20]
water = [21, 19, 13, 22, 15, 22, 15, 22, 20, 12, 24, 24,21, 19, 18, 16, 23, 20]

Теперь я должен сделать тот же метод, но для медианы и распределения δn∗ = MD∗n − mdn и сравните с предыдущим распределением. Вот код для моего среднего примера:

#beer
b = 10000
id = 1111111
np.random.seed(id)
x_star = np.random.choice(beer, (n, b))
x_star_bar = np.mean(x_star, axis=0)
x_bar = np.mean(beer)
delta_star = x_star_bar - x_bar
print('{}% bootstrap confidence interval for mu: ({:5.2f}, {:5.2f})'.format((1 - alpha) * 100,
      x_bar - np.quantile(delta_star, 1 - alpha/2),
      x_bar - np.quantile(delta_star, alpha/2)))
#water
b = 10000
id = 11817339
np.random.seed(id)
x_star = np.random.choice(water, (n, b))
x_star_bar = np.mean(x_star, axis=0)
x_bar = np.mean(water)
delta_star = x_star_bar - x_bar
print('{}% bootstrap confidence interval for mu: ({:5.2f}, {:5.2f})'.format((1 - alpha) * 100,
      x_bar - np.quantile(delta_star, 1 - alpha/2),
      x_bar - np.quantile(delta_star, alpha/2)))

Как мне теперь рассчитать медиану?