How to add lagged values to rugarch-model
Lasse,
Your problem is that you have NAs in your external regressor matrix. Since
you lagged it twice, you have 2 NAs. Try the code below.
Thanks for the reproducible example, it made it really easy.
Bob
require(rugarch)
require(quantmod)
x = runif(n = 1000, min = -5, max = 100)
ext.reg_lagged = cbind(Lag(x, k = 1), Lag(x, k = 2))
ext.reg = cbind(x + runif(n = 1000, min = 5, max = 20), x - runif(n = 1000,
min = 5, max = 15))
*x <- x[-c(1:2)]*
*ext.reg_lagged = ext.reg_lagged[-c(1:2), ]*
fit.spec <- ugarchspec(variance.model = list(model = "sGARCH",
garchOrder = c(1, 1)),
mean.model = list(armaOrder = c(3, 1),
include.mean = TRUE,
external.regressors =
ext.reg_lagged), # with ext.reg.lagged it does
not work.
distribution.model = "sstd")
fit <- ugarchfit(data = x, spec = fit.spec, solver = "hybrid")
On Wed, Nov 5, 2014 at 9:16 AM, Lasse Thorst <thorstlasse at gmail.com> wrote:
Hi - First question here, so please bear with me:
I would like to added lagged values to a rugarch-model, but have been
unable to find answers in the vignette.
Basic test-code:
require(rugarch)
require(quantmod)
x = runif(n = 1000, min = -5, max = 100)
ext.reg_lagged = cbind(Lag(x, k = 1), Lag(x, k = 2))
ext.reg = cbind(x + runif(n = 1000, min = 5, max = 20), x - runif(n = 1000,
min = 5, max = 15))
fit.spec <- ugarchspec(variance.model = list(model = "sGARCH",
garchOrder = c(1, 1)),
mean.model = list(armaOrder = c(3, 1),
include.mean = TRUE,
external.regressors =
ext.reg), # with ext.reg.lagged it does not work.
distribution.model = "sstd")
fit <- ugarchfit(data = x, spec = fit.spec, solver = "hybrid")
Now this works perfectly, but if I exchange the external.regressors =
ext.reg_lagged
it no longer works.
How can you do this in rugarch? Or do I need to try another package?
/Regards Lasse
[[alternative HTML version deleted]]
_______________________________________________ R-SIG-Finance at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-finance -- Subscriber-posting only. If you want to post, subscribe first. -- Also note that this is not the r-help list where general R questions should go.