Python中(复杂)自定义pdf的极大似然估计

2024-03-28 16:30:39 发布

您现在位置:Python中文网/ 问答频道 /正文

我有一个概率密度函数(遗憾的是不直接),它依赖于四个参数p1,p2,p3,p4,并包含诸如erf,sqrt。。。在

有没有使用Python计算MLE估计的通用方法?我知道我可以继承stats.rv_连续(见下文)但在这种情况下,我需要重新定义fit()函数,而且我没有MLE估计的分析公式。。。在

如果可能的话,我宁愿避免基于优化的方法。在

class custom_distrib(stats.rv_continuous):
    def _init_(self, p1, p2, p3, p4):
        super(custom_distrib, self)._init_()
        self.p1 = p1
        self.p2 = p2
        self.p3 = p3
        self.p4 = p4

    def _pdf(self, x):
        res = compute_my_custom_distrib(x, self.p1, self.p2, self.p3, self.p4)
        return res

Tags: 方法selfinitdefstatscustomresp2