garch estimation on impossibly long series

[This article was first published on Portfolio Probe » R language, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

The variability of garch estimates when the series has 100,000 returns.


The post “Variability of garch estimates” showed estimates of 1000 series that were each 2000 observations long.  Here we do the same thing except that the series each have 100,000 observations.

That would be four centuries of daily data.  It’s not presently feasible to think of  a market mechanism being stable for that long.

For 1-minute returns, that is only a year of data.  Hence we could actually gather such data and feed it to a garch estimator.  The problem here is that there is significant seasonality in volatility throughout the trading day.  So to get viable results we would need a more complicated model than a garch(1,1) with t-distributed errors.


The figures show the distributions of the estimates of:

  1. alpha and beta
  2. the half-life
  3. the degrees of freedom of the t distribution
  4. the asymptotic variance

Figure 1: Smoothed scatterplot of the alpha and beta estimates.

Figure 2: Distribution of the estimated half-life.

Figure 3: Distribution of estimated degrees of freedom.

Figure 4: Distribution of estimated asymptotic variance.

The variability in the estimate of the asymptotic variance is perhaps surprisingly high.

Appendix R

The command (see “Variability of garch estimates”  for more details) to do the estimation was:

> system.time(ges.a.100K.07 <- pp.garchEstSim(c(.01, .07, .925),
+        spec=tspec, nobs=1e5, df=7, trials=1000))
    user   system  elapsed 
28445.89  2634.36 32473.76

That is 9 hours of elapsed time.

To leave a comment for the author, please follow the link and comment on their blog: Portfolio Probe » R language. offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)