Page 1 of 1

GARCH-BEKK with initial values vector

Posted: Mon Apr 10, 2017 4:57 am
by myrthevandieijen
Hi,
I would like to estimate a GARCH BEKK without exogenous regressors and use those coefficient estimates as input for another GARCH BEKK estimation with exogenous regressors. I saw that there's an option 'initial', but I am struggling with the way to input previously estimated coefficients. I used %BETA, but the size of that vector is concerning me. It includes the estimates of the VAR(1) model coefficients I've estimated as mean model and the BEKK coefficients, but I would like to use the values in another model via 'initial' where I include dummies in the VAR(1) and external regressors in the GARCH BEKK. I think I would need to adjust the %BETA vector to include extra zeros in order to compute a vector of the same size of the new model I intend to estimate. How would I be able to do that?

I was thinking as well that I could estimate the VAR(1) with dummies separately (just as I did for the VAR(1) without dummies), and include those estimates in the vector and then add the %BETA coefficients of just the GARCH BEKK C, A and B matrices. How would I be able to do that?

I'm doing this because the GARCH model won't converge (even when adjusting the simplex iterations and BFGS iterations), so I want to change the initial values to help it converge.

It concerns this code (the code for the VAR(1) is standard):

Code: Select all

*****************************************************************************************
*
* Estimate GARCH-BEKK without dummies
* with VAR(1) without dummies
*
garch(p=1,q=1,model=var1model,mv=bekk,rvector=rr,hmatrices=hh,pmethod=simplex,piters=20,iters=1000,robusterrors)
*
* need variables with coefficients to use later for initial_values in next GARCH-BEKK
*
compute coefficients = %BETA
display coefficients
*
*****************************************************************************************
*
* Estimate GARCH-BEKK with external event regressors
* with VAR(1) without dummies
*
* use initial values from previous GARCH-BEKK without dummies ---> how??
*
garch(p=1,q=1,model=var1model,mv=bekk,rvector=rr_ev,hmatrices=hh_ev,initial = coefficients, iters=1000,robusterrors,xreg)
# prod_announce strat_moves fin_client neg_events
*
******************************************************************************************
* Estimate GARCH-BEKK with company related events dummies
* with VAR(1) with seasonality dummies
*
* use initial values from previous GARCH-BEKK without dummies ---> how??
*
garch(p=1,q=1,model=var1model_seas,mv=bekk,rvector=rr_ev,hmatrices=hh_ev,initial = coefficients, iters=1000,robusterrors,xreg)
# prod_announce strat_moves fin_client neg_events
*
*****************************************************************************************
Many thanks in advance for your help!

Best,

Myrthe

Re: GARCH-BEKK with initial values vector

Posted: Mon Apr 10, 2017 10:27 am
by TomDoan
Those are big, rather messy, parameter sets and rearranging the information isn't simple. Is there a reason you want to do that?

Re: GARCH-BEKK with initial values vector

Posted: Tue Apr 11, 2017 1:59 am
by myrthevandieijen
My co-authors and I want to investigate volatility spillovers on different type of data than stock returns. I know it's not the standard way multivariate GARCH models are used, but that's one of the novelties of our paper so I'm certain I want to do this, I really hope you can help. As changing the initial values is one of the standard recommendations you receive once the model doesn't converge, I had expected there to be more examples online, but I couldn't find any unfortunately. I hope you might be able to help with creating a vector of coefficients.
I know how to create a vector in various languages, but I seem to get stuck in RATS (still a newbie)..

Many thanks in advance!

Myrthe

Re: GARCH-BEKK with initial values vector

Posted: Thu Apr 13, 2017 3:14 am
by myrthevandieijen
Hi Tom,

Sorry to bother you again, but did you get a chance to look at the issue? I really need some help setting up this vector of initial values, would be great if you could give an example of how to code it.

Many thanks in advance!

Best,

Myrthe