r - Y = log(X) および X ~ R のガンマの場合、Y をシミュレートします。

Question

この質問は、こことここの投稿に非常に密接に関連しています。

確率変数 X = Gamma(k, theta) で、Y = log(X) のシミュレーションに興味があるとしましょう。

それを可能にするRパッケージを提案してもらえますか（またはRコードを提案してください）？事前にどうもありがとうございました！

Answer 1

サンプルを抽出する必要がある場合は、次Yを試してください。

#set the number of samples
n<-10000
Y<-log(rgamma(n,k,theta))

次に、たとえば、hist(Y)分布を見ることができます。