如果 R 中的 Y = log(X) 和 X ~ Gamma,则模拟 Y

simulate Y if Y = log(X) and X ~ Gamma in R

这个问题与帖子 here and also here 密切相关。

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

您能否推荐允许执行此操作的 R 包(或推荐 R 代码)?非常感谢您!

如果您需要提取 Y 个样本,只需尝试:

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

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