如果Y = log（X）并且R中的X~Gamma，则模拟Y.

Question

此问题与帖子here以及here密切相关。

让我们说随机变量X = Gamma（k，theta），我感兴趣的是模拟Y，其中Y = log（X）。

你能否建议R-package允许这样做（或建议一个R代码）？非常感谢你提前！

Answer 1

如果您需要提取Y个样本，请尝试：

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

然后，例如，您可以hist(Y)查看分发。