[0,1] の間の確率変数に対して、正規化されていないノンパラメトリック分布が 2 つあるとします。
unnormalized_pdf_A = abs(sin(linspace(1,10,100)))
と
unnormalized_pdf_B = abs(cos(linspace(1,10,100)))
C = A*B変数の pdf を取得し、その基本的な統計 (平均、分散、ベイジアン間隔など)を取得したいと考えています。
これはstatsmodels、PyMCまたはで実行できscipy.statsますか? どうすればこの問題に取り組み始めることができますか?