等效的MATLAB函数gaminv()

Question

我一直在努力寻找一个与MATLAB的R等价物 gaminv 函数。 据我所知，在R中没有任何一个packagefunctions可以完全做`gaminv的工作。 我也不知道自己该如何编码。 任何指导意见都将是非常感激的。

Answer 1

这将有助于给我们一个 文档链接:

x = gaminv(p,a,b) 返回形状参数a和比例参数b的伽马分布的icdf[逆累积分布函数]，在p中的值进行评估。

反累积分布函数也被称为反累积分布函数。分量函数：这些函数被称为 q<distname> 在R中，特别是： qgamma 是Gamma分布的量子函数：定义从 ?qgamma 州

qgamma(p, shape, rate = 1, scale = 1rate, lower.tail = TRUE, log.p = FALSE)

特别要注意的是，R函数默认是以速率为参数的：如果你指定了 qgamma(0.5, shape=2, 2) 你将得到速率为2的函数值（比例为12）。 如果你想得到比例参数化，你需要 qgamma(0.5, shape=2, scale=2).

我已经用一个任意的例子证实了这一点。gaminv(0.5,2,2) (在Octave中，因为我没有Matlab)给出的答案与 qgamma(0.5,2,scale=2).

Answer 2

你可以尝试从 反累积函数.如果这不是你要找的，我建议先研究一下累积反演的整个理论基础，并尝试在非函数式语言中应用，然后将其整合到R中，或者干脆读一读MatLab原函数中写的内容，并尝试在R中重现。