r - 如果 Y = log(X) 和 X ~ Gamma 在 R 中模拟 Y

Question

这个问题与这里和这里的帖子密切相关。

假设随机变量 X = Gamma(k, theta)，我对模拟 Y 感兴趣，其中 Y = log(X)。

您能否建议允许这样做的 R 包（或建议 R 代码）？非常感谢您！

Answer 1

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

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

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