r - R で変数をガンマ分布に変換する

Question

既知の形状およびレート パラメータを使用してガンマ分布に変換したい変数があります。Rで変数をガンマ分布に変換するにはどうすればよいですか? dgamma、pgamma、および qgamma 関数を調べましたが、希望どおりの機能があるかどうかはわかりません。

以下に小さな例を示します。

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

注: この例では正規分布データ (ガンマ分布に適合しない) を使用していることは理解していますが、変数を元のガンマ分布に再スケーリングする必要があります。

Answer 1

分布関数と分位関数を使用して翻訳します。

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

variableこれは、平均0、sd 1を想定しています（例の場合と同様）。それ以外の場合は、平均と標準偏差を に渡すことができますpnorm。

変換を確認するには:

plot(density(variable))

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