我有一个概率密度函数(遗憾的是不直接),它依赖于四个参数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
目前没有回答
相关问题 更多 >
编程相关推荐