可以使用来自两个分布的抽取来从分布中抽取其产品的密度吗?

机器算法验证 分布 模拟 蒙特卡洛 随机生成
2022-04-10 03:17:26

假设我有一组来自密度为的分布的绘图和一组来自密度为的分布的绘图,其中是未知的。有没有办法使用这些绘制来获得一组(可能更小)绘制,其 pdf 与这两个分布的 pdf 的乘积成正比?fgfg

我知道可以通过内核平滑来估计两个 pdf 的近似分布,但我想知道是否有一种方法可以得到完全正确分布的绘图。

1个回答

从纯模拟的角度来看(警告:链接到我的书),旨在从产品模拟f×g两个(正)函数,使得fg在适当的区域上是可积的,这表明了一些特定的方法:

  1. Sampling-importance-resampling:当来自两者的样本时应用fg可以产生并且当两者的归一化版本的对应数值fg在这些点可以进行数值计算(重新采样对于计算与相关的积分不是必需的fg)。
  2. Accept-Reject : 申请时f或者g与标准概率密度成正比,乘积中的其他项是有界的。
  3. 切片采样:适用于f或者g与标准概率密度成正比,易于模拟,另一个函数使得集合{x;f(x)a}易于构建。或者当两者fg是这样的{x;f(x)a}易于构建。
  4. Metropolis-Hastings 算法:适用于f或者g与标准概率密度成正比,易于模拟,并且可以(数字)计算其他函数。
  5. 桥接采样:两者都适用fg与标准概率密度成正比,当来自两者的样本fg可以在合理的时间内产生,并且当两者的对应数值fg在这些点可以(数字)计算。
  6. 嵌套采样:适用于f或者g与标准概率密度成正比,易于模拟,并且可以(数字)计算其他函数。

这些不同的解决方案不一定考虑到以下事实:f是先验密度和g是一种可能性。实际上,它不一定是产品的最佳分解fg, 因为f可能太平甚至不合适,而g不一定可积在参数中。

注意:去年在交叉验证上提出了类似的问题,但没有吸引任何流量。我在该页面上添加了指向可扩展蒙特卡洛解决方案的链接 。