Documentation

h = waldtest(r,R,EstCov) returns a logical value (h) with the rejection decision from conducting a Wald test of model specification.

waldtest constructs the test statistic using the restriction function and its Jacobian, and the value of the unrestricted model covariance estimator, all evaluated at the unrestricted parameter estimates (r, R, and EstCov, respectively).

If any input argument is a cell vector of length k > 1, then the other input arguments must be cell arrays of length k. waldtest(r,R,EstCov) treats each cell as a separate, independent test, and returns a vector of rejection decisions.
If any input argument is a row vector, then the software returns output arguments as row vectors.

h = waldtest(r,R,EstCov,alpha) returns the rejection decision of the Wald test conducted at significance level alpha.

[h,pValue] = waldtest(___) returns the rejection decision and p-value (pValue) for the hypothesis test, using any of the input arguments in the previous syntaxes.

[h,pValue,stat,cValue] = waldtest(___) additionally returns the test statistic (stat) and critical value (cValue) for the hypothesis test.

Examples

Assess Model Specifications Using the Wald Test

Check for significant lag effects in a time series regression model.

Load the U.S. GDP data set.

load Data_GDP

Plot the GDP against time.

plot(dates,Data)
datetick

The series seems to increase exponentially.

Transform the data using the natural logarithm.

logGDP = log(Data);

logGDP is increasing in time, so assume that there is a significant lag 1 effect. To use the Wald test to check if there is a significant lag 2 effect, you need the:

Estimated coefficients of the unrestricted model
Restriction function evaluated at the unrestricted model coefficient values
Jacobian of the restriction function evaluated at the unrestricted model coefficient values
Estimated, unrestricted parameter covariance matrix.

The unrestricted model is

$y_{t} = β_{0} + β_{1} y_{t - 1} + β_{2} y_{t - 2} + ε_{t} .$

Estimate the coefficients of the unrestricted model.

LagLGDP = lagmatrix(logGDP,1:2);
UMdl = fitlm(table(LagLGDP(:,1),LagLGDP(:,2),logGDP));

UMdl is a fitted LinearModel model. It contains, among other things, the fitted coefficients of the unrestricted model.

The restriction is $β_{2} = 0$ . Therefore, the restriction function (r) and Jacobian (R) are:

$r = β_{2}$
$R = [\begin{array}{c} 0 & 0 & 1 \end{array}]$

Specify r, R, and the estimated, unrestricted parameter covariance matrix.

r = UMdl.Coefficients.Estimate(3);
R = [0 0 1];
EstParamCov = UMdl.CoefficientCovariance;

Test for a significant lag 2 effect using the Wald test.

[h,pValue] = waldtest(r,R,EstParamCov)

h = logical
   1

pValue = 1.2521e-07

h = 1 indicates that the null, restricted hypothesis ( $β_{2} = 0$ ) should be rejected in favor of the alternative, unrestricted hypothesis. pValue is quite small, which suggests that there is strong evidence for this result.

Assess Conditional Heteroscedasticity Using the Wald Test

Test whether there are significant ARCH effects in a simulated response series using waldtest.

Suppose that the model for the simulated data is AR(1) with an ARCH(1) variance. Symbolically, the model is

$y_{t} = 0.9 y_{t - 1} + ε_{t},$

where

$ε_{t} = w_{t} \sqrt{h_{t}}$
$h_{t} = 1 + 0.5 ε_{t - 1}^{2}$
$w_{t}$ is Gaussian with mean 0 and variance 1.

Specify the model for the simulated data.

VarMdl = garch('ARCH',0.5,'Constant',1);
Mdl = arima('Constant',0,'Variance',VarMdl,'AR',0.9);

Mdl is a fully specified AR(1) model with an ARCH(1) variance.

Simulate presample and effective sample responses from Mdl.

T = 100;
rng(1);  % For reproducibility
n = 2;   % Number of presample observations required for the Jacobian
[y,epsilon,condVariance] = simulate(Mdl,T + n);

psI = 1:n;             % Presample indices
esI = (n + 1):(T + n); % Estimation sample indices

epsilon is the random path of innovations from VarMdl. The software filters epsilon through Mdl to yield the random response path y.

Specify the unrestricted model assuming that the conditional mean model is

$y_{t} = c + ϕ_{1} y_{t - 1} + ε_{t},$

where $h_{t} = α_{0} + α_{1} ε_{t - 1}^{2}$ . Fit the simulated data (y) to the unrestricted model using the presample observations.

UVarMdl = garch(0,1);
UMdl = arima('ARLags',1,'Variance',UVarMdl);
[UEstMdl,UEstParamCov] = estimate(UMdl,y(esI),'Y0',y(psI),...
    'E0',epsilon(psI),'V0',condVariance(psI),'Display','off');

UEstMdl is the fitted, unrestricted model, and UEstParamCov is the estimated parameter covariance of the unrestricted model parameters.

The null hypothesis is that $α_{1} = 0$ , i.e., the restricted model is AR(1) with Gaussian innovations that have mean 0 and constant variance. Therefore, the restriction function is $r (θ) = α_{1}$ , where $θ = [c, ϕ_{1}, α_{0}, α_{1}]^{'}$ . The components of the Wald test are:

The restriction function evaluated at the unrestricted parameter estimates is $r = {α_{}^{ˆ}}_{1}$ .
The Jacobian of r evaluated at the unrestricted model parameters is $R = [\begin{array}{c} 0 & 0 & 0 & 1 \end{array}]$ .
The unrestricted model estimated parameter covariance matrix is UEstParamCov.

Specify r and R.

r = UEstMdl.Variance.ARCH{1};
R = [0, 0, 0, 1];

Test the null hypothesis that $α_{1} = 0$ at the 1% significance level using waldtest.

[h,pValue,stat,cValue] = waldtest(r,R,UEstParamCov,0.01)

h = logical
   0

pValue = 0.0549

stat = 3.6846

cValue = 6.6349

h = 0 indicates that the null, restricted model should not be rejected in favor of the alternative, unrestricted model. This result is consistent with the model for the simulated data.

Test Among Multiple Nested Model Specifications

Assess model specifications by testing down among multiple restricted models using simulated data. The true model is the ARMA(2,1)

$y_{t} = 3 + 0.9 y_{t - 1} - 0.5 y_{t - 2} + ε_{t} + 0.7 ε_{t - 1},$

where $ε_{t}$ is Gaussian with mean 0 and variance 1.

Specify the true ARMA(2,1) model, and simulate 100 response values.

TrueMdl = arima('AR',{0.9,-0.5},'MA',0.7,...
    'Constant',3,'Variance',1);
T = 100;
rng(1); % For reproducibility 
y = simulate(TrueMdl,T);

Specify the unrestricted model and the names of the candidate models for testing down.

UMdl = arima(2,0,2);
RMdlNames = {'ARMA(2,1)','AR(2)','ARMA(1,2)','ARMA(1,1)',...
    'AR(1)','MA(2)','MA(1)'};

UMdl is the unrestricted, ARMA(2,2) model. RMdlNames is a cell array of strings containing the names of the restricted models.

Fit the unrestricted model to the simulated data.

[UEstMdl,UEstParamCov] = estimate(UMdl,y,'Display','off');

UEstMdl is the fitted, unrestricted model, and UEstParamCov is the estimated parameter covariance matrix.

The unrestricted model has six parameters. To construct the restriction function and its Jacobian, you must know the order of the parameters in UEstParamCov. For this arima model, the order is $[c, ϕ_{1}, ϕ_{2}, θ_{1} θ_{2}, σ^{2}]$ .

Each candidate model corresponds to a restriction function. Put the restriction function vectors into separate cells of a cell vector.

rf1 = UEstMdl.MA{2};                          % ARMA(2,1)
rf2 = cell2mat(UEstMdl.MA)';                  % AR(2)
rf3 = UEstMdl.AR{2};                          % ARMA(1,2)
rf4 = [UEstMdl.AR{2};UEstMdl.MA{2}]';         % ARMA(1,1)
rf5 = [UEstMdl.AR{2};cell2mat(UEstMdl.MA)'];  % AR(1)
rf6 = cell2mat(UEstMdl.AR)';                  % MA(2)
rf7 = [cell2mat(UEstMdl.AR)';UEstMdl.MA{2}];  % MA(1)
r = {rf1;rf2;rf3;rf4;rf5;rf6;rf7};

r is a 7-by-1 cell vector of vectors corresponding to the restriction function for the candidate models.

Put the Jacobian of each restriction function into separate, corresponding cells of a cell vector. The order of the elements in the Jacobian must correspond to the order of the elements in UEstParamCov.

J1 = [0 0 0 0 1 0];                           % ARMA(2,1)   
J2 = [0 0 0 1 0 0; 0 0 0 0 1 0];              % AR(2)      
J3 = [0 1 0 0 0 0];                           % ARMA(1,2)  
J4 = [0 1 0 0 0 0; 0 0 0 0 1 0];              % ARMA(1,1)  
J5 = [0 1 0 0 0 0; 0 0 0 1 0 0; 0 0 0 0 1 0]; % AR(1)      
J6 = [1 0 0 0 0 0; 0 1 0 0 0 0];              % MA(2)      
J7 = [1 0 0 0 0 0; 0 1 0 0 0 0; 0 0 0 0 1 0]; % MA(1)      
R = {J1;J2;J3;J4;J5;J6;J7};

R is a 7-by-1 cell vector of vectors corresponding to the restriction function for the candidate models.

Put the estimated parameter covariance matrix in each cell of a 7-by-1 cell vector.

EstCov = cell(7,1); % Preallocate
for j = 1:length(EstCov)
    EstCov{j} = UEstParamCov;
end

Apply the Wald test at a 1% significance level to find the appropriate, restricted model specifications.

alpha = .01;
h = waldtest(r,R,EstCov,alpha);
RestrictedModels = RMdlNames(~h)

RestrictedModels = 1x5 cell
    {'ARMA(2,1)'}    {'ARMA(1,2)'}    {'ARMA(1,1)'}    {'MA(2)'}    {'MA(1)'}

RestrictedModels lists the most appropriate restricted models.

You can test down again, but use ARMA(2,1) as the unrestricted model. In this case, you must remove MA(2) from the possible restricted models.

Conduct a Wald Test Using Nonlinear Restriction Functions

Test whether the parameters of a nested model have a nonlinear relationship.

Load the Deutschmark/British Pound bilateral spot exchange rate data set.

load Data_MarkPound

The data set (Data) contains a time series of prices.

Convert the prices to returns, and plot the return series.

returns = price2ret(Data);

figure
plot(returns)
axis tight
ylabel('Returns')
xlabel('Days, 02Jan1984 - 31Dec1991')
title('{\bf Deutschmark/British Pound Bilateral Spot Exchange Rate}')

The returns series shows signs of heteroscedasticity.

Suppose that a GARCH(1,1) model is an appropriate model for the data. Fit a GARCH(1,1) model to the data including a constant.

Mdl = garch(1,1);
[EstMdl,EstParamCov] = estimate(Mdl,returns);

 
    GARCH(1,1) Conditional Variance Model (Gaussian Distribution):
 
                  Value       StandardError    TStatistic      PValue  
                __________    _____________    __________    __________

    Constant    1.0538e-06     3.5053e-07        3.0063       0.0026443
    GARCH{1}       0.80653       0.012913        62.458               0
    ARCH{1}        0.15439       0.011578        13.335      1.4463e-40

g1 = EstMdl.GARCH{1};
a1 = EstMdl.ARCH{1};

g1 is the estimated GARCH effect, and a1 is the estimated ARCH effect.

The following might represent relationships between the GARCH and ARCH coefficients:

$γ_{1} α_{1} = 1$
$γ_{1} + α_{1} = 1$

where $γ_{1}$ is the GARCH effect and $α_{1}$ is the ARCH effect. Specify these relationships as the restriction function $r (θ) = 0$ , evaluated at the unrestricted model parameter estimates. This specification defines a nested, restricted model.

r = [g1*a1; g1+a1] - 1;

Specify the Jacobian of the restriction function vector.

R = [0, a1, g1;0, 1, 1];

Conduct a Wald test to assess whether there is sufficient evidence to reject the restricted model.

[h,pValue,stat,cValue] = waldtest(r,R,EstParamCov)

h = logical
   1

pValue = 0

stat = 1.4589e+04

cValue = 5.9915

h = 1 indicates that there is sufficient evidence to reject the restricted model in favor of the unrestricted model. pValue = 0 indicates that the evidence for rejecting the restricted model is strong.

Input Arguments

`r` — Restriction functions
scalar | vector | cell vector of scalars or vectors

Restriction functions corresponding to restricted models, specified as a scalar, vector, or cell vector of scalars or vectors.

If r is a q-vector or a singleton cell array containing a q-vector, then the software conducts one Wald test. q must be less than the number of unrestricted model parameters.
If r is a cell vector of length k > 1, and cell j contains a q_j-vector, j = 1,...,k, then the software conducts k independent Wald tests. Each q_j must be less than the number of unrestricted model parameters.

Data Types: double | cell

`R` — Restriction function Jacobians
row vector | matrix | cell vector of row vectors or matrices

Restriction function Jacobians, specified as a row vector, matrix, or cell vector of row vectors or matrices.

Suppose r₁,...,r_q are the q restriction functions, and the unrestricted model parameters are θ₁,...,θ_p. Then, the restriction function Jacobian is

$R = (\begin{matrix} \frac{\partial r_{1}}{\partial θ_{1}} & \dots & \frac{\partial r_{1}}{\partial θ_{p}} \\ ⋮ & ⋱ & ⋮ \\ \frac{\partial r_{q}}{\partial θ_{1}} & \dots & \frac{\partial r_{q}}{\partial θ_{p}} \end{matrix}) .$
If R is a q-by-p matrix or a singleton cell array containing a q-by-p matrix, then the software conducts one Wald test. q must be less than p, which is the number of unrestricted model parameters.
If R is a cell vector of length k > 1, and cell j contains a q_j-by-p_j matrix, j = 1,...,k, then the software conducts k independent Wald tests. Each q_j must be less than p_j, which is the number of unrestricted parameters in model j.

Data Types: double | cell

`EstCov` — Unrestricted model parameter covariance estimate
matrix | cell vector of matrices

Unrestricted model parameter covariance estimates, specified as a matrix or cell vector of matrices.

If EstCov is a p-by-p matrix or a singleton cell array containing a p-by-p matrix, then the software conducts one Wald test. p is the number of unrestricted model parameters.
If EstCov is a cell vector of length k > 1, and cell j contains a p_j-by-p_j matrix, j = 1,...,k, then the software conducts k independent Wald tests. Each p_j is the number of unrestricted parameters in model j.

Data Types: double | cell

`alpha` — Nominal significance levels
0.05 (default) | scalar | vector

Nominal significance levels for the hypothesis tests, specified as a scalar or vector.

Each element of alpha must be greater than 0 and less than 1.

When conducting k > 1 tests,

If alpha is a scalar, then the software expands it to a k-by-1 vector.
If alpha is a vector, then it must have length k.

Data Types: double

Output Arguments

`h` — Test rejection decisions
logical | vector of logicals

Test rejection decisions, returned as a logical value or vector of logical values with a length equal to the number of tests that the software conducts.

h = 1 indicates rejection of the null, restricted model in favor of the alternative, unrestricted model.
h = 0 indicates failure to reject the null, restricted model.

`pValue` — Test statistic p-values
scalar | vector

Test statistic p-values, returned as a scalar or vector with a length equal to the number of tests that the software conducts.

`stat` — Test statistics
scalar | vector

Test statistics, returned as a scalar or vector with a length equal to the number of tests that the software conducts.

`cValue` — Critical values
scalar | vector

Critical values determined by alpha, returned as a scalar or vector with a length equal to the number of tests that the software conducts.

More About