## Iterations of dependent random maps and exogeneity in nonlinear dynamics

Zinsou Max Debaly

∗ ∗ †

Lionel Truquet

## Abstract

We discuss the existence and uniqueness of stationary and ergodic nonlinear autoregressive processes when exogenous regressors are incorporated into the dynamic. To this end, we consider the convergence of the backward iterations of dependent random maps. In particular, we give a new result when the classical condition of contraction on average is replaced with a contraction in conditional expectation. Under some conditions, we also discuss the dependence properties of these processes using the functional dependence measure of Wu (2005) that delivers a central limit theorem giving a wide range of applications. Our results are illustrated with CHARN models, GARCH processes, count time series, binary choice models and categorical time series for which we provide many extensions of existing results.

2010 Mathematics Subject Classification: Primary 62M10; secondary 60G05, 60G10.

Keywords and Phrases: time series, random maps, ergodicity, dependence.

## 1 Introduction

Among the various contributions devoted to time series analysis, theoretical results justifying stationarity and ergodicity properties of some standard stochastic processes when exogenous covariates

∗ CREST-ENSAI, UMR CNRS 9194, Campus de Ker-Lann, rue Blaise Pascal, BP 37203, 35172 Bruz cedex, France. † This work was funded by CY Initiative of Excellence (grant 'Investissements d'Avenir' ANR-16-IDEX-0008), Project 'EcoDep' PSI-AAP2020-0000000013.

are incorporated in the dynamic are rather scarce. A notable exception concerns linear models, such as VARMA processes, for which such properties are a consequence of the linearity. See for instance L¨ utkepohl (2005), a standard reference for multivariate time series models. Moreover, linear models represent a very simple setup for discussing various exogeneity notions found in the literature. See for instance Engle et al. (1983). For nonlinear dynamics, a few contributions consider the problem of exogenous regressors. For general GARCH type processes, Francq and Thieu (2019) recently studied stationarity conditions when the noise and the covariate process form a stationary process. Agosto et al. (2016) considered a Poisson autoregressive process with exogenous regressors (PARX models), under a Markov chain assumption for the covariate process. de Jong and Woutersen (2011) consider the case of dynamic binary choice models and provide results about stationarity and mixing properties of a 0 / 1 -valued time series which is autoregressive and defined conditionally on some exogenous regressors. Fokianos and Truquet (2019) studied stationarity and ergodicity of general categorical time series defined conditionally on a strictly exogenous covariate process.

In this paper, we give general results for getting stationarity, ergodicity and stochastic dependence properties for general nonlinear dynamics defined in terms of iterations of random maps. For simplicity, we explain our setup with the following example which represents the basis for studying other processes. Let us consider the following model

$$X _ { t } = F \left ( X _ { t - 1 }, Z _ { t - 1 }, \varepsilon _ { t } \right ), \quad t \in \mathbb { Z },$$

where ( Z t ) t ∈ Z is a covariate process and ( ε t ) t ∈ Z a noise process. One can note that X t = f t ( X t -1 ) for the random function defined by f t ( x ) = F ( x, Z t -1 , ε t ). The sequence ( f t ) t ∈ Z is a sequence of dependent random maps even if the ε ′ t s are i.i.d. because typically the Z ′ t s exhibit temporal dependence. A key point for getting existence of a stationary solution in (1) is to control the behavior of the backward iterations { f t ◦ f t -1 ◦ · · · ◦ f t -n ( x ) : n ≥ } 1 . The convergence of such iterations of random maps has been extensively studied in the independent case. In this case, the process ( X t ) t ∈ Z is a Markov chain. We defer the reader to Letac (1986) and Diaconis and Freedman (1999) for seminal papers on iterated independent random maps and to Wu and Shao (2004) for additional results useful in a time series context. The last contribution is particularly interesting

for getting existence of some moments for the marginal X t and also some dependence properties for the process ( X t ) t ∈ Z that are often needed for statistical applications. All these contributions use average contraction conditions and the interested reader is referred to the interesting survey of Stenflo (2012) for an overview of the available results. There also exist some contributions studying the more general case of iterated stationary random maps ( f t ) t ∈ Z . For instance, Borovkov (1998) gives many results for studying what he calls stochastically recursive sequences, when the independence assumption is removed. See also Iosifescu (2009) for a survey of some available results. The results obtained in the dependent case are based on Lyapunov type exponents and the convergence of the backward iterations is only studied almost surely. We recall the following result which can be found in Elton (1990) (see also Iosifescu (2009), Theorem 6 2) and which generalizes . a widely known result given in Brandt (1986) or Bougerol and Picard (1992) for iterations of affine random maps.

To do so, we introduce some notations and conditions. We assume that f t : E → E are random Lipschitz functions where E denotes a locally compact Polish space endowed with a metric d . We define the Lipschitz constant of a measurable function g : E → E by

$$c ( g ) \coloneqq \sup _ { x \neq y \in E } \frac { d \left ( g ( x ), g ( y ) \right ) } { d ( x, y ) }.$$

/negationslash

Moreover, for any integers s &lt; t , we set f t s = f t ◦ · · · ◦ f s . For a positive real number x , we set log + ( x ) = log( x ) if x ≥ 1 and 0 otherwise.

Theorem 1. Assume that the process (( Z , ε t t )) t ∈ Z in (1) is stationary and ergodic. Assume further that E [ log + c f ( 0 ) ] &lt; ∞ and E [ log + d x , f ( 0 0 ( x 0 )) ] &lt; ∞ for some point x 0 ∈ E .

- 1. There exists a constant χ ∈ R ∪ {-∞} called Lyapunov exponent and such that

$$\lim _ { n \to \infty } \frac { 1 } { n } \log c ( f _ { 1 } ^ { n } ) = \chi \ a. s.$$

Moreover

$$\chi = \lim _ { n \to \infty } \frac { 1 } { n } \mathbb { E } \left [ \log c \left ( f _ { 1 } ^ { n } \right ) \right ] = \inf _ { n \geq 1 } \frac { 1 } { n } \mathbb { E } \left [ \log c \left ( f _ { 1 } ^ { n } \right ) \right ].$$

- 2. If the constant χ is negative, then the almost sure limit f t -∞ = lim k →∞ f t t -k ( x ) exists for any x ∈ E and does not depend on x . Setting X t = f t -∞ , the process ( X t ) t ∈ Z is stationary and ergodic and satisfies the recursions (1). Moreover, ( X t ) t ∈ Z is the unique stationary process satisfying (1).

The affine random maps version of this result has been applied recently by Francq and Thieu (2019) for studying stationarity of asymmetric power GARCH processes. For nonlinear random maps, Theorem 1 is less known in the time series literature. In this paper, we will make use of Theorem 1 for defining a general class of categorical time series with exogenous covariates. In particular, we will see in Section 4.4 how Theorem 1 can be applied to binary time series and lead to an improvement of a result of de Jong and Woutersen (2011).

However, the result presented above has several limitations.

- 1. First, it requires the random maps f t to be almost surely Lipschitz. Such a property is not always valid, for instance for the Poissonian autoregressions discussed in Section 4.3. When there are no exogenous covariates, Davis and Liu (2016) studied integer-valued time series by using a different contraction result, developed by Wu and Shao (2004).
- 2. Existence of some moments for the marginal distributions that are sometimes necessary for statistical applications cannot be obtained directly from this result.
- 3. For autoregressions with several lags, it is not straightforward to get an explicit condition on the parameters of the model to ensure that χ &lt; 0.

To overcome these drawbacks, we will adapt the approach used by Wu and Shao (2004) for independent random maps to the case of dependent random maps. Our main result, see Theorem 2 and its extension Theorem 4, is obtained by replacing the usual contraction on average condition by a contraction in conditional expectation. The assumptions that we use are very simple to check and the proof of our main result is straightforward but its merit is to provide an elegant way for presenting a general approach which encompasses most of the previous attempts to include exogenous regressors in nonlinear dynamics. For strictly exogenous regressors, i.e. the processes ( Z t ) t ∈ Z

and ( ε t ) t ∈ Z are independent, we also provide an additional result, see Theorem 3, with weaker assumptions. In the context of Theorem 2 and Theorem 4, we will then discuss how to control the functional dependence measure of ( X t ) t ∈ Z introduced by Wu (2005), a dependence notion which is an alternative to the standard strong mixing condition and which can be more easily checked for iterations of contracting random maps. Let us mention that even in the independent case, mixing properties of the process ( X t ) t ∈ Z require restrictive assumptions on the noise distribution, otherwise such properties may fail. We refer the reader to the standard textbook of Doukhan (1994), section 2 4 for mixing properties of iterations of independent random functions and to Andrews (1984) for . a famous counterexample of a non strongly mixing sequence defined via iterations of random maps. In the dependent case, as in (1), getting usual strong mixing properties seems to be harder because the process ( X t ) t ∈ Z does not have a Markov structure in general and the criteria for getting mixing properties of Markov chains are useless.

This paper is mainly motivated by dynamics of type (1) with covariates that are not necessarily strictly exogenous, assuming that at any time t , the noise ε t is independent from the past information σ (( Z , ε s s ) : s ≤ t -1). The term predetermindness is sometimes used in the literature. This independence assumption is substantially weaker than the independence between the two processes ( ε t ) t ∈ Z and ( Z t ) t ∈ Z . The latter independence condition implies strict exogeneity, a notion initially defined by Sims (1972) and extended to general models by Chamberlain (1982). Strict exogeneity is useful for deriving the conditional likelihood of the X ′ t s conditionally on the Z ′ t s. However, strict exogeneity is a rather strong assumption. Under additional regularity conditions on the model, Chamberlain (1982) has shown that this assumption is equivalent to the non Granger-causality, i.e. Z t is independent of ( X s ) s ≤ t conditionally on ( Z s ) s ≤ -t 1 . It roughly means that the covariate process ( Z t ) t ∈ Z evolves in a totally autonomous way. In contrast, our exogeneity condition allows general covariates of the form Z t = H η ,η ( t t -1 , . . . ) with H a measurable function and a sequence (( η , ε t t )) t ∈ Z of i.i.d. random vectors, ε t being possibly correlated with η t . The error ε t can then still have an influence on future values of the covariates. For linear models, the two technical independence conditions discussed above between the noise and the covariate processes are often used as a distinction between weak and strict exogeneity. See for instance L¨tkepohl (2005), Section 10 2. u .

Let us mention that there exist additional concepts of exogeneity that are introduced and discussed in Engle et al. (1983), in particular a notion of weak exogeneity. However, this notion is related to the estimation of a specific parameter of the conditional distribution for the bivariate process ( X , Z t t ) and it is necessary to specify the joint dynamic of the process. Since we do not want to consider specific dynamics for the covariate process, we will not use it in this paper. Inclusion of exogenous regressors motivates our approach which is based on conditional average contraction conditions. But our results can be also applied without referring to these concepts of exogeneity, i.e. when (( Z t -1 , ε t )) t ∈ Z is a general stationary and ergodic process in (1). However, in the latter case, a closed form expression for the conditional distribution of X t given F t -1 cannot be obtained directly from the recursions (1). Our contribution is then the first one presenting a general framework for inclusion of covariates in nonlinear dynamics. Our results can be applied to any statistical procedure which require either ergodic properties or the use of some limit theorems developed from the notion of functional dependence introduced by Wu (2005). Since the existing literature already contains many asymptotic results of this type, we do not discuss specific applications.

The rest of this paper is organized as follows. In Section 2, we give our main results for defining stationary and ergodic solutions for recursions of type (1). In Section 3, we study weak dependence properties of the process using the functional dependence measure of Wu (2005). Many examples of nonlinear time series models satisfying our assumptions are given in Section 4 and we revisit some nonlinear dynamics discussed recently in the literature but we also consider new ones. A conclusion is given in Section 5. The proofs of our results are postponed to the last section of the paper.

## 2 General results

In this section, we state several results for controlling the convergence of the backward iterations in some L p spaces. We recall that for a random variable X and a real number p ≥ 1, the quantity ‖ X ‖ p = E 1 /p ( | X | p ) is called the L p -norm of X . Now let ( f t ) t ∈ Z be a sequence of random maps defined on a Polish space ( E,d ) and taking values in the same space. We assume for convenience that f t = F ( · , ζ t ) where ( ζ t ) t ∈ Z is a stochastic process taking values in another Polish space E ′ and F : E × E ′ → E is a measurable map. In connection with our initial example (1), we have

ζ t = ( Z t -1 , ε t ). In the latter case, we will assume throughout the paper that E ′ = E ′ 1 × E ′ 2 where E ′ 1 is a Borel subset of R e and E ′ 2 is another Polish space.

For s &lt; t , we set f t s = f t ◦ f t -1 s with the convention f t t = f t and f t -1 t ( x ) = x . Moreover, we consider a filtration ( F t ) t ∈ Z for which ( ζ t ) t ∈ Z is adapted.

## 2.1 Conditional contraction on average

We first give a general and useful result for getting a convergence in some L p spaces. For some real numbers p ≥ 1, L &gt; 0, κ ∈ (0 1) and an integer , m ≥ 1, we consider the two following assumptions.

- A1 There exists x 0 ∈ E such that sup t ∈ Z E [ d p ( f t ( x 0 ) , x 0 )] &lt; ∞ .
- A2 For every t ∈ Z , almost surely, the following inequalities hold for every ( x, y ) ∈ E 2 .

E [ d p ( f t ( x , f ) t ( y )) |F t -1 ] ≤ L d p p ( x, y ) and E [ d p ( f t + m -1 t ( x , f ) t + m -1 t ( y ) ) |F t -1 ] ≤ κ d p p ( x, y ) .

## Theorem 2. Suppose that Assumptions A1-A2 hold.

- 1. For every ( x, t ) ∈ E × Z , there exists a E -valued random variable X x t ( ) such that

$$\sup _ { t \in \mathbb { Z } } \| d \left ( X _ { t } ( x ), x _ { 0 } \right ) \| _ { p } < \infty, \ \sup _ { t \in \mathbb { Z } } \| d \left ( f _ { t - s } ^ { t } ( x ), X _ { t } ( x ) \right ) \| _ { p } = O \left ( \kappa ^ { s / m } \right ).$$

Moreover the sequence ( f t t -s ( x ) ) s ≥ 0 converges almost surely to X x t ( ) .

- 2. For x = y , we have P ( X x t ( ) = X y t ( )) = 0 . We then set X t = X x t ( ) .

/negationslash

/negationslash

- 3. The process (( X , ζ t t )) t ∈ Z is stationary and also ergodic if the process ( ζ t ) t ∈ Z is itself stationary and ergodic.
- 4. If ( Y t ) t ∈ Z is a non-anticipative process (i.e. Y t ∈ F t ) such that Y t = f t ( Y t -1 ) for t ∈ Z and sup j ∈ Z E [ d p ( Y , x j 0 )] &lt; ∞ , then Y t = X t a.s.

## Notes

- 1. The bounds given in Assumption A2 are required to hold for all ( x, y ) ∈ E 2 at the same time. Since a conditional expectation is only unique up to a set with measure 0 for the probability measure P , the bound given in Assumption A2 has to be understood in term of regular conditional distribution, i.e. there exists a regular version of the conditional distribution of ( ζ , . . . , ζ t t + m -1 ) given F t -1 . On Polish spaces, a regular version always exists. See Kallenberg (2006), Chapter 5.
- 2. When ζ t = ( Z t -1 , ε t ) forms a stationary process and F t = σ (( Z , ε j j ) : j ≤ t ), Theorem 2 guarantees existence and uniqueness of a stationary process possessing a moment of order p and solution of (1). However, stationarity of the covariate/error process is not required for applying this result. In particular, when ( ζ t ) t ∈ Z is non stationary, one can still define solutions of (1) provided that A1-A2 are satisfied. In general, these solutions will be non stationary and the classical law of large numbers is not valid. In this case, studying the asymptotic properties of some classical inferential procedures such as conditional likelihood estimation requires a specific analysis.
- 3. Setting p = 1 and d x, y ( ) = | x - | y o for x, y ∈ E = R and some o ∈ (0 1), one can consider , stochastic recursions (1) with heavy-tailed covariate processes.
- 4. When the process ( ζ t ) t ∈ Z is stationary and the recursions are initialized at time t = 0 with a given state x ∈ E , the probability distribution of the forward iterations f t 1 ( x ) coincides with the probability distribution of the backward iterations f 0 -t +1 ( x ). Since lim t →∞ -f 0 t +1 ( x ) = X 0 a.s., f t 1 ( x ) converges in distribution to X 0 . The same property holds true when the iterations are initialized with a random variable X 0 independent from ( ζ t ) t ≥ 1 . The main interest of the convergence of the backward iterations is to define the good random initialization X 0 = f 0 -∞ ( x ) in order to get a stationary process ( X t ) t ≥ 0 .
- 5. Assumption A2 is a conditional contraction property in L p which is crucial for getting the convergence of the backward iterations in L p norms. Relaxing this assumption by introducing a random coefficient κ t -1 instead of κ can be problematic for getting a similar result. We discuss this point below.

## 2.2 Comments on Assumption A2

Let us consider the dynamic (1), set f t ( x ) = F ( x, Z t -1 , ε t ), with E = R , p = 1 and assume that for every t , ε t is independent from F t -1 = σ (( Z , ε j j ) : j ≤ t -1). Then Assumption A2 is satisfied for m = 1 if and only if there exists κ ∈ (0 1) such that ,

$$\sup _ { z \in E _ { 1 } ^ { \prime } } \mathbb { E } \left | F \left ( x, z, \varepsilon _ { 0 } \right ) - F \left ( y, z, \varepsilon _ { 0 } \right ) \right | \leq \kappa | x - y |.$$

At a first sight, the latter condition is quite strong and it is natural to wonder if the following weaker assumption can be used, i.e. there exists a measurable function κ : G → (0 , ∞ ) such that for every z ∈ G ,

$$\mathbb { E } \left | F \left ( x, z, \varepsilon _ { 0 } \right ) - F \left ( y, z, \varepsilon _ { 0 } \right ) \right | \leq \kappa ( z ) | x - y |.$$

Of course, the challenging question concerns the convergence of the backward iterations when the function κ may take values larger than 1. However, a problem occurs for applying the successive contraction properties to the iterated random maps. Consider the iterations f t ◦ f t -1 . We have for ( x, y ) ∈ E 2 ,

$$\mathbb { E } \left | f _ { t } \circ f _ { t - 1 } ( x ) - f _ { t } \circ f _ { t - 1 } ( y ) \right | \ & = \ \mathbb { E } \left [ \mathbb { E } \left [ \left | f _ { t } \circ f _ { t - 1 } ( x ) - f _ { t } \circ f _ { t - 1 } ( y ) \right | \mathcal { F } _ { t - 1 } \right ] \right ] \\ & \leq \ \mathbb { E } \left [ \kappa ( Z _ { t - 1 } ) \left | f _ { t - 1 } ( x ) - f _ { t - 1 } ( y ) \right | \right ].$$

If the random variable κ Z ( t -1 ) depends on past values of the error ε t -j , j ≥ 1, it is stochastically dependent on the random map f t -1 and also on F t -2 . It is then not possible to use the contraction property of f t -1 unless the function κ can be bounded by a constant. To show that the successive iterations lose the memory with respect to initialization, this constant has to be smaller than 1. Of course, this does not prove that the convergence of the iterations in L 1 is not possible.

To show that the convergence of the backward iterations in L 1 is problematic, we now consider a map f t linear in x , a case for which an explicit solution is available. Accordingly, we assume that

$$f _ { t } ( x ) = \kappa ( Z _ { t - 1 } ) x + \varepsilon _ { t },$$

where the function κ is bounded but not necessarily by 1 and an integrable noise ε 0 . The dynamic is then given by an AR process with a random lag coefficient and it is widely known that the unique solution can be written as

$$X _ { t } = \sum _ { j \geq 1 } \prod _ { i = 1 } ^ { j } \kappa \left ( Z _ { t - i } \right ) \varepsilon _ { t - j } + \varepsilon _ { t },$$

provided that E log κ Z ( 0 ) &lt; 0. The series (2) converges almost surely. The L 1 -convergence is guaranteed from Theorem 2, as soon as κ := ‖ κ Z ( 0 ) ‖ ∞ &lt; 1, where for a random variable X , ‖ X ‖ ∞ denotes its suppremum norm. If κ ≥ 1, convergence in L 1 of the series 2 is much more difficult to get because of the possible stochastic dependence between the coordinates of the process ( Z t ) t ∈ Z . Let us first note that such a problem occurs in the non ergodic case, when Z t = Z 0 a.s. In this case, we have

$$\mathbb { E } \left ( \prod _ { i = 1 } ^ { j } \kappa \left ( Z _ { t - i } \right ) \left | \varepsilon _ { t - j } \right | \right ) = \mathbb { E } \left ( \left | \varepsilon _ { 0 } \right | \right ) \cdot \mathbb { E } \left ( \kappa ( Z _ { 0 } ) ^ { j } \right ) \\ \mathbb { \vee } \, \mathbb { D } \, \left | \right | \left | \mathcal { Z } \, \smallsetminus \, 1 \, \smallsetminus \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \cdot \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \, \longrightarrow \colon$$

and since E ( κ Z ( 0 ) j ) ≥ P ( κ Z ( 0 ) ≥ 1), one cannot get convergence of the series (2) if P ( κ Z ( 0 ) ≥ 1) &gt; 0. In what follows, we also stress that a similar problem of convergence also occurs in the ergodic case. To this end, set φ p ( ) = ‖ κ Z ( 0 ) ‖ p for p ≥ 1. The function φ is non decreasing and φ ( ∞ ) = ‖ κ Z ( 0 ) ‖ ∞ . Assumption A2 is satisfied as soon as ‖ κ Z ( 0 ) ‖ ∞ &lt; 1. It is then tempting to study the convergence of the solution only assuming that ‖ κ Z ( 0 ) ‖ p &lt; 1 but ‖ κ Z ( 0 ) ‖ q ≥ 1 for some 1 ≤ p &lt; q . However for any value of the pair ( p, q ), there always exists an example of a process ( κ Z ( t )) t ∈ Z such that the series (2) is not converging in L 1 . To this end, assuming without loss of generality that q is an integer, we define κ z ( ) = z and Z t -1 = a t -1 · · · a t -q where ( a t ) t ∈ Z is a process of i.i.d. nonnegative random variables, independent of ( ε t ) t ∈ Z , and such that ‖ a 0 ‖ p &lt; 1 and ‖ a 0 ‖ q ≥ 1. Since for j ≥ q ,

$$\prod _ { i = 1 } ^ { j } \kappa \left ( Z _ { t - i } \right ) = \prod _ { i = 1 } ^ { q - 1 } a _ { t - i } ^ { i } \cdot \prod _ { i = q } ^ { j + 1 } a _ { t - i } ^ { q } \prod _ { i = 0 } ^ { q - 2 } a _ { t - j - k + i } ^ { i + 1 },$$

we find

$$\mathbb { E } \left ( \prod _ { i = 1 } ^ { j } \kappa \left ( Z _ { t - i } \right ) \right ) = \prod _ { i = 1 } ^ { q - 1 } \mathbb { E } ^ { 2 } \left ( a _ { 0 } ^ { i } \right ) \cdot \mathbb { E } ^ { j - q + 2 } \left ( a _ { 0 } ^ { q } \right ). \\ \text{$s$ overwattion does not remove to $0$ when $i$} \rightarrow \sim \text{and $t$}.$$

Hence the previous expectation does not converge to 0 when j → ∞ and the series (2) cannot

converge in L 1 . The analysis of this linear case enlightens that a tail condition on κ Z ( t ) is not sufficient for getting this kind of convergence. In particular, the dependence structure of the process ( κ Z ( t )) t ∈ Z is also of major importance. This contrasts with AR processes with i.i.d. random coefficients, since in this case the condition E κ Z ( 0 ) &lt; 1 is necessary and sufficient for the convergence of the series (2) in L 1 . However, imposing an independence assumption on the covariate process is not reasonable.

In the next section, we show that one can investigate a different mode of convergence for the iterations in model (1) and which allow to relax Assumption A2 . However, it is necessary to impose a strict exogeneity assumption on the covariate process.

## 2.3 An additional result for strictly exogenous regressors

In this subsection, we consider specifically equation (1) when the covariate process is independent of the error process. In this case, conditionally on Z , the process is a time-inhomogeneous Markov chain. The terminology Markov chain in random environments is often used in the literature. See for instance Stenflo (2001). The following result will not be central in the rest of the paper because substantial efforts could be needed to derive moment and weak dependence properties for the corresponding solution and it could be also difficult to obtain explicit conditions for dealing with higher-order autoregressive processes. This is why we only provide a result when the f t 's satisfied a one-step contraction (i.e. m = 1 in A2 ). We assume that there exist a real number p ≥ 1 and a state x 0 ∈ E such that the three following conditions are fulfilled.

A0 The process Z := ( Z t ) t ∈ Z is stationary and ergodic, ( ε t ) t ∈ Z is a process of i.i.d. random variables taking values in E ′ 2 and is independent of Z .

- A1' For every z ∈ E ′ 1 , we have E [ d p ( F x , z, ε ( 0 1 ) , x 0 )] &lt; ∞ .

A2' There exists a measurable function κ : E ′ → (0 , ∞ ) satisfying E ( log + κ Z ( 0 ) ) &lt; ∞ , E log κ Z ( 0 ) &lt; 0 and such that for every ( x, y ) ∈ E 2 ,

$$\mathbb { E } \left [ d ^ { p } \left ( F ( x, z, \varepsilon _ { 0 } ), F ( y, z, \varepsilon _ { 0 } ) \right ) \right ] \leq \kappa ^ { p } ( z ) d ^ { p } ( x, y ).$$

Moreover,

$$\mathbb { E } \left [ \log ^ { + } \int d ^ { p } \left ( x _ { 0 }, F ( x _ { 0 }, Z _ { 0 }, u ) \right ) d \mathbb { P } _ { \epsilon _ { 1 } } ( u ) \right ] < \infty.$$

We remind the notation f t ( x ) = F ( x, Z t -1 , ε t ). Here, we set F t = σ (( Z , ε j j ) : j ≤ t ) and E [ X Z | ] will denote the expectation of a random variable X conditionally on the covariate process ( Z t ) t ∈ Z .

## Theorem 3. Suppose that Assumptions A0,A1'-A2' hold.

- 1. For every ( x, t ) ∈ E × Z , there exists a random variable X x t ( ) such that E [ d p ( X x ,x t ( ) 0 ) | Z &lt; ] ∞ a.s. and lim s →∞ E [ d p ( f t t -s ( x , X ) t ( x ) ) | Z ] = 0 a.s. The sequence ( f t t -s ( x ) ) s ≥ 0 also converges almost surely to X x t ( ) .
- 2. For x = y , we have P ( X x t ( ) = X y t ( )) = 0 . We then set X t = X x t ( ) .

/negationslash

/negationslash

- 3. The process (( X , Z t t )) t ∈ Z is stationary and ergodic.
- 4. If ( Y t ) t ∈ Z is a non-anticipative process (i.e. Y t ∈ F t ) such that (( Y , Z t t )) t ∈ Z is stationary, for every t ∈ Z , Y t = f t ( Y t -1 ) and E [ d p ( Y , x 0 0 ) | Z , Z 0 -1 , . . . ] &lt; ∞ a.s., then Y t = X t a.s.

## Notes

- 1. The contraction inequality in Assumption A2' can be restated as

$$\mathbb { E } \left [ d ^ { p } \left ( f _ { t } ( x ), f _ { t } ( y ) \right ) | Z \right ] \leq \kappa ^ { p } ( Z _ { t - 1 } ) d ^ { p } ( x, y ) \text{ a.s.}$$

It is then another example of contraction in conditional average.

- 2. In our context, our result can be seen as an improvement of Theorem 1 given in Stenflo (2001) for Markov chains in random environments. In particular, we do not assume a uniform

contraction with respect to the environment which is given by the exogenous process Z in our random maps f t .

## 2.4 Example

We compare the contraction conditions necessary to apply Theorem 1, Theorem 2 or Theorem 3 on a specific example. Let (( ε , Z t t )) t ∈ Z be a stationary sequence of pair of random variables and for t ∈ Z , set F t = σ (( Z , ε j j ) : j ≤ t ). Assume that E ( ε 2 1 |F 0 ) = 1 and E ( ε 1 |F 0 ) = 0. For i = 1 2, , let a i : E ′ 1 → R and b i : E ′ 1 → R + be some measurable maps such that E log + a i ( Z 0 ) &lt; ∞ and E log + b i ( Z 0 ) &lt; ∞ . We consider the following AR-ARCH model with functional coefficients

$$X _ { t } = a _ { 0 } ( Z _ { t - 1 } ) + a _ { 1 } ( Z _ { t - 1 } ) X _ { t - 1 } + \varepsilon _ { t } \sqrt { b _ { 0 } ( Z _ { t - 1 } ) + b _ { 1 } ( Z _ { t - 1 } ) X _ { t - 1 } ^ { 2 } }.$$

Here we set E = R and d x, y ( ) = | x - | y and p = 2. Setting σ z, x ( ) = √ b 0 ( z ) + b 1 ( z x ) 2 , one can note that | σ z, x ( ) -σ z, y ( ) | ≤ √ b 1 ( z ) | x - | y .

- 1. To apply Theorem 1, we compute the Lipschitz constant c f ( 1 ) of the random map f 1 . We

have

∣ Making the change of variable v = sign ε ( 1 ) sign ( a 1 ( Z 0 )) v , we have

$$c ( f _ { 1 } ) = \sup _ { v \in \mathbb { R } } | f _ { 1 } ^ { \prime } ( v ) | = \sup _ { v \in \mathbb { R } } \left | a _ { 1 } ( Z _ { 0 } ) + \frac { \varepsilon _ { 1 } b _ { 1 } ( Z _ { 0 } ) v } { \sqrt { b _ { 0 } ( Z _ { 0 } ) + b _ { 1 } ( Z _ { 0 } ) v ^ { 2 } } } \right |. \\ ; \text{the change of variable $\overline{v}=sign(\varepsilon_{1})sign\left(a_{1}(Z_{0})\right)v,\text{ we have}$$

$$c ( f _ { 1 } ) = \sup _ { \overline { v } \in \mathbb { R } } \left [ | a _ { 1 } ( Z _ { 0 } ) | + \frac { | \varepsilon _ { 1 } | \, b _ { 1 } ( Z _ { 0 } ) \overline { v } } { \sqrt { b _ { 0 } ( Z _ { 0 } ) + b _ { 1 } ( Z _ { 0 } ) \overline { v } ^ { 2 } } } \right |. \\ \intertext { 4 c ( f _ { 1 } ) } = | a _ { 1 } ( Z _ { 0 } ) | + \sqrt { b _ { 1 } ( Z _ { 0 } ) } | \varepsilon _ { 1 } | \text{ and Theorem 1 annlies as soo}$$

∣ We then obtain c f ( 1 ) = | a 1 ( Z 0 ) | + √ b 1 ( Z 0 ) | ε 1 | and Theorem 1 applies as soon as

$$\mathbb { E } \log c ( f _ { 1 } ) < 0.$$

## 2. To apply Theorem 2, note that

$$\mathbb { E } \left [ | f _ { t } ( x ) - f _ { t } ( y ) | ^ { 2 } \, | \mathcal { F } _ { t - 1 } \right ] \ & = \ \ a _ { 1 } ( Z _ { t - 1 } ) ^ { 2 } ( x - y ) ^ { 2 } + ( \sigma \, ( Z _ { t - 1 }, x ) - \sigma \, ( Z _ { t - 1 }, x ) ) ^ { 2 } \\ & \leq \ \left ( a _ { 1 } ( Z _ { t - 1 } ) ^ { 2 } + b _ { 1 } ( Z _ { t - 1 } ) \right ) \cdot | x - y | ^ { 2 }.$$

One can then show that Theorem 2 applies as soon as

$$\kappa ^ { 2 } \coloneqq \sup _ { z } \left ( a _ { 1 } ( z ) ^ { 2 } + b _ { 1 } ( z ) \right ) < 1.$$

- 3. If the two processes ( ε t ) t ∈ Z and ( Z t ) t ∈ Z are independent and the ε ′ t s are i.i.d., Theorem 3 applies as soon as

$$\mathbb { E } \log \left ( a _ { 1 } ( Z _ { 0 } ) ^ { 2 } + b _ { 1 } ( Z _ { 0 } ) \right ) < 0.$$

Under the strict exogeneity assumption, we note that (5) is weaker than (4). However, (4) ensures the existence of a second order moment for the solution whereas (5) only guarantees that E ( X Z 2 t | ) &lt; ∞ a.s. On the other hand, (3), which only ensures existence of a stationary solution, is not necessarily weaker than (4) or (5). For instance, if the noise process has a Rademacher distribution, P ( ε 1 = 1) = P ( ε 1 = -1) = 1 / 2 and the functional coefficients are constant, (3) writes | a 1 | + √ b 1 &lt; 1 which is more restrictive than (4) or (5). But if a 1 is identically equal to 0, (3) writes as 1 2 E log b 1 ( Z 0 ) + E log | ε 0 | &lt; 0 which is weaker than (5), since from Jensen's inequality, we have

$$\mathbb { E } \log | \varepsilon _ { 0 } | \leq \log \mathbb { E } | \varepsilon _ { 0 } | \leq \log \mathbb { E } ^ { 1 / 2 } \left ( \varepsilon _ { 0 } ^ { 2 } \right ) = 0. \\ \dots \, \cdot \quad \cdot \quad \cdot$$

## 2.5 A result for higher-order autoregressions

In this subsection, we extend Theorem 2 to higher-order autoregressive processes. We only consider stationary processes in this part. The main result, Theorem 4, is particularly interesting for multivariate autoregressions for which Lipschitz type properties can be obtained equation by equation. See Section 4 for an application of Theorem 4 to various examples.

For a given real number 0 &lt; o ≤ 1, we define the distance ∆ on R , by ∆( u, v ) = | u - | v o for u, v ∈ R . Let E be a subset of R k and ‖ · ‖ an arbitrary norm on R k . Our aim is to study existence of solutions for the following recursive equations:

$$X _ { t } = F \left ( X _ { t - 1 }, \dots, X _ { t - q }, \zeta _ { t } \right ), \quad t \in \mathbb { Z },$$

where F : E q × E ′ → E is a measurable function. Note that one can always associate a random map f t on E q to the dynamic (6). To this end, for t ∈ Z and x , . . . , x 1 q ∈ E , we set

$$f _ { t } ( x _ { 1 }, \dots, x _ { q } ) = \left ( F \left ( x _ { 1 }, \dots, x _ { q }, \zeta _ { t } \right ), x _ { 1 }, \dots, x _ { q - 1 } \right ).$$

We first introduce additional notations. We denote by M k the set of square matrices with real coefficients and k rows and if A ∈ M k , ρ A ( ) the spectral radius of the matrix A . Moreover, for x, y ∈ R k and p ≥ 1, the vector (∆ p ( x , y 1 1 ) , . . . , ∆ ( p x , y k k )) ′ will be denoted by ∆ p vec ( x, y ). Finally, we introduce a partial order relation /precedesequal on R k and such that x /precedesequal y means x i ≤ y i for i = 1 , . . . , k .

The following assumptions will be needed.

B1 The process ( ζ t ) t ∈ Z is stationary and ergodic adapted to a filtration ( F t ) t ∈ Z .

$$\text{B2 For any } y \in E ^ { q }, \, \mathbb { E } \left [ \left \| F \left ( y, \zeta _ { 1 } \right ) \left \| ^ { o p } \right ] < \infty.$$

B3 There exist some matrices A , . . . , A 1 q ∈ M k with nonnegative elements, satisfying ρ ( A 1 + + · A q ) &lt; 1 and such that for y, y ′ ∈ E q ,

$$\mathbb { E } \left [ \Delta _ { v e c } ^ { p } \left ( F ( y, \zeta _ { t } ), F ( y ^ { \prime }, \zeta _ { t } ) \right ) | \mathcal { F } _ { t - 1 } \right ] \preceq \sum _ { i = 1 } ^ { q } A _ { i } \Delta _ { v e c } ^ { p } ( y _ { i }, y _ { i } ^ { \prime } ).$$

Though the following result is stated for an arbitrary pair o ∈ (0 1) , , p ≥ 1, the two interesting cases are o ∈ (0 1) , , p = 1 and o = 1 , p &gt; 1.

Theorem 4. Suppose that Assumptions B1-B3 hold. There then exists a unique stationary and non-anticipative process ( X t ) t ∈ Z solution of (6) and such that E [ ‖ X t ‖ op ] &lt; ∞ . Moreover, this process is ergodic.

## Notes

- 1. If ζ t = ( Z t -1 , ε t ) with ε t independent of F t -1 = σ (( Z , ε j j ) : j ≤ t ) and k = o = 1, Assumption B3 writes

$$\sup _ { z } \mathbb { E } \left [ \left | F ( y _ { 1 }, \dots, y _ { q }, z, \varepsilon _ { 1 } ) - F ( y ^ { \prime } _ { 1 }, \dots, y ^ { \prime } _ { q }, z, \varepsilon _ { 1 } ) \right | ^ { p } \right ] \leq \sum _ { i = 1 } ^ { q } A _ { i } | y _ { i } - y ^ { \prime } _ { i } | ^ { p }, \\ \intertext { w h e m } \lambda / \mathcal { M } ^ { q } \quad \Lambda \, \cdot \, - \, \mathcal { M } ^ { q } \quad \Lambda \, \cdot \, 1 \, \text{ this modification } \lambda \, \text{ andation } \lambda \, \text{ andation } \lambda \, \text{ andation } \lambda \, \text{ andation } \lambda \, \text{ } \lambda \,.$$

with ρ ( ∑ q i =1 A i ) = ∑ q i =1 A &lt; i 1. This provides a quite simple criterion for application to autoregressive processes.

- 2. In the spirit of Section 2.4, the previous criterion can be checked for models with varying parameters, directly constructed from smooth parametric autoregressive processes. Consider the model

$$Y _ { t } = \overline { F } _ { \theta } \left ( Y _ { t - 1 }, \dots, Y _ { t - q }, \varepsilon _ { t } \right ), \ \ t \in \mathbb { Z }, \theta \in \Theta \subset \mathbb { R } ^ { \overline { e } }.$$

If the model

$$\mathbb { E } \left [ | \overline { F } _ { \theta } ( y _ { 1 }, \dots, y _ { q }, \varepsilon _ { 1 } ) - \overline { F } _ { \theta } ( y _ { 1$$

$$Y _ { t } = \overline { F } _ { \theta ( Z _ { t - 1 } ) } \left ( Y _ { t - 1 }, \dots, Y _ { t - q }, \varepsilon _ { t } \right )$$

satisfies B3 as soon as ∑ q i =1 ‖ A i ‖ ∞ &lt; 1. We then obtain a model with exogenous covariates by replacing parameter θ with a varying parameter θ ( Z t -1 ) where θ : R e → Θ is a measurable map. See also the note after Proposition 3 for a discussion.

## 3 Functional dependence measure

The functional dependence measure has been introduced by Wu (2005) and is particularly interesting for autoregressive processes which are not necessarily strong mixing or for which getting strong mixing conditions requires additional regularity conditions on the noise distribution. The single requirement is to get a Bernoulli shift representation of the stochastic process of interest, i.e. X t = H ξ ,ξ ( t t -1 , . . . ) where ( ξ t ) t ∈ Z is a sequence of i.i.d. random variables taking values in a mea-

surable space ( G, G ). The functional dependence measure is expressed in terms of some coefficients which evaluate for t ≥ 0 the L p -distance between X t and a copy X t , obtained by replacing ξ 0 with ξ ′ 0 , ξ ′ 0 following the same distribution as ξ 0 and being independent from the sequence ( ξ t ) t ∈ Z . See below for the definition of these coefficients. Most of the limit theorems and deviation inequalities have been derived under such dependence measures. See for instance Wu (2007) and Wu and Wu (2016). Such asymptotic results have been applied to various statistical problems. See for instance Wu et al. (2010) for kernel estimation for time series, Xiao and Wu (2012) for covariance estimation or Liu and Wu (2010a) for spectral density estimation. The notion of functional dependence is then an attractive alternative to the usual strong mixing when the process is defined by stochastic recursions. Our aim in this section is to show that under the assumptions of Theorem 2 or Theorem 4, when the process ( ζ t ) t ∈ Z or more specifically the covariate process ( Z t ) t ∈ Z satisfies this kind of dependence, the functional dependence measure of the solution ( X t ) t ∈ Z can be controlled. We will then provide a new wide class of examples for which the aforementioned references provide an important number of statistical applications. In this section, we assume that the state space E is a subspace of R k and the distance d is given by d x, y ( ) = | x - | y o where | · | is an arbitrary norm on R k and 0 &lt; o ≤ 1.

## 3.1 Dependence coefficients for general iterated random functions

Assume that the process ζ has a Bernoulli shift representation, i.e. ζ t = H ξ ,ξ ( t t -1 , . . . ) for some measurable map H : G N → E ′ and ( ξ t ) t ∈ Z is a sequence of i.i.d., G -valued random variables. We then have for every ( x, t ) ∈ E × Z ,

$$f _ { t } ( x ) = F \left ( x, H \left ( \xi _ { t }, \xi _ { t - 1 }, \dots \right ) \right )$$

/negationslash and the map f t has itself a Bernoulli shift representation. To define the functional measure coefficients, we then define a new sequence ( ξ t ) t ∈ Z such that ξ 0 = ξ ′ 0 and ξ t = ξ t for t = 0. Here ξ ′ 0 is a copy of ξ 0 which is assumed to be independent from ( ξ t ) t ∈ Z . Moreover, for t &gt; 0, let

$$\overline { f } _ { t } ( x ) = F ( x, H \left ( \overline { \xi } _ { t }, \overline { \xi } _ { t - 1 }, \dots \right ) \right ). \\ \quad \quad \quad$$

we define for t ≥ 0 and p ≥ 1,

$$\theta _ { p, t } & = \mathbb { E } ^ { 1 / p } \left [ d \left ( f _ { - \infty } ^ { t }, \overline { f } _ { - \infty } ^ { t } \right ) ^ { p }$$

Moreover, for h ∈ N , let Θ p,h = ∑ t ≥ h θ p,t . Two cases of interest are p = 1 and o = 1 , p &gt; 1.

Our aim is to get an upper bound for the functional dependence coefficients Θ p,h . To this end, we add other assumptions. Here we set for t ∈ Z , F t = σ ( ξ t -j : j ≥ 0).

A3 There exists a measurable function S : E → R + and r, s ≥ p such that r -1 + s -1 = p -1 and E [ S f ( 0 -∞ ) s ] &lt; ∞ and for all x ∈ E and t ≥ 1,

$$\mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( \overline { f } _ { t } ( x ), f _ { t } ( x ) \right ) ^ { p } | \sigma ( \xi _ { 0 } ^ { \prime } ) \vee \mathcal { F } _ { t - 1 } \right ] \leq S ( x ) H _ { t - 1 },$$

where H t -1 is a random variable measurable with respect to σ ξ ( ′ 0 ) ∨F t -1 and such that E | H t -1 | r &lt; ∞ .

An immediate consequence of Assumption A4 is that for any random variable V t -1 , measurable with respect to σ ξ ( ′ 0 ) ∨ F t -1 , we have E [ d p ( f t ( V t -1 ) , f t ( V t -1 ) )] ≤ E [ S V ( t -1 ) p H p t -1 ] and from H¨ older's inequality, we get

$$\| d \left ( \overline { f } _ { t } ( V _ { t - 1 } ), f _ { t } ( V _ { t - 1 } ) \right ) \| _ { p } \leq \| S ( V _ { t - 1 } ) \| _ { s } \| H _ { t - 1 } \| _ { r }.$$

When ζ t = ( Z t -1 , ε t ), with Z t taking values in a Borel subset E ′ 1 of R e and ε t taking values in a Polish space E ′ 2 , we will still denote by |·| an arbitrary norm on R e and we also set d z, z ( ′ ) = | z -z ′ | o for z, z ′ ∈ R e to avoid additional notations. We will use two specific assumptions.

A3' There exists a measurable function S : E → R + and r, s ≥ p such that r -1 + s -1 = p -1 and E [ S f ( 0 -∞ ) s ] &lt; ∞ , E | Z 0 | ro &lt; ∞ and for all x ∈ E and t ≥ 1,

$$\mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( F \left ( x, z, \varepsilon _ { 0 } \right ), F \left ( x, z ^ { \prime }, \varepsilon _ { 0 } \right ) \right ) \right ] \leq S ( x ) d ( z, z ^ { \prime } ),$$

A4 Let ( η t ) t ∈ Z be a sequence of random variables taking values in a measurable space ( G , 1 G 1 ) and such that Z t = H ′ ( η , η t t -1 , . . . ) for a measurable function H ′ . Moreover, setting ξ t = ( ε , η t t ),

we assume that ( ξ t ) t ∈ Z is a sequence of i.i.d. of random variables taking values in G = E ′ 2 × G 1 .

From A4 , we have the specific representation

$$\zeta _ { t } = H \left ( \xi _ { t }, \xi _ { t - 1 }, \dots \right ) \coloneqq \left ( H ^ { \prime } ( \eta _ { t - 1 }, \eta _ { t - 2 }, \dots ), \varepsilon _ { t } \right ), \ \ t \in \mathbb { Z }.$$

The map H takes values in E ′ = E ′ 1 × E ′ 2 . Note that our formulation allows the covariate process to have a general form, including a VARMA or GARCH process among others.

Proposition 1. 1. Suppose that Assumptions A1-A4 hold. For any h ≥ 2 , there then exists

C 1 &gt; 0 not depending on h such that

$$\Theta _ { p, h } & \leq C _ { 1 } \left [ \kappa ^ { h / m } + \sum _ { i = 0 } ^ { h - 1 } \kappa ^ { i / m } \eta _ { r, h - i } + \sum _ { i \geq h } \kappa ^ { i / m } \eta _ { r, 1 } \right ], \\ \dots \colon \iota.$$

with

$$\eta _ { r, j } = \sum _ { t \geq j } \| H _ { t - 1 } \| _ { r }, \quad j \geq 1. \\. \quad \. \quad \. \quad \.$$

In particular, if η r, 1 &lt; ∞ , there exists C &gt; 0 , not depending on h , such that

$$\Theta _ { p, h } \leq C \left [ \kappa ^ { h / m } + \sum _ { i = 0 } ^ { h - 1 } \kappa ^ { i / m } \eta _ { r, h - i } \right ].$$

- 2. Suppose that Assumptions A1-A2-A4 and A3' hold. Then we get the bound (8) with η r,j = Θ r,j -1 ( Z ) .

## Notes

- 1. Let us comment on Assumption A4 . Under this assumption, the η ′ t s are i.i.d. as well as the ε ′ t s and for any t ∈ Z , ε t is independent from F t -1 = σ ( ξ s : s ≤ t - } 1 . Note that we allow simultaneous dependence between ε t and η t . For instance, we can set η t = K ε , U ( t t ) where K is a measurable function and U is a sequence of i.i.d. random variables, independent from the sequence ε . This assumption is then more flexible than the complete independence between the two error processes ε and η , which implies strict exogeneity.

- 2. It can happen that our assumptions are satisfied with some p , leading to an upper bound for the functional dependence coefficients θ p,t , while it is required a condition on θ q,t or Θ q,h for q &gt; p for applying some limit theorems or statistical results. This is still possible if one can prove finiteness of higher-order moments for the solution, e.g. if E [ | f 0 -∞ | q o ′ ] &lt; ∞ for some q ′ &gt; q . Indeed, from H¨lder's inequality, we have o

$$\theta _ { q, t } \leq \theta _ { p, t } ^ { \frac { p ( q ^ { \prime } - q ) } { q ( q ^ { \prime } - p ) } } \theta _ { q ^ { \prime }, t } ^ { \frac { q ^ { \prime } ( q - p ) } { q ( q ^ { \prime } - p ) } }.$$

Moreover, θ q ,t ′ ≤ ‖ 2 ∣ ∣ f 0 -∞ ∣ ∣ o ‖ q ′ .

## 3.2 Dependence coefficients for higher-order autoregressions with exogenous covariates

Here, we revert to higher-order autoregressions considered in Section 2.5. We consider directly the case ζ t = ( Z t -1 , ε t ) with Assumption A4 being satisfied. Additionally to Assumptions B1-B2-B3 and A4 , the following assumption, which is the analogue of A3' , will be needed.

B4 If r &gt; 0 and s ∈ R + ∪ {∞} are such that r -1 + s -1 = p -1 , there exists a measurable function S : E q → R such that E [ S X ,... , X ( q 1 ) s ] &lt; ∞ , E | Z 0 | or &lt; ∞ and for all z, z ′ ∈ E ′ 1 and x , . . . , x 1 q ∈ E q ,

$$\mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( F ( x _ { 1 }, \dots, x _ { q }, z, \varepsilon _ { 0 } ), F ( x _ { 1 }, \dots, x _ { q }, z ^ { \prime }, \varepsilon _ { 0 } ) \right ) \right ] \leq S ( x _ { 1 }, \dots, x _ { q } ) d ( z, z ^ { \prime } ).$$

The result is the following.

Proposition 2. Suppose that Assumptions B2 -B4 and A4 hold. There then exists C &gt; 0 and ρ ∈ (0 1) , , such that for all h ≥ 1 , Θ p,h ( X ) ≤ C ρ [ h + ∑ h i =1 ρ i Θ r,h -i ( Z ) ] .

Note. From the upper bound given in Proposition 2, we note that the decay of Θ p,h ( X ) is polynomial (respectively geometric) in h if the decay of Θ r,h is polynomial (respectively geometric) in h .

## 3.3 A central limit theorem

To illustrate the usefulness of our results, we give below a central limit theorem for partial sums

$$S _ { n } = \sum _ { i = 1 } ^ { n } f \left ( X _ { t }, X _ { t - 1 }, \dots, X _ { t - k } \right ),$$

where f is some real-valued measurable function and ( X t ) t ∈ Z is a stochastic process solution of X t = F ( X t -1 , . . . , X t -q , ζ t ) and the assumptions of either Proposition 1 or Proposition 2 are satisfied. The following result, which is a straightforward corollary of the invariance principle given in Wu (2005), is not the most general as possible. In particular, when a moment of order greater than p is available for the stationary solution, different assumptions on the function f could be used.

Theorem 5. Suppose that either Assumptions A1-A4 or Assumptions B2-B3-B4-A4 hold true for some p &gt; 2 and Θ p, 0 ( X &lt; ) ∞ . If there exists C &gt; 0 and 0 ≤ /lscript ≤ p -2 2 such that for x , x i ′ i ∈ E , 0 ≤ i ≤ k ,

$$\left | f \left ( x _ { 0 }, \dots, x _ { k } \right ) - f \left ( x _ { 0 } ^ { \prime }, \dots, x _ { k } ^ { \prime } \right ) \right | \leq C \left [ 1 + \sum _ { i = 0 } ^ { k } \left ( | x _ { i } | ^ { o \ell } + | x _ { i } ^ { \prime } | ^ { o \ell } \right ) \right ] \cdot \sum _ { i = 0 } ^ { k } d ( x _ { i }, x _ { i } ^ { \prime } ). \\ \intertext { T h e n n e h a n e t h e n e c e }$$

Then we have the weak convergence

$$\frac { 1 } { \sqrt { n } } \left ( S _ { n } - \mathbb { E } S _ { n } \right ) \Rightarrow \mathcal { N } \left ( 0, \sigma ^ { 2 } \right ),$$

with σ 2 = ∑ j ∈ Z Cov ( Y , Y j 0 ) , Y t = f ( X , . . . , X t t -k ) .

## 4 Examples

## 4.1 CHARN models

In this section, we consider conditional heteroscedastic autoregressive nonlinear (CHARN) models such as in H¨rdle and Tsybakov (1997) or H¨rdle et al. (1998) but that can encompass exogenous a a

regressors. More precisely, we consider the dynamic

$$Y _ { t } = f _ { 1 } \left ( Y _ { t - 1 }, \dots, Y _ { t - q }, Z _ { t - 1 } \right ) + \varepsilon _ { t } f _ { 2 } \left ( Y _ { t - 1 }, \dots, Y _ { t - q }, Z _ { t - 1 } \right ),$$

where q is a positive integer and f , f 1 2 : R q × E ′ 1 → R are measurable functions. In order to study stationary solutions of the recursions (9), the following assumptions will be needed.

CH1 The process (( Z , ε t t )) t ∈ Z is stationary and ergodic.

CH2 For j = 1 2, there exist measurable functions , a i,j : R d × R + , 1 ≤ i ≤ q such that

$$\left | f _ { j } ( y _ { 1 }, \dots, y _ { q }, z ) - f _ { j } ( y _ { 1 } ^ { \prime }, \dots, y _ { q } ^ { \prime }, z ) \right | \leq \sum _ { i = 1 } ^ { q } a _ { i, j } ( z ) \left | y _ { i } - y _ { i } ^ { \prime } \right |. \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad \quad \quad$$

CH3 There exist a real number p ≥ 1 such that ‖ ε 1 ‖ p &lt; ∞ and r, s ≥ p such that r -1 + s -1 = p -1 , s can be infinite, E | Z 0 | r &lt; ∞ and two functions L , L 1 2 defined on R q and such that for j = 1 2, , y , . . . , y 1 q ∈ R and z, z ′ ∈ E ′ 1 ,

$$& | f _ { j } ( y _ { 1 }, \dots, y _ { q }, z ) - f _ { j } ( y _ { 1 }, \dots, y _ { q }, z ^ { \prime } ) | \leq L _ { j } ( y _ { 1 }, \dots, y _ { q } ) | z - z ^ { \prime } |. \\ & \\ - \quad. & \quad. - \quad. - \quad. - \quad. - \quad. - \quad. - \quad. -$$

For t ∈ Z and i = 1 , . . . , q , we set c i,t = a i, 1 ( Z t -1 ) + a i, 2 ( Z t -1 ) | ε t | , We then define a sequence of random matrices A = ( A t ) t ∈ Z by

$$A _ { t } = \begin{pmatrix} c _ { 1, t } & c _ { 2, t } & \cdots & c _ { q, t } \\ & & 0 \\ I _ { q - 1 } & & \vdots \\ & & 0 \end{pmatrix}. \\ \text{the Lyapunov exponent of the sequence } A$$

Finally, we denote by χ ( A ) the Lyapunov exponent of the sequence A , i.e.

$$\chi ( \mathbf A ) = \lim _ { n \to \infty } \frac { \mathbb { E } \left [ \log \left \| A _ { n } \cdots A _ { 1 } \right \| \right ] } { n },$$

where ‖ · ‖ is an arbitrary norm on the space of square matrices of size q × q .

## Proposition 3. Suppose that Assumptions CH1-CH2 hold.

- 1. Suppose that χ ( A ) &lt; 0 . There then exists a unique stationary process ( Y t ) t ∈ Z solution of (9) which is also ergodic.
- 2. Assume additionally that for every t ∈ Z , ε t is independent from F t -1 . If there exist x ∈ R q , p ≥ 1 such that f 1 ( x, Z 0 ) + ε f 1 2 ( x, Z 0 ) ∈ L p and

$$\sum _ { i = 1 } ^ { q } \sup _ { z } \| a _ { i, 1 } ( z ) + a _ { i, 2 } ( z ) | \varepsilon _ { 1 } | \| _ { p } < 1,$$

there then exists a unique stationary and non-anticipative process solution of (9) which is also ergodic and such that E | Y 1 | p &lt; ∞ .

- 3. Assume furthermore that Assumptions ( CH3 ) and ( A4 ) hold true with E L s j ( Y , . . . , Y q 1 ) &lt; ∞ for j = 1 2 , . There then exists C &gt; 0 and ρ ∈ (0 1) , , such that for all h ≥ 1 , Θ p,h ( X ) ≤ C [ ρ h + ∑ h i =1 ρ i Θ r,h -i ( Z ) ] .

## Notes

- 1. Our results can be useful for dealing with models with functional coefficients in the spirit of the example given in Section 2.4. See also the notes after the statement of Theorem 4. For j = 1 2, let , m ( j ) θ : R q → R be some functions depending on some parameters θ ∈ R e and such that

for some nonnegative real numbers d i,j ( θ ), 1 ≤ i ≤ q . If θ is replaced by a function θ ( ) · : R d → R e and f j ( y , . . . , y 1 q , z ) = m ( j ) θ z ( ) ( y , . . . , y 1 q ) for j = 1 2, one can then consider some , standard autoregressive processes and obtain a version with functional parameters depending on exogenous covariates. For instance, threshold autoregressions or power-ARCH volatility,

$$\left | m _ { \theta } ^ { ( j ) } ( y _ { 1 }, \dots, y _ { q } ) - m _ { \theta } ^ { ( j ) } ( y _ { 1 } ^ { \prime }, \dots, y _ { q } ^ { \prime } ) \right | \leq \sum _ { i = 1 } ^ { q } d _ { i, j } ( \theta ) \left | y _ { i } - y _ { i } ^ { \prime } \right |, \\ \text{e nonnegative real numbers $d_{i,j}(\theta), 1 \leq i \leq q$.  If $\theta$ is replaced by a functio}$$

$$m _ { \theta } ^ { ( 1 ) } ( y _ { 1 }, \dots, y _ { q } ) = \theta _ { 0 } + \sum _ { i = 1 } ^ { q } \left ( \theta _ { i } y _ { i } ^ { + } + \theta _ { i + q } y _ { i } ^ { - } \right ), \quad m _ { \theta } ^ { ( 2 ) } ( y _ { 1 }, \dots, y _ { q } ) = \left ( \theta _ { 0 } + \sum _ { i = 1 } ^ { q } \theta _ { i } | y _ { i } | ^ { \delta } \right ) ^ { 1 / \delta },$$

where x + and x -denotes respectively the positive part and the negative part of a real number x and δ ≥ 1.

- 2. When ε t is not necessarily independent from F t -1 , (10) can be replaced with the following more abstract condition. There exists η ∈ (0 1) such that ,

$$\sum _ { i = 1 } ^ { q } \mathbb { E } ^ { 1 / p } \left [ ( a _ { i, 1 } ( Z _ { 0 } ) + a _ { i, 2 } ( Z _ { 0 } ) | \varepsilon _ { 1 } | ) ^ { p } \left | \mathcal { F } _ { 0 } \right ] \leq 1 - \eta \text{ a.s.}$$

For q = 1, let us compare (11) with the condition χ ( A ) &lt; 0, which reduces to

$$\mathbb { E } \left [ \log \left ( a _ { 1, 1 } ( Z _ { 0 } ) + a _ { 1, 2 } ( Z _ { 0 } ) | \varepsilon _ { 1 } | \right ) \right ] < 0.$$

This latter condition is much weaker than (11). Indeed, (11) entails that

$$\mathbb { E } \left [ a _ { 1, 1 } ( Z _ { 0 } ) + a _ { 1, 2 } ( Z _ { 0 } ) | \varepsilon _ { 1 } | \right ] \leq \| a _ { 1, 1 } ( Z _ { 0 } ) + a _ { 1, 2 } ( Z _ { 0 } ) | \varepsilon _ { 1 } | \| _ { p } < 1$$

and from Jensen's inequality, log E [ a 1 1 , ( Z 0 ) + a 1 2 , ( Z 0 ) | ε 1 | ] ≤ χ ( A ).

For q ≥ 2, it is more difficult to obtain explicit conditions which guaranty that χ ( A ) &lt; 0.

- 3. Using the results of Liu and Wu (2010b), a nonparametric kernel estimation of the functions f and g is possible. Proposition 3 gives precise assumptions under which it is possible to control the functional dependence measure of some CHARN models when the regressors include lag values of the response as well as exogenous covariates. We then obtain additional examples of time series models for which standard nonparametric estimators of the regression function are still consistent.

## 4.2 GARCH processes

GARCHprocesses with exogenous regressors have been considered recently by Pedersen and Rahbek (2018) or Francq and Thieu (2019). We consider here the asymmetric power GARCH studied by

Francq and Thieu (2019). The model is defined as follows.

$$Y _ { t } = \varepsilon _ { t } h _ { t } ^ { 1 / \delta }, \ \ h _ { t } = \pi ^ { \prime } Z _ { t - 1 } + \sum _ { i = 1 } ^ { q } \left \{ \beta _ { i } h _ { t - i } + \alpha _ { i + } ( Y _ { t - i } ^ { + } ) ^ { \delta } + \alpha _ { i - } ( Y _ { t - i } ^ { - } ) ^ { \delta } \right \},$$

where ( ε t ) t ∈ Z and ( Z t ) t ∈ Z are two sequences of random variables taking values in R and R d + respectively, δ &gt; 0, π ∈ R d + and the β ′ i s, α ′ i + s and α ′ i -s are nonnegative real numbers. Optimal stationarity properties of time series models defined by (12) have been obtained by Francq and Thieu (2019), using a version of Theorem 1 for affine random maps. In contrast, we use our results to get existence of a moment of order δ for the unique stationary solution. The following assumptions will be needed.

G1 The process (( Z , ε t t )) t ∈ Z is stationary and ergodic and E | Z 0 | &lt; ∞ .

G2 There exist s -, s + such that E [ ( ε + t ) δ |F t -1 ] ≤ s + and E [ ( ε -t ) δ |F t -1 ] ≤ s -a.s. and γ := ∑ p i =1 ( β i + s + α i + + s -α i -) &lt; 1.

## Proposition 4. Suppose that Assumptions G1-G2 hold.

- 1. There then exists a unique stationary and non-anticipative solution ( Y t ) t ∈ Z of (12). This solution is ergodic and satisfies E | Y 0 | δ &lt; ∞ .
- 2. Additionally, assume that Assumption A4 holds true. Let H t = ( ( Y + t ) δ , ( Y + t ) δ , h t ) . There then exists C &gt; 0 and ρ ∈ (0 1) , , such that for all h ≥ 1 , Θ 1 ,h ( H ) ≤ C ρ [ h + ∑ h i =1 ρ i Θ 1 ,h -i ( Z ) ] . Moreover, if δ ≥ 1 , we have the bound

$$\theta _ { \delta, t } ( Y ) \leq \theta _ { 1, t } ^ { 1 / \delta } \left ( Y ^ { + } \right ) + \theta _ { 1, t } ^ { 1 / \delta } \left ( Y ^ { - } \right ), \ \ t \in \mathbb { N }.$$

Note. Let us consider the example of a GARCH process. We then set δ = 2, α j + = α j -= α j and we assume that ( ε t ) t ∈ Z is a martingale difference, adapted to the filtration ( F t ) t ∈ Z with F t = σ (( ε , Z s s ) : s ≤ t ). Set v t -1 = E [ ε 2 t |F t -1 ] . If there exists a positive real number v + such that v t -1 ≤ v + a.s., the contraction condition in G2 reduces to v + ∑ q j =1 ( α j + β j ) &lt; 1. For GARCH

processes with i.i.d. innovations ε t , we recover a standard condition ensuring the existence of a solution with a finite second moment.

## 4.3 Poisson autoregressions

We consider the PARX model introduced in Agosto et al. (2016). The idea is to model the conditional distribution of Y t given F t -1 by a Poisson distribution with a random intensity λ t depending on past values and a covariate process. More precisely, we assume that

$$Y _ { t } = N _ { \lambda _ { t } } ^ { ( t ) }, \quad \lambda _ { t } = \beta _ { 0 } + \sum _ { j = 1 } ^ { q } \beta _ { j } \lambda _ { t - j } + \sum _ { j = 1 } ^ { q } \alpha _ { j } Y _ { t - j } + \pi ^ { \prime } Z _ { t - 1 },$$

where ( N ( ) t ) t ∈ Z is a sequence of i.i.d. Poisson processes with intensity 1, β , . . . , β 0 q , α , . . . , α 1 q are nonnegative real numbers and π is a vector of R d with nonnegative coordinates.

$$\text{PA1} & \text{ We have } \gamma \coloneqq \sum _ { j = 1 } ^ { q } \alpha _ { j } + \sum _ { j = 1 } ^ { q } \beta _ { j } < 1. \\ \text{PA2} & \text{ The process } ( ( Z _ { * } N ^ { ( t ) } ) ) \quad _ { - } \text{ is stationary.  ergo}$$

PA2 The process ( ( Z , N t ( ) t ) ) t ∈ Z is stationary, ergodic and adapted to a filtration ( F t ) t ∈ Z such that for all t ∈ Z , N ( ) t is independent from F t -1 . Moreover, E | Z 1 | &lt; ∞ .

Proposition 5. 1. Suppose that Assumptions PA1-PA2 hold. There then exists a unique nonanticipative, stationary and ergodic process ( Y t ) t ∈ Z solution of (13).

- 2. Additionally, if Assumption A4 is also satisfied with ε t = N ( ) t , there then exists C &gt; 0 and ρ ∈ (0 1) , , such that for all h ≥ 1 , Θ 1 ,h (( Y , λ t t ) ) t ≤ C [ ρ h + ∑ h i =1 ρ i Θ 1 ,h -i ( Z ) ] .

Note. Our result extends substantially that of Agosto et al. (2016). First, we prove ergodicity properties in PARX models without assuming that the covariate process ( Z t ) t ∈ Z is a Markov chain defined by a random map contracting in average. Secondly, for the stochastic dependence properties, we control the coefficient of functional dependence measure only assuming a general Bernoulli shift representation for ( Z t ) t ∈ Z . For instance, ( Z t ) t ∈ Z can be defined by an infinite moving average process and is not necessarily Markovian.

## 4.4 Dynamic binary choice model

We consider the dynamic

$$Y _ { t } = \mathbb { 1 } _ { g ( Y _ { t - 1 }, \dots, Y _ { t - q }, \zeta _ { t } ) > 0 },$$

where ( ζ t ) t ∈ Z is a stationary process taking values in a measurable space E ′ and g : { 0 1 , } q × E ′ → R is a measurable function. This kind of binary model is popular in econometrics for studying the dynamics of recessions. See de Jong and Woutersen (2011) who studied the case g linear and Kauppi and Saikkonen (2008) for a study of US recessions.

Proposition 6. 1. Assume that ( ζ t ) t ∈ Z is a stationary and ergodic process such that

$$\mathbb { P } \left ( \min _ { y \in \{ 0, 1 \} ^ { q }, 1 \leq t \leq q } g ( y, \zeta _ { t } ) > 0 \right ) + \mathbb { P } \left ( \max _ { y \in \{ 0, 1 \} ^ { q }, 1 \leq t \leq q } g ( y, \zeta _ { t } ) \leq 0 \right ) > 0.$$

There then exists a unique stationary and ergodic solution ( Y t ) t ∈ Z for the recursions (14).

- 2. Assume that for some real numbers a , . . . , a 1 q and π ∈ R e , g y, ζ ( t ) = ∑ q i =1 a y i i + π Z ′ t -1 + ε t , with ζ t = ( Z t -1 , ε t ) satisfying A4 and the c.d.f. F ε of ε t being Lipschitz and taking values in (0 1) , . Moreover, setting υ t = π Z ′ t -1 + ε t , we assume that there exists δ &gt; 0 and a positive integer K such that

$$\mathbb { P } \left ( \phi _ { - } + \min _ { 1 \leq t \leq q } v _ { t } > 0 | \mathcal { F } _ { - K } \right ) + \mathbb { P } \left ( \phi _ { + } + \max _ { 1 \leq t \leq q } v _ { t } \leq 0 | \mathcal { F } _ { - K } \right ) \geq \delta \ a. s.,$$

where

$$\phi _ { + } = \max \left \{ \sum _ { i = 1 } ^ { q } a _ { i } y _ { i } \colon ( y _ { 1 }, \dots, y _ { n } ) \in \{ 0, 1 \} ^ { n } \right \}, \quad \phi _ { - } = \min \left \{ \sum _ { i = 1 } ^ { q } a _ { i } y _ { i } \colon ( y _ { 1 }, \dots, y _ { n } ) \in \{ 0, 1 \} ^ { n } \right \}.$$

There then exists C &gt; 0 and ρ ∈ (0 1) , , such that for all h ≥ 1 ,

$$\Theta _ { 1, h } ( Y ) \leq C \left [ \rho ^ { h } + \sum _ { i = 1 } ^ { h } \rho ^ { i } \Theta _ { 1, h - i } ( Z ) \right ].$$

## Notes

- 1. Consider the case of g linear as in the second point of Proposition 6. In this case, de Jong and Woutersen (2011) derived existence of a unique stationary and ergodic solution for when Condition (16) holds true. As shown in de Jong and Woutersen (2011), Condition (16) holds in particular when the process ( υ t ) t ∈ Z is m -dependent or for some infinite moving averages. Condition (15) is much weaker since it holds as soon as

$$\mathbb { P } \left ( \phi _ { + } + \max _ { 1 \leq t \leq q } v _ { t } \leq 0 \right ) + \mathbb { P } \left ( \phi _ { - } + \min _ { 1 \leq t \leq q } v _ { t } > 0 \right ) > 0.$$

Condition (17) holds true as soon as the random vector ( υ , . . . , υ 1 q ) has full support. Another sufficient condition for (17) is the following. If ( υ t ) t ∈ Z is adapted to a filtration ( F t ) t ∈ Z , we assume that for any x ∈ R and t ∈ Z , P ( υ t ≤ x |F t -1 ) &gt; 0 a.s. or for any t ∈ Z and x ∈ R P , ( υ t &gt; x |F t -1 ) &gt; 0 a.s. Recall that υ t = π Z ′ t -1 + ε t . The latter condition is valid in particular when ε t has full support and is independent from F t -1 = σ (( ε , Z s s ) : s ≤ t -1).

- 2. As the proof of Proposition 6 will show, the condition (16) implies Assumption A2 . Condition (15) is only used for applying Theorem 1. However, (15) does not entail mixing properties. In contrast, Condition (16) does. See de Jong and Woutersen (2011), Theorem 2. Our results (see point 2 . of Proposition 6) give a complement when the covariate process is not necessarily strongly mixing and has a Bernoulli shift representation.
- 3. When ζ t = ( Z t -1 , ε t ) ∈ R d +1 in (14), one can allow interactions between lag values of the response and the covariates. For example,

$$g ( y, \zeta _ { t } ) = \sum _ { i = 1 } ^ { d } c _ { i } y _ { i } + \sum _ { i = 1 } ^ { q } \sum _ { j = 1 } ^ { d } \left [ a _ { i, j } y _ { i } + b _ { i, j } ( 1 - y _ { i } ) \right ] Z _ { j, t - i } + \varepsilon _ { t }.$$

When ε t is independent of F t -1 = σ (( ε t -j , Z t -j ) : j ≥ 1), one can show, using the same arguments as in the previous point, that condition (15) is satisfied as soon as the distribution of ε t has support equal to the whole real line. We will not give a control of the functional dependence measure for this model because we were not able to check A2 when the covariate

process ( Z t ) t ∈ Z is not bounded. However when the cdf of ε t is known (e.g. for the logistic or the probit model), it is widely known that ergodicity of the process is sufficient for showing consistency and asymptotic normality of conditional pseudo likelihood estimators of the parameters.

## 4.5 Categorical time series with covariates

We consider a finite set E = { 1 2 , , . . . , N } , an integer q ≥ 1, a process ( Z t ) t ∈ Z taking values in Z ⊂ R d and a family { K z ( ·|· ) : z ∈ Z} of probability kernels from E q to E . Our aim is to construct a process ( Y t ) t ∈ Z , taking values in E and such that

$$\mathbb { P } \left ( Y _ { t } = i | Y _ { t - 1 } ^ { - }, Z _ { t - 1 } ^ { - } \right ) = K _ { Z _ { t - 1 } } \left ( i | Y _ { t - 1 }, \dots, Y _ { t - q } \right ).$$

A particular example is given by the multinomial autoregression, i.e.

$$K _ { z } ( i, y _ { 1 }, \dots, y _ { q } ) = \frac { \exp \left ( \sum _ { j = 1 } ^ { q } a _ { i, j } y _ { j } + \gamma _ { i } ^ { \prime } z \right ) } { \sum _ { k = 1 } ^ { N } \exp \left ( \sum _ { j = 1 } ^ { q } a _ { k, j } y _ { j } + \gamma _ { k } ^ { \prime } z \right ) } \\ \intertext { a model for categorical time series. See Fokianos and Kedem (2003 }$$

and is a classical model for categorical time series. See Fokianos and Kedem (2003). In econometrics, Russell and Engle (2005) studied the dynamic of price changes using this kind of model but with a more general observation-driven form such as in GARCH models and that will not fall into our framework.

For applying our results, we now define some random maps. For t ∈ Z , let ε t be a random variable uniformly distributed over [0 , 1]. For u ∈ [0 , 1], z ∈ E ′ 1 , y ∈ E q and u ∈ [0 , 1], we set and

$$K _ { z } ^ { - } ( u | y ) = \inf \left \{ i = 1, \dots, N \colon \sum _ { j = 1 } ^ { i } K _ { z } ( j | y ) \geq u \right \}$$

$$f _ { t } ( y _ { 1 }, \dots, y _ { q } ) = \left ( K _ { Z _ { t - 1 } } ^ { - } \left ( \varepsilon _ { t } | y _ { 1 }, \dots, y _ { q } \right ), y _ { 1 }, \dots, y _ { q - 1 } \right ) ^ { \prime }. \\. \quad. \quad. \quad.. \quad..$$

We introduce the following assumptions.

- C1 The probability kernels K z are lower bounded by a positive constant, i.e. for any z ∈ E ′ , η z ( ) := min ( i,y ) ∈ E q +1 K z ( i y | ) &gt; 0.
- C2 The process (( Z , ε t t )) t ∈ Z is stationary and ergodic. Moreover, for t ∈ Z , ε t is independent from F t -1 = σ (( Z , ε s s ) : s ≤ t -1).
- C3 There exists a constant C &gt; 0 such that for all y , . . . , y 1 q ∈ E ,

$$\sum _ { i = 1 } ^ { N } | K _ { z } ( i | y _ { 1 }, \dots, y _ { q } ) - K _ { \overline { z } } ( i | y _ { 1 }, \dots, y _ { q } ) | \leq C | z - \overline { z } |.$$

Proposition 7. Suppose that Assumptions C1-C2 hold.

- 1. There exists a unique stationary process satisfying the recursions

$$Y _ { t } = K _ { Z _ { t - 1 } } ^ { - } \left ( \varepsilon _ { t } | Y _ { t - 1 }, \dots, Y _ { t - q } \right ), \quad t \in \mathbb { Z }.$$

Moreover, the process (( Y , Z t t )) t ∈ Z is ergodic.

- 2. Additionally, assume that Assumption A4 and C3 hold true and that η -= inf z ∈ E ′ 1 η z ( ) &gt; 0 in C1 . There then exist C &gt; 0 and ρ ∈ (0 1) , , such that for all h ≥ 1 , Θ 1 ,h ( Y ) ≤ C [ ρ h + ∑ h i =1 ρ i Θ 1 ,h -i ( Z ) ] .

Note. A proof of the first point of Proposition 7 is based on Theorem 1 and provides a general result for existence of stationary categorical time series with covariates. In particular, probit, logistic and multinomial autoregressions can be considered without restriction for the covariate process ( Z t ) t ∈ Z . However, the derivation of the dependence properties in the second point imposes a more restrictive assumption on the transition kernel K because it is necessary to check Assumption A2 . Recently, Fokianos and Truquet (2019) studied categorical time series under the strict exogeneity assumption for the covariate process. For the recursions (18), strict exogeneity holds true as soon as the two processes ( Z t ) t ∈ Z and ( ε t ) t ∈ Z are independent. Assumption C2 is weaker than strict exogeneity in general.

## 4.6 Categorical time series and coalescence of the paths

In this section, we give another interpretation of the convergence of the backward iterations for categorical time series. This interpretation has a link with some perfect simulation schemes that are widely known for Markov chains. See Propp and Wilson (1996). Since the state space is discrete, the iterations should be automatically constant after some steps. Figure 1 illustrates the convergence when q = 1 and N = 3. In this case, f t ( j ) = K -Z t -1 ( ε t | j ) for j = 1 2 3. , , Setting

$$T = \inf \left \{ k \geq 1 \, \colon \varepsilon _ { t - k } \leq \min _ { i, j } K _ { Z _ { t - k - 1 } } ( i | j ) \right \},$$

we know that from the ergodicity assumption C2 and the positivity assumption C1 , T is finite almost surely. In this case, f t -T ( j ) = 1 a.s., f t -∞ := lim n →∞ f t t -n ( j ) = f t t -T ( j ) a.s. and the limit does not depend on the state j . All the paths corresponding to f t t -n ( j ) for n ≥ T coalesce through State 1.

Figure 1: Illustration of the convergence for N = 3 modalities and q = 1 lag

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00845v3-figure-46123073-9fbd-4e34-b0f0-5fc9bde542c0.jpg)

When q ≥ 2, our assumptions guarantee that it is possible to get q times successively the value 1 for the time series whatever the previous values. A coalescence property for the paths will then also occur in this case.

This interpretation is also relevant for getting an interpretation of Condition (15) for binary choice models. When (15) is satisfied, it is possible to get, starting at any time t , either q times the value 1 or q -times the value 0 whatever the previous values of the binary time series. Running the backward iterations, we have coalescence of the paths at the first (random) time t -T such that such an event occurs.

## 5 Conclusion

A general theoretical analysis of nonlinear autoregressive time series models with exogenous covariates is absent from the present time series literature and only a few references consider such a problem, mainly for specific examples. The aim of this paper was to provide some results for a reasonable class of nonlinear time series models for which the required assumptions can be checked. In particular, we provide two results, Theorem 2 and Theorem 4, which justify existence and uniqueness of stationary and ergodic solutions possessing some moments. The crucial assumption to check, A2 or B3 , involves a uniform conditional contraction condition. Assumption B3 is the main assumption to check for autoregressive models with several lags. For some nonlinear models already considered in the literature, such as GARCH or autoregressive Poisson processes, this contraction condition is easily checked because the exogenous covariates have an additive contribution in the expression of the latent process and play the role of a random intercept which does not modify the usual stability conditions. However, our results can also be applied to autoregressive processes for which lag parameters depend on the covariates (see Section 2.4, the notes after Theorem 4, Proposition 3 and Proposition 6 for some examples). In this case, a uniform control of the random lag parameters is necessary to check our assumptions which shows the limit of our approach.

It may be possible to weaken our uniform contraction condition, as shown in Theorem 3, at least under a strict exogeneity assumption. However, getting additional general results to ensure existence of some unconditional moments, to control dependence coefficients and to consider higher-order autoregressive processes would require substantial effort. The proposed framework is also useful for deriving weak dependence properties of the solution, leading to the possibility to apply many existing statistical inference procedures, the central limit theorem of Section 5 providing

an illustration. To this end, the functional dependence measure discussed in Section 3 is of primary importance. Note that a general result for getting weak dependence properties of autoregressive processes with exogenous covariates is also new and it is another contribution of this paper. Finally, we also derived results for categorical time series in Sections 4.4-4.5. Apart from the weak dependence properties which can be derived from our general results, we also obtained stationarity conditions with weaker assumptions, applying Theorem 1. Note that whatever the results used in the paper (Theorems 1, 2 or 3), the convergence of the backward iterations of random maps appears to be a central point of view for considering many different types of autoregressive systems with exogenous regressors.

## 6 Proofs of the results

## 6.1 Proof of Theorem 2

We use the convention f t -1 t ( x ) = x for ( x, t ) ∈ E × Z . From Assumption A2 , if ( t, s, s ′ ) ∈ Z 3 are such that s ′ ≤ s ≤ t and x, y ∈ E ,

$$\mathbb { E } \left [ d ^ { p } \left ( f _ { s } ^ { t } ( x ), f _ { s ^ { \prime } } ^ { t } ( y ) \right ) | \mathcal { F } _ { t - 1 } \right ] \leq L ^ { p } d ^ { p } \left ( f _ { s } ^ { t - 1 } ( x ), f _ { s ^ { \prime } } ^ { t - 1 } ( y ) \right ),$$

and then

$$\| d \left ( f _ { s } ^ { t } ( x ), f _ { s ^ { \prime } } ^ { t } ( y ) \right ) \| _ { p } \leq L \| d \left ( f _ { s } ^ { t - 1 } ( x ), f _ { s ^ { \prime } } ^ { t - 1 } ( y ) \right ) \| _ { p }.$$

Applying (19) with y = x 0 , s = s ′ = , we get sup t t ∈ Z ‖ d ( f t ( x , f ) t ( x 0 )) ‖ p &lt; ∞ . Next, using A1 and the triangular inequality, we get

$$\sup _ { t \in \mathbb { Z } } \| d \left ( f _ { t } ( x ), y \right ) \| _ { p } < \infty \text{ for every } ( x, y ) \in E ^ { 2 }.$$

With the same kind of arguments, we get for s ≤ t -m ,

$$\| d \left ( f _ { s } ^ { t } ( x ), f _ { s ^ { \prime } } ^ { t } ( y ) \right ) \| _ { p } & \leq \kappa \| d \left ( f _ { s } ^ { t - m } ( x ), f _ { s ^ { \prime } } ^ { t - m } ( y ) \right ) \| _ { p }.$$

- 1. We denote by [ z ] the integer part of a real number z . To apply recursively the previous bounds, we note that for any integer i , i +1 = r m 1 + r 2 with r 1 = [( i +1) /m ] and r 2 = i +1 -r m 1 . We then get from (19) and (21),

$$\sum _ { i \geq 0 } \| d \left ( f _ { t - i } ^ { t } ( x ), f _ { t - i - 1 } ^ { t } ( x ) \right ) \| _ { p } & \ \leq \ \sum _ { i \geq 0 } \kappa ^ { [ ( i + 1 ) / m ] } L ^ { i + 1 - m [ ( i + 1 ) / m ] } \| d \left ( x, f _ { t - i - 1 } ( x ) \right ) \| _ { p } \\ & \ \leq \ \frac { ( L + 1 ) ^ { m } \kappa ^ { ( 1 - m ) / m } } { 1 - \kappa ^ { 1 } / m } \sup _ { j \in \mathbb { Z } } \| d \left ( x, f _ { j } ( x ) \right ) \| _ { p }.$$

This latter bound entails that the series ∑ i ≥ 0 d ( f t t -i ( x , f ) t t - -i 1 ( x ) ) is almost surely finite. By the Cauchy criterion, there exists a random variable X x t ( ) such that lim i →∞ d ( f t t -i ( x , X ) t ( x ) ) = 0 a.s. Moreover, from the previous bound, (20) and the triangular inequality, we deduce that

$$\sum _ { s, t \in \mathbb { Z }, s \leq t } \| d \left ( y, f _ { s } ^ { t } ( x ) \right ) \| _ { p } < \infty \text{ for every } ( x, y ) \in E ^ { 2 }.$$

Next we note that the convergence also holds in L p , since from Fatou's lemma,

$$\| d \left ( X _ { t } ( x ), f _ { t - s } ^ { t } ( x ) \right ) \| _ { p } & \ \leq \ \liminf _ { j \to \infty } \| d \left ( f _ { t - j } ^ { t } ( x ), f _ { t - s } ^ { t } ( x ) \right ) \| _ { p } \\ & \ \leq \ \sum _ { i \geq s } \| d \left ( f _ { t - i } ^ { t } ( x ), f _ { t - i - 1 } ^ { t } ( x ) \right ) \| _ { p } \\ & \ \leq \ \kappa ^ { s / m } \frac { ( L + 1 ) ^ { m } \kappa ^ { ( 1 - m ) / m } } { 1 - \kappa ^ { 1 / m } } \sup _ { j \in \mathbb { Z } } \| d \left ( f _ { j } ( x ), x \right ) \| _ { p } \to 0 \text{ as } s \to \infty.$$

Finally, using (22) and the triangular inequality, we get the last assertion sup t ∈ Z ‖ d ( X x ,x t ( ) 0 ) ‖ p &lt; ∞ .

- 2. If x = y , we have from the almost sure convergence and Fatou's lemma

/negationslash

$$\| d \left ( X _ { t } ( x ), X _ { t } ( y ) \right ) \| _ { p } & \ \leq \ \liminf _ { s \to \infty } \| d \left ( f _ { t - s } ^ { t } ( x ), f _ { t - s } ^ { t } ( y ) \right ) \| _ { p } \\ & \ \leq \ \liminf _ { s \to \infty } \kappa ^ { [ ( s + 1 ) / m ] } L ^ { s + 1 - [ ( s + 1 ) / m ] m } d ( x, y ) = 0.$$

This shows the second point.

- 3. For the third point, we observe that for any j ≥ 1, there exists a measurable function H ( x ) j : E ′ j +1 → E such that f t t -j ( x ) = H ( x ) j ( ζ , . . . , ζ t t -j ). Since lim j →∞ H ( x ) j exists P ( ζ t -j ) j ≥ 0 a.s., it is then possible to define a measurable function H : ( E ′ ) N → E such that X t = H (( ζ t -j ) j ≥ 0 ) a.s. The process ( X t ) t ∈ Z has a Bernoulli shift structure with dependent entries and is then stationary and ergodic provided that the process ( ζ t ) t ∈ Z satisfies the same properties.
- 4. The last property follows from the following bounds which hold for any j ≥ 1:

$$\| d ( X _ { t }, Y _ { t } ) \| _ { p } \ & = \ \| d \left ( f _ { t - m j + 1 } ^ { t } ( X _ { t - m j } ), f _ { t - m j + 1 } ^ { t } ( Y _ { t - m j } \right ) \| _ { p } \\ & \leq \ \kappa ^ { j } \left [ \sup _ { j \in \mathbb { Z } } \| d ( x _ { 0 }, X _ { j } ) \| _ { p } + \sup _ { j \in \mathbb { Z } } \| d ( x _ { 0 }, Y _ { j } ) \| _ { p } \right ]. \quad \Box$$

## 6.2 Proof of Theorem 3

- 1. From Assumption A0 , f t t -s ( x ) is, conditionally on Z , an iteration of s +1 independent random maps. Using Assumption A2' , we get

$$\sum _ { s \geq 0 } \mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( f _ { t - s } ^ { t } ( x ), f _ { t - s - 1 } ^ { t } ( x ) \right ) | Z \right ] \leq \sum _ { s \geq 0 } \prod _ { i = 1 } ^ { s + 1 } \kappa ( Z _ { t - i } ) b _ { t - s - 2 } ( x ),$$

with

$$b _ { t } ^ { p } ( x ) = \int d ^ { p } \left ( x, F ( x, Z _ { t }, u ) \right ) d \mathbb { P } _ { \varepsilon _ { 0 } } ( u ).$$

From Assumption A2' and the triangular inequality, we have E log + b t ( x ) &lt; ∞ for any x ∈ E . We are going to show that

$$\sum _ { s \geq 0 } \mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( f _ { t - s } ^ { t } ( x ), f _ { t - s - 1 } ^ { t } ( x ) \right ) | Z \right ] < \infty \text{ a.s.}$$

This follows from the assumptions on the logarithmic moments. Indeed, (( κ Z ( t ) , b t ( x ))) t ∈ Z is a stationary process and it is widely known that the stochastic recursions

$$Y _ { t } = \kappa ( Z _ { t - 1 } ) Y _ { t - 1 } + b _ { t - 1 } ( x )$$

have a unique stationary solution given by

$$Y _ { t } = b _ { t - 1 } ( x ) + \sum _ { s \geq 0 } \prod _ { i = 1 } ^ { s + 1 } \kappa ( Z _ { t - i } ) b _ { t - s - 2 } ( x ),$$

the latter series being convergent almost surely. See for instance Brandt (1986), Theorem 1. This shows (23). Using Minkowski's inequality and Fatou's lemma for conditional expectations, see for instance Kallenberg (2006), Chapter 5, we then deduce that

$$\mathbb { E } ^ { 1 / p } \left [ S ( x ) ^ { p } | Z \right ] < \infty \text{ a.s. } S ( x ) \coloneqq \sum _ { s \geq 0 } d \left ( f _ { t - s } ^ { t } ( x ), f _ { t - s - 1 } ^ { t } ( x ) \right ).$$

As a consequence, we have P ( S x ( ) &lt; ∞| Z ) = 1 a.s. and then P ( S x ( ) &lt; ∞ ) = 1. From the Cauchy criterion, we then conclude the existence of a random variable X x t ( ) such that lim s →∞ f t t -s ( x ) = X x t ( ) a.s. Note that from Assumption A1' and (23), we have E [ d p ( f t t -s ( x , x ) 0 ) | Z ] &lt; ∞ a.s. for every positive integer s . The convergence lim s →∞ E [ d p ( f t t -s ( x , X ) t ( x ) ) | Z ] = 0 a.s. and E [ d p ( X x ,x t ( ) 0 ) | Z ] &lt; ∞ a.s. follow as in the proof of Theorem 2, using Fatou's lemma for conditional expectation.

- 2. For a positive integer s , we have from A2 ,

$$\mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( f _ { t - s } ^ { t } ( x ), f _ { t - s } ^ { t } ( y ) \right ) | Z \right ] \leq \prod _ { i = 1 } ^ { s } \kappa ( Z _ { t - i } ) d ( x, y ) \to 0 \text{ a.s.}$$

/negationslash

/negationslash

Letting s → ∞ , we have E [ d p ( X x ,X t ( ) t ( y )) | Z ] = 0 a.s. and then P ( X x t ( ) = X y t ( ) | Z ) = 0 a.s. Taking the expectation, we conclude that P ( X x t ( ) = X y t ( )) = 0.

- 3. From the almost sure convergence of the sequence ( f t t -s ( x ) ) s ≥ 0 , stationarity and ergodicity of the process (( X , Z t t )) t ∈ Z follows exactly as in the proof of point 3 of Theorem 2.
- 4. Let ( Y t ) t ∈ Z be a stochastic process satisfying the proposed conditions. If the process is non-

anticipative, we have from A2' ,

$$\mathbb { E } \left [ d ^ { p } ( X _ { t }, Y _ { t } ) | Z \right ] \ & = \ \mathbb { E } \left [ d ^ { p } \left ( f _ { t - s } ^ { t } ( X _ { t - s - 1 } ), f _ { t - s } ^ { t } ( Y _ { t - s - 1 } ) \right ) | Z \right ] \\ & \leq \ \prod _ { i = 0 } ^ { s } \kappa ^ { p } ( Z _ { t - i - 1 } ) \cdot \mathbb { E } \left [ d ^ { p } \left ( X _ { t - s - 1 }, Y _ { t - s - 1 } \right ) | Z \right ].$$

Note that, from A2' , ∏ s i =0 κ p ( Z t - -i 1 ) = o P (1). Moreover

$$\mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( X _ { t - s - 1 }, Y _ { t - s - 1 } \right ) | Z \right ] \leq \mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( X _ { t - s - 1 }, x _ { 0 } \right ) | Z \right ] + \mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( x _ { 0 }, Y _ { t - s - 1 } \right ) | Z \right ].$$

Note that

$$\mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( x _ { 0 }, Y _ { t - s - 1 } \right ) | Z \right ] = \mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( x _ { 0 }, Y _ { t - s - 1 } \right ) | Z _ { t - s - 1 }, Z _ { t - s - 2 }, \dots \right ]$$

and if the process (( Z , Y t t )) t ∈ Z is stationary, then the process ( V t ) t ∈ Z defined by

$$V _ { t } = \mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( x _ { 0 }, Y _ { t } \right ) | Z _ { t }, Z _ { t - 1 }, \dots \right ]$$

is also stationary and takes finite values from our assumptions. Then V t - -s 1 = O P (1) where s is the index of the sequence. The same property holds if Y t is replaced with X t . As a consequence

$$\mathbb { E } \left [ d ^ { p } \left ( X _ { t - s - 1 }, Y _ { t - s - 1 } \right ) | Z \right ] = O _ { \mathbb { P } } ( 1 ).$$

/negationslash

We then conclude that E [ d p ( X , Y t t ) | Z ] = 0 a.s. Then P [ X t = Y t | Z ] = 0 a.s. and by integration, we get the conclusion.

## 6.3 Proof of Proposition 1

## 1. We use the decomposition

$$\overline { f } _ { - \infty } ^ { t } - f _ { - \infty } ^ { t } \ & = \ \sum _ { i = 0 } ^ { t - 1 } \left [ f _ { t - i } ^ { t } \circ \overline { f } _ { - \infty } ^ { t - i - 1 } - f _ { t - i - 1 } ^ { t } \circ \overline { f } _ { - \infty } ^ { t - i - 2 } \right ] \\ & + \ \overline { f } _ { t } \circ \overline { f } _ { - \infty } ^ { t - 1 } - f _ { t } \circ \overline { f } _ { - \infty } ^ { t - 1 }.$$

From Assumption A2 and Assumption A4 , we have, for i = 0 , . . . , t -2,

$$\| d \left ( f _ { t - i } ^ { t } \circ \overline { f } _ { - \infty } ^ { t - i - 1 }, f _ { t - i - 1 } ^ { t } \circ \overline { f } _ { - \infty } ^ { t - i - 2 } \right ) \| _ { p } & \ \leq \ \kappa ^ { \frac { i + 1 } { m } - 1 } L ^ { m } \| d \left ( \overline { f } _ { t - i - 1 } \circ \overline { f } _ { - \infty } ^ { t - i - 2 }, f _ { t - i - 1 } \circ \overline { f } _ { - \infty } ^ { t - i - 2 } \right ) \| _ { p } \\ & \ \leq \ \kappa ^ { \frac { i + 1 } { m } - 1 } L ^ { m } \| S ( \overline { f } _ { - \infty } ^ { t - i - 2 } ) H _ { t - i - 2 } \| _ { p } \\ & \ \leq \ \kappa ^ { \frac { i + 1 } { m } - 1 } L ^ { m } \| S ( f _ { - \infty } ^ { 0 } ) \| _ { s } \| H _ { t - i - 2 } \| _ { r }.$$

If i = t -1, we have ‖ d ( f t t -i ◦ f t - -i 1 -∞ , f t t - -i 1 ◦ f t - -i 2 -∞ ) ‖ p ≤ ‖ 2 d ( 0 , f 0 -∞ ) ‖ p κ t/m -1 . Using the triangular inequality, we get for t ≥ 2,

$$\theta _ { p, t } \leq \kappa ^ { - 1 } L ^ { m } \| S ( f _ { - \infty } ^ { 0 } ) \| _ { s } \sum _ { i = 0 } ^ { t - 2 } \kappa ^ { ( i + 1 ) / m } \| H _ { t - i - 2 } \| _ { r } + \| S ( f _ { - \infty } ^ { 0 } ) \| _ { s } \| H _ { t - 1 } \| _ { r } + 2 \kappa ^ { t / m - 1 } \| d \left ( 0, f _ { - \infty } ^ { 0 } \right ) \| _ { p }.$$

The bound (7) is obtained by summation and entails the simpler bound (8).

- 2. From A4' , we have A4 with H t -1 = d ( Z t -1 , Z t -1 ) with

$$\overline { Z } _ { t } = H ^ { \prime } \left ( \eta _ { t }, \dots, \eta _ { 1 }, \eta _ { 0 } ^ { \prime }, \eta _ { - 1 }, \dots \right ).$$

We then deduce the result from the previous point, noticing that η r,j = Θ r,j -1 ( Z ).

## 6.4 Proof of Theorem 4

Define the following random map

$$f _ { t } ( u _ { 1 }, \dots, u _ { q } ) = \left ( F \left ( u _ { 1 }, \dots, u _ { q }, \zeta _ { t } \right ) ^ { \prime }, u _ { 1 } ^ { \prime }, \dots, u _ { q - 1 } ^ { \prime } \right ) ^ { \prime }.$$

We set x = ( u , . . . , u 1 q ) ∈ E q and for 1 ≤ t ≤ q , U t ( x ) = u q -t +1 . Next for t ≥ q + 1, we define U t ( x ) recursively by

$$U _ { t } ( x ) = F \left ( U _ { t - 1 } ( x ), \dots, U _ { t - q } ( x ), \varepsilon _ { t } \right ).$$

We then have for t ≥ q +1,

$$( U _ { t } ( x ), \dots, U _ { t - q + 1 } ( x ) ) = f _ { q + 1 } ^ { t } ( x ).$$

Using our assumptions, we have for t ≥ q +1,

$$\mathbb { E } \left [ \Delta _ { v e c } ^ { p } \left ( U _ { t } ( x ), U _ { t } ( x ^ { \prime } ) \right ) | \mathcal { F } _ { t - 1 } \right ] \preceq \sum _ { i = 1 } ^ { q } A _ { i } \Delta _ { v e c } ^ { p } \left ( U _ { t - i } ( x ), U _ { t - i } ( x ^ { \prime } ) \right ).$$

We introduce the matrix

The condition ρ A ( 1 + · · · + A q ) &lt; 1 entails that ρ B ( ) &lt; 1. Indeed, if v = ( v , . . . , v ′ 1 ′ q ) ′ ∈ R kq \ { 0 } is such that Bv = λv for | λ | ≥ 1, we get the equality v 1 = [ λ -1 A 1 + · · · + λ -q A q ] v 1 . Since the coefficients of the A ′ i s are nonnegative, we get

$$B & = \begin{pmatrix} A _ { 1 } & \cdots & A _ { q - 1 } & A _ { q } \\ I _ { k ( q - 1 ) } & & 0 _ { k ( q - 1 ), 1 } \end{pmatrix}. \\ \cdots + A _ { q } ) & < 1 \text{ entails that } \rho ( B ) < 1. \text{ Indeed, if } v = ($$

$$| v _ { 1 } | _ { v e c } \preceq \sum _ { i = 1 } ^ { q } | \lambda | ^ { - i } A _ { i } | v _ { 1 } | _ { v e c } \preceq \sum _ { i = 1 } ^ { q } A _ { i } | v _ { 1 } | _ { v e c },$$

where | v 1 | vec denotes the vector of the absolute values of the coordinates of v 1 . We then get | v 1 | vec /precedesequal ( ∑ q j =1 A j ) k | v 1 | vec for any positive integer k . Letting k → ∞ , we obtain v 1 = 0. Since v i = λv i +1 for i = 1 , . . . , q -1, we get v = 0 which is a contradiction. Then | λ &lt; | 1 and ρ B ( ) &lt; 1.

Next, we set

$$V _ { t } ( x ) = \left ( U _ { t } ( x ) ^ { \prime }, \dots, U _ { t - q + 1 } ( x ) ^ { \prime } \right ), \quad t \geq q + 1. \\ \intertext { e t } \tilde { t } \quad \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \text{$\mathbb{ }$} \end{$$

Note that V t ( x ) = f t q +1 ( x ). We then have

$$\mathbb { E } \left [ \Delta _ { v e c } ^ { p } \left ( V _ { t } ( x ), V _ { t } ( x ^ { \prime } ) \right ) | \mathcal { F } _ { t - 1 } \right ] \preceq B \Delta _ { v e c } ^ { p } \left ( V _ { t - 1 } ( x ), V _ { t - 1 } ( x ^ { \prime } ) \right ) \preceq \dots \preceq B ^ { t - q } \Delta ^ { p } \left ( x, x ^ { \prime } \right ).$$

$$\text{Setting for } v, v ^ { \prime } \in \mathbb { R } ^ { k q }, \, \overline { d } ( v, v ^ { \prime } ) & = \left ( \sum _ { i = 1 } ^ { k q } \Delta ^ { p } ( v _ { i }, v _ { i } ^ { \prime } ) \right ) ^ { 1 / p }, \, \text{we get} \\ & \quad \Gamma _ { m } \,. \,. \quad. \quad 1 \quad. \,. \,. \quad - n$$

$$\mathbb { E } \left [ \overline { d } ^ { p } \left ( f _ { q + 1 } ^ { t } ( x ), f _ { q + 1 } ^ { t } ( x ^ { \prime } ) \right ) | \mathcal { F } _ { q } \right ] \leq \left | 1 ^ { \prime } B ^ { t - q } \right | _ { \infty } \overline { d } ^ { p } ( x, x ^ { \prime } ), \\ \intertext { \in t h a v e t o r o f $ \mathbb { P } k q $ h a r i n e $ all $ i t e $ common $ n o n t e $ a n a l $ t $ a n d $ l $ t $ h $ }$$

where 1 denotes the vector of R kq having all its components equal to 1, and | · | ∞ the infinite norm in R kq . Since ρ B ( ) &lt; 1, if t is large enough, we have ∣ ∣ 1 ′ B t -q ∣ ∣ ∞ &lt; 1. We then conclude that A2 is satisfied with d = d , m = inf { j ≥ 1 : ∣ ∣ 1 ′ B j ∣ ∣ ∞ &lt; 1 } and κ p = | 1 ′ B m | ∞ . Moreover, A1 is a direct consequence of B2 . The result then follows from Theorem 2. /square

## 6.5 Proof of Proposition 2

Defining

$$f _ { t } \left ( x _ { 1 }, \dots, x _ { q } \right ) = \left ( F \left ( x _ { 1 }, \dots, x _ { q }, Z _ { t - 1 }, \varepsilon _ { t } \right ), x _ { 1 }, \dots, x _ { q - 1 } \right ),$$

we set for t ≥ 0, Z t = H η ,.. . , η ( t 1 , η ′ 0 , η -1 , . . . ). Using Assumptions A5-B5 , we have, for t ≥ 1,

$$\mathbb { E } ^ { 1 / p } \left [ d ^ { p } \left ( f _ { t } ( x _ { 1 }, \dots, x _ { q } ), \overline { f } _ { t } ( x _ { 1 }, \dots, x _ { q } ) \right ) | \mathcal { F } _ { t - 1 } \vee \sigma ( \xi _ { 0 } ^ { \prime } ) \right ] \leq S ( x _ { 1 }, \dots, x _ { q } ) d \left ( Z _ { t - 1 }, \overline { Z } _ { t - 1 } \right ).$$

We then check Assumption A4 , setting H t -1 = d ( Z t -1 , Z t -1 ) . One can then apply Proposition 1 using for instance the /lscript 1 -distance on E q defined by d u, v ( ) = ∑ q i =1 d u , v ( i i ) which is equivalent to the distance d used in the proof of Theorem 4. Since E 1 /r [ d r ( Z t -1 , Z t -1 )] = θ r,t -1 ( Z ), the proof of the proposition is now complete.

## 6.6 Proof of Theorem 5

The proof is a consequence of Theorem 3 in Wu (2005). It is simply necessary to prove that

$$\Gamma ( Y ) = \sum _ { t = 0 } ^ { \infty } \| Y _ { t } - \overline { Y } _ { t } \| _ { 2 } < \infty.$$

To this end, for t ≥ 0, let X t be the random variables obtained by replacing ξ 0 by an independent copy ξ ′ 0 in its Bernoulli shift representation. We then have Y t = f ( X ,. . . , X t t -k ) for t ≥ k . From the assumption on the function f , stationarity and H¨ older's inequality, we have, setting

$$\ell _ { 1 } = 2 p / ( p - 2 ),$$

$$\theta _ { 2, t } ( Y ) \leq C \left ( 1 + 2 k \| | X _ { 0 } | ^ { o \ell } \| _ { \ell _ { 1 } } \right ) \cdot \sum _ { i = 0 } ^ { k } \theta _ { p, t - i } ( X ). \\ \intertext { \text{ition coefficient for common} }$$

The required condition easily follows by summation.

## 6.7 Proof of Proposition 3

- 1. Let d be the distance induced by the /lscript 1 -norm on R q , i.e. d x, y ( ) = ∑ q i =1 | x i -y i | . We use the notation | x - | y instead of d x, y ( ). For a square matrix A of size q × q , we denote by ‖ A ‖ the corresponding operator norm of A . We define the sequence of random maps as follows:

$$g _ { t } ( x ) = ( f _ { 1 } ( x _ { 1 }, \dots, x _ { q }, Z _ { t - 1 } ) + \varepsilon _ { t } f _ { 2 } ( x _ { 1 }, \dots, x _ { q }, Z _ { t - 1 } ), x _ { 1 }, \dots, x _ { q - 1 } ) ^ { \prime } \,.$$

We then have

$$| g _ { t } ( x ) - g _ { t } ( y ) | _ { v e c } \preceq \left ( \sum _ { i = 1 } ^ { q } c _ { i, t } | x _ { i } - y _ { i } |, | x _ { 1 } - y _ { 1 } |, \dots, | x _ { q - 1 } - y _ { q - 1 } | \right ) ^ { \prime } = A _ { t } \cdot | x - y | _ { v e c }.$$

Iterating the previous bound, we get for any positive integer t ,

$$\left | g _ { 1 } ^ { t } ( x ) - g _ { 1 } ^ { t } ( y ) \right | _ { v e c } \preceq A _ { t } \cdots A _ { 1 } \cdot | x - y | _ { v e c }. \\ \intertext { a t a t } \lambda ^ { t } \, \cdot \, < \, \| \, \Delta \, \dots \, \Delta \, \| \, \text{Tho result is then } a \, \text{conconc}.$$

We then deduce that c g ( t 1 ) ≤ ‖ A t · · · A 1 ‖ . The result is then a consequence of Theorem 1, using the condition χ ( A ) &lt; 0.

- 2. We check the assumptions of Theorem 4. First note that from Assumption CH1 , the process ζ t = ( Z t -1 , ε t ) is ergodic. This entails B1 . Next we set

$$F ( x _ { 1 }, \dots, x _ { q }, \zeta _ { t } ) = f _ { 1 } ( x _ { 1 }, \dots, x _ { q }, Z _ { t - 1 } ) + \varepsilon _ { t } f _ { 2 } ( x _ { 1 }, \dots, x _ { q }, Z _ { t - 1 } ).$$

Our assumptions guarantee that F ( x, ζ 1 ) ∈ L p and, using CH2 , we deduce that F x,ζ ( 1 ) ∈ L p for any x ∈ R q . This shows B2 . Finally, we check B3 . To this end, for i = 1 , . . . , q , we set δ i = sup z ‖ a i, 1 ( z ) + a i, 2 ( z ) | ε 1 |‖ p .

Using Minkowski's inequality for conditional expectations (see for instance Doob (2012), Chapter XI, Section 3), we have

$$\mathbb { E } ^ { 1 / p } \left [ | F ( x _ { 1 }, \dots, x _ { q }, \zeta _ { t } ) - F ( y _ { 1 }, \dots, y _ { q }, \zeta _ { t } ) | ^ { p } \left | \mathcal { F } _ { t - 1 } \right ] & \ \leq \ \mathbb { E } ^ { 1 / p } \left [ \left ( \sum _ { i = 1 } ^ { q } c _ { i, t } | x _ { i } - y _ { i } | \right ) ^ { p } | \mathcal { F } _ { t - 1 } \right ] \\ & \ \leq \ \sum _ { i = 1 } ^ { q } \mathbb { E } ^ { 1 / p } \left [ c _ { i, t } ^ { p } | \mathcal { F } _ { t - 1 } \right ] \cdot | x _ { i } - y _ { i } | \\ & \ \leq \ \sum _ { i = 1 } ^ { q } \delta _ { i } | x _ { i } - y _ { i } |. \\ \text{Next, using convexitv.} \text{ we set}$$

Next using convexity, we get

$$\mathbb { E } \left [ | F ( x _ { 1 }, \dots, x _ { q }, \zeta _ { t } ) - F ( y _ { 1 }, \dots, y _ { q }, \zeta _ { t } ) | ^ { p } \, | \mathcal { F } _ { t - 1 } \right ] \leq \left ( \sum _ { i = 1 } ^ { q } \delta _ { i } \right ) ^ { p - 1 } \sum _ { i = 1 } ^ { q } \delta _ { i } | x _ { i } - y _ { i } |.$$

B3 is then a consequence of CH3 .

- 3. We apply Proposition 2. From the previous points, it is only required to check B4 which is a consequence of Assumption CH3 .

## 6.8 Proof of Proposition 4

For the first part, we apply Theorem 4. To this end, we set E = R 3 + , F = ( F , F 1 2 , F 3 ), F 2 ( y , . . . , y 1 q , ζ t ) = ( ε + t ) δ y 1 1 , , F 3 ( y , . . . , y 1 q , ζ t ) = ( ε -t ) δ y 1 1 , and

$$F _ { 1 } \left ( y _ { 1 }, \dots, y _ { q }, \zeta _ { t } \right ) = \pi ^ { \prime } Z _ { t - 1 } + \sum _ { j = 1 } ^ { q } \beta _ { j } y _ { 1, j } + \left ( \alpha _ { 1 + } ( \varepsilon _ { t } ^ { + } ) ^ { \delta } + \alpha _ { 1 - } ( \varepsilon _ { t } ^ { - } ) ^ { \delta } \right ) y _ { 1, 1 } + \sum _ { j = 2 } ^ { q } \left ( \alpha _ { j + } y _ { 2, j } + \alpha _ { j - } y _ { 3, j } \right ).$$

We then deduce that Assumption B3 holds true with

$$A _ { 1 } = \begin{pmatrix} \beta _ { 1 } + \alpha _ { 1 + } s _ { + } + \alpha _ { 1 - } s _ { - } & 0 & 0 \\ s _ { + } & 0 & 0 \\ s _ { - } & 0 & 0 \end{pmatrix}, \ \ A _ { j } = \begin{pmatrix} \beta _ { j } & \alpha _ { j + } & \alpha _ { j - } \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \ \ j \geq 2.$$

It is straightforward to show that the matrix Γ := ∑ q j =1 A j has eigenvalues 0 and a ± √ a 2 +4( bd + ce ) 2 with a = ∑ q j =1 β j + α 1+ s + + α 1 - -s , b = ∑ q j =2 α j + , c = ∑ q j =2 α j -, d = s + and e = s -. Condition ρ (Γ) &lt; 1 is equivalent to γ &lt; 1. It is then clear that B1-B3 follow from G1-G2 .

For the second part, it is easily seen that B4 is satisfied for a constant function S , p = 1 and s = ∞ . This gives the bound for Θ 1 ,h ( H ). If δ ≥ 1, the last bound for θ δ,t ( Y ) can be obtained from the inequalities

$$| x - y | ^ { \delta } \leq \left | x ^ { \delta } - y ^ { \delta } \right |, \ \ x, y \geq 0. \\ \intertext { \i on } 5$$

## 6.9 Proof of Proposition 5

- 1. To show the first point, we check the assumptions of Theorem 4. We set ζ t = ( N ( ) t , Z t -1 ) , E = N × R + and the state space E is endowed with the /lscript 1 -norm. We first note that ( Y t ) t ∈ Z is a stationary solution of (13) if and only if X t = ( Y , λ t t ) ′ is the solution of X t = F ( X t -1 , . . . , X t -q , ζ t ) with

$$F \left ( x _ { 1 }, \dots, x _ { q }, \zeta _ { t } \right ) = \left ( N _ { f ( x _ { 1 }, \dots, x _ { q }, Z _ { t - 1 } ) } ^ { ( t ) }, f ( x _ { 1 }, \dots, x _ { q }, Z _ { t - 1 } ) \right ) ^ { \prime }$$

and f ( x , . . . , x 1 q , Z t -1 ) = β 0 + ∑ q j =1 β s j j + ∑ q j =1 α y j j + π Z ′ t -1 , x i = ( y , s i i ), 1 ≤ i ≤ q . For x ∈ ( N × R + ) q ,

$$\mathbb { E } \left [ \left | F \left ( x, \zeta _ { 1 } \right ) \right | \right ] = 2 \left ( \beta _ { 0 } + \sum _ { j = 1 } ^ { q } \beta _ { j } s _ { j } + \sum _ { j = 1 } ^ { q } \alpha _ { j } y _ { j } + \pi ^ { \prime } \mathbb { E } ( Z _ { 1 } ) \right ) < \infty \\ \mathbb { m } / \left | \mathcal { C } \right | \right | \right.$$

since E ( | Z 1 | ) &lt; ∞ .

We then have, for ( x, x ′ ) ∈ (( N × R + ) q ) 2 with x = ( x , . . . , x 1 q ) , x ′ = ( x , . . . , x ′ 1 ′ q ) ∀ j = 1 , . . . , q, x j = ( y , s j j ) , x j = ( y , s ′ j ′ j ),

$$\mathbb { E } \left [ \left | F ( x, \zeta _ { t } ) - F ( x ^ { \prime }, \zeta _ { t } ) \right | _ { v e c } | \mathcal { F } _ { t - 1 } \right ] \preceq \sum _ { j = 1 } ^ { q } \begin{pmatrix} \alpha _ { j } & \beta _ { j } \\ \alpha _ { j } & \beta _ { j } \end{pmatrix} \left | x _ { j } - x _ { j } ^ { \prime } \right | _ { v e c }. \\ \Gamma _ { \cdots } \quad \omega \quad \omega \. \quad 1 \.$$

In the previous bounds, we have used the identity E [ | N ( ) t h t -1 -N ( ) t g t -1 | ∣ ∣ F t -1 ] = | h t -1 -g t -1 |

which is valid for two nonnegative random variables h t -1 , g t -1 measurable with respect to F t -1 . The previous equality follows from the properties of the Poisson process. Letting

$$\Gamma = \sum _ { j = 1 } ^ { q } \begin{pmatrix} \alpha _ { j } & \beta _ { j } \\ \alpha _ { j } & \beta _ { j } \end{pmatrix}, \\. \circ \infty \infty \lambda \sim \Lambda \cos \max \min \mathfrak { D } \lambda$$

the matrix Γ has two eigenvalues: 0 and γ . Assumption PA1 then guarantees that ρ (Γ) &lt; 1. Assumptions B1-B3 of Theorem 4 are satisfied. Hence, according to Theorem 4, there exists a unique stationary and non-anticipative process ( X t ) t ∈ Z solution of (13) and such that E [ | X t | ] &lt; ∞ . This process is ergodic. This completes the proof of the first point.

- 2. For the second point, we use Proposition 2. To this end, it is only necessary to check B4 for p = r = 1 and s = ∞ . This is straightforward since we have the equality

$$\mathbb { E } \left [ | F ( x _ { 1 }, \dots, x _ { q }, z, N ^ { ( t ) } ) - F ( x _ { 1 }, \dots, x _ { q }, \overline { z }, N ^ { ( t ) } ) | \right ] = 2 | \pi ^ { \prime } ( z - \overline { z } ) |.$$

The proof of the second point is now complete.

## 6.10 Proof of Proposition 6

- 1. We apply Theorem 1. To this end, we define the random map from E = 0 1 { , } q to E by

$$f _ { t } ( x ) = \left ( \mathbf 1 _ { g ( x, \zeta _ { t } ) > 0 }, x _ { 1 }, \dots, x _ { q - 1 } \right ) ^ { \prime }.$$

We set

$$\delta _ { t } \coloneqq \max _ { y, y ^ { \prime } \in \{ 0, 1 \} ^ { q } } \left | 1 _ { g ( y, \zeta _ { t } ) > 0 } - 1 _ { g ( y ^ { \prime }, \zeta _ { t } ) > 0 } \right | \leq 1 _ { \max _ { y \in \{ 0, 1 \} ^ { q } } g ( y, \zeta _ { t } ) > 0 } - 1 _ { \min _ { y \in \{ 0, 1 \} ^ { q } } g ( y, \zeta _ { t } ) > 0 }. \\ \intertext { c u t i m a t i n u n d } \Delta \colon = \Delta \colon \sum _ { t \in \omega } \Delta \colon = \Delta \colon \sum _ { t \in \omega \lambda } \Delta \colon = \Delta \colon \sum _ { t \in \omega \lambda } \Delta \colon = \Delta \colon \sum _ { t \in \omega \lambda } \Delta \colon = \Delta \colon \sum _ { t \in \omega \lambda \lambda } \Delta \colon = \Delta \colon \sum _ { t \in \omega \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \alpha \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \xi \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \eta \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \skip \} } }$$

Setting ( y t ( x , . . . , y ) t -q +1 ( x )) ′ = f t 1 ( x ) for t ≥ q , we have

$$y _ { t } ( x ) = 1 \text{ if and only if } g \left ( y _ { t - 1 } ( x ), \dots, y _ { t - q } ( x ), \zeta _ { t } \right ) > 0.$$

We have, setting c t = | y t ( x ) -y t ( x ′ ) , |

$$c _ { t } \leq \delta _ { t } \max _ { 1 \leq j \leq q } c _ { t - j }.$$

Using the fact that δ t ≤ 1, a straightforward induction on i = 0 , . . . , q -1 shows that

$$c _ { t + i } \leq \delta _ { t + i } \max _ { 1 \leq j \leq q } c _ { t - j }.$$

Setting d y, y ( ′ ) = max 1 ≤ ≤ i q | y i -y ′ i | ,

$$d \left ( f _ { 1 } ^ { t + q - 1 } ( x ), f _ { 1 } ^ { t + q - 1 } ( x ^ { \prime } ) \right ) \leq \max _ { 0 \leq i \leq q - 1 } \delta _ { t + i } d \left ( f _ { 1 } ^ { t - 1 } ( x ), f _ { 1 } ^ { t - 1 } ( x ^ { \prime } ) \right ).$$

Setting t = 1, this shows in particular that

$$c \left ( f _ { 1 } ^ { q } \right ) \leq \max _ { 1 \leq i \leq q } \delta _ { i } \leq 1 _ { \max _ { y, i } g ( y, \zeta _ { i } ) > 0 } - 1 _ { \min _ { y, i } g ( y, \zeta _ { i } ) > 0 }.$$

From our assumptions, the last upper bound can vanish with positive probability, and then E [log c f ( q 1 )] = -∞ = χ . Theorem 1 leads to the result.

- 2. The result will follow from Proposition 1. To this end, we check Assumptions A1-A3 . A1 is automatic. We use the metric d on { 0 1 , } q which is bounded. We set h t = max t -q +1 ≤ ≤ i t δ i . Note that in the linear case, we have

$$h _ { t } \leq 1 _ { \phi _ { + } + \max _ { t - q + 1 \leq i \leq t } \upsilon _ { i } > 0 } - 1 _ { \phi _ { - } + \min _ { t - q + 1 \leq i \leq t } \upsilon _ { i } > 0 }.$$

To check A2 , we use the bound (24) and the inequality h s ≤ 1 for all s ∈ Z to get

$$c \left ( f _ { t - J q + 1 } ^ { t } \right ) \leq \prod _ { j = 0 } ^ { J - 1 } h _ { t - j q } \leq h _ { t }.$$

Moreover using (16), we have

$$\mathbb { E } \left ( h _ { t } | \mathcal { F } _ { t - J q } \right ) & \ \leq \ 1 - \mathbb { P } \left ( \phi _ { + } + \max _ { t - q + 1 \leq i \leq t } v _ { i } \leq 0 | \mathcal { F } _ { t - J q } \right ) - \mathbb { P } \left ( \phi _ { - } + \min _ { t - q + 1 \leq i \leq t } v _ { i } > 0 | \mathcal { F } _ { t - J q } \right ) \\ & \ \leq \ 1 - \delta,$$

provided that Jq ≥ q + K . This guarantees A2 , with κ = 1 -δ and m = Jq .

Finally, let us check A3 . Setting Z t = H η ,.. . , η ( t 1 , η ′ 0 , η -1 , . . . ). We have

$$\mathbb { E } \left [ d \left ( f _ { t } ( y ), \overline { f } _ { t } ( y ) \right ) | \mathcal { F } _ { t - 1 } \vee \sigma ( \xi _ { 0 } ^ { \prime } ) \right ] \ & \ \leq \ \left | F _ { \varepsilon } \left ( - \sum _ { i = 1 } ^ { q } a _ { i } y _ { i } - \pi ^ { \prime } Z _ { t - 1 } \right ) - F _ { \varepsilon } \left ( - \sum _ { i = 1 } ^ { q } a _ { i } y _ { i } - \pi ^ { \prime } \overline { Z } _ { t - 1 } \right ) \right | \\ & \ \leq \ \ L _ { \varepsilon } \cdot \max _ { j = 1 } ^ { d } | \pi _ { j } | \cdot \sum _ { j = 1 } ^ { d } | Z _ { j, t - 1 } - \overline { Z } _ { j, t - 1 } |, \\ \\ \text{where $L_{\varepsilon}$ denotes the Linschitz constant of $F_{\varepsilon}$.}$$

where L ε denotes the Lipschitz constant of F ε .

## 6.11 Proof of Proposition 7

/negationslash

- 1. For the first point, we will apply Theorem 1 for the discrete metric d x, y ( ) = 1 x = y . Setting for x ∈ E q and t ∈ Z , f t t -q +1 ( x ) = f t ◦ · · · ◦ f t -q +1 ( x ), we have for x, y ∈ E q ,

$$\left \{ f _ { t - q + 1 } ^ { t } ( x ) = f _ { t - q + 1 } ^ { t } ( y ) \right \} & \ \supset \ \left \{ f _ { t - q + 1 } ^ { t } ( x ) = f _ { t - q + 1 } ^ { t } ( y ) = ( 1, \dots, 1 ) \right \} \\ & = \ \left \{ \varepsilon _ { j } \in [ 0, \eta ( Z _ { j - 1 } ) ] \, \colon j = t - q + 1, \dots t \right \}.$$

We then obtain

$$d \left ( f _ { t - q + 1 } ^ { t } ( x ), f _ { t - q + 1 } ^ { t } ( y ) \right ) \leq \left ( 1 - \prod _ { j = t - q + 1 } ^ { t } \mathbf 1 _ { \varepsilon _ { j } \in [ 0, \eta ( Z _ { j - 1 } ) ] } \right ) d ( x, y ). \\ \dots \dots \dots \dots$$

Let us show that

$$p \coloneqq \mathbb { P } \left ( \varepsilon _ { j } \in [ 0, \eta ( Z _ { j - 1 } ) ] \, \colon j = t - q + 1, \dots t \right ) > 0.$$

Note that by stationarity, p does not depend on t . Assume that p = 0. From Assumption C2 and the properties of the conditional expectations, we have

$$p = \mathbb { E } \left [ \eta ( Z _ { t - 1 } ) \prod _ { j = t - q + 1 } ^ { t - 1 } 1 _ { \{ \varepsilon _ { j } \leq \eta ( Z _ { j - 1 } ) \} } \right ] = 0. \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Since η is positive, we get P ( ε j ∈ [0 , η ( Z j -1 )] : j = t -q +1 , . . . t -1) = 0. By finite induction, we deduce that P ( ε t -q +1 ≤ η Z ( t -q )) = 0. Since this latter probability equals to E ( η Z ( t -q )) &gt; 0, we obtain a contradiction. The property (25) is then valid and Theorem 1 applies, which leads to the conclusion.

- 2. For the second point, we will use Theorem 2 and Proposition 1. Since d is a bounded metric, A1 is automatically satisfied. Next, observe that

$$\mathbb { P } \left ( f _ { t - q + 1 } ^ { t } ( x ) = f _ { t - q + 1 } ^ { t } ( y ) | \mathcal { F } _ { t - q } \right ) & \ \geq \ \mathbb { P } \left ( f _ { t - q + 1 } ^ { t } ( x ) = f _ { t - q + 1 } ^ { t } ( y ) = ( 1, \dots, 1 ) | \mathcal { F } _ { t - q } \right ) \\ & \ \geq \ \mathbb { P } \left ( \varepsilon _ { t }, \dots, \varepsilon _ { t - q + 1 } \in [ 0, \eta _ { - } ] \mathcal { F } _ { t - q } \right ) \\ & = \ \mathbb { P } \left ( \varepsilon _ { t }, \dots, \varepsilon _ { t - q + 1 } \in [ 0, \eta _ { - } ] \right ) \\ & \ \geq \ \eta _ { - } ^ { q }.$$

This yields to the bound E [ d ( f t t -q +1 ( x , f ) t t -q +1 ( y ) |F t -q )] ≤ 1 -η q -, which shows the second part of A2 . The first part is automatic.

/negationslash

It remains to check A3' . Note that for i, j ∈ E , we have ( N -1) -1 | i - | ≤ j 1 i = j ≤ | i - | j . Using the /lscript 1 -metric on E q which is equivalent to the discrete metric, we have

$$\mathbb { E } \left [ | K _ { z } ^ { - } ( \varepsilon _ { 1 } | y _ { 1 }, \dots, y _ { q } ) - K _ { z ^ { \prime } } ^ { - } ( \varepsilon _ { 1 } | y _ { 1 }, \dots, y _ { q } ) | \right ] & \ \leq \ \int _ { 0 } ^ { 1 } | K _ { z } ^ { - } ( u | y _ { 1 }, \dots, y _ { q } ) - K _ { z ^ { \prime } } ^ { - } ( u | y _ { 1 }, \dots, y _ { q } ) | \, d u \\ & \ \leq \ \sum _ { j = 1 } ^ { N } \left | \sum _ { i = 1 } ^ { j } K _ { z } ( i | y _ { 1 }, \dots, y _ { q } ) - \sum _ { i = 1 } ^ { j } K _ { z ^ { \prime } } ( i | y _ { 1 }, \dots, y _ { q } ) \right | \\ & \ \leq \ N C \, | z - z ^ { \prime } | \,.$$

Assumption A3' then follows with s = ∞ and r = p = 1. The bound for the functional dependence coefficients is then a direct consequence of Proposition 1. /square

## References

- A. Agosto, G. Cavaliere, D. Kristensen, and A. Rahbek. Modeling corporate defaults: Poisson autoregressions with exogenous covariates (PARX). Journal of Empirical Finance , 38:640-663, 2016.
- D.W.K. Andrews. Non-strong mixing autoregressive processes. Journal of Applied Probability , 21 (4):930-934, 1984.
- A.A. Borovkov. Ergodicity and stability of stochastic processes . J. Wiley, 1998.
- P. Bougerol and N. Picard. Strict stationarity of generalized autoregressive processes. The Annals of Probability , pages 1714-1730, 1992.
- A. Brandt. The stochastic equation Y n +1 = A Y n n + B n with stationary coefficients. Advances in Applied Probability , 18(1):211-220, 1986.
- G. Chamberlain. The general equivalence of Granger and Sims causality. Econometrica: Journal of the Econometric Society , pages 569-581, 1982.
- R. A. Davis and H. Liu. Theory and inference for a class of nonlinear models with application to time series of counts. Statistica Sinica , pages 1673-1707, 2016.
- R.M. de Jong and T. Woutersen. Dynamic time series binary choice. Econometric Theory , 27(4): 673-702, 2011.
- P. Diaconis and D. Freedman. Iterated random functions. SIAM review , 41(1):45-76, 1999.
- J. L Doob. Measure theory , volume 143. Springer Science &amp; Business Media, 2012.
- P. Doukhan. Mixing: properties and examples . Number 85 in Lecture Notes in Statistics. SpringerVerlag, New York, 1994.

- J.H. Elton. A multiplicative ergodic theorem for lipschitz maps. Stochastic Processes and their Applications , 34(1):39-47, 1990.
- R.F. Engle, D.F. Hendry, and J-F. Richard. Exogeneity. Econometrica: Journal of the Econometric Society , pages 277-304, 1983.
- K. Fokianos and B. Kedem. Regression theory for categorical time series. Statist. Sci. , 18:357-376, 2003. ISSN 0883-4237. doi: 10.1214/ss/1076102425. URL http://dx.doi.org/10.1214/ss/1076102425 .
- K. Fokianos and L. Truquet. On categorical time series models with covariates. Stochastic processes and their applications , 129(9):3446-3462, 2019.
- C. Francq and L.Q. Thieu. QML inference for volatility models with covariates. Econometric Theory , 35(1):37-72, 2019.
- W. H¨ ardle and A. Tsybakov. Local polynomial estimators of the volatility function in nonparametric autoregression. Journal of Econometrics , 81(1):223-242, 1997.
- W. H¨ ardle, A. Tsybakov, and L. Yang. Nonparametric vector autoregression. Journal of Statistical Planning and Inference , 68(2):221-245, 1998.
- M. Iosifescu. Iterated function systems: A critical survey . Univ., 2009.
- O. Kallenberg. Foundations of modern probability . Springer Science &amp; Business Media, 2006.
- H. Kauppi and P. Saikkonen. Predicting US recessions with dynamic binary response models. The Review of Economics and Statistics , 90(4):777-791, 2008.
- G. Letac. A contraction principle for certain Markov chains and its applications. Contemp. Math , 50:263-273, 1986.
- W. Liu and W.B. Wu. Asymptotics of spectral density estimates. Econometric Theory , 26(4): 1218-1245, 2010a.

- W. Liu and W.B. Wu. Simultaneous nonparametric inference of time series. The Annals of Statistics , 38(4):2388-2421, 2010b.
- H. L¨tkepohl. u New introduction to multiple time series analysis . Springer Science &amp; Business Media, 2005.
- R.S. Pedersen and A. Rahbek. Testing GARCH-X type models. Econometric Theory , pages 1-36, 2018.
- J.G. Propp and D.B. Wilson. Exact sampling with coupled Markov chains and applications to Statistical Mechanics. Random Structures &amp; Algorithms , 9(1-2):223-252, 1996.
- J.R. Russell and R.F. Engle. A discrete-state continuous-time model of financial transactions prices and times: The autoregressive conditional multinomial-autoregressive conditional duration model. Journal of Business &amp; Economic Statistics , 23(2):166-180, 2005.
- C.A. Sims. Money, income, and causality. The American Economic Review , 62(4):540-552, 1972.
- O. Stenflo. Markov chains in random environments and random iterated function systems. Transactions of the American Mathematical Society , 353(9):3547-3562, 2001.
- O. Stenflo. A survey of average contractive iterated function systems. Journal of Difference Equations and Applications , 18(8):1355-1380, 2012.
- W.B. Wu. Nonlinear system theory: Another look at dependence. Proceedings of the National Academy of Sciences , 102(40):14150-14154, 2005.
- W.B. Wu. Strong invariance principles for dependent random variables. The Annals of Probability , 35(6):2294-2320, 2007.
- W.B. Wu and X. Shao. Limit theorems for iterated random functions. Journal of Applied Probability , 41(2):425-436, 2004.
- W.B. Wu and Y.N. Wu. Performance bounds for parameter estimates of high-dimensional linear models with correlated errors. Electronic Journal of Statistics , 10(1):352-379, 2016.

W.B. Wu, Y. Huang, and Y. Huang. Kernel estimation for time series: An asymptotic theory. Stochastic Processes and their Applications , 120(12):2412-2431, 2010.

- H. Xiao and W.B. Wu. Covariance matrix estimation for stationary time series. The Annals of Statistics , 40(1):466-493, 2012.
Max Zinsou Debaly Lionel Truquet
Iterations of dependent random maps and exogeneity in nonlinear dynamics
We discuss existence and uniqueness of stationary and ergodic nonlinear
autoregressive processes when exogenous regressors are incorporated in the
dynamic. To this end, we consider the convergence of the backward iterations of
dependent random maps. In particular, we give a new result when the classical
condition of contraction on average is replaced with a contraction in
conditional expectation. Under some conditions, we also derive an explicit
control of the functional dependence of Wu (2005) which guarantees a wide range
of statistical applications. Our results are illustrated with CHARN models,
GARCH processes, count time series, binary choice models and categorical time
series for which we provide many extensions of existing results.