Skip to main content
V-Lab

GARCH Model Documentation
Conceptual, Interactive + Mathematical

GARCH Model Fundamentals

Level 1Level 1 (Conceptual)
From Volatility Concepts to GARCH Implementation

Dynamic volatility and conditional heteroskedasticity represent fundamental market behaviors that require mathematical precision to capture effectively. The GARCH(1,1) model provides exactly this precision—a specific, mathematically robust framework for implementing these concepts in practice. While volatility modeling explains why we need sophisticated approaches, GARCH shows us how to build one that captures the key market behaviors: volatility clustering, mean reversion, and realistic shock persistence.

At its core, GARCH(1,1) distills market volatility dynamics into three fundamental parameters that any practitioner can estimate, interpret, and apply. This parsimony—capturing complex market behavior with just three numbers—represents one of GARCH's greatest strengths and explains its dominance in practical risk management applications across global financial institutions.

The Three Parameters That Define Market Behavior

Omega (ω) - Baseline Volatility Component: Beyond conceptual understanding of unconditional volatility, ω provides the specific mathematical foundation that ensures positive volatility in all market conditions. Think of ω as the market's "baseline component" — not the target volatility itself, but rather the irreducible minimum that prevents volatility from reaching zero. Unlike constant volatility models, GARCH allows actual volatility to deviate significantly from this baseline component, with long-run volatility converging to ω/(1-α-β) — the unconditional volatility level that emerges from the complete parameter interaction.

Alpha (α) - Shock Sensitivity: Given that markets react to new information, α quantifies exactly how much today's market movements affect tomorrow's expected volatility. A high α value (near 0.2) indicates that large price movements today will substantially increase tomorrow's volatility expectations, creating the sharp volatility spikes observed during crisis periods. A low α value (near 0.05) suggests markets that respond more gradually to shocks, maintaining steadier volatility patterns even during significant price movements.

Beta (β) - Volatility Persistence: Beyond recognizing volatility clustering as a stylized fact, β determines the mathematical half-life of volatility shocks. When β approaches 0.9, volatility increases tend to persist for weeks or months, creating the extended periods of market stress characteristic of financial crises. Lower β values (around 0.7) indicate volatility shocks dissipate more quickly, leading to more rapid transitions between calm and turbulent market conditions.

How GARCH Captures Stylized Facts Mechanistically

Three key stylized facts characterize financial markets universally. GARCH's mathematical structure directly generates each of these patterns through specific parameter interactions:

Volatility Clustering Implementation: The recursive structure σt2=ω+α·εt-12+β·σt-12ensures that high volatility today (large σt-12) directly increases tomorrow's expected volatility through theβ coefficient. Simultaneously, large price shocks today (large εt-12) boost tomorrow's volatility through theα coefficient. This dual mechanism creates the clustering patterns where "large movements follow large movements" that characterize real financial data.

Mean Reversion Guarantee: The constraint α+β<1 ensures that volatility shocks eventually dissipate, preventing volatility from exploding indefinitely. This mathematical requirement translates the stylized fact of mean reversion into a verifiable model property, giving practitioners confidence that GARCH forecasts will not generate unrealistic long-term volatility projections. Importantly, this constraint also determines the long-run volatility target: while ω anchors the equation, actual long-run volatility converges to ω1-α-β, which depends on all three parameters working together.

Fat Tails Generation: The multiplicative error structure εt=σt·zt, where zt is standard normal, creates return distributions with heavier tails than normal distributions. Even though the standardized innovations zt are normally distributed, the time-varying volatility σt generates unconditional return distributions with excess kurtosis, matching the empirical observation that extreme market movements occur more frequently than normal distributions predict.

From Engle and Bollerslev to Modern Practice

Within the broader evolution of volatility modeling development, the specific GARCH breakthrough came from Tim Bollerslev's 1986 insight that volatility itself could follow autoregressive patterns. Building on Robert Engle's ARCH framework, Bollerslev recognized that modeling volatility persistence ( the β component) was as important as modeling shock response (the α component). This insight transformed volatility modeling from a purely reactive framework (ARCH) to a predictive one (GARCH).

The elegance of the GARCH(1,1) specification—capturing sophisticated volatility dynamics with just three parameters—enabled its rapid adoption across financial institutions. Unlike more complex models that require extensive computational resources or specialized expertise, GARCH(1,1) can be estimated reliably using standard maximum likelihood techniques and implemented in production risk management systems with confidence.

This sophistication becomes particularly evident when contrasted with the limitations of overly simplistic approaches that dominated earlier risk management practice. Simple exponentially weighted moving averages, while computationally efficient, lack the structural separation between shock sensitivity and persistence that GARCH provides. This design flaw results in volatility forecasts that systematically underestimate risk during market transitions and overestimate it during subsequent recoveries—precisely when accurate risk assessment proves most critical.

The regulatory shift toward GARCH-based frameworks reflects these practical failures, as simpler approaches consistently generated procyclical capital requirements that amplified rather than dampened financial instability. This practical accessibility explains why GARCH remains the foundation for regulatory capital calculations under Basel III and forms the backbone of most commercial risk management platforms.

Professional Applications and Use Cases

Risk Management Implementation: Moving from conceptual understanding to practical implementation, GARCH provides the specific mathematical framework that translates volatility forecasts into actionable risk limits. Portfolio managers use GARCH parameter estimates to set position sizes that adapt to changing market conditions: reducing exposure when GARCH forecasts indicate rising volatility and accepting larger positions during predicted calm periods. The three-parameter structure enables systematic, rule-based risk management that removes emotional decision-making from volatility timing.

Regulatory Capital Calculations: Banking regulators specifically mandate GARCH-type models for calculating Value-at-Risk metrics that determine required capital reserves. The model's mathematical properties—particularly the mean reversion constraint and finite unconditional variance—provide the stability that regulators require for systemic risk assessment. Unlike more exotic volatility models, GARCH's well-understood parameter space enables regulatory validation and cross-institutional comparison.

Derivatives Pricing Enhancement: Options traders use GARCH volatility forecasts to identify mispriced derivatives in markets that still rely on constant volatility assumptions. The model's ability to predict changes in volatility regimes provides systematic advantages in volatility arbitrage strategies. Professional derivatives desks integrate GARCH forecasts into their pricing models to capture volatility risk premiums that constant volatility models miss entirely.

Preparation for Advanced Exploration

This foundation prepares you for hands-on parameter manipulation in the Advanced Parameter Laboratory, where you'll observe how changing α, β, and ω values affects simulated market data in real-time. The mathematical precision introduced here becomes crucial when you estimate GARCH models on actual financial data and interpret the economic significance of parameter estimates.

Understanding these fundamentals also provides the conceptual framework needed for the complete mathematical development that follows, where we'll derive the theoretical properties that make GARCH both practically useful and theoretically sound.

Advanced Parameter Laboratory

Level 2Level 2 (Interactive)

Building upon the Interactive Parameter Exploration introduced in the volatility overview, this advanced laboratory provides deeper tools for understanding GARCH parameter behavior. While the application-level exploration focused on basic parameter effects, here you can explore preset scenarios based on actual maximum likelihood estimates from real financial markets, observe advanced statistics like half-life and unconditional variance, and understand how parameter combinations relate to specific asset classes and market regimes. Each preset configuration represents a different market structure or economic regime, helping you understand the economic interpretation of parameter values.

Enhanced Interactive Exploration

The enhanced parameter explorer below includes real-time calculation of advanced statistics that help interpret the economic significance of parameter combinations. The explorer starts with S&P 500 typical parameter values. Rather than exploring parameters in isolation, you can select from market-based presets above or adjust parameters manually to observe how the half-life of volatility shocks, unconditional variance target, and excess kurtosis change in response to your modifications.

GARCH(1,1) Parameter Explorer
Model Parameters

ω (omega) - Long-run variance

0.005

0.001Step: 0.001 (normal)0.150
Base level of volatility (in %²)

α (alpha) - ARCH effect

0.080

0.009Step: 0.010 (normal)0.300
Reaction to recent shocks

β (beta) - GARCH effect

0.900

0.690Step: 0.010 (normal)0.990
Persistence of volatility

Volatility Comparison

GARCH vs Historical Measures (252 days)


Impulse Response Function

Response to a 1% shock over time

Parameter Interpretation
Note: The simulation uses the same random shock sequence for all parameter combinations, allowing you to see the pure effect of parameter changes on volatility patterns.

ω = 0.0050
Base volatility level. Higher values increase the floor of volatility.

α = 0.080
Shock sensitivity. Higher values mean stronger reaction to recent returns.

β = 0.900
Volatility persistence. Higher values mean volatility changes last longer.

Understanding Parameter Economic Significance
Half-Life Interpretation

Formula: ln(0.5)/ln(α+β)

• < 5 days: Rapid volatility decay (rare)
• 5-20 days: Typical equity markets
• 20-50 days: Persistent markets (FX)
• > 50 days: Near-IGARCH behavior

Unconditional Variance

Formula: ω/(1-α-β)

• Represents long-run volatility target
• Must be positive (requires α+β<1)
• Typical values: 15-25% annual volatility
• Crisis periods: 30-60% annual volatility

Excess Kurtosis

Condition: 3α2+2αβ+β2<1

• Normal distribution: kurtosis = 3
• GARCH typically generates: 4-8
• Higher persistence → fatter tails
• Crisis periods: extreme tail behavior

Experimental Scenarios for Learning

Try these parameter experiments to deepen your understanding of GARCH behavior:

Experiment 1: Shock Sensitivity

Start with α = 0.05, β = 0.90. Gradually increase α to 0.15 while reducing β to 0.83. Observe how volatility spikes become more pronounced and the unconditional variance increases. This simulates moving from stable to crisis market conditions while maintaining stationarity.

Experiment 2: Persistence Effects

Fix α = 0.04 and vary β from 0.70 to 0.94. Notice how higher β values create longer-lasting volatility episodes and increase the half-life of shocks. This demonstrates why FX markets with central bank intervention show high persistence while maintaining stationarity.

Experiment 3: IGARCH Boundary

Set α+β = 0.99 (very close to 1) and observe how volatility shocks become extremely persistent. Compare this with α+β = 0.85 to see the dramatic difference in shock persistence and mean reversion speed near the stationarity boundary.

Experiment 4: Baseline Volatility

With fixed α = 0.08, β = 0.90, vary ω from 0.001 to 0.100 STD. This changes the unconditional variance target while keeping the dynamic response structure constant. Notice how this affects the overall volatility level without changing clustering patterns.

Preparation for Estimation and Forecasting

The parameter relationships you've explored here become crucial when you move to the Estimation Playground, where you'll fit GARCH models to real data and interpret the economic significance of the estimated parameters. Understanding realistic parameter ranges helps distinguish between economically meaningful estimates and numerical artifacts.

Similarly, the half-life and persistence concepts explored here directly inform volatility forecasting accuracy and confidence interval construction in the Forecasting Interface that follows.

Building upon the mathematical framework introduced in the volatility overview, this section develops the complete theoretical foundation for GARCH(1,1) models. We examine the stochastic properties that ensure model consistency, derive moment conditions for parameter identification, and establish the asymptotic theory underlying maximum likelihood estimation. This rigorous development connects the practical applications explored in previous sections to their theoretical foundations.

Stochastic Properties and Existence Conditions

The GARCH(1,1) process must satisfy specific regularity conditions to ensure well-defined stochastic behavior. Consider the complete probability space (Ω, ℱ, P) with natural filtration{ℱₜ} generated by the return process. The existence of a strictly stationary solution requires careful analysis of the recursive variance equation.

🔄 Loading mathematical content... (1 equations)
Moment Existence and Unconditional Distribution

The existence of unconditional moments for the GARCH process requires increasingly restrictive conditions as the moment order increases. These conditions directly impact the applicability of asymptotic theory and the validity of standard econometric inference procedures.

🔄 Loading mathematical content... (1 equations)
💡 Economic Intuition: Why These Conditions Matter

The mathematical conditions derived above ensure that volatility does not spiral out of control. The strict stationarity condition E[ln(αzt2+β)]<0 and the simpler covariance stationarity condition α+β<1 both prevent volatility from building on itself indefinitely. When these conditions hold, past shocks eventually die out, and the process reverts toward a long-run average variance — a key feature of financial markets where risk fluctuates but does not explode indefinitely.

The higher-order moment conditions (particularly the fourth moment condition 3α2+2αβ+β2<1) capture how “fat” the tails of the return distribution are. GARCH models naturally produce distributions with heavier tails than the normal distribution, reflecting the fact that markets exhibit more extreme events than Gaussian theory would predict. As the kurtosis formula shows, when persistence ( α+β) rises toward unity, the tails become dramatically fatter, meaning the model assigns more probability to crises or turbulence. These conditions thus provide the bridge between statistical rigor and the intuitive market story: volatility clusters, shocks fade slowly, and risk is punctuated by rare but severe events.

Autocorrelation Structure and Temporal Dependencies

While GARCH innovations {εt} form a martingale difference sequence, the squared innovations exhibit rich autocorrelation structure that captures volatility clustering. This section derives the exact autocorrelation functions and their implications for model identification.

🔄 Loading mathematical content... (1 equations)
💡 Economic Intuition: Volatility Clustering and Predictability

The autocorrelation structure derived above reveals why volatility clustering creates useful patterns in financial risk measurement. While returns themselves are unpredictable (the martingale property), their magnitudes show strong persistence—today's volatility predicts tomorrow's volatility with meaningful accuracy. This enables professional risk managers to anticipate periods of elevated uncertainty and adjust position sizes accordingly, though these forecasts remain inherently uncertain and require careful validation.

The ARMA(1,1) representation with autoregressive coefficient ( α+β) shows that volatility behaves like a predictable time series with its own momentum and mean reversion. This mathematical structure enables forecasting models that help institutional investors time their risk exposure: reducing positions when volatility is predicted to rise and increasing exposure during forecasted calm periods. The exponential decay rate (α+β)k shown in the autocorrelation formula directly determines how far ahead meaningful volatility predictions can be made.

Asymptotic Theory for Maximum Likelihood Estimation

The asymptotic properties of GARCH parameter estimators depend on regularity conditions that ensure the applicability of standard likelihood theory. The quasi-maximum likelihood estimator (QMLE) possesses desirable properties even under distributional misspecification.

🔄 Loading mathematical content... (1 equations)
💡 Economic Intuition: Robust Estimation During Market Stress

The quasi-maximum likelihood estimator (QMLE) derived above provides important protection when market assumptions break down. Real financial data rarely follows the perfect Gaussian distributions that textbooks assume—during crises, markets exhibit extreme skewness, changing correlations, and structural breaks. The QMLE's robustness means that GARCH parameter estimates remain consistent under milder conditions than full maximum likelihood, though this protection has limits when innovations become extremely heavy-tailed or when structural breaks occur.

The sandwich estimator A-1BA-1 shown in the information matrix components captures the extra uncertainty that comes from model misspecification. When markets behave as expected, A=B and standard errors simplify to the familiar A-1. But when market behavior deviates from assumptions, the sandwich form automatically inflates uncertainty estimates, providing more conservative confidence intervals. This mathematical machinery helps practitioners get more honest assessments of parameter uncertainty, though it cannot eliminate all risks from model misspecification.

Multi-Step Ahead Forecasting Theory

The theoretical development of GARCH forecasting formulas provides the foundation for the practical forecasting tools explored in previous sections. We derive the exact forecast formulas and characterize their asymptotic properties.

🔄 Loading mathematical content... (1 equations)
💡 Economic Intuition: Volatility Forecasting and Investment Horizons

The h-step ahead variance forecast formula derived above provides the mathematical foundation for matching risk models to investment decision horizons. The exponential convergence at rate (α+β)h toward the unconditional variance ω1-α-β captures how market turbulence gradually gives way to historical norms. This enables informed decisions about when temporary volatility spikes will fade versus when structural shifts may be occurring, though these forecasts become increasingly uncertain as the horizon extends.

The linear growth of forecast error variance ( VarC·h) reflects a fundamental trade-off in risk management: longer forecasts become less precise but remain centered on economically sensible long-run averages. This mathematical property guides practical decisions about forecast confidence intervals and helps determine optimal rebalancing frequencies for different investment strategies. The constant C depends on fourth moments and parameter values, connecting forecast uncertainty directly to the underlying volatility dynamics captured by the GARCH process.

Extensions to GARCH(p,q) and Theoretical Connections

The theoretical framework developed for GARCH(1,1) extends naturally to higher-order specifications and provides the foundation for understanding model extensions like EGARCH, GJR-GARCH, and multivariate GARCH models.

🔄 Loading mathematical content... (1 equations)
Connection to Practical Applications

The theoretical properties developed here provide the rigorous foundation for all practical GARCH applications explored in previous sections. Parameter constraints ensure model coherence, moment conditions validate asymptotic inference procedures, and forecasting theory supports risk management applications.

Understanding these mathematical foundations enables proper model specification, accurate uncertainty quantification, and valid statistical inference—essential skills for professional econometric modeling and advanced research applications.

From Stylized Facts to GARCH Implementation

Level 1Level 1 (Conceptual)
Translating Market Observations into Mathematical Structure

The empirical stylized facts described in the volatility overview represent universal patterns observed across financial markets. GARCH's innovation lies in providing a specific mathematical framework that not only captures these patterns but generates them naturally through parameter interactions. This section shows exactly how GARCH's three-parameter structure implements each stylized fact and why certain parameter ranges produce realistic market behavior while others generate unrealistic dynamics.

Stylized Fact Implementation Matrix
1. Volatility Clustering
Market Observation:

Large price movements tend to be followed by large movements, creating distinct periods of market turbulence and calm.

GARCH Implementation:

The equation σt2=ω+α·εt-12+β·σt-12 creates clustering through both α (recent shock impact) and β (volatility persistence).

Parameter Impact:

• High α (≈0.15): Sharp volatility spikes after shocks
• High β (≈0.83): Long-lasting volatility periods
α+β ≈ 0.98: Strong clustering, realistic persistence

2. Fat Tails
Market Observation:

Extreme returns occur much more frequently than normal distributions predict, with excess kurtosis typically 3-10 times normal levels.

GARCH Implementation:

The multiplicative structure εt=σt·zt creates fat tails even when zt is normally distributed, through time-varying volatility.

Parameter Impact:

• Higher α: More extreme tail events
• Higher persistence (α+β): Fatter tails
• Typical kurtosis: 4-8 (vs. 3 for normal)

3. Mean Reversion
Market Observation:

While volatility clusters short-term, it exhibits long-run stability, reverting to historical averages over months or years.

GARCH Implementation:

The constraint α+β<1 ensures convergence to unconditional variance ω/ (1-α-β), providing mathematical mean reversion.

Parameter Impact:

α+β = 0.95: Half-life ≈ 14 days
α+β = 0.99: Half-life ≈ 69 days
α+β = 1.00: No mean reversion (IGARCH)

Why Naive Models Fail: Quantitative Evidence

Understanding GARCH's sophistication requires comparing it against simpler alternatives that fail to capture stylized facts. Consider four approaches to modeling daily S&P 500 returns during the 2008 financial crisis:

Constant Volatility Model

Assumes σt=σ=20% annually

Problems:

• Misses volatility spikes during crisis
• Underestimates tail risk by factor of 10
• Cannot predict volatility changes
• VaR violations: 15% vs. expected 1%

Rolling Window Model

Uses σt=std(returns,30-day window)

Problems:

• Reacts slowly to regime changes
• Artificial jumps at window boundaries
• Equal weights for all observations
• Poor forecasting performance

EWMA Model

Uses σt2=λ·σt-12+(1-λ)·rt-12

Problems:

• IGARCH special case (λ + (1-λ) = 1)
• No mean reversion property
• Shocks persist indefinitely
• Fixed λ = 0.94 limits flexibility

GARCH(1,1) Model

Uses σt2=ω+α·εt-12+β·σt-12

Advantages:

• Captures volatility clustering naturally
• Adapts quickly to regime changes
• Provides forward-looking forecasts
• VaR violations: 1.2% vs. expected 1%

Parameter Constraints and Economic Interpretation

GARCH's parameter constraints aren't arbitrary mathematical requirements — they ensure the model generates economically sensible behavior. Each constraint translates directly into market realism:

Mathematical Constraints

ω>0: Ensures positive volatility always

α0: Shocks increase (never decrease) volatility

β0: Past volatility persists positively

α+β<1: Volatility shocks eventually dissipate

Economic Interpretation

Realistic Parameter Ranges:

ω: 0.001 - 0.1 (scaled by 1000 for display)

α: 0.05 - 0.20 (shock sensitivity)

β: 0.75 - 0.95 (persistence)

α+β: 0.85 - 0.99 (overall persistence)

Real-World Parameter Examples Across Asset Classes

Different financial markets exhibit characteristic GARCH parameter patterns that reflect their underlying economic structures and participant behavior:

Equity Markets

S&P 500 typical:
ω ≈ 0.005
α ≈ 0.08
β ≈ 0.90

Moderate shock response, high persistence reflecting institutional investor behavior

Foreign Exchange

EUR/USD typical:
ω ≈ 0.002
α ≈ 0.04
β ≈ 0.94

Lower shock response, very high persistence due to central bank interventions

Commodity Markets

Oil futures typical:
ω ≈ 0.03
α ≈ 0.15
β ≈ 0.82

High shock response, moderate persistence reflecting supply/demand shocks

Bond Markets

10Y Treasury typical:
ω ≈ 0.001
α ≈ 0.03
β ≈ 0.96

Very low shock response, extreme persistence reflecting monetary policy stability

Bridge to Interactive Exploration

Understanding how parameters generate stylized facts prepares you for the Advanced Parameter Laboratory, where you can manipulate α, β, and ω values and observe their real-time effects on simulated market data. The parameter ranges and asset class examples provided here will guide your exploration of realistic versus unrealistic parameter combinations.

This implementation knowledge also provides the foundation for understanding why maximum likelihood estimation typically produces parameter values within these ranges and what economic forces drive the parameter differences observed across asset classes.

While the volatility overview introduced model comparison conceptually, this dashboard provides quantitative evidence of GARCH's superiority across different market conditions. By examining performance metrics during normal markets, financial crises, and low volatility periods, you'll understand when and why GARCH outperforms simpler alternatives like constant volatility assumptions, rolling window estimates, and exponentially weighted moving averages.

Market Period Selection
Normal Market (2004-2006)

Low volatility, steady growth period

Gradual volatility changes, few extreme events

Financial Crisis (2007-2009)

High volatility, regime changes

Volatility clustering, frequent extreme movements

Post-Crisis (2010-2012)

Declining volatility, policy intervention

Mean reversion, central bank influence

Low Volatility (2014-2017)

Persistently low volatility regime

Extended calm period, compressed risk premiums

Performance Comparison: Financial Crisis (2007-2009)
GARCH(1,1)

Loading real data...

Historical Average

Loading real data...

EWMA

Loading real data...

When Different Models Excel
GARCH Advantages

• Crisis periods with volatility clustering
• Markets with regime changes
• Medium to long-term forecasting
• Regulatory capital calculations
• Options pricing with volatility smile

Alternative Model Use Cases

EWMA: Real-time risk monitoring (RiskMetrics)
Rolling Window: Very short-term tactical decisions
Constant: Baseline comparison, simple backtesting
Regime Models: Markets with clear structural breaks

Foundation for Advanced Model Selection

The quantitative evidence presented here demonstrates why GARCH has become the industry standard for volatility modeling across financial institutions. The consistent outperformance across different market conditions validates the theoretical insights from the parameter exploration.

This comparison framework prepares you for understanding GARCH extensions (EGARCH, GJR-GARCH) that address specific limitations revealed by diagnostic testing, particularly the asymmetric response to positive and negative shocks not captured by the basic GARCH specification.

Volatility Forecast Visualization

Interactive comparison of GARCH forecasting performance against simpler alternatives. Explore how different models adapt to market volatility across various time periods.

Model Diagnostic Analysis
GARCH Model Diagnostics

Ljung-Box (10 lags): p-value = 1.00

ARCH-LM (5 lags): p-value = 1.00

Jarque-Bera: p-value = 1.00

? No data available
? No data available
? No data available

Constant Volatility Diagnostics

Ljung-Box (10 lags): p-value = 1.00

ARCH-LM (5 lags): p-value = 1.00

Jarque-Bera: p-value = 1.00

? No data available
? No data available
? No data available

Rolling Window Diagnostics

Ljung-Box (10 lags): p-value = 1.00

ARCH-LM (5 lags): p-value = 1.00

Jarque-Bera: p-value = 1.00

? No data available
? No data available
? No data available

Diagnostic Test Interpretation

GARCH models should show no remaining autocorrelation in standardized residuals (Ljung-Box test) and no remaining ARCH effects (ARCH-LM test). P-values above 0.05 indicate the model has successfully captured the volatility clustering. Some non-normality may remain, which can be addressed with alternative error distributions (Student-t, GED) in advanced implementations.

Educational Insights
Crisis Adaptation Speed

GARCH adapts 40% faster to volatility spikes than alternatives

Forecast Accuracy

GARCH typically achieves 20-40% better Q-LIKE score than historical average

Volatility Clustering

Notice how GARCH captures volatility persistence during market stress

GARCH Adaptation Speed

Notice how GARCH (blue line) responds quickly to volatility spikes and gradually returns to normal levels.

Model Performance Differences

Compare how different models respond to the same volatility events. GARCH typically shows 20-40% better accuracy.

Volatility Clustering Effect

High volatility periods tend to cluster together. GARCH captures this pattern better than simple models.

Moving beyond parameter exploration to practical analysis, this lab lets you examine real GARCH parameter estimates from financial data. You'll analyze how different assets exhibit varying volatility dynamics through their parameter values, understand statistical significance, and interpret the economic meaning of these estimates. This bridges theory with empirical evidence from financial markets.

Dataset Selection and Characteristics
S&P 500 (1990-2020)

Daily returns including dot-com crash, financial crisis, and COVID-19

7812 observations

Moderate volatility clustering, multiple crisis periods

Parameter Analysis Results
ParameterEstimateStd. Errort-statisticp-value

ω (omega)

0.01760.00453.940.000

α (alpha)

0.1026920.00536419.140.000

β (beta)

0.8828020.004821183.150.000
Model Statistics

Log-Likelihood: -10358.4

AIC: 20722.8

BIC: 20749.4

Persistence (α+β): 0.985

Half-life: 47.4 days

Uncond. Variance: 0.001211

Economic Interpretation of Results
Parameter Significance

All parameters show t-statistics well above 2.0 and p-values below 0.05, indicating strong statistical significance. This confirms that both shock response (α) and volatility persistence (β) are essential for modeling this time series.

Persistence Analysis

The persistence measure (α+β = 0.985) indicates highly persistent volatility shocks with a half-life of 47.4 days. This level is typical for financial markets and confirms strong volatility clustering.

Model Selection

The AIC (20723) and BIC (20749) values can be compared across different model specifications (GARCH(1,1) vs GARCH(1,2) vs EGARCH) to select the optimal model. Lower values indicate better fit while penalizing model complexity. Note that absolute values depend on sample size.

Understanding Parameter Analysis
Parameter Interpretation

These parameter estimates reveal the volatility dynamics of each asset. ω determines the long-run volatility level, α captures immediate shock response, and β measures volatility persistence. Together, they define how each market responds to and maintains volatility over time.

Statistical Reliability

These parameter estimates come from rigorous maximum likelihood estimation with robust standard errors. The t-statistics and p-values confirm statistical significance, while the model selection criteria help compare across specifications.

Common Interpretation Challenges
Boundary Solutions

When α or β approach constraint boundaries (especially β ≈ 1), this suggests extremely persistent volatility that may require more sophisticated models or indicate structural breaks in the data.

Extreme Persistence

Very high persistence (α+β > 0.99) suggests volatility shocks decay extremely slowly, which is common in FX markets but may indicate model misspecification.

Outliers and Breaks

Extreme outliers or structural breaks can distort parameter estimates. Consider outlier-robust estimation or break-point detection before estimation.

Small Sample Issues

GARCH estimation requires substantial data (500+ observations) for reliable parameter estimates. Small samples may yield unstable or biased estimates.

Bridge to Forecasting and Application

The parameter estimates analyzed here feed directly into the Forecasting Interface, where you'll see how these values translate into volatility predictions. Understanding parameter significance and economic interpretation helps assess forecast reliability and model appropriateness for different applications.

The statistical properties explored here (standard errors, likelihood values, information criteria) connect directly to diagnostic testing and model validation, helping you understand when GARCH models are appropriate and when extensions may be needed.

Successful GARCH modeling requires systematic validation that the fitted model adequately captures the data characteristics. This comprehensive diagnostic framework tests the key assumptions underlying GARCH models and identifies potential specification problems. Proper diagnostic testing is essential for reliable inference, accurate forecasting, and robust risk management applications.

Diagnostic Testing Framework

GARCH model validation follows a systematic approach, testing whether the standardized residuals exhibit the properties expected under correct specification. The testing sequence proceeds from basic autocorrelation tests through advanced specification diagnostics.

Autocorrelation Tests for Standardized Residuals

The Ljung-Box test examines whether standardized residuals z^t=εt/σt exhibit remaining autocorrelation. Under correct GARCH specification, these residuals should be serially uncorrelated, indicating that the conditional mean and variance models have captured all temporal dependencies.

TestLagsStatisticp-valueResult
Ljung-Box53.240.663
Pass
Ljung-Box108.910.540
Pass
Ljung-Box2018.450.558
Pass
Interpretation Guidelines

p-value > 0.05: Fail to reject null hypothesis of no autocorrelation. This indicates the GARCH model has successfully captured temporal dependencies in the conditional variance.

p-value ≤ 0.05: Evidence of remaining autocorrelation, suggesting model misspecification or need for higher-order GARCH terms.

Common Remedies

Higher-order GARCH: Try GARCH(2,1) or GARCH(1,2)
Mean equation: Include AR/MA terms if needed
Structural breaks: Test for parameter instability
Alternative specifications: Consider EGARCH or GJR-GARCH

Comprehensive Diagnostic Summary
Model Adequacy Assessment

✓ No remaining autocorrelation in residuals
✓ No remaining ARCH effects detected
✓ Independence tests passed
⚠ Non-normality in residuals

Overall: Model captures volatility dynamics well, consider alternative distributions.

Recommended Actions

1. Re-estimate with Student-t errors to address non-normality
2. Use robust standard errors for inference
3. Investigate outliers for data quality
4. Consider asymmetric models if leverage effects present
5. Validate with out-of-sample forecasting

Next Steps for Research

Model comparison: Compare AIC/BIC across specifications
Forecast evaluation: Out-of-sample performance testing
Structural stability: Parameter constancy tests
Alternative specifications: EGARCH, GJR-GARCH evaluation
Risk model validation: VaR backtesting

Integration with Model Development Process

Diagnostic testing represents the crucial validation step in the iterative model development process. Failed diagnostics indicate specific areas for model improvement, whether through alternative specifications (EGARCH, GJR-GARCH), different error distributions, or structural break accommodation.

The diagnostic framework developed here provides the foundation for rigorous empirical research and practical model implementation. Understanding these validation procedures ensures that GARCH models meet the statistical requirements for reliable inference and accurate risk assessment.

The ultimate test of any volatility model lies in its forecasting performance. This interface transforms the theoretical GARCH forecasting formulas into a practical tool for generating multi-step ahead volatility predictions. You'll experience how GARCH forecasts evolve over time, how parameter uncertainty affects confidence intervals, and how these predictions connect to real-world risk management applications like Value-at-Risk calculations and portfolio optimization.

Key Forecasting Concepts
Mean Reversion

GARCH forecasts always converge to the unconditional (long-run) volatility over time. High volatility periods eventually return to normal levels.

Persistence

Persistence (α+β) measures how long volatility shocks last. Values near 1.0 indicate very persistent volatility that decays slowly.

Forecast Uncertainty

Confidence intervals widen over longer horizons, reflecting greater uncertainty about distant future volatility. Parameter estimation adds additional uncertainty.

Forecast Configuration

Select from real GARCH estimation results across different asset classes to explore how model parameters affect volatility forecasting behavior.

Forecast Horizon: 180 days

Select forecast horizon from 1 month to 1 year
Current Model: S&P 500 (1990-2020)

Estimated from 30 years of S&P 500 daily returns

ω = 0.0176
α = 0.103
β = 0.883

Persistence:0.985494
Half-life: 47.4 days
Unconditional Vol: 17.5% annual

Volatility Forecast Visualization
Forecast Performance Metrics
Forecast Accuracy

MSE: 0.0015
MAE: 0.025
QLIKE: 2.34

vs. Historical Average

GARCH achieves 25%better
forecast accuracy than simple
historical volatility estimation

vs. EWMA

GARCH provides 15%improvement
over exponentially weighted
moving average forecasts

Forecast Horizon Analysis
1-Day Ahead Forecast

Daily Vol: NaN%
Annual Vol: NaN%
95% CI: [NaN%, NaN%]

30-Day Ahead Forecast

Daily Vol: NaN%
Annual Vol: NaN%
95% CI: [NaN%, NaN%]

Long-term Convergence

Converges to: 17.5% annual
95% Converged: 206 days
Forecast Certainty: Decreases with horizon

Risk Management Applications
Position Sizing

With tomorrow's forecasted volatility of NaN%, a $100M portfolio targeting 15% annual volatility should maintain a position size of approximately NaNx leverage. This adjusts automatically as volatility forecasts change.

Options Pricing Impact

The forecasted volatility term structure affects options pricing across different expiration dates. Near-term options use the NaN% forecast, while 30-day options incorporate the NaN% prediction and longer-term convergence to 17.5%.

Understanding Forecast Properties
Mean Reversion

GARCH forecasts always converge to the unconditional volatility over time. The speed of convergence depends on persistence (α+β): higher persistence means slower convergence and longer-lasting volatility shocks.

Forecast Uncertainty

Confidence intervals widen as the forecast horizon increases, reflecting greater uncertainty about distant future volatility. Parameter estimation uncertainty and model specification uncertainty both contribute to forecast errors.

Forecast Limitations

GARCH cannot predict volatility direction changes or structural breaks. Forecasts assume model stability and may perform poorly during regime changes or extreme market stress not captured in the estimation sample.

Forecast Evaluation Techniques
Out-of-Sample Testing

Evaluate forecast accuracy using volatility-specific metrics: QLIKE (Quasi-Likelihood) for distributional fit, Diebold-Mariano tests for statistical significance, and Model Confidence Sets for ranking multiple models. These address volatility's unobservable nature and heteroskedasticity better than basic MSE/MAE.

Model Comparison

Compare GARCH forecasts against simpler alternatives (EWMA, rolling volatility) and more sophisticated models (EGARCH, stochastic volatility). Diebold-Mariano tests provide statistical significance for forecast accuracy differences.

Parameter Sensitivity Analysis
Model Parameters

ω (Constant Term): 0.0000

Long-run variance component. Higher ω increases unconditional volatility.

α (ARCH Effect): 0.080

Shock impact. Higher α means stronger immediate response to innovations.

β (Persistence): 0.900

Volatility persistence. Higher β means longer-lasting volatility effects.
Parameter Impact on Forecasts

Parameter Summary

Persistence (α + β): 0.980
Unconditional Vol: 25.1%
Half-life: 35 days
Status: ✓ Stationary
ω (Omega) Impact

Controls the long-run volatility level. Changes in ω affect all forecast horizons equally, shifting the entire forecast path up or down without changing the dynamics. Higher ω increases the level where forecasts converge.

α (Alpha) Impact

Controls immediate shock response. Higher α increases the initial forecast jump after volatility innovations but does not affect long-term persistence. Creates more reactive but less persistent volatility forecasts.

β (Beta) Impact

Controls volatility persistence. Higher β makes volatility shocks last longer, creating slower mean reversion and wider confidence intervals. Most critical parameter for forecast accuracy and financial risk management.

Advanced Forecast Evaluation
Quasi-Maximum Likelihood (QLIKE) Analysis

QLIKE provides robust evaluation of volatility forecasts using the log-likelihood ratio. This metric is particularly valuable because it is robust to the fourth moment assumption and provides scale-invariant comparisons across different volatility regimes.

QLIKE=[ln(σt2)+rt2σt2]

2.340

GARCH QLIKE

2.761

Historical QLIKE

Performance: GARCH shows 18.0% better QLIKE score than historical average. Lower scores indicate superior forecast performance.

Why QLIKE Over MSE?

Unlike MSE, QLIKE is derived from likelihood principles and remains robust under fat tails and asymmetric distributions common in financial data.

Risk Management Impact: The 18% QLIKE improvement translates to approximately 10-15% better VaR accuracy and more efficient portfolio allocation.

Technical Note: QLIKE evaluation uses realized volatility as proxy for true unobservable volatility, requiring substantial sample sizes for reliable inference.

Diebold-Mariano Forecast Comparison

Statistical test for forecast accuracy differences. Null hypothesis: equal predictive accuracy. Tests whether GARCH forecasts are significantly better than benchmarks.

DM=d¯Var(d¯)N(0,1)

GARCH vs Historical Average

DM Statistic: -3.47 | p-value: 0.0005

Statistically Significant

GARCH vs EWMA

DM Statistic: -2.13 | p-value: 0.033

Significant at 5%

GARCH vs t-GARCH

DM Statistic: -0.87 | p-value: 0.384

Not Significant

Statistical Interpretation: Negative DM statistics indicate GARCH superiority. The -3.47 statistic against historical volatility provides strong evidence (p < 0.001) that GARCH forecasts are significantly more accurate.

Practical Significance: The non-significant result versus t-GARCH (p = 0.384) suggests both models provide similar forecast accuracy, validating GARCH(1,1) as a parsimonious choice.

Model Confidence Set (MCS)

Hansen et al. (2011) procedure for identifying the set of models with superior forecast performance.

✓ In MCS at 90%
GARCH(1,1) included in Model Confidence Set
Reality Check Bootstrap

White (2000) test controlling for data-snooping bias when comparing multiple forecasting models.

0.021
p-value (Bootstrap)
Loss Differential Analysis

Asymmetric loss functions for directional accuracy and tail risk evaluation.

Directional Accuracy: 67.8%

Hit Rate (VaR): 98.2%

Economic Significance Framework

Statistical significance does not always translate to economic significance. Evaluating the practical value of superior forecast accuracy requires considering transaction costs, portfolio constraints, and implementation feasibility in real-world trading environments.

Portfolio Applications

Superior volatility forecasts can enhance portfolio optimization through better risk budgeting, more efficient diversification, and improved timing of portfolio rebalancing decisions. The magnitude of benefits depends on the specific strategy and market conditions.

Implementation Considerations

Economic value realization requires careful consideration of transaction costs, liquidity constraints, and model parameter stability. Academic studies suggest potential benefits, but actual outcomes vary significantly across different market environments and implementation approaches.

Integration with Advanced Topics

The forecasting experience here connects directly to the mathematical framework that follows, where you will explore the theoretical foundations behind these forecasting formulas. Understanding how forecasts are generated and their limitations prepares you for advanced topics like forecast combination, regime-switching models, and alternative volatility specifications.

The practical applications demonstrated here—position sizing, VaR calculation, options pricing—represent the ultimate purpose of volatility modeling. These connections help bridge the gap between academic theory and professional risk management practice.

While GARCH(1,1) represents a significant advancement over constant volatility models, empirical research has identified several systematic limitations that motivate more sophisticated specifications. This section analyzes these limitations and shows how they provide the theoretical foundation for the rich family of GARCH extensions available in V-Lab, including EGARCH, GJR-GARCH, and multivariate models.

Systematic Limitations of Basic GARCH

Despite its success in capturing volatility clustering and fat tails, GARCH(1,1) exhibits several empirical deficiencies that become apparent through careful diagnostic analysis and comparison with observed market behavior. Understanding these limitations provides the motivation and theoretical foundation for advanced volatility models.

1. Symmetric Response to Shocks

GARCH treats positive and negative shocks of equal magnitude identically, failing to capture the empirically observed leverage effect where negative returns tend to increase volatility more than positive returns of the same size.

🔄 Loading mathematical content... (1 equations)

This symmetry contradicts empirical evidence from equity markets

2. Exponential Decay Structure

The fixed exponential decay rate (α+β)ʰ may be too restrictive for some markets that exhibit more complex volatility persistence patterns, including hyperbolic decay or regime-dependent persistence.

🔄 Loading mathematical content... (1 equations)

Real markets may show slower or faster decay depending on conditions

3. Parameter Constancy Assumption

GARCH assumes constant parameters over the entire sample period, which may be violated during structural breaks, policy regime changes, or crisis periods that fundamentally alter market volatility dynamics.

This limitation motivates regime-switching and time-varying parameter models

4. Univariate Framework

Standard GARCH cannot capture volatility spillovers, correlation dynamics, or contagion effects that are crucial for portfolio risk management and understanding financial system interconnections.

This limitation drives development of multivariate GARCH specifications

Asymmetric Volatility Models: Theoretical Development

The most empirically significant limitation of basic GARCH is its symmetric treatment of positive and negative shocks. This has motivated several asymmetric extensions that allow differential volatility responses to good and bad news.

🔄 Loading mathematical content... (1 equations)
Multivariate Extensions: From Univariate to System Modeling

Modern financial analysis requires understanding volatility interactions across multiple assets, markets, and time zones. This has driven the development of multivariate GARCH models that capture volatility spillovers and time-varying correlations.

🔄 Loading mathematical content... (1 equations)
Advanced Topics and Research Frontiers

The GARCH framework continues to evolve, with current research addressing high-frequency data modeling, regime switching, machine learning integration, and alternative risk measures. Understanding these developments helps practitioners stay current with evolving best practices.

High-Frequency Volatility Modeling

Realized GARCH: Incorporates realized volatility measures from high-frequency data into GARCH frameworks, improving forecasting accuracy.

MIDAS-GARCH: Mixed-frequency models that combine daily returns with intraday volatility measures for enhanced prediction.

Regime-Switching Extensions

Markov-Switching GARCH: Allows GARCH parameters to change across different volatility regimes (crisis vs. normal periods).

Smooth Transition GARCH: Enables gradual parameter changes based on threshold variables or time trends.

Machine Learning Integration

Neural Network GARCH: Uses neural networks to model complex non-linear relationships in volatility dynamics.

Support Vector GARCH: Applies machine learning classification techniques to volatility regime identification and forecasting.

Alternative Risk Measures

Expected Shortfall GARCH: Models conditional expected shortfall rather than just volatility for enhanced tail risk assessment.

Quantile GARCH: Directly models conditional quantiles of return distributions for comprehensive risk analysis.

Model Selection and Specification Strategy

Choosing among GARCH variants requires systematic evaluation based on data characteristics, intended application, and computational constraints. This framework guides practitioners through the model selection process.

Decision Framework for GARCH Extensions
Step 1: Diagnostic Analysis

• Test basic GARCH(1,1) specification
• Check for asymmetric responses
• Examine parameter stability
• Assess forecasting performance

Step 2: Extension Selection

Leverage effects: EGARCH or GJR-GARCH
Multivariate needs: DCC or BEKK
Regime changes: Markov-switching
High-frequency data: Realized GARCH

Step 3: Validation and Comparison

• Information criteria (AIC, BIC)
• Out-of-sample forecasting tests
• Risk model backtesting
• Computational efficiency assessment

Cross-References to V-Lab Model Family

The theoretical foundations developed here directly support the other volatility models available in V-Lab. Each extension addresses specific limitations identified in this analysis while building upon the core GARCH framework.

Application-Specific Guidance

Portfolio Risk Management: Start with DCC for multiple assets
Options Pricing: Consider asymmetric models (EGARCH/GJR)
Regulatory Capital: Basic GARCH often sufficient
High-Frequency Trading: Realized GARCH specifications
Emerging Markets: Regime-switching extensions

This comprehensive bibliography provides scholarly references for deepening your understanding of GARCH models, their theoretical foundations, empirical applications, and extensions. The references are organized by topic to support different learning objectives, from foundational theory through advanced research applications.

Foundational Papers
Engle, R. F.

1982

"Autoregressive Conditional Heteroskedasticity with Estimates of the Variance of United Kingdom Inflation"

The seminal paper introducing ARCH models and conditional heteroskedasticity. Essential reading for understanding the historical development of volatility modeling.

Econometrica, 50(4), 987-1007

Bollerslev, T.

1986

"Generalized Autoregressive Conditional Heteroskedasticity"

The original GARCH paper introducing the generalized framework that forms the basis for all subsequent GARCH model development.

Journal of Econometrics, 31(3), 307-327

Bollerslev, T., Chou, R. Y., & Kroner, K. F.

1992

"ARCH modeling in finance: A review of the theory and empirical evidence"

Comprehensive early survey of ARCH/GARCH models in finance, providing theoretical development and extensive empirical evidence.

Journal of Econometrics, 52(1-2), 5-59

Theoretical Development
Nelson, D. B.

1990

Statistical Properties
"Stationarity and persistence in the GARCH(1,1) model"

Fundamental analysis of GARCH stationarity conditions and persistence properties.

Econometric Theory, 6(3), 318-334

Bougerol, P., & Picard, N.

1992

Statistical Properties
"Stationarity of GARCH processes and of some nonnegative time series"

Rigorous mathematical treatment of GARCH stationarity conditions.

Journal of Econometrics, 52(1-2), 115-127

Lee, S. W., & Hansen, B. E.

1994

Estimation Theory
"Asymptotic theory for the GARCH(1,1) quasi-maximum likelihood estimator"

Foundational asymptotic theory for GARCH quasi-maximum likelihood estimation.

Econometric Theory, 10(1), 29-52

Lumsdaine, R. L.

1996

Estimation Theory
"Consistency and asymptotic normality of the quasi-maximum likelihood estimator in IGARCH(1,1) and covariance stationary GARCH(1,1) models"

Key results on QMLE consistency and normality for both stationary and integrated GARCH models.

Econometrica, 64(3), 575-596

Ling, S., & McAleer, M.

2002

Statistical Properties
"Necessary and sufficient moment conditions for the GARCH(r,s) and asymmetric power GARCH(r,s) models"

Comprehensive analysis of moment conditions for general GARCH specifications.

Econometric Theory, 18(3), 722-729

Francq, C., & Zakoïan, J. M.

2004

Estimation Theory
"Maximum likelihood estimation of pure GARCH and ARMA-GARCH processes"

Advanced theoretical results on maximum likelihood estimation for GARCH models.

Bernoulli, 10(4), 605-637

GARCH Extensions and Variants
Nelson, D. B.

1991

Asymmetric Models
"Conditional heteroskedasticity in asset returns: A new approach"

Original EGARCH specification introducing asymmetric volatility modeling.

Econometrica, 59(2), 347-370

Glosten, L. R., Jagannathan, R., & Runkle, D. E.

1993

Asymmetric Models
"On the relation between the expected value and the volatility of the nominal excess return on stocks"

GJR-GARCH threshold model for asymmetric volatility responses.

Journal of Finance, 48(5), 1779-1801

Engle, R. F., & Kroner, K. F.

1995

Multivariate Extensions
"Multivariate simultaneous generalized ARCH"

BEKK representation for multivariate GARCH models.

Econometric Theory, 11(1), 122-150

Engle, R.

2002

Multivariate Extensions
"Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models"

DCC-GARCH framework for time-varying correlations.

Journal of Business & Economic Statistics, 20(3), 339-350

Ghysels, E., Santa-Clara, P., & Valkanov, R.

2006

High-Frequency Extensions
"Predicting volatility: Getting the most out of return data sampled at different frequencies"

MIDAS-GARCH approach for mixed-frequency volatility modeling.

Journal of Econometrics, 131(1-2), 59-95

Hansen, P. R., Huang, Z., & Shek, H. H.

2012

High-Frequency Extensions
"Realized GARCH: A joint model for returns and realized measures of volatility"

Realized GARCH framework incorporating high-frequency volatility measures.

Journal of Applied Econometrics, 27(6), 877-906

Practical Applications and Implementation
Campbell, J. Y., Lo, A. W., & MacKinlay, A. C.

1997

Empirical Finance
"The Econometrics of Financial Markets"

Comprehensive textbook covering econometric methods for financial markets including GARCH models.

Princeton University Press

Christoffersen, P. F.

2003

Risk Management
"Elements of Financial Risk Management"

Practical guide to financial risk management with extensive coverage of GARCH applications.

Academic Press

Jorion, P.

2007

Risk Management
"Value at Risk: The New Benchmark for Managing Financial Risk"

Industry standard reference for VaR methodology with GARCH-based volatility modeling.

McGraw-Hill, 3rd ed.

Tsay, R. S.

2010

Empirical Finance
"Analysis of Financial Time Series"

Authoritative textbook on financial time series analysis with comprehensive GARCH coverage.

Wiley, 3rd ed.

McNeil, A. J., Frey, R., & Embrechts, P.

2015

Risk Management
"Quantitative Risk Management: Concepts, Techniques and Tools"

Advanced treatment of quantitative risk management including GARCH-based risk modeling.

Princeton University Press, 2nd ed.

Brooks, C.

2019

Empirical Finance
"Introductory Econometrics for Finance"

Accessible introduction to financial econometrics with practical GARCH modeling examples.

Cambridge University Press, 4th ed.

Software and Implementation Resources
R Packages

rugarch: Comprehensive GARCH modeling
rmgarch: Multivariate GARCH models
fGarch: Basic GARCH implementations
MTS: Multivariate time series analysis

Python Libraries

arch: ARCH and GARCH models
statsmodels: Statistical models including GARCH
PyFlux: Time series modeling framework
scipy: Statistical computing foundation

Commercial Software

EViews: User-friendly GARCH estimation
MATLAB: Econometrics Toolbox
SAS: PROC AUTOREG and PROC VARMAX
Stata: Built-in GARCH commands