Scipy Stats rv_continuous 没有归一化常数答案

【问题标题】：Scipy Stats rv_continuous without normalization constantScipy Stats rv_continuous 没有归一化常数
【发布时间】：2016-04-29 05:18:26
【问题描述】：

我想在不知道归一化常数的分布中绘制随机偏差。该分布是形状未知且尺度 = 1 的 gamma 似然的共轭先验。 pdf 在这里的维基百科上给出https://en.wikipedia.org/wiki/Conjugate_prior#Table_of_conjugate_distributions as

此 pdf 的积分没有封闭形式。在这里使用 scipy.stats.rv_continuous 是一个选项吗？如果没有，还有什么可以推荐用于绘制偏离任意 pdf 的内容？

【问题讨论】：

【解决方案1】：

以数值方式进行归一化积分 (scipy.integrate.quad) 始终是一种选择。

【讨论】：

感谢您的提示！有点担心那个速度。根据维基百科，实际上不需要知道归一化常数来使用某些算法进行采样。我想避免以数字方式计算积分会更快。 en.wikipedia.org/wiki/Rejection_sampling
是的，这也是一个选项 :-)。不过你需要一个好的少校。