MH算法的编程问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

题目说要使用MH算法。

我的程序是：

1. f <- function(n = n, x = x){
   
2.   chain = matrix(NA, n, 1)
   
3.   chain[1,] = x
   
4.   for(i in 2:n){
   
5.     prop = sample(exp(-y),size = 1)
   
6.     logratio = -0.5*prop^2 + 0.5*x^2
   
7.     logalpha = min(0, logratio)
   
8.     if(log(runif(1)) < logalpha){
   
9.       x = prop
   
10.     }
    
11.     x = x
    
12.     chain[i,] = x
    
13.   }
    
14.   return(chain)
    
15. 
    
16. }

复制代码

我使用：
p <- f(n = 1000, x = 2)
但是不收敛。请问懂MH算法的朋友帮帮忙，看看哪出错了。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：程序

大家有没有会的呢，不要沉贴嘛