拒绝采样(Rejection sampling)

已知分布p(x),如何从该分布中采样样本呢?今天介绍一种称为拒绝采样的方法。具体如下:

假设已知一个容易采样的候选分布q(x)p(x)q(x)之间满足:

                        ​​​​​​​        Mq(x)\geq p(x),

M为某一个常数,上式的意思即选择一个候选分布q(x)和某一常数M,使得Mq(x)包含p(x)

接下来,拒绝采样的步骤可以总结如下:

1. 确定需采样的样本个数N;

2. 从候选分布q(x)中采样样本x^{(i)}

3. 从均匀分布U_{(0,Mq(x^{(i)}))}​​​​​​​中采样系数u

4. 计算\frac{p({x^{(i)}})}{Mq(x^{(i)})},若u<\frac{p({x^{(i)}})}{Mq(x^{(i)})},则接受当前采样样本x^{(i)},否则,拒绝该次采样;

5. 重复步骤2-4,直到采样样本数量等于N

  • 2
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值