---
title: Asymmetric Portmanteau Test
url: https://www.emergentmind.com/topics/asymmetric-portmanteau-test
type: topic
---

# Asymmetric Portmanteau Test

Searching arXiv for recent papers on asymmetric portmanteau tests and related diagnostics.
An asymmetric portmanteau test is a residual-based omnibus diagnostic that departs from the classical Box–Pierce/Ljung–Box paradigm by targeting asymmetric, nonlinear, or directional dependence structures that are not exhausted by linear residual autocorrelation. In the literature represented here, the label does not denote a single universal statistic. Instead, it covers several non-equivalent constructions: standard portmanteau statistics applied to squared residuals to detect conditional heteroskedasticity; mixed tests that jointly stack residual and absolute-residual autocorrelations; generalized-covariance tests based on odd or sign-sensitive transforms; model-specific adequacy tests for asymmetric power GARCH classes; and a one-sided cross-correlation test that removes inverse-causality effects when testing weak exogeneity [2005.00931] [2007.09649] [2312.05373] [1811.08769] [2404.12685] [2606.07715].

## 1. Terminological scope and conceptual distinctions

In this literature, “asymmetric” has several precise meanings. In diagnostic checking for ARMA- and VARMA-type models, it can mean replacing residuals $\hat e_t$ by squared residuals $\hat e_t^2$ so that the test targets dependence in volatility magnitudes rather than linear mean dynamics. In asymmetric volatility models, it can refer to separate positive and negative shock effects entering the conditional scale recursion. In generalized-covariance frameworks, it means choosing odd or sign-sensitive transforms such as $x$, $x^3$, or $\operatorname{sign}(x)|x|^\alpha$ to interrogate skewness-related or leverage-type serial dependence. In the weak-exogeneity setting, asymmetry is directional: the test isolates dependence from past omitted variables to current shocks and excludes the reverse direction from its variance [2005.00931] [1811.08769] [2312.05373] [2606.07715].

| Variant | Core object | Null calibration |
|---|---|---|
| Squared-residual portmanteau in `portes` | Classical BP/LB/Hosking/Li–McLeod/Mahdi–McLeod statistics applied to $\hat e_t^2$ | Approximate $\chi^2$ or Monte Carlo [2005.00931] |
| ALDAR mixed portmanteau | Stacked ACFs of $\hat\eta_t$ and $|\hat\eta_t|$ | $\chi^2_{2M}$ [2007.09649] |
| GCov asymmetric subtest | Generalized autocovariances of odd/sign-sensitive transforms | $\chi^2$ with projector-rank df [2312.05373] |
| Weak-exogeneity asymmetric Portmanteau Test | Corrected one-sided cross-covariances $\widehat{\Gamma}_{uz}(j)$ | Asymptotically $\mathcal N(0,1)$ [2606.07715] |
| APGARCH and CCC-APGARCH adequacy tests | Squared standardized residual autocovariances or $\|z_t\|^2-d$ autocovariances | $\chi^2_m$ [1811.08769] [2404.12685] |

A recurrent misconception is that there exists a single canonical “Asymmetric Portmanteau Test.” The supplied sources instead show a family of diagnostics with different null hypotheses, different objects of aggregation, and different reference laws. This is explicit in the `portes` paper, which states that it does not define a separate “Asymmetric Portmanteau Test” statistic; the asymmetry-oriented diagnostic there is simply the standard portmanteau machinery run on squared residuals via `squared.residuals = TRUE` [2005.00931].

## 2. Squared-residual and mixed formulations

The classical univariate Box–Pierce and Ljung–Box statistics are
\[
Q_m = n\sum_{\ell=1}^{m}\hat r_\ell^2,
\qquad
\hat Q_m = n(n+2)\sum_{\ell=1}^{m}\frac{\hat r_\ell^2}{n-\ell},
\]
with $\hat r_\ell$ the sample autocorrelation at lag $\ell$. In `portes`, the nonlinear diagnostic is obtained by replacing $\hat e_t$ with $\hat e_t^2$ throughout, so that the same statistic tests dependence in magnitudes rather than linear autocorrelation. This operationalizes the McLeod–Li idea for detecting ARCH-type effects and volatility clustering. The package applies the same principle to multivariate Hosking, Li–McLeod, and generalized-variance statistics, and supports asymptotic $\chi^2$ calibration or Monte Carlo p-values with `innov.dist = c("Gaussian","t","stable","bootstrap")` [2005.00931].

The null hypothesis changes when squared residuals are used. For residual-based tests, adequacy means no remaining linear autocorrelation. For squared-residual tests, adequacy means no remaining nonlinear dependence in conditional variance. In the `portes` framework, the relevant degrees of freedom for Box–Pierce, Ljung–Box, Hosking, and Li–McLeod are
\[
\mathrm{df} = k^2(m-o),
\]
where $k$ is the series dimension and $o=p+q+ps+qs$ counts autoregressive and moving-average parameters, including seasonal components. For the generalized variance statistic $\mathfrak D_m$,
\[
\mathfrak D_m \overset{approx}{\sim} \chi^2_{\nu},
\qquad
\nu = k^2\left(\frac{1.5\,m(m+1)}{2m+1} - o\right).
\]
The paper emphasizes Monte Carlo p-values when $n$ is modest, $m$ is large relative to $n$, innovations are non-Gaussian or stable, or residuals exhibit conditional heteroskedasticity [2005.00931].

A distinct mixed construction appears in the asymmetric linear double autoregressive model, ALDAR($p$). After QMLE, the standardized residuals are
\[
\widehat{\eta}_t = \frac{y_t-\widehat{\alpha}'Y_{t-1}}{\widehat{\beta}'X_{t-1}}.
\]
For lags $k=1,\ldots,M$, the paper defines the residual ACF
\[
\widehat{\rho}_k
=
\frac{\sum_{t=p+k+1}^{n}(\widehat{\eta}_t-\bar\eta_1)(\widehat{\eta}_{t-k}-\bar\eta_1)}
{\sum_{t=p+1}^{n}(\widehat{\eta}_t-\bar\eta_1)^2}
\]
and the absolute-residual ACF
\[
\widehat{\gamma}_k
=
\frac{\sum_{t=p+k+1}^{n}(|\widehat{\eta}_t|-\bar\eta_2)(|\widehat{\eta}_{t-k}|-\bar\eta_2)}
{\sum_{t=p+1}^{n}(|\widehat{\eta}_t|-\bar\eta_2)^2}.
\]
The mixed portmanteau statistic stacks both vectors and uses a full covariance correction:
\[
Q(M)
=
n
\begin{pmatrix}
\widehat{\rho}\\
\widehat{\gamma}
\end{pmatrix}^{\!\prime}
(\widehat V\widehat G\widehat V')^{-1}
\begin{pmatrix}
\widehat{\rho}\\
\widehat{\gamma}
\end{pmatrix}
\Rightarrow \chi^2_{2M}.
\]
Here “mixed” means joint testing of serial dependence in the conditional mean and in the conditional volatility/asymmetric components. The paper reports that for S\&P 500 weekly returns, BIC-based selection chose $p=4$, asymmetry tests strongly rejected equality of positive and negative volatility coefficients, while mixed portmanteau p-values at $M=6,12,18$ were $0.41$, $0.18$, and $0.27$, indicating no lack of fit [2007.09649].

## 3. Directional asymmetry: weak exogeneity versus inverse causality

The most literal use of the title “asymmetric Portmanteau test” appears in the weak-exogeneity problem studied in "Causality versus Serial Correlation: an Asymmetric Portmanteau Test" [2606.07715]. The setup involves a shock process $u_t \in \mathbb R^{d_1}$ and an omitted-variable process $z_t \in \mathbb R^{d_2}$, with joint past
\[
\mathcal F_{t-1} = \sigma(\{u_s,z_s; s<t\}).
\]
The null is weak exogeneity:
\[
\mathcal H_0:\ \mathbb E[u_t\mid \mathcal F_{t-1}] = 0.
\]
Under stationarity and square integrability, this implies
\[
\Gamma_{uz}(j)=\mathbb E[u_tz_{t-j}']=0,\qquad j\ge 1.
\]

The paper’s central point is that conventional one-sided Hong-type tests based on
\[
\mathcal T_\omega=\sum_{j=1}^{T-1}\omega(j)\|\widehat{\Gamma}_{uz}(j)\|_F^2,
\qquad
\widehat{\Gamma}_{uz}(j)=T^{-1}\sum_{t=j+1}^{T}u_tz_{t-j}',
\]
are still “symmetric” in their variance. Their variance incorporates not only the direction relevant for weak exogeneity, from past $z$ to current $u$, but also “inverse causality,” from past $u$ to current $z$. The decomposition
\[
\mathcal T_\omega = \mathcal T_{1\omega} + \mathcal T_{2\omega}
\]
isolates a sum-of-squares term $\mathcal T_{1\omega}$ preserving the desired ordering and a sum-of-cross-products term $\mathcal T_{2\omega}$ whose variance is inflated by inverse causality [2606.07715].

The proposed correction subtracts the unordered part:
\[
\mathcal T_\omega^c = \mathcal T_\omega - \mathcal C_\omega
= \mathcal T_{1\omega} + \mathcal T_{2\omega}^c.
\]
This is the asymmetry of the test: the retained cross-products respect the martingale ordering that is consistent with $u_t$ being a martingale difference under $\mathcal H_0$. After centering and separate variance normalization, the implemented statistic
\[
\mathrm{Hete}
=
\frac{\mathcal T_{1\omega}^{(f)}-\mu_{\omega,T}^{(f)}}{\sqrt{D_{\omega,T}^{(f)}}}
+
\frac{\mathcal T_{2\omega}^{(f),c}}{\sqrt{D_{\omega,T}^{(f),(\mathrm{Hete})}}}
\]
is asymptotically standard normal under $\mathcal H_0$. The regularity conditions include conditional homoskedasticity and conditional homokurtosis of $u_t$ given $\mathcal F_{t-1}$, strict stationarity and finite $(8+\delta)$-order moments for $z_t$, strong mixing of $z_t$, and a kernel bandwidth satisfying $M/T\to 0$ or $M^2/T\to 0$ under weaker dependence [2606.07715].

The interpretive contrast is explicit. If a symmetric Hong test rejects but the asymmetric test does not, the paper interprets this as likely inverse causality rather than failure of weak exogeneity. If the asymmetric test rejects, past omitted variables help predict current shocks, so weak exogeneity fails. In the empirical application to Economic Policy Uncertainty shocks, both the benchmark and asymmetric tests reject, but at different horizons; augmenting local projections with lagged macro factors changes the estimated inflation response from negative to positive [2606.07715].

## 4. Adequacy testing in asymmetric volatility models

For asymmetric power GARCH models, the portmanteau problem is not formulated on raw residual autocorrelations but on squared standardized residual autocovariances. In the APARCH/APGARCH($p,q$) model with unknown power $\delta>0$,
\[
\sigma_t^{\delta}
=
\omega
+\sum_{i=1}^{q}\alpha_i\big(|\varepsilon_{t-i}|-\gamma_i\varepsilon_{t-i}\big)^\delta
+\sum_{j=1}^{p}\beta_j\sigma_{t-j}^\delta,
\]
or equivalently, in split-coefficient form,
\[
\varepsilon_t=\zeta_t\eta_t,\qquad
\zeta_t^\delta=\omega_0+\sum_{i=1}^{q}\alpha_{0i}^+(\varepsilon_{t-i}^+)^\delta+\alpha_{0i}^-(-\varepsilon_{t-i}^-)^\delta+\sum_{j=1}^{p}\beta_{0j}\zeta_{t-j}^\delta.
\]
The paper defines
\[
\hat\gamma(h)=\frac{1}{n}\sum_{t=h+1}^{n}(\hat\eta_t^2-1)(\hat\eta_{t-h}^2-1),
\qquad
\hat{\boldsymbol\gamma}_m=(\hat\gamma(1),\ldots,\hat\gamma(m))',
\]
and proves
\[
\sqrt n\,\hat{\boldsymbol\gamma}_m \xrightarrow{d}\mathcal N(0,D),
\qquad
D=(\kappa_4-1)^2I_m-(\kappa_4-1)C_mJ^{-1}C_m'.
\]
The portmanteau statistic is
\[
Q_m = n\,\hat{\boldsymbol\gamma}_m' \hat D^{-1}\hat{\boldsymbol\gamma}_m
\Rightarrow \chi^2_m.
\]
Here the estimation effect of all APARCH parameters, including the unknown power, is absorbed into $J$ and $C_m$, so no ad hoc degree-of-freedom reduction is required [1811.08769].

The multivariate analogue in the CCC-APGARCH class uses
\[
\varepsilon_t=D_tR^{1/2}u_t,\qquad
H_t=D_tRD_t,\qquad
z_t=H_t^{-1/2}\varepsilon_t,
\]
with asymmetric positive and negative matrices $A_i^+$ and $A_i^-$ entering the power recursion. The paper defines
\[
u_t=\|z_t\|^2-d,\qquad
\hat r_h = \frac{1}{n}\sum_{t=h+1}^{n}\hat u_t\hat u_{t-h},
\qquad
\hat r_m=(\hat r_1,\ldots,\hat r_m)'.
\]
Under the null,
\[
\sqrt n\,\hat r_m \xrightarrow{d}\mathcal N(0,D),
\qquad
Q_m = n\,\hat r_m'\hat D^{-1}\hat r_m \Rightarrow \chi^2_m,
\]
and similarly for the autocorrelation-based version
\[
Q_m^{(\rho)} = n\,\hat\rho_m' \hat D_\rho^{-1}\hat\rho_m \Rightarrow \chi^2_m.
\]
The covariance correction now involves $\hat J^{-1}\hat I\hat J^{-1}$, $\hat C_m$, and $\hat L_{rm}$, reflecting the QMLE effect in the multivariate asymmetric power recursion [2404.12685].

Both papers explicitly contrast these diagnostics with classical Ljung–Box testing. In GARCH-type settings, raw residual autocorrelations are not the primary adequacy criterion; the relevant misspecification appears in second-order dependence of standardized residual squares or their multivariate analogue. This is why the test is built from $\hat\eta_t^2-1$ in the univariate APGARCH case and from $\|z_t\|^2-d$ in the CCC-APGARCH case [1811.08769] [2404.12685].

## 5. Generalized-covariance, rank-sign, and high-dimensional extensions

The GCov-based framework generalizes portmanteau testing from linear autocovariances to nonlinear generalized covariances
\[
C_{k,\ell}(h)=E[g_k(X_t)g_\ell(X_{t-h})],\qquad h\ge 1.
\]
The full GCov portmanteau statistic is
\[
\hat\xi_T(H)=T\sum_{h=1}^{H}\operatorname{Tr}\big[\hat R^2(h)\big],
\qquad
\hat R^2(h)=\hat\Gamma(h)\hat\Gamma(0)^{-1}\hat\Gamma(h)'\hat\Gamma(0)^{-1}.
\]
For asymmetry, the paper restricts the transformation set to odd or sign-sensitive functions, producing the asymmetric-only statistic
\[
Q_{\mathrm{asym}}
=
T\sum_{h=1}^{H}
\operatorname{Tr}\!\left[
\hat\Gamma_A(h)\hat\Gamma_A(0)^{-1}\hat\Gamma_A(h)'\hat\Gamma_A(0)^{-1}
\right].
\]
Under GCov estimation and regularity conditions, the df is the rank of the projector $\Pi_A(H;\theta_0,f_0)$; in the common case of full identification using the same moments for estimation, $\mathrm{df}_{\mathrm{asym}}=K_A^2H-\dim(\theta)$. If estimation is performed by ML, QML, or AML instead of GCov, the paper recommends a bootstrap test [2312.05373].

A closely related high-dimensional extension is the shrinkage-regularized NLSD test. It starts from the augmented vector
\[
X_t^a=(a_1(X_t),\ldots,a_K(X_t))',\qquad a_1(X_t)=X_t,
\]
and replaces the potentially ill-conditioned $\hat\Gamma_T^a(0)$ by the Ledoit–Wolf shrinkage estimator
\[
\hat\Gamma_T^{a*}(0)=\hat\rho_{1,T}I+\hat\rho_{2,T}\hat\Gamma_T^a(0).
\]
The resulting statistic
\[
\hat\xi_{SR}^a(H)
=
T\sum_{h=1}^{H}
\operatorname{Tr}\!\left[
\hat\Gamma_T^a(h)\{\hat\Gamma_T^{a*}(0)\}^{-1}
\hat\Gamma_T^a(h)'
\{\hat\Gamma_T^{a*}(0)\}^{-1}
\right]
\]
has asymptotic $\chi^2(p^2H)$ null law, with $p=NK$, when $p/T\to 0$. The paper describes asymmetry as being operationalized by including sign- or state-sensitive transforms alongside linear and magnitude-sensitive ones [2603.10152].

A different extension, although not labeled an asymmetric portmanteau test in the title, is the center-outward rank- and sign-based VARMA portmanteau procedure. It constructs center-outward ranks and signs from residuals and forms score-based lagged cross-covariances. The resulting statistic is asymptotically $\chi^2_{d^2(m-p-q)}$ for sufficiently large $m$ and is distribution-free under broad innovation densities. Its relevance here is that the paper explicitly emphasizes robustness to skewness, heavy tails, and non-elliptical innovation densities; the method therefore functions as a portmanteau diagnostic robust to innovation asymmetry rather than a test of asymmetric dependence per se [2208.12143].

## 6. Calibration, implementation, and interpretation

Across these constructions, the reference distribution depends on the null being tested. Squared-residual diagnostics in `portes` are approximately $\chi^2$ with df adjusted for model order, but the paper repeatedly recommends Monte Carlo calibration for heavy tails, stable innovations, small $n$, or large $m$ [2005.00931]. The ALDAR mixed statistic is $\chi^2_{2M}$ after covariance correction [2007.09649]. GCov-based asymmetric subtests are $\chi^2$ with projector-rank df under GCov estimation and use bootstrap critical values when parameters are estimated by another method [2312.05373]. The weak-exogeneity asymmetric statistic is asymptotically standard normal, not chi-square [2606.07715]. APGARCH and CCC-APGARCH adequacy tests return to $\chi^2_m$, but only after nontrivial parameter-estimation corrections in the covariance matrix [1811.08769] [2404.12685].

Lag or bandwidth selection is likewise model-dependent. `portes` recommends using a range of lags such as $5,10,15,\ldots$ and notes that larger $m$ interrogates longer-range dependence but increases df and small-sample distortions [2005.00931]. The ALDAR paper uses modest values such as $M=6,12,18$ and states that no explicit data-driven rule is imposed [2007.09649]. The GCov paper commonly uses small $H$, for example $H=3$ [2312.05373]. The weak-exogeneity paper recommends kernel weighting $\omega(j)=k^2(j/M)$ with practical choices $M\in[0.75\cdot T^{1/3},\sqrt T-1]$ [2606.07715].

Interpretation also differs sharply across variants. In squared-residual diagnostics, rejection indicates remaining variance dynamics such as ARCH/GARCH or volatility clustering, even when linear residual autocorrelations are negligible [2005.00931]. In ALDAR, inspecting $\widehat{\rho}_k$ and $\widehat{\gamma}_k$ separately locates whether misfit is in the conditional mean or in the volatility/asymmetric component [2007.09649]. In GCov and SR-NLSD settings, significant odd-transform blocks indicate sign-related nonlinear serial dependence, whereas even-transform blocks indicate volatility-type nonlinear dependence [2312.05373] [2603.10152]. In the weak-exogeneity setting, rejection of the asymmetric test indicates that omitted variables’ past predicts current shocks; failure to reject when Hong rejects suggests inverse causality rather than a violation of weak exogeneity [2606.07715].

A plausible implication of this literature is that “asymmetry” in portmanteau testing is best understood as a design principle rather than a unique formula. The design principle may target magnitudes instead of levels, positive and negative shocks separately, odd versus even transforms, or directional predictability rather than bidirectional dependence. The common feature is that each statistic augments or restructures the portmanteau aggregation so that the null distribution remains tractable while the test becomes sensitive to dependence patterns that the classical symmetric residual autocorrelation check was not constructed to detect.

Source: https://www.emergentmind.com/topics/asymmetric-portmanteau-test