我正在按照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'中的初始值不是有限的
谢谢!
答案 0 :(得分:1)
我刚刚意识到我的错误在哪里。通过计算收益,就缺少了第一个观测值,因此,在第一个观测值中使用NaN会导致错误。