[R-SIG-Finance] [R-sig-finance] Garch problem

Patrick Burns patrick at burns-stat.com
Wed Mar 18 10:19:20 CET 2009


I was hoping to leave the "Doubtless" [1] as
an exercise for the reader -- mainly as I'm not
at all well versed in what is available in R for
garch these days.

One idea would be to try a components model
(may not be available).

Another idea would be to try a Bayesian estimate
(may not be available).

A method that certainly is available is to pick a
"reasonable" set of parameters (no estimation).

The course of action may well depend on the use
to which the model is to be put.

[1] Stephen Crane "The Wayfarer"

Pat


RON70 wrote:
> Dear Patrick, thank you so much for this reply. You said one solution is to
> increase the data point. However at this point I can not get more. Therefore
> if you please tell more about "doubtless other paths" I will be truly
> grateful.
>
> Regards,
>
>
> Patrick Burns-2 wrote:
>   
>> The fit is essentially saying that the half-life
>> of a shock is infinite.  This generally occurs
>> when the in-sample volatility has a general
>> trend.  One solution is more data.  There are
>> doubtless other paths as well.
>>
>> RON70 wrote:
>>     
>>> I have following dataset as monthly percentage return for a stock :
>>>
>>> 0.173741362
>>> -0.062237174
>>>
>>>   
>>>       
>> [ ... ]
>>     
>>> -0.001652893
>>> -0.092301325
>>>
>>> Now I fit a GARCH (1,1) model on that :
>>>
>>>   
>>>       
>>>> garch(Delt(dat)[-1], c(1,1))
>>>>     
>>>>         
>>>  ***** ESTIMATION WITH ANALYTICAL GRADIENT ***** 
>>>
>>>
>>>      I     INITIAL X(I)        D(I)
>>>
>>>      1     4.331103e-03     1.000e+00
>>>      2     5.000000e-02     1.000e+00
>>>      3     5.000000e-02     1.000e+00
>>>
>>>     IT   NF      F         RELDF    PRELDF    RELDX   STPPAR   D*STEP  
>>> NPRELDF
>>>      0    1 -4.507e+02
>>>      1    6 -4.508e+02  2.00e-04  3.20e-04  1.5e-03  6.3e+06  1.5e-04 
>>> 1.01e+03
>>>      2    7 -4.508e+02  1.57e-05  1.69e-05  1.4e-03  2.0e+00  1.5e-04 
>>> 3.19e-01
>>>      3   13 -4.521e+02  2.85e-03  4.72e-03  5.6e-01  2.0e+00  1.3e-01 
>>> 3.16e-01
>>>      4   16 -4.602e+02  1.76e-02  4.41e-03  8.1e-01  6.7e-01  5.1e-01 
>>> 1.99e-02
>>>      5   23 -4.607e+02  1.13e-03  2.77e-03  1.6e-04  7.4e+00  1.8e-04 
>>> 8.48e+00
>>>      6   24 -4.607e+02  4.81e-05  4.37e-05  1.6e-04  2.0e+00  1.8e-04 
>>> 1.77e+01
>>>      7   30 -4.638e+02  6.60e-03  8.81e-03  9.8e-02  2.0e+00  1.2e-01 
>>> 1.84e+01
>>>      8   31 -4.645e+02  1.52e-03  7.73e-03  8.2e-02  1.3e+00  1.2e-01 
>>> 1.39e-02
>>>      9   33 -4.688e+02  9.18e-03  6.28e-03  6.8e-02  0.0e+00  1.2e-01 
>>> 6.94e-03
>>>     10   35 -4.693e+02  9.32e-04  9.33e-04  8.9e-03  1.9e+00  1.8e-02 
>>> 2.86e-02
>>>     11   37 -4.699e+02  1.34e-03  1.59e-03  1.6e-02  1.8e+00  3.5e-02 
>>> 5.99e-02
>>>     12   38 -4.704e+02  1.05e-03  1.43e-03  1.6e-02  1.6e+00  3.5e-02 
>>> 9.10e-03
>>>     13   40 -4.705e+02  1.84e-04  2.85e-04  5.3e-03  1.2e+00  1.3e-02 
>>> 7.52e-04
>>>     14   42 -4.705e+02  3.71e-05  5.18e-05  2.4e-03  8.1e-01  5.0e-03 
>>> 7.09e-05
>>>     15   44 -4.705e+02  8.51e-07  3.04e-06  4.9e-04  8.2e-01  9.5e-04 
>>> 5.29e-06
>>>     16   57 -4.705e+02 -7.73e-15  1.09e-15  5.0e-15  4.4e+06  9.1e-15 
>>> 2.87e-07
>>>
>>>  ***** FALSE CONVERGENCE *****
>>>
>>>  FUNCTION    -4.704848e+02   RELDX        4.961e-15
>>>  FUNC. EVALS      57         GRAD. EVALS      16
>>>  PRELDF       1.088e-15      NPRELDF      2.867e-07
>>>
>>>      I      FINAL X(I)        D(I)          G(I)
>>>
>>>      1    2.824235e-05     1.000e+00     5.619e+01
>>>      2    8.649332e-02     1.000e+00    -5.899e-01
>>>      3    9.175397e-01     1.000e+00    -6.866e-01
>>>
>>>
>>> Call:
>>> garch(x = Delt(dat)[-1], order = c(1, 1))
>>>
>>> Coefficient(s):
>>>        a0         a1         b1  
>>> 2.824e-05  8.649e-02  9.175e-01  
>>>
>>> Warning message:
>>> In sqrt(pred$e) : NaNs produced
>>>
>>> What we see that sum of alpha and beta coef is more than 1. Therefore
>>> probably I choose a wrong model on my dataset. Can anyone please guide me
>>> how to modify that model?
>>>
>>> Regards,
>>>
>>>       
>> _______________________________________________
>> R-SIG-Finance at stat.math.ethz.ch mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-sig-finance
>> -- Subscriber-posting only.
>> -- If you want to post, subscribe first.
>>
>>
>>     
>
>



More information about the R-SIG-Finance mailing list