带有参数的自定义PDF，并带有Scipy的支持

Question

我想使用RV_Continouul类构建自定义PDF；这样的PDF应取决于几个参数，l1，l2，l3在下面的示例中。这是我的尝试：

k=1
class estimated_pdf(scipy.stats.rv_continuous):


    def _get_support(self, l1, l2, l3):
         self.a=0
         self.b=np.inf
    

    def _pdf(self, x, l1, l2, l3):
        return np.exp(-l1*x/k-l2*x**2/k-l3*x**3/k)


    def Z(self, l1, l2, l3):
        return quad(self._pdf, 0, np.inf, args=(l1, l2, l3))

    # edited
    def my_freeze(self, l1, l2, l3):
        self.l1 = l1
        self.l2 = l2
        self.l3 = l3


p = estimated_pdf()
# p.__call__(l1=1, l2=2, l3=3)
p.my_freeze(l1=1, l2=2, l3=3)

我会遇到以下问题：

1. “ _get_support（）采用1个位置参数，但给出了4个。 我不确定如何处理它
2. 支持不应取决于参数，因此，如果我评论_get_support部分， ：“ _parse_args（）缺少3所必需的位置参数：'l1'，'l2'和''和' l3'“

我不知道发生了什么，1。不开心，因为我给它l1，l2，l3参数看起来2。不开心，因为我没有给它。

编辑：我意识到使用__调用__或freeze会将我的对象再次变为rv_continuul对象，而不是estectated_pdf对象，因此我添加了我自己的my_freeze 解决该问题的功能。 我还将l1，l2，l3参数添加到_Get_Support。目前，没有提出错误，但是通过检查支持，我知道这是完整的真实行，即dafault值，因此_get> _get_support仍然不起作用。

Answer 1

_get_support有望接收与PDF相同的形状参数，此处l1 l2，l3。

IIUC，如果支撑是实际的半轴，并且不依赖于L1，L2，L3，则您无需明确定义get_support。只是省略它。

PS：您的片段还有很多其他问题。 k在PDF中不确定； Z方法不会返回任何内容。

编辑：喜欢：

In [8]: from numpy import exp

In [9]: class my_distr_gen(rv_continuous):
   ...:     def _get_support(self, l1, l2, l3):
   ...:         return 0, float('inf')
   ...:     def _pdf(self, x, l1, l2, l3):
   ...:         return exp(-l1*x - l2*x**2 - l3*x**3)  # NB: normalize!
   ...: 

In [10]: my_distr = my_distr_gen()

In [11]: my_distr.pdf(1, l1=1, l2=2, l3=3)
Out[11]: 0.0024787521766663585

In [12]: my_distr.pdf(-1, l1=1, l2=2, l3=3)
Out[12]: 0.0