我创建了一个代表统计模拟的类。其中一个方面是p值的分布,并且该类包含表征该分布的方法:
class Simulation:
...
def pdf_p(self, p):
...
def cdf_p(self, p):
...
def ppf_p(self, P):
...
def rvs_p(self, size):
...
我现在想将该分布也公开为scipy.stats
-style 分布对象。为此__init__
,我的类的方法包含一个语句
self.p = PValueDist(self.pdf_p, self.cdf_p, self.ppf_p, self.rvs_p)
其中类PValueDist
定义为
from scipy.stats import rv_continuous
class PValueDist (rv_continuous):
def __init__(self, pdf, cdf, ppf, rvs):
self._pdf = pdf
self._cdf = cdf
self._ppf = ppf
self._rvs = rvs
super().__init__(self)
这似乎可行,但我想知道这是否是正确或规范的方式?
当然,我也可以创建Simulation
一个子类rv_continuous
并重命名我的方法。然而,模拟包含几个不同的分布,用其中一个来识别模拟在语义上似乎不正确。