---
title: Edgeworth–Cornish–Fisher Expansions
url: https://www.emergentmind.com/topics/edgeworth-cornish-fisher-expansions
type: topic
---

# Edgeworth–Cornish–Fisher Expansions

Edgeworth–Cornish–Fisher Expansions are higher-order asymptotic series that provide refined approximations to the distribution and quantiles of sample statistics, surpassing the accuracy of the classical central limit theorem (CLT). These expansions systematically incorporate standardized cumulants beyond variance (notably, skewness and kurtosis), leading to corrections which enhance the normal approximation for finite-sample inferential procedures. The Edgeworth expansion addresses the distribution function, while the Cornish–Fisher expansion is its quantile inversion, producing expansions for percentiles or confidence bounds. The formalism applies broadly: i.i.d., non-i.i.d., weighted, dependent, or even randomly-sized samples, provided certain regularity conditions on cumulant growth and smoothness are met [1210.4052, 1002.4338, 1604.00539, 1810.00289].

## 1. Edgeworth Expansion: Fundamentals and Formulation

The Edgeworth expansion approximates the CDF (and, if desired, density) of a standardized statistic whose cumulants admit power series in the inverse sample size. Let $\hat\theta$ be a standard estimate so that
\[
\kappa_r(\hat\theta) = a_{r,0} n^{1-r} + a_{r,1} n^{-r} + \cdots,
\quad r \geq 1, \quad a_{1,0}=0, \quad a_{2,0}>0
\]
and define the normalized variable
\[
Y_n = (n/a_{2,0})^{1/2}(\hat\theta - \theta_0)
\]
whose scaled cumulants $\ell_r$ satisfy $\kappa_r(Y_n) = r!\,\ell_r + O(n^{-1/2})$, $\ell_r = O(n^{1-r/2})$ [1210.4052]. The $R$-term one-sided Edgeworth expansion is
\[
P_n(x) = \mathbb{P}\{Y_n \leq x\}
        = P(x) - p(x) \sum_{r=1}^R n^{-r/2} h_r(x) + O(n^{-(R+1)/2})
\]
where $P(x), p(x)$ are the CDF and PDF of a selected reference law (usually $N(0,1)$), and each $h_r(x)$ is a polynomial in $x$ and the scaled cumulants, with a representation in terms of generalized Hermite polynomials:
\[
h_r(x) = \sum_{k=r,\,r+2,\dots,3r} C_{rk}\, H_{k-1}(x)
\]
with $C_{rk}$ combinatorially defined by Bell polynomials in $\ell_1,\dots,\ell_r$ [1210.4052].

The expansion increases the accuracy of finite-sample approximations, with the order of the error term decreasing as $R$ increases.

## 2. Cornish–Fisher Expansion for Quantiles

The Cornish–Fisher expansion inverts the Edgeworth series to obtain asymptotic expansions of quantiles (percentiles) for the statistic. Given $z_p = P^{-1}(p)$, the $p$-quantile of the reference law, the expansion for the $p$-quantile $Q_n(p)$ of $\hat\theta$ is
\[
Q_n(p) = (a_{2,0}/n)^{1/2}\left[z_p + \sum_{r=1}^R n^{-r/2} g_r(z_p)\right] + \theta_0 + O(n^{-(R+1)/2})
\]
where $g_r(z_p)$, like $h_r$, are polynomials in $z_p$ and the scaled cumulants, and can be calculated using formulas involving Hermite polynomials [1210.4052, 1810.00289]. For example, to $O(n^{-1})$ with normal reference:
\[
g_1(x) = \ell_1 + \frac{1}{3} \ell_3(x^2-1)
\]
\[
g_2(x) = \frac{1}{4}\ell_4(x^3-3x) + \frac{1}{6} \ell_3^2(2x^3-5x) + \frac{1}{2} \ell_2 \ell_1(x^2-1) + \frac{1}{2}\ell_1^2 x
\]
These expansions serve to correct coverage of confidence intervals and significance thresholds under non-negligible skewness or kurtosis [1210.4052, 1810.00289].

## 3. Generalizations: Weighted, Dependent, and Non-i.i.d. Cases

Edgeworth–Cornish–Fisher expansions have been generalized for weighted empirical distributions and for samples with non-identically distributed components. For independent $(X_{1n},\ldots,X_{nn})$ with predetermined weights $w_{1n},\ldots,w_{nn}$, cumulants of smooth functionals $T(\hat F_n)$ are expanded using higher-order von Mises derivatives [1002.4338]. This yields
\[
\kappa_r\bigl(T(\hat F_n)\bigr)
= \sum_{j=r-1}^\infty a_{rj} n^{-j}
\]
with explicit formulas for cumulant coefficients up to $O(n^{-3/2})$, enabling third-order expansion for both the distribution and quantiles. The expansions are valid provided weights are uniformly bounded and the power sums $n^{-1}\sum_i w_{i,n}^r$ remain $O(1)$. Applications include nonparametric estimates, rank statistics, and regression functionals [1002.4338].

For dependent sequences, valid Edgeworth and Cornish–Fisher expansions require additional terms reflecting the bias and stochastic variability of variance estimators. For studentized statistics under strong mixing, the expansion involves three series: powers of $n^{-1/2}$, $(n/\ell)^{-1/2}$ (where $\ell$ is the number of lag-covariances estimated), and the bias $a_n$ of the studentizing factor. Each enters into the distributional and quantile expansions as specific correction polynomials [1001.1858].

## 4. Distributional Reference Laws and Matched-Skewness Expansions

While the standard normal is a frequent reference law, expansions about alternative distributions yield improved efficiency for skewed statistics. Expanding about a gamma law with matched skewness (i.e., choosing the gamma parameter $m$ so $\kappa_3$ of the standardized statistic matches $\kappa_3$ of the gamma) nullifies leading-order skewness corrections, drastically reducing the polynomial complexity required for a fixed accuracy [1210.4052]. This "matched-skewness gamma" approach can cut the number of necessary nonzero correction terms by up to 90% for $O(n^{-2})$ accuracy relative to the standard-normal expansion.

## 5. Explicit Formulas, Implementation, and Non-Asymptotic Error

Edgeworth correction polynomials are built from cumulant-derived coefficients and generalized Hermite polynomials. For the standard normal, $H_0=1$, $H_1=x$, $H_2=x^2-1$, $H_3=x^3-3x$, $H_4=x^4-6x^2+3$. For a regular smooth function of the mean, cumulants up to order 4 can be computed by symbolic differentiation and substitution of population moments, with implementation facilitated via automated tools (e.g., Maple sheets and R code) for complex parametric or nonparametric estimators [1810.00289].

The Cornish–Fisher coefficients for quantiles are derived via formal Taylor expansion or via explicit back-substitution, with non-asymptotic error bounds available when the Edgeworth remainder is controlled. The minimum required for second-order accuracy is typically explicit control of the Edgeworth remainder and the possibility of a monotone Bartlett correction transformation [1604.00539]. Bounds hold uniformly only on central probability intervals.

## 6. Applications and Impact

Edgeworth–Cornish–Fisher expansions have broad applicability in contemporary statistics:

- Confidence intervals and hypothesis tests for sample means, variances, medians (including under random sample size), U-statistics, M-estimators, and network moments [2004.06615, 1905.07765].
- Correcting critical values for heavy-tailed, skew, or dependent samples, with applications in finance (VaR), biostatistics, insurance, and network inference [1604.00539, 2004.06615].
- Improved approximations for non-normal limit laws (Student's $t$, Laplace, Gamma) in scenarios with random or overdispersed sample sizes [1905.07765, 1210.4052].
- Accurate quantile estimation for high-dimensional or weighted empirical functionals [1002.4338].
- Enhanced coverage accuracy for bootstrap confidence intervals and bias-corrected and accelerated (BCA) resampling [1810.00289].

The expansions provide significant reductions in coverage error compared to first-order (CLT-based) inference, with quantitative benefits persisting in moderate $n$ regimes when higher cumulants are non-negligible. The practical implementation is enabled by explicit coefficient expressions and algorithmic differentiation, with code templates available for symbolic and numerical computation [1810.00289].

## 7. Limitations, Regularity, and Error Rates

The validity of these expansions is contingent on regularity conditions: existence of a cumulant expansion (i.i.d. or weakly dependent data), nonlattice structure or sufficient "self-smoothing" (particularly for U-statistics and network moments), boundedness of higher moments, and adequate smoothness of the statistic's functional representation [1210.4052, 2004.06615]. For non-i.i.d. or weighted cases, von Mises derivatives must remain uniformly bounded. In certain settings (e.g., extreme tails or excessively small sample size), the expansions do not converge uniformly, and explicit error bounds require knowledge of specific Edgeworth remainder terms [1604.00539].

The accuracy is characterized by $O(n^{-R/2})$ error for $R$-term expansions; for non-i.i.d. or weighted cases, achieving $O(n^{-3/2})$ error is standard. For dependent or studentized statistics, the error in quantile estimation accumulates from normal, estimation-variance, and bias contributions, whose coefficients must be balanced according to bias-variance tradeoffs (e.g., choice of lag-length in long-run variance estimation) [1001.1858].

In applied contexts, these expansions permit computation of accurate, higher order-corrected critical values and confidence intervals for a wide variety of statistics, over a flexible range of sample regimes and underlying data structures [1210.4052, 1002.4338, 1810.00289].

Source: https://www.emergentmind.com/topics/edgeworth-cornish-fisher-expansions