Я пытаюсь создать пользовательский pdf в python с помощью scipy stats.rv_continuous ().
Поскольку я пытаюсь узнать об этом классе, я пытаюсь воспроизвести нормальное распределение.
import numpy as np
import scipy.stats as sps
class my_pdf(sps.rv_continuous):
def _pdf(self,x, mu, sigma):
return np.exp(-((x-mu)/sigma)**2 / 2) / (sigma*np.sqrt(2 * 3.1416))
my_normal = my_pdf( name='custom_normal')
values = my_normal.rvs(mu=10, sigma=2, size=200)
С маленьким му у меня не возникнет никаких проблем. Однако для mu = 10, например, я получаю следующую ошибку: OverflowError: (34, 'Result too large')
Может ли кто-нибудь исправить мой код, чтобы я мог использовать любые mu и sigma, которые я хочу? Спасибо