【发布时间】:2012-01-26 11:40:57
【问题描述】:
我一直在使用 scipy.stats.invgamma.rvs 函数从逆伽马分布中随机选择值。
我正在努力寻找有关 Python/Scipy 用于生成这些值的方法的任何信息(例如,使用 CDF 或拒绝采样、从另一个分布转换等)。
有人知道这是如何进行的吗?非常感谢您的帮助!抱歉,如果我重复一个问题,但我无法在任何地方找到此信息!
【问题讨论】:
我一直在使用 scipy.stats.invgamma.rvs 函数从逆伽马分布中随机选择值。
我正在努力寻找有关 Python/Scipy 用于生成这些值的方法的任何信息(例如,使用 CDF 或拒绝采样、从另一个分布转换等)。
有人知道这是如何进行的吗?非常感谢您的帮助!抱歉,如果我重复一个问题,但我无法在任何地方找到此信息!
【问题讨论】:
Scipy source code 说:“默认使用基本的逆 cdf 算法生成 RV。”:
def _rvs(self, *args):
## Use basic inverse cdf algorithm for RV generation as default.
U = mtrand.sample(self._size)
Y = self._ppf(U,*args)
return Y
【讨论】: