R中的广义逆Gamma分布

时间:2012-06-29 11:09:22

标签: r distribution

Mathematica有一个四参数广义逆伽马分布:

http://reference.wolfram.com/mathematica/ref/InverseGammaDistribution.html

并在该页面上提供PDF。有人在R?

中实现了密度,分布,分位数和采样函数

我确实快速启动(PDF只是将该页面上的公式转换为R),但如果已经完成,我将不会为执行CDF和分位数函数而烦恼。

是否存在计算CDF(通过PDF的集成)和给定PDF的任何分布的分位数(通过CDF的反转)的一般函数?

[注意这是广义逆高斯]

另请注意Mathematica页面上的“属性和关系”下拉列表,这似乎暗示它不是特殊情况或任何事物的概括(除了反伽玛)。

2 个答案:

答案 0 :(得分:3)

我启动了一个包来实现这个:

https://github.com/barryrowlingson/geninvgamma

它只使用简单的反转和密度的集成,所以没有什么聪明的。目前来自分布的随机样本是通过生成U(0,1)并获得分位数来完成的,这似乎不是很有效或非常准确。

无论如何,这是一个开始。

答案 1 :(得分:1)

根据this vignette(附录C2),反伽马分布是广义双曲分布的一个特例,由ghyp包实现。