quasi MLE standard errors
Moderators: EViews Gareth, EViews Jason, EViews Moderator, EViews Matt
quasi MLE standard errors
Hey. I'm using Eviews7.
I am estimating a GARCH model using a likelihood approach, i.e. using a LoglObject.
I want to construct robust tstats by using the BollerslevWoolridge estimator for the standard errors. Since I didn't find it already implemented, I'd like to calculate those "by hand".
Herefore, I need an estimator of the Hessian of the likelihood function and the outer product of the gradient vector with itself, both evaluated at the parameter estimates.
I'm not sure in what way these are readily available in the Logl Object.
I'd be very grateful for any hints...
I am estimating a GARCH model using a likelihood approach, i.e. using a LoglObject.
I want to construct robust tstats by using the BollerslevWoolridge estimator for the standard errors. Since I didn't find it already implemented, I'd like to calculate those "by hand".
Herefore, I need an estimator of the Hessian of the likelihood function and the outer product of the gradient vector with itself, both evaluated at the parameter estimates.
I'm not sure in what way these are readily available in the Logl Object.
I'd be very grateful for any hints...

 EViews Developer
 Posts: 2620
 Joined: Wed Oct 15, 2008 9:17 am
Re: quasi MLE standard errors
The gradients are available using the Make Gradients proc. You can then use the stom command to save them in a matrix. Alternately, you can get the OPG by inverting the default coefficient covariance matrix.
Hessians are not available in the Logl. If I recall correctly, however, the BollerslevWooldrige covariances don't use the Hessian as the estimate of the expected Information. The BollerslevWooldrige [edit]estimate of the Information matrix can be computed from the gradients, but I don't recall off of the top of my head the formula. It's certainly described in the original paper
Hessians are not available in the Logl. If I recall correctly, however, the BollerslevWooldrige covariances don't use the Hessian as the estimate of the expected Information. The BollerslevWooldrige [edit]estimate of the Information matrix can be computed from the gradients, but I don't recall off of the top of my head the formula. It's certainly described in the original paper
Last edited by EViews Glenn on Thu Apr 14, 2011 10:23 am, edited 1 time in total.
Re: quasi MLE standard errors
Thank you.
Well, the formula I found for calculating the QMLE covariance estimator is: T^{1} ( I_1 * I_2^{1} * I_1 )^{1},
where I_1 is the outer product of the gradient vector with itself,
and I_2 is the negative Hessin Matrix evaluated at the parameters.
Is there any way to compute the Hessian?
Well, the formula I found for calculating the QMLE covariance estimator is: T^{1} ( I_1 * I_2^{1} * I_1 )^{1},
where I_1 is the outer product of the gradient vector with itself,
and I_2 is the negative Hessin Matrix evaluated at the parameters.
Is there any way to compute the Hessian?

 EViews Developer
 Posts: 2620
 Joined: Wed Oct 15, 2008 9:17 am
Re: quasi MLE standard errors
That's the standard form for QMLE but not the one typically used in ARCH.
Computing the Hessian using the LogL is not going to be easy. There's nothing builtin. You might be able to do it by putting analytic equations for the Hessian into the LogL or by putting equations into the LogL that allow you to do the numeric derivatives. But that's about it...
Computing the Hessian using the LogL is not going to be easy. There's nothing builtin. You might be able to do it by putting analytic equations for the Hessian into the LogL or by putting equations into the LogL that allow you to do the numeric derivatives. But that's about it...
Re: quasi MLE standard errors
Thank you for your replies.
I didn't look at the original paper yet, but I found the formula in hamiltons time series analysis.
Could you tell me, how EViews calculates the robust standard errors when estimating GARCH in an equation object? That would help me a lot.
Basically, I just want to replicate the computation of robust errors using my Logl Object. I didn't find an explicit formula in the User Guides.
I didn't look at the original paper yet, but I found the formula in hamiltons time series analysis.
Could you tell me, how EViews calculates the robust standard errors when estimating GARCH in an equation object? That would help me a lot.
Basically, I just want to replicate the computation of robust errors using my Logl Object. I didn't find an explicit formula in the User Guides.
Re: quasi MLE standard errors
Hello.
I had a look at the B&W paper. They derive an expresson for an estimator for the B&W errors which requires no second derivatives. Could you tell me if Eviews internally uses this formula with first derivatives normally calculate the B&W errors, that is when estimating an GARCH equation? Otherwise, please tell me if I should post this question to another forum.
I'd really appreciate your help.
I had a look at the B&W paper. They derive an expresson for an estimator for the B&W errors which requires no second derivatives. Could you tell me if Eviews internally uses this formula with first derivatives normally calculate the B&W errors, that is when estimating an GARCH equation? Otherwise, please tell me if I should post this question to another forum.
I'd really appreciate your help.

 EViews Developer
 Posts: 2620
 Joined: Wed Oct 15, 2008 9:17 am
Re: quasi MLE standard errors
Yes. EViews computes BollerslevWooldridge standard errors as described in the paper...
Re: quasi MLE standard errors
Unfortunately, I am still struggeling with computing robust standard errors.
The formula given in B&W uses first derivatives of the conditional variance wrt to the parameters. In my case, analytical derivatives are very tricky to program, so I was wondering whether there is anything builtin,
to get numerical derivatives?
Does anyone have any suggestins how to compute robust errors when working with a Logl object?
I'm sure I am not the first one esitmating a GARCH model with likelihood methods in Eviews and needing robust errors...
The formula given in B&W uses first derivatives of the conditional variance wrt to the parameters. In my case, analytical derivatives are very tricky to program, so I was wondering whether there is anything builtin,
to get numerical derivatives?
Does anyone have any suggestins how to compute robust errors when working with a Logl object?
I'm sure I am not the first one esitmating a GARCH model with likelihood methods in Eviews and needing robust errors...

 EViews Developer
 Posts: 2620
 Joined: Wed Oct 15, 2008 9:17 am
Re: quasi MLE standard errors
The Make Gradients proc will return the matrix of first derivatives.
Re: quasi MLE standard errors
thats for the gradients of the likelihood, right?
However, the formula to calculate the B&W errors uses also gradients of the conditional variance process with respect to the parameters.
If Eviews internally uses this formula, how does it get those derivatives?
However, the formula to calculate the B&W errors uses also gradients of the conditional variance process with respect to the parameters.
If Eviews internally uses this formula, how does it get those derivatives?

 EViews Developer
 Posts: 2620
 Joined: Wed Oct 15, 2008 9:17 am
Re: quasi MLE standard errors
Those are for the likelihood with respect to the parameters.
Numeric derivatives. Looks like you'll have to take your expressions which compute the conditional variance and use them to compute numeric derivatives.
Numeric derivatives. Looks like you'll have to take your expressions which compute the conditional variance and use them to compute numeric derivatives.
Re: quasi MLE standard errors
How can I get numeric derivatives of my expression for the conditional variance? Is there any builtin function? I didn't find anything in the manual.

 EViews Developer
 Posts: 2620
 Joined: Wed Oct 15, 2008 9:17 am
Re: quasi MLE standard errors
What you are doing with the loglikelihood object is to produce a custom estimator. The calculation of the conditional variance is specific to the model you are trying to estimate. That's why there's nothing in the manual... Note also that what you want to do isn't really in the realm of what the object is currently designed to support. I'm afraid that you are not going to find anything specific to what you want to calculate in our manual.
There will, however, be general discussions that will be of use. For example, bear in mind that the expressions in your logl are series expressions. They can be evaluated inside of the logl or outside of the logl. Putting them in the logl is a convenient way of evaulating for each update of the parameters.
Since you are computing an arch model you already have equations for the conditional variance. Given that, what you can do is to put additional equations in your logl which evaluate the conditional variance at the coefficients plus a small value and (if you wish) minus a small value of what you are trying to differentiate with respect to. This is going to be messy, but can be done. Then put additional equations for the derivatives based on these small changes...
Or you could do it outside the logl using the series expressions.
There will, however, be general discussions that will be of use. For example, bear in mind that the expressions in your logl are series expressions. They can be evaluated inside of the logl or outside of the logl. Putting them in the logl is a convenient way of evaulating for each update of the parameters.
Since you are computing an arch model you already have equations for the conditional variance. Given that, what you can do is to put additional equations in your logl which evaluate the conditional variance at the coefficients plus a small value and (if you wish) minus a small value of what you are trying to differentiate with respect to. This is going to be messy, but can be done. Then put additional equations for the derivatives based on these small changes...
Or you could do it outside the logl using the series expressions.
Who is online
Users browsing this forum: No registered users and 6 guests