python - 一般分布的产品

Question

假设我对 [0,1] 之间的随机变量有两个未归一化的非参数分布，例如：

unnormalized_pdf_A = abs(sin(linspace(1,10,100)))

和

unnormalized_pdf_B = abs(cos(linspace(1,10,100)))

我想获得变量的 pdfC = A*B并获得它的一些基本统计数据（例如均值、方差和/或贝叶斯区间）

这对statsmodels,PyMC或是否可行scipy.stats？我怎样才能开始解决这个问题？