Main Content

egarch

EGARCH conditional variance time series model

Description

Use egarch to specify a univariate EGARCH (exponential generalized autoregressive conditional heteroscedastic) model. The egarch function returns an egarch object specifying the functional form of an EGARCH(P,Q) model, and stores its parameter values.

The key components of an egarch model include the:

  • GARCH polynomial, which is composed of lagged, logged conditional variances. The degree is denoted by P.

  • ARCH polynomial, which is composed of the magnitudes of lagged standardized innovations.

  • Leverage polynomial, which is composed of lagged standardized innovations.

  • Maximum of the ARCH and leverage polynomial degrees, denoted by Q.

P is the maximum nonzero lag in the GARCH polynomial, and Q is the maximum nonzero lag in the ARCH and leverage polynomials. Other model components include an innovation mean model offset, a conditional variance model constant, and the innovations distribution.

All coefficients are unknown (NaN values) and estimable unless you specify their values using name-value pair argument syntax. To estimate models containing all or partially unknown parameter values given data, use estimate. For completely specified models (models in which all parameter values are known), simulate or forecast responses using simulate or forecast, respectively.

Creation

Description

example

Mdl = egarch creates a zero-degree conditional variance egarch object.

example

Mdl = egarch(P,Q) creates an EGARCH conditional variance model object (Mdl) with a GARCH polynomial with a degree of P, and ARCH and leverage polynomials each with a degree of Q. All polynomials contain all consecutive lags from 1 through their degrees, and all coefficients are NaN values.

This shorthand syntax enables you to create a template in which you specify the polynomial degrees explicitly. The model template is suited for unrestricted parameter estimation, that is, estimation without any parameter equality constraints. However, after you create a model, you can alter property values using dot notation.

example

Mdl = egarch(Name,Value) sets properties or additional options using name-value pair arguments. Enclose each name in quotes. For example, 'ARCHLags',[1 4],'ARCH',{0.2 0.3} specifies the two ARCH coefficients in ARCH at lags 1 and 4.

This longhand syntax enables you to create more flexible models.

Input Arguments

expand all

The shorthand syntax provides an easy way for you to create model templates that are suitable for unrestricted parameter estimation. For example, to create an EGARCH(1,2) model containing unknown parameter values, enter:

Mdl = egarch(1,2);
To impose equality constraints on parameter values during estimation, set the appropriate property values using dot notation.

GARCH polynomial degree, specified as a nonnegative integer. In the GARCH polynomial and at time t, MATLAB® includes all consecutive logged conditional variance terms from lag t – 1 through lag tP.

You can specify this argument using the egarch(P,Q) shorthand syntax only.

If P > 0, then you must specify Q as a positive integer.

Example: egarch(1,1)

Data Types: double

ARCH polynomial degree, specified as a nonnegative integer. In the ARCH polynomial and at time t, MATLAB includes all consecutive magnitudes of standardized innovation terms (for the ARCH polynomial) and all standardized innovation terms (for the leverage polynomial) from lag t – 1 through lag tQ.

You can specify this argument using the egarch(P,Q) shorthand syntax only.

If P > 0, then you must specify Q as a positive integer.

Example: egarch(1,1)

Data Types: double

Name-Value Arguments

Specify optional pairs of arguments as Name1=Value1,...,NameN=ValueN, where Name is the argument name and Value is the corresponding value. Name-value arguments must appear after other arguments, but the order of the pairs does not matter.

Before R2021a, use commas to separate each name and value, and enclose Name in quotes.

The longhand syntax enables you to create models in which some or all coefficients are known. During estimation, estimate imposes equality constraints on any known parameters.

Example: 'ARCHLags',[1 4],'ARCH',{NaN NaN} specifies an EGARCH(0,4) model and unknown, but nonzero, ARCH coefficient matrices at lags 1 and 4.

GARCH polynomial lags, specified as the comma-separated pair consisting of 'GARCHLags' and a numeric vector of unique positive integers.

GARCHLags(j) is the lag corresponding to the coefficient GARCH{j}. The lengths of GARCHLags and GARCH must be equal.

Assuming all GARCH coefficients (specified by the GARCH property) are positive or NaN values, max(GARCHLags) determines the value of the P property.

Example: 'GARCHLags',[1 4]

Data Types: double

ARCH polynomial lags, specified as the comma-separated pair consisting of 'ARCHLags' and a numeric vector of unique positive integers.

ARCHLags(j) is the lag corresponding to the coefficient ARCH{j}. The lengths of ARCHLags and ARCH must be equal.

Assuming all ARCH and leverage coefficients (specified by the ARCH and Leverage properties) are positive or NaN values, max([ARCHLags LeverageLags]) determines the value of the Q property.

Example: 'ARCHLags',[1 4]

Data Types: double

Leverage polynomial lags, specified as the comma-separated pair consisting of 'LeverageLags' and a numeric vector of unique positive integers.

LeverageLags(j) is the lag corresponding to the coefficient Leverage{j}. The lengths of LeverageLags and Leverage must be equal.

Assuming all ARCH and leverage coefficients (specified by the ARCH and Leverage properties) are positive or NaN values, max([ARCHLags LeverageLags]) determines the value of the Q property.

Example: 'LeverageLags',1:4

Data Types: double

Properties

expand all

You can set writable property values when you create the model object by using name-value argument syntax, or after you create the model object by using dot notation. For example, to create an EGARCH(1,1) model with unknown coefficients, and then specify a t innovation distribution with unknown degrees of freedom, enter:

Mdl = egarch('GARCHLags',1,'ARCHLags',1);
Mdl.Distribution = "t";

This property is read-only.

GARCH polynomial degree, specified as a nonnegative integer. P is the maximum lag in the GARCH polynomial with a coefficient that is positive or NaN. Lags that are less than P can have coefficients equal to 0.

P specifies the minimum number of presample conditional variances required to initialize the model.

If you use name-value pair arguments to create the model, then MATLAB implements one of these alternatives (assuming the coefficient of the largest lag is positive or NaN):

  • If you specify GARCHLags, then P is the largest specified lag.

  • If you specify GARCH, then P is the number of elements of the specified value. If you also specify GARCHLags, then egarch uses GARCHLags to determine P instead.

  • Otherwise, P is 0.

Data Types: double

This property is read-only.

Maximum degree of ARCH and leverage polynomials, specified as a nonnegative integer. Q is the maximum lag in the ARCH and leverage polynomials in the model. In either type of polynomial, lags that are less than Q can have coefficients equal to 0.

Q specifies the minimum number of presample innovations required to initiate the model.

If you use name-value pair arguments to create the model, then MATLAB implements one of these alternatives (assuming the coefficients of the largest lags in the ARCH and leverage polynomials are positive or NaN):

  • If you specify ARCHLags or LeverageLags, then Q is the maximum between the two specifications.

  • If you specify ARCH or Leverage, then Q is the maximum number of elements between the two specifications. If you also specify ARCHLags or LeverageLags, then egarch uses their values to determine Q instead.

  • Otherwise, Q is 0.

Data Types: double

Conditional variance model constant, specified as a numeric scalar or NaN value.

Data Types: double

GARCH polynomial coefficients, specified as a cell vector of positive scalars or NaN values.

  • If you specify GARCHLags, then the following conditions apply.

    • The lengths of GARCH and GARCHLags are equal.

    • GARCH{j} is the coefficient of lag GARCHLags(j).

    • By default, GARCH is a numel(GARCHLags)-by-1 cell vector of NaN values.

  • Otherwise, the following conditions apply.

    • The length of GARCH is P.

    • GARCH{j} is the coefficient of lag j.

    • By default, GARCH is a P-by-1 cell vector of NaN values.

The coefficients in GARCH correspond to coefficients in an underlying LagOp lag operator polynomial, and are subject to a near-zero tolerance exclusion test. If you set a coefficient to 1e–12 or below, egarch excludes that coefficient and its corresponding lag in GARCHLags from the model.

Data Types: cell

ARCH polynomial coefficients, specified as a cell vector of positive scalars or NaN values.

  • If you specify ARCHLags, then the following conditions apply.

    • The lengths of ARCH and ARCHLags are equal.

    • ARCH{j} is the coefficient of lag ARCHLags(j).

    • By default, ARCH is a Q-by-1 cell vector of NaN values. For more details, see the Q property.

  • Otherwise, the following conditions apply.

    • The length of ARCH is Q.

    • ARCH{j} is the coefficient of lag j.

    • By default, ARCH is a Q-by-1 cell vector of NaN values.

The coefficients in ARCH correspond to coefficients in an underlying LagOp lag operator polynomial, and are subject to a near-zero tolerance exclusion test. If you set a coefficient to 1e–12 or below, egarch excludes that coefficient and its corresponding lag in ARCHLags from the model.

Data Types: cell

Leverage polynomial coefficients, specified as a cell vector of numeric scalars or NaN values.

  • If you specify LeverageLags, then the following conditions apply.

    • The lengths of Leverage and LeverageLags are equal.

    • Leverage{j} is the coefficient of lag LeverageLags(j).

    • By default, Leverage is a Q-by-1 cell vector of NaN values. For more details, see the Q property.

  • Otherwise, the following conditions apply.

    • The length of Leverage is Q.

    • Leverage{j} is the coefficient of lag j.

    • By default, Leverage is a Q-by-1 cell vector of NaN values.

The coefficients in Leverage correspond to coefficients in an underlying LagOp lag operator polynomial, and are subject to a near-zero tolerance exclusion test. If you set a coefficient to 1e–12 or below, egarch excludes that coefficient and its corresponding lag in LeverageLags from the model.

Data Types: cell

This property is read-only.

The model unconditional variance, specified as a positive scalar.

The unconditional variance is

σε2=exp{κ(1i=1Pγi)}.

κ is the conditional variance model constant (Constant).

Data Types: double

Innovation mean model offset, or additive constant, specified as a numeric scalar or NaN value.

Data Types: double

Conditional probability distribution of the innovation process εt, specified as a string or structure array. egarch stores the value as a structure array.

DistributionStringStructure Array
Gaussian"Gaussian"struct('Name',"Gaussian")
Student’s t"t"struct('Name',"t",'DoF',DoF)

The 'DoF' field specifies the t distribution degrees of freedom parameter.

  • DoF > 2 or DoF = NaN.

  • DoF is estimable.

  • If you specify "t", DoF is NaN by default. You can change its value by using dot notation after you create the model. For example, Mdl.Distribution.DoF = 3.

  • If you supply a structure array to specify the Student's t distribution, then you must specify both the 'Name' and the 'DoF' fields.

Example: Distribution=struct('Name',"t",'DoF',10)

Model description, specified as a string scalar or character vector. egarch stores the value as a string scalar. The default value describes the parametric form of the model, for example "EGARCH(1,1) Conditional Variance Model (Gaussian Distribution)".

Data Types: string | char

Since R2023a

Response series name, specified as a string scalar or character vector. egarch stores the value as a string scalar.

Example: "StockReturn"

Data Types: string | char

Note

  • All NaN-valued model parameters, which include coefficients and the t-innovation-distribution degrees of freedom (if present), are estimable. When you pass the resulting egarch object and data to estimate, MATLAB estimates all NaN-valued parameters. During estimation, estimate treats known parameters as equality constraints, that is,estimate holds any known parameters fixed at their values.

  • Typically, the lags in the ARCH and leverage polynomials are the same, but their equality is not a requirement. Differing polynomials occur when:

    • Either ARCH{Q} or Leverage{Q} meets the near-zero exclusion tolerance. In this case, MATLAB excludes the corresponding lag from the polynomial.

    • You specify polynomials of differing lengths by specifying ARCHLags or LeverageLags, or by setting the ARCH or Leverage property.

    In either case, Q is the maximum lag between the two polynomials.

Object Functions

estimateFit conditional variance model to data
filterFilter disturbances through conditional variance model
forecastForecast conditional variances from conditional variance models
inferInfer conditional variances of conditional variance models
simulateMonte Carlo simulation of conditional variance models
summarizeDisplay estimation results of conditional variance model

Examples

collapse all

Create a default egarch model object and specify its parameter values using dot notation.

Create an EGARCH(0,0) model.

Mdl = egarch
Mdl = 
  egarch with properties:

     Description: "EGARCH(0,0) Conditional Variance Model (Gaussian Distribution)"
      SeriesName: "Y"
    Distribution: Name = "Gaussian"
               P: 0
               Q: 0
        Constant: NaN
           GARCH: {}
            ARCH: {}
        Leverage: {}
          Offset: 0

Mdl is an egarch model. It contains an unknown constant, its offset is 0, and the innovation distribution is 'Gaussian'. The model does not have GARCH, ARCH, or leverage polynomials.

Specify two unknown ARCH and leverage coefficients for lags one and two using dot notation.

Mdl.ARCH = {NaN NaN};
Mdl.Leverage = {NaN NaN};
Mdl
Mdl = 
  egarch with properties:

     Description: "EGARCH(0,2) Conditional Variance Model (Gaussian Distribution)"
      SeriesName: "Y"
    Distribution: Name = "Gaussian"
               P: 0
               Q: 2
        Constant: NaN
           GARCH: {}
            ARCH: {NaN NaN} at lags [1 2]
        Leverage: {NaN NaN} at lags [1 2]
          Offset: 0

The Q, ARCH, and Leverage properties update to 2, {NaN NaN}, {NaN NaN}, respectively. The two ARCH and leverage coefficients are associated with lags 1 and 2.

Create an egarch model object using the shorthand notation egarch(P,Q), where P is the degree of the GARCH polynomial and Q is the degree of the ARCH and leverage polynomial.

Create an EGARCH(3,2) model.

Mdl = egarch(3,2)
Mdl = 
  egarch with properties:

     Description: "EGARCH(3,2) Conditional Variance Model (Gaussian Distribution)"
      SeriesName: "Y"
    Distribution: Name = "Gaussian"
               P: 3
               Q: 2
        Constant: NaN
           GARCH: {NaN NaN NaN} at lags [1 2 3]
            ARCH: {NaN NaN} at lags [1 2]
        Leverage: {NaN NaN} at lags [1 2]
          Offset: 0

Mdl is an egarch model object. All properties of Mdl, except P, Q, and Distribution, are NaN values. By default, the software:

  • Includes a conditional variance model constant

  • Excludes a conditional mean model offset (i.e., the offset is 0)

  • Includes all lag terms in the GARCH polynomial up to lag P

  • Includes all lag terms in the ARCH and leverage polynomials up to lag Q

Mdl specifies only the functional form of an EGARCH model. Because it contains unknown parameter values, you can pass Mdl and time-series data to estimate to estimate the parameters.

Create an egarch model object using name-value pair arguments.

Specify an EGARCH(1,1) model. By default, the conditional mean model offset is zero. Specify that the offset is NaN. Include a leverage term.

Mdl = egarch('GARCHLags',1,'ARCHLags',1,'LeverageLags',1,'Offset',NaN)
Mdl = 
  egarch with properties:

     Description: "EGARCH(1,1) Conditional Variance Model with Offset (Gaussian Distribution)"
      SeriesName: "Y"
    Distribution: Name = "Gaussian"
               P: 1
               Q: 1
        Constant: NaN
           GARCH: {NaN} at lag [1]
            ARCH: {NaN} at lag [1]
        Leverage: {NaN} at lag [1]
          Offset: NaN

Mdl is an egarch model object. The software sets all parameters to NaN, except P, Q, and Distribution.

Since Mdl contains NaN values, Mdl is appropriate for estimation only. Pass Mdl and time-series data to estimate.

Create an EGARCH(1,1) model with mean offset,

yt=0.5+εt,

where εt=σtzt,

σt2=0.0001+0.75logσt-12+0.1(|εt-1|σt-1-2π)-0.3εt-1σt-1+0.01εt-3σt-3,

and zt is an independent and identically distributed standard Gaussian process.

Mdl = egarch('Constant',0.0001,'GARCH',0.75,...
    'ARCH',0.1,'Offset',0.5,'Leverage',{-0.3 0 0.01})
Mdl = 
  egarch with properties:

     Description: "EGARCH(1,3) Conditional Variance Model with Offset (Gaussian Distribution)"
      SeriesName: "Y"
    Distribution: Name = "Gaussian"
               P: 1
               Q: 3
        Constant: 0.0001
           GARCH: {0.75} at lag [1]
            ARCH: {0.1} at lag [1]
        Leverage: {-0.3 0.01} at lags [1 3]
          Offset: 0.5

egarch assigns default values to any properties you do not specify with name-value pair arguments. An alternative way to specify the leverage component is 'Leverage',{-0.3 0.01},'LeverageLags',[1 3].

Access the properties of a created egarch model object using dot notation.

Create an egarch model object.

Mdl = egarch(3,2)
Mdl = 
  egarch with properties:

     Description: "EGARCH(3,2) Conditional Variance Model (Gaussian Distribution)"
      SeriesName: "Y"
    Distribution: Name = "Gaussian"
               P: 3
               Q: 2
        Constant: NaN
           GARCH: {NaN NaN NaN} at lags [1 2 3]
            ARCH: {NaN NaN} at lags [1 2]
        Leverage: {NaN NaN} at lags [1 2]
          Offset: 0

Remove the second GARCH term from the model. That is, specify that the GARCH coefficient of the second lagged conditional variance is 0.

Mdl.GARCH{2} = 0
Mdl = 
  egarch with properties:

     Description: "EGARCH(3,2) Conditional Variance Model (Gaussian Distribution)"
      SeriesName: "Y"
    Distribution: Name = "Gaussian"
               P: 3
               Q: 2
        Constant: NaN
           GARCH: {NaN NaN} at lags [1 3]
            ARCH: {NaN NaN} at lags [1 2]
        Leverage: {NaN NaN} at lags [1 2]
          Offset: 0

The GARCH polynomial has two unknown parameters corresponding to lags 1 and 3.

Display the distribution of the disturbances.

Mdl.Distribution
ans = struct with fields:
    Name: "Gaussian"

The disturbances are Gaussian with mean 0 and variance 1.

Specify that the underlying disturbances have a t distribution with five degrees of freedom.

Mdl.Distribution = struct('Name','t','DoF',5)
Mdl = 
  egarch with properties:

     Description: "EGARCH(3,2) Conditional Variance Model (t Distribution)"
      SeriesName: "Y"
    Distribution: Name = "t", DoF = 5
               P: 3
               Q: 2
        Constant: NaN
           GARCH: {NaN NaN} at lags [1 3]
            ARCH: {NaN NaN} at lags [1 2]
        Leverage: {NaN NaN} at lags [1 2]
          Offset: 0

Specify that the ARCH coefficients are 0.2 for the first lag and 0.1 for the second lag.

Mdl.ARCH = {0.2 0.1}
Mdl = 
  egarch with properties:

     Description: "EGARCH(3,2) Conditional Variance Model (t Distribution)"
      SeriesName: "Y"
    Distribution: Name = "t", DoF = 5
               P: 3
               Q: 2
        Constant: NaN
           GARCH: {NaN NaN} at lags [1 3]
            ARCH: {0.2 0.1} at lags [1 2]
        Leverage: {NaN NaN} at lags [1 2]
          Offset: 0

To estimate the remaining parameters, you can pass Mdl and your data to estimate and use the specified parameters as equality constraints. Or, you can specify the rest of the parameter values, and then simulate or forecast conditional variances from the GARCH model by passing the fully specified model to simulate or forecast, respectively.

Fit an EGARCH model to an annual time series of Danish nominal stock returns from 1922-1999.

Load the Data_Danish data set. Plot the nominal returns (RN).

load Data_Danish;
nr = DataTable.RN;

figure;
plot(dates,nr);
hold on;
plot([dates(1) dates(end)],[0 0],'r:'); % Plot y = 0
hold off;
title('Danish Nominal Stock Returns');
ylabel('Nominal return (%)');
xlabel('Year');

Figure contains an axes object. The axes object with title Danish Nominal Stock Returns, xlabel Year, ylabel Nominal return (%) contains 2 objects of type line.

The nominal return series seems to have a nonzero conditional mean offset and seems to exhibit volatility clustering. That is, the variability is smaller for earlier years than it is for later years. For this example, assume that an EGARCH(1,1) model is appropriate for this series.

Create an EGARCH(1,1) model. The conditional mean offset is zero by default. To estimate the offset, specify that it is NaN. Include a leverage lag.

Mdl = egarch('GARCHLags',1,'ARCHLags',1,'LeverageLags',1,'Offset',NaN);

Fit the EGARCH(1,1) model to the data.

EstMdl = estimate(Mdl,nr);
 
    EGARCH(1,1) Conditional Variance Model with Offset (Gaussian Distribution):
 
                     Value      StandardError    TStatistic     PValue  
                   _________    _____________    __________    _________

    Constant        -0.62723       0.74401        -0.84304        0.3992
    GARCH{1}         0.77419       0.23628          3.2766     0.0010507
    ARCH{1}          0.38636       0.37361          1.0341       0.30107
    Leverage{1}    -0.002499       0.19222       -0.013001       0.98963
    Offset           0.10325      0.037727          2.7368     0.0062047

EstMdl is a fully specified egarch model object. That is, it does not contain NaN values. You can assess the adequacy of the model by generating residuals using infer, and then analyzing them.

To simulate conditional variances or responses, pass EstMdl to simulate.

To forecast innovations, pass EstMdl to forecast.

Simulate conditional variance or response paths from a fully specified egarch model object. That is, simulate from an estimated egarch model or a known egarch model in which you specify all parameter values.

Load the Data_Danish data set.

load Data_Danish;
rn = DataTable.RN;

Create an EGARCH(1,1) model with an unknown conditional mean offset. Fit the model to the annual, nominal return series. Include a leverage term.

Mdl = egarch('GARCHLags',1,'ARCHLags',1,'LeverageLags',1,'Offset',NaN);
EstMdl = estimate(Mdl,rn);
 
    EGARCH(1,1) Conditional Variance Model with Offset (Gaussian Distribution):
 
                     Value       StandardError    TStatistic     PValue  
                   __________    _____________    __________    _________

    Constant         -0.62723       0.74401        -0.84304        0.3992
    GARCH{1}          0.77419       0.23628          3.2766     0.0010508
    ARCH{1}           0.38636       0.37361          1.0341       0.30107
    Leverage{1}    -0.0024989       0.19222          -0.013       0.98963
    Offset            0.10325      0.037727          2.7368     0.0062047

Simulate 100 paths of conditional variances and responses from the estimated EGARCH model.

numObs = numel(rn); % Sample size (T)
numPaths = 100;     % Number of paths to simulate
rng(1);             % For reproducibility
[VSim,YSim] = simulate(EstMdl,numObs,'NumPaths',numPaths);

VSim and YSim are T-by- numPaths matrices. Rows correspond to a sample period, and columns correspond to a simulated path.

Plot the average and the 97.5% and 2.5% percentiles of the simulate paths. Compare the simulation statistics to the original data.

VSimBar = mean(VSim,2);
VSimCI = quantile(VSim,[0.025 0.975],2);
YSimBar = mean(YSim,2);
YSimCI = quantile(YSim,[0.025 0.975],2);

figure;
subplot(2,1,1);
h1 = plot(dates,VSim,'Color',0.8*ones(1,3));
hold on;
h2 = plot(dates,VSimBar,'k--','LineWidth',2);
h3 = plot(dates,VSimCI,'r--','LineWidth',2);
hold off;
title('Simulated Conditional Variances');
ylabel('Cond. var.');
xlabel('Year');

subplot(2,1,2);
h1 = plot(dates,YSim,'Color',0.8*ones(1,3));
hold on;
h2 = plot(dates,YSimBar,'k--','LineWidth',2);
h3 = plot(dates,YSimCI,'r--','LineWidth',2);
hold off;
title('Simulated Nominal Returns');
ylabel('Nominal return (%)');
xlabel('Year');
legend([h1(1) h2 h3(1)],{'Simulated path' 'Mean' 'Confidence bounds'},...
    'FontSize',7,'Location','NorthWest');

Figure contains 2 axes objects. Axes object 1 with title Simulated Conditional Variances, xlabel Year, ylabel Cond. var. contains 103 objects of type line. Axes object 2 with title Simulated Nominal Returns, xlabel Year, ylabel Nominal return (%) contains 103 objects of type line. These objects represent Simulated path, Mean, Confidence bounds.

Forecast conditional variances from a fully specified egarch model object. That is, forecast from an estimated egarch model or a known egarch model in which you specify all parameter values. The example follows from Estimate EGARCH Model.

Load the Data_Danish data set.

load Data_Danish;
nr = DataTable.RN;

Create an EGARCH(1,1) model with an unknown conditional mean offset and include a leverage term. Fit the model to the annual nominal return series.

Mdl = egarch('GARCHLags',1,'ARCHLags',1,'LeverageLags',1,'Offset',NaN);
EstMdl = estimate(Mdl,nr);
 
    EGARCH(1,1) Conditional Variance Model with Offset (Gaussian Distribution):
 
                     Value       StandardError    TStatistic     PValue  
                   __________    _____________    __________    _________

    Constant         -0.62723       0.74401        -0.84304        0.3992
    GARCH{1}          0.77419       0.23628          3.2766     0.0010508
    ARCH{1}           0.38636       0.37361          1.0341       0.30107
    Leverage{1}    -0.0024989       0.19222          -0.013       0.98963
    Offset            0.10325      0.037727          2.7368     0.0062047

Forecast the conditional variance of the nominal return series 10 years into the future using the estimated EGARCH model. Specify the entire returns series as presample observations. The software infers presample conditional variances using the presample observations and the model.

numPeriods = 10;
vF = forecast(EstMdl,numPeriods,nr);

Plot the forecasted conditional variances of the nominal returns. Compare the forecasts to the observed conditional variances.

v = infer(EstMdl,nr);

figure;
plot(dates,v,'k:','LineWidth',2);
hold on;
plot(dates(end):dates(end) + 10,[v(end);vF],'r','LineWidth',2);
title('Forecasted Conditional Variances of Nominal Returns');
ylabel('Conditional variances');
xlabel('Year');
legend({'Estimation sample cond. var.','Forecasted cond. var.'},...
    'Location','Best');

Figure contains an axes object. The axes object with title Forecasted Conditional Variances of Nominal Returns, xlabel Year, ylabel Conditional variances contains 2 objects of type line. These objects represent Estimation sample cond. var., Forecasted cond. var..

More About

expand all

Tips

  • You can specify an egarch model as part of a composition of conditional mean and variance models. For details, see arima.

  • An EGARCH(1,1) specification is complex enough for most applications. Typically in these models, the GARCH and ARCH coefficients are positive, and the leverage coefficients are negative. If you get these signs, then large unanticipated downward shocks increase the variance. If you get signs opposite to those signs that are expected, you can encounter difficulties inferring volatility sequences and forecasting. A negative ARCH coefficient is problematic. In this case, an EGARCH model might not be the best choice for your application.

References

[1] Tsay, R. S. Analysis of Financial Time Series. 3rd ed. Hoboken, NJ: John Wiley & Sons, Inc., 2010.

Version History

Introduced in R2012a

expand all