假设我有两个随机变量:
X ~ Beta(α1,β1)
Y ~ Beta(α2,β2)
我想计算 Z = XY 的分布(随机变量的乘积)
使用scipy
,我可以通过以下方式获取单个 Beta 的 pdf:
from scipy.stats import beta
rv = beta(a, b)
x = np.linspace(start=0, stop=1, num=200)
my_pdf = rv.pdf(x)
但是两个 Beta 的乘积呢?我可以分析地做到这一点吗?(Python/Julia/R 解决方案很好)。