推 thomas83:同蛇 ,你來錯版了 ! 01/25 09:36
ㄟ 想說班上很多統研的 想問一個生亂數的問題
就是 我只是需要 從一個截尾分配抽樣而已
像是 Z~N(0,1) k為一常數
需要一些 random number from f(z|z>k)
那現在我只知道 Accept-Reject Method
一種就是 從 Z~N(0,1) 抽 一直抽到 z>k 為止
一種就是 X~ g(x) (如用 Exp(β,k) β是scale,k是location )
1. Generate X~ Exp(β,k) , U ~uniform[0,1]
2. if U < f(x)/cg(x)
return X
else
go to step 1
所以,對於截尾分配有常用的生亂數方法嗎?麻煩同學告知ㄧ下
像是 Slice Sampler, Gibbs Sampler可以做嗎 謝謝!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.128.124
※ 編輯: kaishi 來自: 140.113.128.124 (01/18 15:15)