将变量转换为R中的伽玛分布

Question

我有一个要转换为具有已知形状和速率参数的伽马分布的变量。 如何在R中将变量转换为伽马分布？ 我已经研究了dgamma，pgamma和qgamma函数，但无法确定是否有任何功能可以满足我的要求。

这是一个小例子：

variable <- rnorm(100)
shape <- .83
rate <- .01

注意：我意识到本示例使用正态分布的数据（不适合伽玛分布），但是我需要将变量重新缩放为原始伽玛分布。

Answer 1

使用分布和分位数功能来翻译：

qgamma(pnorm(variable), shape=.83, rate=.01)

假设variable平均值为0，标准差为1（如您的示例所示）。 否则，您可以传递均值并将sd转换为pnorm 。

要查看转换：

plot(density(variable))

plot(density(qgamma(pnorm(variable), shape=.83, rate=.01)))