Выборка из экспоненциального Бернулли

Question

Бернулли — это распределение вероятностей. Мне нужно сделать выборку из экспоненциального бернулли и вернуть двоичное значение (т.е. либо 0, либо 1). Я нашел этот алгоритм экспоненциальной выборки Бернулли , и я хочу его реализовать, но я не понимаю шаг 3 алгоритма, где: r1 = r1 & (2^h - 1 ). Может ли кто-нибудь помочь?

Answer 1

Вы можете использовать библиотеку, которая реализует выборку из распределения Бернулли, например, np.random.binomial(поскольку биномиальное распределение с n = 1 является распределением Бернулли ).

      import numpy as np
np.random.binomial(n=1, p=.2, size=20)
# output: array([0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0])