Я пытаюсь перевести некоторый код Джулии на Python.Это код для полиномиального распределения, и я застрял в последней части этого.Я не знаю, как написать это на Python, потому что я хочу знать, есть ли пакет, который будет делать то, что я хочу.Я не знаю, смогу ли я сделать это, используя SciPy.stats, потому что документация кажется ограниченной.
Здесь вы найдете часть кода Джулии, где я застрял: x1s - это массив, x - это также массив, OrthoNNDist - это имя struct:
Base.length(d::OrthoNNDist) = length(d.x0)
Distributions.rand(d::OrthoNNDist) = rand(d.x1s)
Distributions.pdf(d::OrthoNNDist, x::Vector) = x in d.x1s ? d.prob : 0.0
Distributions.pdf(d::OrthoNNDist) = fill(d.prob, size(d.x1s))
Distributions.logpdf(d::OrthoNNDist, x::Vector) = log(pdf(d, x))