tvgarch | R Documentation |
Quasi Maximum Likelihood (ML) estimation of a univariate multiplicative TV(s)-GARCH(p,q,r)-X model, where s indicates the number and the shape of the transition functions, r is the asymmetry order, p is the ARCH order, q is the GARCH order, and 'X' indicates that covariates can be included. Any transition variable, deterministic or stochastic, can be used to drive the transitions between the variance states. The TV long-term component introduces non-stationarity in the variance process, where the GARCH-X short-term component describes conditional heteroscedasticity. Maximization by parts leads to consistent and asymptotically normal estimates.
tvgarch(y, order.g = 1, order.h = c(1,1,0), xtv = NULL, xreg = NULL,
initial.values = list(), opt = 2, upper.speed = NULL, tvgarch = FALSE,
turbo = FALSE, trace = FALSE)
y |
numeric |
order.g |
|
order.h |
|
initial.values |
a list containing the initial parameter values passed on
to the optimisation routines (constrOptim for the TV component and
nlminb for the GARCH-X component). If list(), the default, then the
values are chosen automatically. TV component: |
xtv |
|
xreg |
numeric |
opt |
|
upper.speed |
|
tvgarch |
|
turbo |
|
trace |
|
An object of class 'tvgarch'.
Susana Campos-Martins
Cristina Amado and Timo Terasvirta (2013) Modelling volatility by variance decomposition, Journal of Econometrics 175, 142-153.
Cristina Amado and Timo Terasvirta (2014) Modelling changes in the unconditional variance of long stock return series, Journal of Empirical Finance 25, 15-35.
Susana Campos-Martins and Genaro Sucarrat (2024) Modeling Nonstationary Financial Volatility with the R Package tvgarch, Journal of Statistical Software 108, 1-38.
garchx
,
tvgarchSim
,
nlminb
,
constrOptim
set.seed(123)
## Simulate from a TV(1)-GARCH(1,1) model (default):
ySim <- tvgarchSim(n = 1500)
## Estimate a TV(1)-GARCH(1,1) model:
yEst <- tvgarch(y = ySim)
## Print estimation results:
print(yEst)
## Extract coefficients:
coef(yEst)
## Plot conditional volatilities:
plot(yEst)
## Extract log-likelihood:
logLik(yEst)
## Extract and store standardised residuals:
etaEst <- residuals(yEst)
## Generate predictions:
predict(yEst)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.