Hodrick t-statistics in Persistent Time Series

• Hodrick t-statistics are variants of the classical t-test tailored for time series models with persistent (nearly unit root) autoregression.
• They rely on exact and numerical distributional characterizations and employ whitening to enable standard t-distribution inference.
• CVF tests use simulation and linear programming to achieve robust size control and higher power compared to conventional methods.

Hodrick t-statistics refer to variants and generalizations of the classical Student t-statistic, specifically in the context of time series models with persistent or autocorrelated regressors, such as AR(1) processes with coefficients near unity. The central issue addressed by Hodrick t-statistics is the breakdown of the standard t-distribution’s finite-sample and asymptotic critical values in the presence of strong regressor persistence, resulting in nonpivotal distributions that depend on nuisance parameters. Recent research provides both exact distributional characterizations for the AR(1) case and algorithms for testing with approximately similar size, even when conventional simulation methods fail.

1. Failure of the Classical t-statistic under Persistent Autoregression

The classical Student t-statistic,

$$T_n = \frac{\sqrt{n}\,(\bar X_n - \mu_0)}{s_n}, \quad s_n^2 = \frac{1}{n-1} \sum_{i=1}^n (X_i - \bar X_n)^2,$$

assumes that the sequence $X_1,\dots,X_n$ is i.i.d. Gaussian. Under these conditions, $T_n$ follows the $t_{n-1}$ distribution, because the numerator and denominator are independent (by Cochran’s theorem).

Under an AR(1) process,

$$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$

the sequence $\{X_t\}$ exhibits autocorrelation, never i.i.d. The covariance structure

$\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$

invalidates Cochran’s factorization, so the distribution of $T_n$ is a ratio of dependent variables and is not $t_{n-1}$ for any $n$. This failure is especially pronounced for $X_1,\dots,X_n$0 near unity, causing the asymptotic distribution of $X_1,\dots,X_n$1 to become nonpivotal, that is, dependent on $X_1,\dots,X_n$2 (Benhamou, 2018).

2. Exact and Asymptotic Distributions of the AR(1) t-statistic

For finite samples, the AR(1) t-statistic $X_1,\dots,X_n$3 admits a formal expression as the ratio of a zero-mean normal and the Euclidean norm of a correlated normal vector: $X_1,\dots,X_n$4 where

$X_1,\dots,X_n$5

with $X_1,\dots,X_n$6 and $X_1,\dots,X_n$7, $X_1,\dots,X_n$8. $X_1,\dots,X_n$9 and $T_n$0 are no longer independent.

The explicit joint law of $T_n$1 is multivariate normal for $T_n$2, and $T_n$3 involves the norm of a projected Gaussian vector. Consequently, $T_n$4’s marginal law is given by integrating over the joint density of $T_n$5, as derived in (Benhamou, 2018). No closed-form in elementary functions exists, but distributions can be simulated or evaluated numerically.

As $T_n$6 with $T_n$7 fixed, $T_n$8 converges in law to standard normal, reflecting the central limit theorem. For $T_n$9, classical results are recovered. However, for $t_{n-1}$0 near 1, dependence on $t_{n-1}$1 remains significant in finite and moderate samples.

3. Modified t-statistics and Whitened Inference under AR(1)

To bypass nonpivotality for inference under AR(1), a modified statistic $t_{n-1}$2 can be constructed by linear transformation: $t_{n-1}$3 yielding

$t_{n-1}$4

where $t_{n-1}$5 are i.i.d. normal, making $t_{n-1}$6 exactly noncentral t-distributed, allowing for exact confidence intervals and tests (Benhamou, 2018). This whitening can also be effected via diagonalization of the circulant Toeplitz covariance using the discrete Fourier transform basis.

A plausible implication is that practitioners with sufficient knowledge of $t_{n-1}$7 (i.e., $t_{n-1}$8) and access to the data can always transform the sample so that standard t-distribution inference holds. This requires full knowledge of the covariance, which is often not feasible with highly persistent regressor processes.

4. Critical Value Function (CVF) Testing in Predictive Regressions

For inference about predictive slope coefficients in models with persistent regressors,

$t_{n-1}$9

the limiting null distribution of the t-statistic

$$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$0

depends on the autoregressive parameter $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$1 and is nonpivotal for $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$2 (Moreira et al., 2016).

The CVF test rejects $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$3 when $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$4, where

$$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$5

with $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$6 chosen to ensure exact size $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$7 at each $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$8 in a grid $$X_t = \mu + \epsilon_t, \quad \epsilon_t = \rho\,\epsilon_{t-1} + \sigma v_t, \quad v_t \sim \mathrm{i.i.d.}\ N(0,1),\quad |\rho|<1,$$9. The grid can be refined to achieve uniform size control (to within $\{X_t\}$0) over the nuisance parameter space.

This approach implements similarity directly using linear programming over simulated draws, so it does not rely on consistent estimation of $\{X_t\}$1, splitting statistics, or randomization.

CVF Algorithm (summary)

Step	Description	Notes
1	Choose $\{X_t\}$2 grid $\{X_t\}$3 covering region of interest	e.g., $\{X_t\}$4
2	Simulate $\{X_t\}$5 draws $\{X_t\}$6 from baseline mixture $\{X_t\}$7	Large $\{X_t\}$8 for precision
3	Compute score matrix $\{X_t\}$9
4	Solve linear program to find $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$0	Imposes size $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$1
5	Evaluate and refine as necessary	Add grid points as needed

This method controls empirical size at the nominal level, and for $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$2, $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$3, achieves size within 1 percentage point for $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$4 even with autocorrelated errors and trend, outperforming bootstrap and L2/Wright similar tests (Moreira et al., 2016).

5. Comparison to Conventional Methods and Power Properties

Traditional methods, including constant critical value t-tests, bootstrap, and subsampling, do not achieve uniform size over $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$5 in highly persistent settings. Bootstraps may under- or overestimate the critical value, and subsampling can be overly conservative or anti-conservative. The Wright L2 and Jansson & Moreira UMPCU tests use various conditioning and similarity-enforcing constructions, but can have non-monotonic power or intricate implementation (multiple integrals).

The CVF approach, by construction, yields approximately similar tests over all persistence levels considered, bridging the gap between stationary and nearly integrated regressors without user tuning of nuisance parameter estimators or block sizes. In power comparisons, CVF tests outperform L2 and UMPCU, especially in the nearly unit root regime; for example, with $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$6 and $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$7, the CVF achieves $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$8 power, compared to $\Cov(X_t, X_u) = \frac{\sigma^2}{1-\rho^2}\rho^{|t-u|}$9 for alternatives (Moreira et al., 2016).

6. Applications: Sharpe Ratio, Confidence Bounds, and Financial Time Series

In financial econometrics, the inference on the mean of an autocorrelated process directly determines the properties of Sharpe ratio estimators: $T_n$0 Under AR(1) assumptions, its sampling distribution follows $T_n$1 or the whitened $T_n$2 law, requiring exact or numerically integrated confidence intervals, rather than the classical t-distribution. Monte Carlo or transformation-based methods are necessary to reflect autocorrelation in the uncertainty quantification (Benhamou, 2018).

A plausible implication is that practitioners should avoid using off-the-shelf t-tables for Sharpe ratio inference when data exhibit serial correlation, even for large $T_n$3 or moderate autocorrelation.

7. Broader Implications and Limitations

The study of t-statistics under autocorrelation and near-unit-root persistence highlights the need for nonstandard inference procedures in time series. While whitening and CVF-based methods restore valid inferential properties, they entail computational cost and may require knowledge (or accurate estimation) of nuisance parameters (e.g., $T_n$4, $T_n$5) or fine grid specification.

Uniformly similar tests constructed by CVF ensure robust inference across model regimes, and their linear programming-based implementation is numerically feasible for moderate $T_n$6 and grids. The explicit breakdown of classical t-distribution results in persistent environments demonstrates the importance of context-aware critical value derivation in high-persistence time series analysis.