Why a Robust Weak-Instrument Diagnostic Is Needed
A first-stage statistic is used to assess how much information the excluded instruments provide about the endogenous regressor or regressors. Under homoskedastic errors, conventional measures such as the Cragg-Donald statistic can be used in weak-identification analysis. If first-stage errors are heteroskedastic, however, the conventional covariance calculation may not be appropriate. The Kleibergen-Paap rk Wald F statistic (often abbreviated KP rk Wald F) is a heteroskedasticity-robust weak-identification diagnostic intended for this setting.
Its robust covariance estimate changes how first-stage uncertainty is measured; it does not make instruments valid, prove that identification is strong, or supply universal critical values. The introductory material on heteroskedasticity-robust first-stage tests establishes how a robust covariance estimate can be used for a joint test. Here the focus is how that idea relates to the KP rk Wald F and how far its interpretation can safely go.
Model, Assumptions, and Interpretation
Write the first stage for one endogenous regressor as
where \(D_i\) is the endogenous regressor, \(X_i\) contains included exogenous covariates, \(Z_i\) contains the \(L\) excluded instruments, and \(\boldsymbol{\pi}\) is the vector of first-stage coefficients on those instruments. The robust joint Wald test concerns \(H_0:\boldsymbol{\pi}=\boldsymbol{0}\). The robust covariance matrix allows the variance of \(u_i\) to differ across observations; it does not remove the need for a correctly specified first-stage conditional mean, suitable sampling conditions, and a full-rank covariance estimate.
For the usual heteroskedasticity-robust large-sample interpretation, observations must satisfy appropriate independence or sampling assumptions, the relevant moments must exist, and the first-stage regressors must not be perfectly collinear. The included covariates and excluded instruments must be handled consistently in the first stage. If observations are clustered or serially dependent, an ordinary heteroskedasticity-robust covariance estimate may be insufficient; the covariance estimator must reflect the dependence structure. A statistic computed with the wrong covariance estimator is not made robust merely by its name.
Under \(H_0\), the robust Wald statistic \(W\) has an asymptotic chi-square reference distribution with \(L\) degrees of freedom when the usual regularity conditions hold. Dividing by \(L\) gives the scalar-case KP rk Wald F normalization. That normalization is useful for instrument-strength diagnostics, but with robust covariance it should not automatically be interpreted as an exact finite-sample \(F\) random variable with the usual numerator and denominator degrees of freedom. Moreover, a test of the null that all first-stage coefficients are zero is not the same thing as a calibrated test of whether instruments are “strong enough” for reliable conventional IV inference.
Why the Scalar Normalization Is \(W/L\)
The division by \(L\) has a familiar connection to the classical first-stage F statistic. In the homoskedastic linear model, a joint test of \(L\) restrictions can be written as a Wald quadratic form divided by the number of restrictions. With the conventional homoskedastic covariance estimate and the usual linear-model assumptions, this gives the ordinary finite-sample F statistic for testing the excluded-instrument coefficients jointly.
The same scaling is used in the scalar-endogenous-regressor KP rk Wald F. To see its form, suppose the robust Wald statistic is calculated from the coefficient vector and its covariance estimate. Then the KP normalization is
This follows directly from the definition: \(W\) measures the joint discrepancy from the zero-coefficient null in covariance-scaled units, and \(L\) is the number of restrictions. The scaling makes the statistic comparable in form to a first-stage F measure. It does not restore the conventional finite-sample F distribution after replacing the homoskedastic covariance with a robust one. In particular, dividing an asymptotic chi-square statistic by \(L\) is not a proof that ordinary Stock-Yogo critical values apply.
Worked Example: Two Excluded Instruments
Consider a hypothetical first stage with one endogenous regressor and \(L=2\) excluded instruments. Suppose the estimated coefficients and their heteroskedasticity-robust covariance matrix are
The determinant of the covariance matrix is \(0.010(0.008)-0.002^2=0.000076\). Its inverse is therefore
The robust Wald statistic for \(H_0:\pi_1=\pi_2=0\) is
Since \(L=2\), the scalar-case KP rk Wald F is \(11.579/2\approx5.789\). This is the robust, normalized diagnostic for this hypothetical first stage. The calculation does not establish that the instruments are sufficiently strong for a particular IV inference procedure: no universally valid cutoff follows from the value alone. Under the joint zero-coefficient null, \(W\) has an asymptotic chi-square reference distribution with 2 degrees of freedom under the stated regularity conditions; the normalized value \(5.789\) should not be mistaken for a finite-sample F test with automatically applicable weak-instrument critical values.
Reading the Statistic Without Overclaiming
A larger KP rk Wald F generally indicates a stronger first-stage signal relative to its robust sampling uncertainty. A small value is a warning that the data may provide limited information about the endogenous regressor through the excluded instruments. It is a diagnostic, not a verdict: the consequences of a particular value depend on the number of instruments and endogenous regressors, the design, the covariance estimator, and the inferential procedure being used.
The statistic also does not test the exclusion restriction or the independence of the instruments from the structural error. Those are substantive conditions for IV validity, not conclusions that can be drawn from a first-stage strength statistic. Nor does a nonzero first stage alone establish that the instrument-outcome relationship is identified in a useful way. Weak identification means the instruments contain limited identifying information; underidentification concerns failure of the relevant rank condition. The KP rk Wald F is used to diagnose weak identification, not to certify instrument validity.
Worked Example: One Excluded Instrument
Suppose a hypothetical first stage has one excluded instrument, so \(L=1\). The estimated coefficient is \(0.18\), and its heteroskedasticity-robust standard error is \(0.06\). For a single restriction, the robust Wald statistic is the squared robust \(t\) statistic:
Dividing by \(L=1\) gives \(F_{\mathrm{KP}}=9/1=9\). The robust Wald test of the zero-coefficient null has an asymptotic \(\chi^2_1\) statistic of 9, equivalent to a two-sided standard-normal test with \(|z|=3\) in this large-sample reference. But a small p-value against the zero-first-stage null is not a guarantee of strong identification: with a sufficiently large sample, even a limited association can be statistically distinguishable from zero. Strength diagnostics concern the amount of identifying information and the reliability of IV inference, not merely rejection of \(H_0\).
Worked Example: Three Instruments and Robust Uncertainty
Suppose a hypothetical first stage has one endogenous regressor and three excluded instruments. Let the estimated coefficient vector be \((0.12,0.09,0.06)'\). For a transparent illustration, suppose the robust covariance matrix is diagonal with entries \(0.0004\), \(0.0009\), and \(0.0004\). These entries are variances, so the robust standard errors are \(0.020\), \(0.030\), and \(0.020\).
The inverse covariance matrix has diagonal entries \(1/0.0004=2500\), \(1/0.0009\approx1111.111\), and \(2500\). Thus the joint robust Wald statistic is
With \(L=3\), the scalar-case KP rk Wald F is \(54/3=18\). The arithmetic shows how each coefficient contributes according to its robust uncertainty. If the covariance matrix had nonzero off-diagonal entries, the inverse covariance matrix would also account for dependence among coefficient estimates; adding three squared individual \(t\) statistics would not generally produce the same joint statistic. The value 18 is a robust strength diagnostic, not an automatic pass mark, and this illustration does not show that the instruments satisfy the structural IV assumptions.
Critical Values, Alternatives, and Failure Modes
A common mistake is to read “robust” as “all standard weak-instrument cutoffs remain valid.” That reasoning is incorrect. Stock-Yogo critical values are tied to particular homoskedastic weak-identification statistics and specified instrument configurations. They are not automatically justified for the heteroskedasticity-robust KP rk Wald F. Likewise, the conventional rule of ten is not a universal guarantee. A KP value above ten does not prove that conventional IV confidence intervals have acceptable coverage, and a value below ten is not, by itself, a complete assessment of the analysis.
The main comparisons help clarify the statistic’s role:
| Diagnostic | What it uses | Important limitation |
|---|---|---|
| Conventional first-stage F or Cragg-Donald statistic | Homoskedastic covariance structure | May be inappropriate when heteroskedasticity matters |
| Heteroskedasticity-robust first-stage Wald test | Robust covariance for testing zero excluded-instrument coefficients | Tests a null; rejection alone does not establish adequate strength |
| KP rk Wald F | Robust rank-based strength diagnostic; in the scalar case it is the robust joint Wald statistic divided by \(L\) | Conventional weak-ID critical values are not automatically valid |
For multiple endogenous regressors, one must take particular care not to treat a collection of separate first-stage F statistics as if it fully described system-wide strength. Identification depends on whether the instruments provide distinct information across the endogenous regressors, and the rank-based statistic accounts for that system structure. Report the statistic together with the number of endogenous regressors, excluded instruments, covariance estimator, and model specification. When software produces a KP rk Wald F, use the definition and inference conventions documented for that implementation rather than reconstructing a multi-equation statistic from unrelated tests.
Diagnostics should also prompt checks of the analysis itself. Confirm that included covariates and excluded instruments were specified as intended, that the robust covariance estimator matches the dependence structure, and that the reported first-stage specification is the one used for IV estimation. If the data are clustered, a non-clustered heteroskedasticity-robust estimate may understate or overstate uncertainty. If instrument strength is doubtful, a robust KP statistic does not itself repair the resulting weaknesses of conventional 2SLS inference; it is a warning to choose and report inference with those limitations in view.
The KP rk Wald F is valuable because it replaces a homoskedastic covariance calculation with a heteroskedasticity-robust one in a widely used weak-identification diagnostic. Its limits are equally important: robust covariance does not provide universal critical values, establish instrument validity, or guarantee reliable conventional IV inference. A careful conclusion is therefore conditional: report the value, state the assumptions and covariance estimator, and explain why the chosen benchmark is appropriate rather than treating a familiar cutoff as decisive.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- With one endogenous regressor and \(L=4\) excluded instruments, a robust joint Wald statistic is \(W=20\). What is the scalar-case KP rk Wald F?
- Why does using a heteroskedasticity-robust covariance matrix not automatically justify Stock-Yogo critical values?
- A robust first-stage test rejects the null that all excluded-instrument coefficients are zero. Does this prove that the instruments are valid and strong enough for conventional IV inference? Explain.
- Why may summing individual first-stage \(t\) statistics be wrong when the coefficient estimates are correlated?
Solutions. (1) \(F_{\mathrm{KP}}=W/L=20/4=5\). This is the normalized robust diagnostic, not necessarily a statistic with an ordinary finite-sample F distribution. (2) Those critical values are calibrated for particular homoskedastic weak-identification statistics and configurations; changing the covariance estimator changes the statistic’s sampling behavior, so the comparison needs separate justification. (3) No. Rejection indicates evidence against the joint zero-coefficient null, but does not establish instrument exclusion or exogeneity, and does not alone guarantee adequate strength for conventional IV inference. (4) Correlation affects the joint covariance matrix. The joint Wald statistic uses the inverse of the full covariance matrix, including off-diagonal terms; adding individual squared \(t\) statistics generally ignores those terms and need not equal the joint test.