Самозагрузка биномиального распределения

Question

Я работаю над простой проблемой, я хочу написать функцию, которая будет случайным образом рисовать n выборок размера q из биномиального распределения (скажем, 60% 1 и 40% 0), а затем сохранять долю 1 из каждый образец в массиве (чтобы я мог проанализировать этот массив позже). Например, я хотел бы смоделировать, если бы я нарисовал 1000 выборок каждого размера 30 из распределения с 60% 1 с и 40% 0 с, и я хотел бы в качестве вывода получить массив длиной 1000, где каждая строка представляет долю 1 из этого образца.

Answer 1

Вы имеете в виду распределение Бернулли?

from scipy.stats import bernoulli
data = [bernoulli.rvs(0.6, size=30) for _ in range(1000)]

Answer 2

Чтобы получить число 1 с, это будет

import numpy as np
sim = np.random.binomial(n=30,p=0.6,size=1000)

Для фактического получения доктора aws:

sim = np.random.binomial(n=1,p=0.6,size=(1000,30))
sim.shape
import matplotlib.pyplot as plt
plt.hist(np.mean(sim,axis=1))