R语言用极大似然估计方法估计GARCH模型参数报错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我想估计GARCH模型的一个变型模型的参数。用Optim函数进行极大似然估计，会报错“L-BFGS-B needs finite values of 'fn'”，我自己通过browser（），发现其中第二个参数取负值时会出现inf值，但是我明明设置了下界lower的。不知道咋办了，求助大神们！！
这是我的代码:
r_t = Garch_data$return
rate = Garch_data$rf
T = nrow(Garch_data)
sigma2_t = matrix(rep(NA,T*1),nrow=T)
z_t = matrix(rep(NA,T*1),nrow=T)
sigma2_t[1] = r_t[1]^2

loglike = function(param){
        lamda = param[1]
        omiga = param[2]
        b = param[3]
        a =param[4]
        c = param[5]
       
        z_t[1] = (r_t[1]-rate[1]-(lamda-0.5)*sigma2_t[1])/sqrt(sigma2_t[1])
        term = log(sigma2_t[1]) + z_t[1]^2
        for (t in 2:T){
                sigma2_t[t] = omiga + b*sigma2_t[t-1] + a*(z_t[t-1]-c*sqrt(sigma2_t[t-1]))^2
                z_t[t] = (r_t[t]-rate[1]-(lamda-0.5)*sigma2_t[t])/sqrt(sigma2_t[t])
                term = term + log(sigma2_t[t]) + z_t[t]^2
        }
        smallest = -1/2*term
        return(-smallest)       
}
x = c(1,0.01,0.1,0.1)
result = optim(x,loglike,method = "L-BFGS-B",lower=c(0,0.001,0.001,0.001),hessian = T)



下面是报错信息：

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：GARCH模型 极大似然估计 ARCH模型 GARCH 似然估计

第一次用这个论坛，不知道发帖格式是否正确，求大家帮助啊

发现其中第二个参数取负值时会出现inf值，但是我明明设置了下界lower的

没看错的话，需要优化的函数有5个参数，初始值和lower bounds都只给了4个。

zhouhao211314 发表于 2018-1-12 17:29
没看错的话，需要优化的函数有5个参数，初始值和lower bounds都只给了4个。

哎呀，我贴错图片了。四个参数是我在查找bug时做的一种尝试。原始的代码和错误是下面这张图。第二个参数我设置了下界为0.001，但是却取到-0.001，而且这时我的程序报错了。

zhouhao211314 发表于 2018-1-12 17:29
没看错的话，需要优化的函数有5个参数，初始值和lower bounds都只给了4个。

欲宇内 发表于 2018-1-12 16:42
第一次用这个论坛，不知道发帖格式是否正确，求大家帮助啊

你好！请问你这个问题最后怎么解决的？

zhouhao211314 发表于 2018-1-12 17:29
没看错的话，需要优化的函数有5个参数，初始值和lower bounds都只给了4个。

你好！请问你这个问题最后怎么解决的？

哪位大侠好心帮帮忙吧！在线等

问题已解决！错在初值上。