“ vamin”中的初始值不确定

时间:2018-07-22 09:50:21

标签: r optimization quantitative-finance volatility

我正在按照http://www.unstarched.net/2013/03/20/high-frequency-garch-the-multiplicative-component-garch-mcsgarch-model/中的说明尝试拟合可乘garch模型。虽然我每个月使用15分钟间隔。当我运行最后一条命令时,出现以下错误:“ optim(init [mask],armaCSS,method = optim.method,hessian = TRUE,   'vmmin'中的初始值不是有限的。 “对于每日变化,我使用的样本始于2010年,否则我无法估算每日变化。艰难的日内数据,我只有一个月的数据。我需要估算110只股票的波动率指标。

我非常感谢您的帮助。我不知道该如何处理。

我提供了代码,并提供了上一条命令的错误。

 sub <- subset(ITCH_volat, ITCH_volat$tickerid == 2,                                     
 select=c(returnmidend, time))
 sub_t <- xts(sub$returnmidend, sub$time)
 C = quantmod::getSymbols('AAPL', from = '2010-03-02',auto.assign=FALSE)
 C = quantmod::adjustOHLC(C, use.Adjusted = TRUE)
 R_d = TTR::ROC(Cl(C), na.pad = FALSE)
 spec_d = ugarchspec(mean.model = list(armaOrder = c(1, 1)), variance.model   
  = list(model = 'eGARCH', garchOrder = c(2, 1)), distribution = 'nig')
 roll = ugarchroll(spec_d, data = R_d['/2013-03-28'], forecast.length = n,     
 refit.every = 5, refit.window = 'moving', moving.size = 5, calculate.VaR =    
 FALSE)# extract the sigma forecast
 df = as.data.frame(roll)
 f_sigma = as.xts(df[, 'Sigma', drop = FALSE])
 spec = ugarchspec(mean.model = list(armaOrder = c(1, 1), include.mean =               
 TRUE), variance.model = list(model = 'mcsGARCH'), distribution = 'nig')
 fit = ugarchfit(data = sub_t, spec = spec, DailyVar = f_sigma^2)

最后一条命令给我以下错误:

optim(init [mask],armaCSS,method = optim.method,hessian = TRUE,中的错误:   'vmmin'中的初始值不是有限的

谢谢!

1 个答案:

答案 0 :(得分:1)

我刚刚意识到我的错误在哪里。通过计算收益,就缺少了第一个观测值,因此,在第一个观测值中使用NaN会导致错误。