R中的广义逆Gamma分布

Mathematica具有四参数广义逆伽马分布：

并在该页面上提供其PDF。有人在R中实现了密度、分布、分位数和函数采样吗

我确实做了一个快速的开始（PDF只是那个页面上的方程式翻译成R），但是如果已经做了，我就不必费心实现CDF和分位数函数了

是否存在计算给定PDF的任何分布的CDF（通过PDF积分）和分位数（通过CDF反转）的通用函数

[注意，这不是广义逆高斯分布]

---

另请注意Mathematica页面上的“属性和关系”下拉列表，这似乎意味着它不是任何事物的特例或概括（除了逆伽马）。

根据（附录C2），逆伽马分布是广义双曲线分布的一种特殊情况，由

ghyp

软件包实现。

我启动了一个软件包来实现：

它只使用简单的密度反演和积分，所以并没有什么高明之处。目前，分布中的随机样本是通过生成U（0,1）并获得分位数来完成的，这似乎不是非常有效或非常准确

---

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

是的，但是广义逆伽马还有另外两个参数，我不知道它们如何适合这个框架。但是，相同分布的几个不同的参数化对清晰性没有帮助……好吧，您可以将这些参数转换为它们的固定值，可能只需将

ghyp

中的函数包装在一组新的d/q/p/r函数中。或者，您可以使用该软件包中的函数作为起点，以分析/编程方式删除固定的参数。问题不是删除固定的参数，而是添加可变的新参数……如果我们知道您需要哪种参数化，也许会有所帮助？@Spacedman啊，对不起。我的思路过于密集，没有抓住你的重点，尽管是斜体字。我现在明白为什么这不起作用了。