Генерация случайных значений с пользовательским PDF

Question

У меня большой набор данных (более 300000 значений), и я создал PDF для него со stats.gaussian_kde следующим образом:

kernel=scipy.stats.gaussian_kde(vals)
ind=np.linspace(np.min(vals), np.max(vals),1000)
PDE=kernel.evaluate(ind)

Теперь я хочу сгенерировать новые случайные значения из этого PDF-файла, но я не уверен, как это сделать. У кого-нибудь есть примеры?

Answer 1

Используйте метод resample .

vals = [1,2,3,4,5]
kernel = scipy.stats.gaussian_kde(vals)    
my_random_vals = kernel.resample(size=10)

print(my_random_vals)
[[ 2.18467151 -1.22069121  2.64074108  0.32268845  2.78022227  6.29483273
   1.94848359  4.79865993  1.30363424  3.71441989]]