Tutorials › Econometrics (Graduate) › Wald Tests And Confidence Regions In GMM

Generalized Method of Moments · Tutorial 212 of 952

Wald Tests And Confidence Regions In GMM

Use the GMM sandwich covariance to test smooth restrictions and construct asymptotic confidence regions, including joint regions for multiple parameters.

Advanced 11 min read

What You'll Learn

  • Define a Wald test for smooth restrictions on a GMM parameter vector
  • Distinguish the asymptotic covariance of the estimator from that of its scaled limit
  • Calculate a Wald test for a linear contrast using the sandwich covariance
  • Test joint restrictions while accounting for covariance between estimates
  • Apply the delta method to a nonlinear restriction
  • Construct and interpret an asymptotic confidence region for a parameter vector

Why Wald Inference Uses the GMM Sandwich Covariance

A GMM estimate is often used to answer questions about combinations of parameters, not just about individual coefficients. For example, an analyst may want to test whether two coefficients sum to a specified value, whether several coefficients jointly equal a benchmark vector, or whether a nonlinear function of the coefficients satisfies a restriction. Wald inference turns these questions into comparisons between an estimated restriction and its hypothesized value, scaled by its estimated uncertainty.

The uncertainty calculation must account for the covariance among the GMM estimates. Two estimates may be correlated, so the uncertainty of a contrast or a joint hypothesis cannot generally be found by treating each coefficient separately. The GMM sandwich covariance provides the required joint covariance matrix.

Setup and Wald Statistic

Let \(\theta_0\) be a \(k\)-dimensional parameter vector and let \(\hat{\theta}\) be a GMM estimator based on \(n\) observations. Suppose the estimator is consistent and asymptotically normal:

$$ \sqrt{n}(\hat{\theta}-\theta_0)\xrightarrow{d}N(0,V). $$

Here \(V\) is the asymptotic covariance matrix of the scaled estimator. In standard GMM notation, if \(G\) is the Jacobian of the population moments, \(S\) is their covariance matrix, and \(W\) is the limiting GMM weight matrix, the sandwich covariance is

$$ V=(G'WG)^{-1}G'WSWG(G'WG)^{-1}, \qquad \widehat{\operatorname{Var}}(\hat{\theta})=\frac{\hat V}{n}. $$

The estimate \(\hat V\) estimates the covariance of \(\sqrt n(\hat{\theta}-\theta_0)\); the estimated covariance of \(\hat{\theta}\) itself is \(\hat V/n\). This distinction determines where the factor \(n\) belongs in a Wald statistic. We assume the GMM estimator satisfies the regularity and identification conditions for the stated asymptotic normality, \(\hat V\) consistently estimates \(V\), and the covariance matrices used below are nonsingular.

Definition: Let \(r(\theta)\) be a differentiable \(q\)-dimensional vector of restrictions, and consider \(H_0:r(\theta_0)=r_0\). Write \(R(\theta)=\partial r(\theta)/\partial\theta'\). If \(R(\theta_0)\) has rank \(q\), the Wald statistic is $$ W_n=n\{r(\hat{\theta})-r_0\}'\{R(\hat{\theta})\hat V R(\hat{\theta})'\}^{-1}\{r(\hat{\theta})-r_0\}. $$ Under \(H_0\) and the stated regularity conditions, \(W_n\) converges in distribution to a chi-square random variable with \(q\) degrees of freedom.

The rank condition means the \(q\) restrictions provide \(q\) locally independent pieces of information. The matrix \(R\hat V R'\) estimates the covariance of the scaled restriction error \(\sqrt n\{r(\hat{\theta})-r(\theta_0)\}\). A large Wald statistic indicates that the estimated restrictions are far from their hypothesized values relative to their uncertainty.

Linear Contrasts

For a linear restriction \(r(\theta)=C\theta\), where \(C\) has \(q\) independent rows, the Jacobian is simply \(C\). The test of \(H_0:C\theta_0=r_0\) becomes

$$ W_n=n(C\hat{\theta}-r_0)'(C\hat V C')^{-1}(C\hat{\theta}-r_0). $$

For one restriction, this is the square of a standard-normal Wald statistic. For several restrictions, it is a quadratic form that uses the full covariance among the estimated contrasts.

Worked Example: Testing a Linear Contrast

Suppose a hypothetical GMM analysis estimates \(\theta=(\theta_1,\theta_2)'\) as \(\hat{\theta}=(0.8,1.3)'\) from \(n=100\) observations. Its estimated sandwich covariance for the scaled estimator is

$$ \hat V= \begin{pmatrix} 0.16&0.04\\ 0.04&0.09 \end{pmatrix}. $$

Test \(H_0:\theta_1+\theta_2=2\) against the two-sided alternative. The contrast vector is \(C=(1,1)\), and the estimated contrast is \(C\hat{\theta}=0.8+1.3=2.1\). Its estimated asymptotic variance is

$$ C\hat V C'=0.16+0.09+2(0.04)=0.33. $$

Therefore its estimated standard error is \(\sqrt{0.33/100}=0.05745\), and the Wald statistic is

$$ W_n=\frac{(2.1-2)^2}{0.33/100} =\frac{0.01}{0.0033} =3.030. $$

Against a chi-square distribution with one degree of freedom, the p-value is approximately \(0.082\). At the \(5\%\) significance level, the statistic is below the critical value \(3.841\), so we do not reject \(H_0\). This result is not evidence that the null is true; it says the observed contrast is not sufficiently far from 2 relative to its estimated uncertainty at that level.

As a check, the 95% Wald confidence interval for the contrast is \(2.1\pm1.96(0.05745)\), or approximately \((1.987,2.213)\). Since 2 is in this interval, the interval and two-sided test give the same decision.

Joint Hypotheses and Confidence Regions

A joint hypothesis tests several restrictions together. For example, \(H_0:\theta_1=a_1,\theta_2=a_2\) has two restrictions and uses a chi-square reference distribution with two degrees of freedom. The off-diagonal entries of the covariance matrix matter: they describe how the estimates move together and affect the distance from the null vector.

The corresponding \(100(1-\alpha)\%\) Wald confidence region for the full parameter vector is the set of values \(\theta\) satisfying

$$ n(\hat{\theta}-\theta)'\hat V^{-1}(\hat{\theta}-\theta) \leq \chi^2_{k,1-\alpha}, $$

where \(\chi^2_{k,1-\alpha}\) is the \(1-\alpha\) quantile of a chi-square distribution with \(k\) degrees of freedom. The region is typically an ellipsoid centered at \(\hat{\theta}\). Its shape and orientation reflect the estimated variances and covariances. It is a joint region, not a collection of separate intervals with the same nominal coverage.

Worked Example: A Joint Test and Its Confidence Ellipse

Use the estimates and \(\hat V\) from the previous hypothetical calculation to test \(H_0:\theta_1=0.5,\theta_2=1.0\). The difference between the estimate and the null vector is \(d=(0.3,0.3)'\). The inverse sandwich matrix is

$$ \hat V^{-1}= \begin{pmatrix} 7.03125&-3.125\\ -3.125&12.5 \end{pmatrix}. $$

The joint Wald statistic is

$$ W_n=100d'\hat V^{-1}d =100(0.3,0.3) \begin{pmatrix} 7.03125&-3.125\\ -3.125&12.5 \end{pmatrix} \begin{pmatrix} 0.3\\0.3 \end{pmatrix} =119.5. $$

The reference distribution has two degrees of freedom. The p-value is approximately \(e^{-119.5/2}=1.1\times10^{-26}\), so the joint null is strongly rejected at conventional significance levels. This conclusion concerns the pair of restrictions together; it does not, by itself, identify a cause for either discrepancy.

The 95% confidence ellipse for \((\theta_1,\theta_2)'\) consists of values satisfying

$$ 100(\hat{\theta}-\theta)'\hat V^{-1}(\hat{\theta}-\theta)\leq\chi^2_{2,0.95}=5.991. $$

The tested null point \((0.5,1.0)'\) is outside this ellipse because its quadratic distance is \(119.5\), well above \(5.991\). For comparison, the marginal 95% Wald interval for \(\theta_1\) is \(0.8\pm1.96\sqrt{0.16/100}\), or approximately \((0.722,0.878)\). That one-coordinate interval is not a substitute for the joint ellipse.

Smooth Nonlinear Restrictions

Restrictions need not be linear. If \(r(\theta)\) is differentiable, a first-order expansion around the true value gives

$$ r(\hat{\theta})-r(\theta_0) =R(\theta_0)(\hat{\theta}-\theta_0)+o_p(n^{-1/2}). $$

Consequently, asymptotic normality of \(\hat{\theta}\) implies that \(\sqrt n\{r(\hat{\theta})-r(\theta_0)\}\) is asymptotically normal with covariance \(R(\theta_0)VR(\theta_0)'\). Replacing the unknown quantities by consistent estimates and forming the corresponding squared, covariance-scaled distance yields the chi-square Wald statistic. This is the delta-method justification for testing smooth restrictions.

Worked Example: Testing a Product Restriction

Suppose a hypothetical GMM estimate is \(\hat{\theta}=(1.10,0.95)'\), with \(n=100\) and

$$ \hat V= \begin{pmatrix} 0.04&0.01\\ 0.01&0.09 \end{pmatrix}. $$

Test the nonlinear restriction \(H_0:\theta_1\theta_2=1\). Define \(r(\theta)=\theta_1\theta_2\), so \(r(\hat{\theta})=1.10(0.95)=1.045\) and \(R(\hat{\theta})=(0.95,1.10)\). The estimated variance of the scaled restriction is

$$ R(\hat{\theta})\hat V R(\hat{\theta})' =(0.95)^2(0.04)+(1.10)^2(0.09)+2(0.95)(1.10)(0.01) =0.1659. $$

The Wald statistic for a single restriction is

$$ W_n=\frac{100(1.045-1)^2}{0.1659} =\frac{0.2025}{0.1659} =1.221. $$

The chi-square reference distribution has one degree of freedom, giving a p-value of approximately \(0.269\). At the \(5\%\) level, we do not reject the product restriction. The calculation uses the gradient of the product at the estimated parameter and includes the covariance term; omitting that term would give a different, generally incorrect standard error.

Intervals and Regions for Smooth Functions

A scalar smooth function \(\psi(\theta)\) can also be reported with a delta-method Wald interval. Let \(a(\theta)=\partial\psi(\theta)/\partial\theta\), written as a column vector. Its estimated standard error is

$$ \widehat{SE}\{\psi(\hat{\theta})\} =\sqrt{\frac{a(\hat{\theta})'\hat V a(\hat{\theta})}{n}}. $$

An approximate 95% interval is \(\psi(\hat{\theta})\pm1.96\,\widehat{SE}\{\psi(\hat{\theta})\}\). For a vector-valued smooth function with a full-row-rank Jacobian, the analogous Wald confidence region is obtained by replacing the scalar standard error with its estimated covariance matrix and using a chi-square critical value with degrees of freedom equal to the number of function components.

These confidence regions are asymptotic: their coverage approaches the stated level as the sample size grows under the model and regularity conditions. They are not guarantees that a particular realized region contains the fixed parameter, nor do they measure the probability that a null hypothesis is true.

Common Pitfalls and Practical Checks

  • Using the wrong covariance scale. If \(\hat V\) estimates the covariance of the scaled estimator, the covariance of \(\hat{\theta}\) is \(\hat V/n\). Equivalently, use \(n\) in the quadratic Wald formula. Do not divide by \(n\) and then multiply by \(n\) a second time.
  • Ignoring covariance. For a contrast or joint test, using only the diagonal variances discards dependence between estimates. Use the full matrix in the quadratic form.
  • Replacing a joint test with separate tests. Testing each coefficient separately does not generally test the joint null at the same error rate or account for their dependence. Use the joint Wald statistic with the appropriate number of degrees of freedom.
  • Using a singular restriction covariance. The inverse in the statistic requires \(R\hat V R'\) to be nonsingular. Redundant restrictions or a rank-deficient Jacobian violate the stated setup; the restriction system must be expressed in terms of independent restrictions.
  • Reading a p-value as the probability the null is true. A p-value is calculated under the null model and describes how extreme the statistic is relative to its reference distribution. It is not the probability that the null holds.

Wald inference is useful because the same principle applies to linear contrasts, nonlinear functions, and joint restrictions. Its reliability depends on a suitable asymptotic approximation, a consistently estimated sandwich covariance, and restrictions that are locally independent. In finite samples, the resulting tests and regions can be only approximate, so the assumptions and the scale of the covariance estimate should be checked carefully.

Check Your Understanding

Use the definitions and calculations above to answer these questions.

  1. In the GMM sandwich notation, what does \(\hat V\) estimate, and how does it differ from the estimated covariance of \(\hat{\theta}\)?
  2. A researcher tests one linear contrast using a Wald statistic of \(4.2\). Which chi-square degrees of freedom apply, and what is the approximate \(5\%\) decision?
  3. Why does a joint Wald test use the off-diagonal entries of the covariance matrix?
  4. For a differentiable scalar restriction \(r(\theta)=0\), what role does its Jacobian play in the Wald statistic?
  5. Explain why a 95% joint confidence ellipse is not the same object as two separate 95% marginal confidence intervals.
Clinical Biostats