如果我有一个像这样生成的随机变量的分析pdf
y= (x**(d-1)*e**(-x/g))/(math.gamma(d)*g**d)
其中 d 和 g 是常数
如何使用蒙特卡罗模拟生成此 pdf
我用 x 代替解析方程并生成 PDF，但是对于模拟我不知道该怎么做
这是伽玛分布。
形状为 d，比例为 g。
取样
import numpy as np
shape=d
scale=g
s = np.random.gamma(shape, scale, 1000) # draw 1000 RVs from gamma distribution