Mathematica有一个四参数广义逆伽马分布:
http://reference.wolfram.com/mathematica/ref/InverseGammaDistribution.html
并在该页面上提供PDF。有人在R?
中实现了密度,分布,分位数和采样函数我确实快速启动(PDF只是将该页面上的公式转换为R),但如果已经完成,我将不会为执行CDF和分位数函数而烦恼。
是否存在计算CDF(通过PDF的集成)和给定PDF的任何分布的分位数(通过CDF的反转)的一般函数?
[注意这是不广义逆高斯]
另请注意Mathematica页面上的“属性和关系”下拉列表,这似乎暗示它不是特殊情况或任何事物的概括(除了反伽玛)。
答案 0 :(得分:3)
我启动了一个包来实现这个:
https://github.com/barryrowlingson/geninvgamma
它只使用简单的反转和密度的集成,所以没有什么聪明的。目前来自分布的随机样本是通过生成U(0,1)并获得分位数来完成的,这似乎不是很有效或非常准确。
无论如何,这是一个开始。
答案 1 :(得分:1)
根据this vignette(附录C2),反伽马分布是广义双曲分布的一个特例,由ghyp
包实现。