我有一个变量,我想转换为具有已知形状和速率参数的伽玛分布.如何将变量转换为R中的gamma分布?我已经研究了dgamma,pgamma和qgamma函数,但我不知道是否会有我想做的事情.
这是一个小例子:
variable <- rnorm(100) shape <- .83 rate <- .01
注意:我意识到这个例子使用的是正态分布式数据(不适合伽玛分布),但我需要将变量重新调整为原始的伽玛分布.
使用分布和分位数函数进行翻译:
qgamma(pnorm(variable), shape=.83, rate=.01)
假设它的variable
平均值为0,sd 1(就像你的例子一样).否则你可以将均值和sd传递给pnorm
.
要看转型:
plot(density(variable))
plot(density(qgamma(pnorm(variable), shape=.83, rate=.01)))