作者luckyb (New * for U)
站内Statistics
标题[问题] Metropolis-Hastings
时间Thu Dec 25 15:35:48 2008
想请问一下各位大大有关於MH的问题
就是如果我现在要生成A和B两个参数
ex--------------
A[1]<-0.6
B[1]<-0.27
for(k in 2:q){
a<-rtnorm(1,A[q-1],sqrt(0.01),0,1-B[q-1])
接受率a
b<-rtnorm(1,B[q-1],sqrt(0.01),0,1-A[q])
接受率b}
--------------------
其中因为是用truncated normal当proposal density
且A跟B相加要小於1
那算接受率的时候范围的地方应该要怎麽办呢?
谢谢大家
& Merry X'mas
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.115.45.158
※ 编辑: luckyb 来自: 140.115.45.158 (12/25 15:39)
1F：→ clickhere:你的 a, b 是 A[k], B[k] 吧? 12/26 03:29
2F：→ clickhere:如果不是, 那 A[2:q] 和 B[2:q] 为何? 12/26 03:30
3F：→ clickhere:如果是, 那得改用 repeat 去取代 for 12/26 03:31
4F：→ clickhere:MH要的是正确的sample,而非全部的sample. 12/26 03:32
5F：→ clickhere:用 A[k]+B[k] <= 1 && A[k]+B[k] >=0 判断是否 k <- k+1 12/26 03:33