我正在尝试将一些Julia代码转换为Python。这是用于多项式分布的代码，我被困在它的最后一部分。我不知道如何用Python编写它，因为我想知道是否有可以满足我需要的软件包。我不知道我是否可以使用SciPy.stats来完成此操作，因为文档似乎有限。
在这里，您会发现Julia卡代码的一部分被卡住了：
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))