Table of Contents
Fetching ...

Estimation of time series by Maximum Mean Discrepancy

Pierre Alquier, Jean-David Fermanian, Benjamin Poignard

TL;DR

This work develops two maximum mean discrepancy (MMD)–based minimum-distance estimators for dependent time series, enabling estimation when likelihoods are intractable by simulating from the model. It provides non-asymptotic bounds, consistency, and asymptotic normality under beta-mixing, and derives the asymptotic behavior for both the ideal and simulation-based estimators, including mixed-rate regimes. The authors implement gradient-based algorithms (ISMMD-sgd, PSMMD-sgd) with a Gaussian kernel and validate the approach through extensive simulations on SV, GARCH, ARMA, nonlinear MA, and Ricker models, highlighting robustness to misspecification and competitive or superior performance to standard likelihood-based methods. A data-driven lag selection method based on an MMD criterion is also proposed to adaptively choose the transform $y_t$. Overall, the paper provides a practical, theoretically grounded framework for simulation-based MMD inference in dependent settings with latent structure.

Abstract

We define two minimum distance estimators for dependent data by minimizing some approximated Maximum Mean Discrepancy distances between the true empirical distribution of observations and their assumed (parametric) model distribution. When the latter one is intractable, it is approximated by simulation, allowing to accommodate most dynamic processes with latent variables. We derive the non-asymptotic and the large sample properties of our estimators in the context of absolutely regular/beta-mixing random elements. Our simulation experiments illustrate the robustness of our procedures to model misspecification, particularly in comparison with alternative standard estimation methods.

Estimation of time series by Maximum Mean Discrepancy

TL;DR

This work develops two maximum mean discrepancy (MMD)–based minimum-distance estimators for dependent time series, enabling estimation when likelihoods are intractable by simulating from the model. It provides non-asymptotic bounds, consistency, and asymptotic normality under beta-mixing, and derives the asymptotic behavior for both the ideal and simulation-based estimators, including mixed-rate regimes. The authors implement gradient-based algorithms (ISMMD-sgd, PSMMD-sgd) with a Gaussian kernel and validate the approach through extensive simulations on SV, GARCH, ARMA, nonlinear MA, and Ricker models, highlighting robustness to misspecification and competitive or superior performance to standard likelihood-based methods. A data-driven lag selection method based on an MMD criterion is also proposed to adaptively choose the transform . Overall, the paper provides a practical, theoretically grounded framework for simulation-based MMD inference in dependent settings with latent structure.

Abstract

We define two minimum distance estimators for dependent data by minimizing some approximated Maximum Mean Discrepancy distances between the true empirical distribution of observations and their assumed (parametric) model distribution. When the latter one is intractable, it is approximated by simulation, allowing to accommodate most dynamic processes with latent variables. We derive the non-asymptotic and the large sample properties of our estimators in the context of absolutely regular/beta-mixing random elements. Our simulation experiments illustrate the robustness of our procedures to model misspecification, particularly in comparison with alternative standard estimation methods.
Paper Structure (26 sections, 8 theorems, 111 equations, 48 figures)

This paper contains 26 sections, 8 theorems, 111 equations, 48 figures.

Key Result

Proposition 3.1

Assume that the kernel $k$ can be written as $k(y,y') = F(\|y-y'\|)$, where $F(x) = \int_{x}^\infty f(u) \,du$ for some nonnegative continuous map $f$, with $\int_0^{+\infty} f(t)\, dt=1$. Then $\varrho_t \leq \beta(t)$.

Figures (48)

  • Figure 1: Case 1: $\eta_t \sim \mathcal{N}(0,1)$. $T=300$. Estimation under correct specification.
  • Figure 2: Case 2: $\eta_t = t(3)/\sqrt{3}$. $T=300$. Estimation under misspecification.
  • Figure 3: Case 1: $\eta_t \sim \mathcal{N}(0,1)$. $T=1000$. Estimation under correct specification.
  • Figure 4: Case 2: $\eta_t = t(3)/\sqrt{3}$. $T=1000$. Estimation under misspecification.
  • Figure 6: Case 1: $u_t \sim \mathcal{N}(0,1)$. $T=300$. Estimation under correct specification.
  • ...and 43 more figures

Theorems & Definitions (25)

  • Remark 1
  • Remark 2
  • Proposition 3.1: Proposition 4.4 in cherief2022finite
  • Proposition 3.2
  • Example 1: Markov chains
  • Example 2: ARMA process
  • Example 3: GARCH and stochastic volatility models
  • Example 4: An example with non-exponential decay
  • Example 5: HMM
  • Example 6
  • ...and 15 more