julia - 如何在 julia0.6.2 中将 2 个 Beta 函数相乘？

Question

我有两个随机变量：

B1 ~ Beta(α1,β1)

B2 ~ Beta(α2,β2)

我想计算 Z = B1*B2 的分布（随机变量的乘积）

使用Julia，我怎样才能获得单个 Beta 的 pdf：

x = collect(0:0.001:1);
b1 = Beta(11,41);
b2 = Beta(41,11);
posterior = pdf.(b1.*b2,x);

我收到这个错误MethodError: no method matching *(::Distributions.Beta{Float64}, ::Distributions.Beta{Float64})