2

我正在尝试按照http://www.unstarched.net/2013/03/20/high-frequency-garch-the-multiplicative-component-garch-mcsgarch-model/中的说明来拟合乘法 garch 模型。虽然我使用 15 分钟的间隔 1 个月。当我运行最后一个命令时，我收到以下错误：“优化错误（init [mask]，armaCSS，方法 = 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' 中的初始值不是有限的

谢谢！

4

1 回答 1

2

我刚刚意识到我的错误在哪里。通过计算返回值，第一个观测值丢失了，因此，错误来自于第一个观测值中的 NaN。

于 2018-07-22T10:21:50.123 回答