*Article* **Cointegration and Error Correction Mechanisms for Singular Stochastic Vectors**

#### **Matteo Barigozzi 1, Marco Lippi 2,\* and Matteo Luciani <sup>3</sup>**


Received: 28 March 2018; Accepted: 7 January 2020; Published: 4 February 2020

**Abstract:** Large-dimensional dynamic factor models and dynamic stochastic general equilibrium models, both widely used in empirical macroeconomics, deal with singular stochastic vectors, i.e., vectors of dimension *r* which are driven by a *q*-dimensional white noise, with *q* < *r*. The present paper studies cointegration and error correction representations for an *I*(1) singular stochastic vector **y***t*. It is easily seen that **y***<sup>t</sup>* is necessarily cointegrated with cointegrating rank *c* ≥ *r* − *q*. Our contributions are: (i) we generalize Johansen's proof of the Granger representation theorem to *I*(1) singular vectors under the assumption that **y***<sup>t</sup>* has rational spectral density; (ii) using recent results on singular vectors by Anderson and Deistler, we prove that for generic values of the parameters the autoregressive representation of **y***t* has a finite-degree polynomial. The relationship between the cointegration of the factors and the cointegration of the observable variables in a large-dimensional factor model is also discussed.

**Keywords:** singular stochastic vectors; cointegration for singular vectors; Granger representation theorem; large-dimensional dynamic factor models)

**JEL Classification:** C0; C01; E0

#### **1. Introduction**

An *r*-dimensional stochastic vector **y***<sup>t</sup>* such that **y***<sup>t</sup>* = **A**0**u***<sup>t</sup>* + **A**1**u***t*−<sup>1</sup> + ··· , where the matrices **A***<sup>j</sup>* are *r* × *q* and **u***<sup>t</sup>* is a *q*-dimensional white noise, with *q* < *r*, is said to be singular. Singular stochastic vectors have been systematically analyzed in a number of papers starting with (Anderson and Deistler 2008a, 2008b). A motivation for studying the consequences of singularity, as argued by these authors, is that the factors' vector in large-dimensional dynamic factor models (DFM), such as those introduced in Forni et al. (2000); Forni and Lippi (2001), (Stock and Watson 2002a, 2002b), is typically singular. Singularity is also an important feature of dynamic stochastic general equilibrium models (DSGE), see e.g., Sargent (1989), Canova (2007), pp. 230–2. Singularity as it arises in DFMs is presented in some detail below.

DFMs are based on the idea that all the observed variables in an economic system are driven by a few common (macroeconomic) shocks and by idiosyncratic components which may result from measurement errors and sectoral or regional shocks. Formally, each variable in the *n*-dimensional dataset *xit*, *i* = 1, 2, ... , *n*, *t* = 1, 2, ... , *T*, is decomposed into the sum of a common component *χit*, and an idiosyncratic component *it*: *xit* = *χit* + *it*, where *χit* and *js* are orthogonal for all *i*, *j*, *t*,*s*. In the standard version of the DFM the common components are linear combinations of an *r*-dimensional vector of common factors **F***<sup>t</sup>* = (*F*1*<sup>t</sup> F*2*<sup>t</sup>* ··· *Frt*) ,

$$
\lambda \chi\_{it} = \lambda\_{i1} F\_{1t} + \lambda\_{i2} F\_{2t} + \dots + \lambda\_{ir} F\_{rt} = \lambda\_i \mathbf{F}\_t. \tag{1}
$$

Now suppose that the observable variables *xit* and the common factors **F***<sup>t</sup>* are *I*(1) and that

$$(1 - L)\mathbf{F}\_l = \mathbf{C}(L)\mathbf{u}\_{l\prime} \tag{2}$$

where **u***<sup>t</sup>* is a nonsingular *q*-dimensional white-noise vector1, the common shocks. A number of papers analyzing macroeconomic databases find strong empirical support for the assumption that the vector **F***<sup>t</sup>* is singular, i.e., that *q* < *r*. See, for US datasets, Giannone et al. (2005); Amengual and Watson (2007); Forni and Gambetti (2010), Luciani (2015). For a Euro-area dataset, see Barigozzi et al. (2014).

Such results can be easily understood observing that usually the static Equation (1) is just a convenient representation derived from a "primitive" set of dynamic equations linking the common components *χit* to the common shocks **u***t*. As a simple example, suppose that the variables *xit* are driven by a common one-dimensional cyclical process *ft*, such that (1 − *αL*)*ft* = *ut*, where *ut* is scalar white noise, and that the variables *xit* load *ft* dynamically:

$$x\_{it} = a\_{i0}f\_t + a\_{i1}f\_{t-1} + \epsilon\_{it}.\tag{3}$$

In this case we can set *F*1*<sup>t</sup>* = *ft*, *F*2*<sup>t</sup>* = *ft*−<sup>1</sup> = *F*1,*t*−1, *λi*<sup>1</sup> = *ai*0, *λi*<sup>2</sup> = *ai*1, so that Equations (1) and (2) take the form

$$\mathbf{x}\_{\mathrm{it}} = \lambda\_{i1} F\_{1t} + \lambda\_{i2} F\_{2t} + \varepsilon\_{\mathrm{it}} \quad \text{and} \quad \begin{pmatrix} F\_{1t} \\ F\_{2t} \end{pmatrix} = \begin{pmatrix} (1 - \alpha L)^{-1} \\ L(1 - \alpha L)^{-1} \end{pmatrix} u\_{t\prime}.$$

respectively. Here *r* = 2 and *q* = 1 so that **F***<sup>t</sup>* is singular. For a general analysis of the relationship between representation (1) and "deeper" dynamic representations like (3), see e.g., Forni et al. (2009); Stock and Watson (2016).

Now suppose that the factors **F***<sup>t</sup>* have been estimated. Obtaining **u***<sup>t</sup>* and the impulse-response functions of the variables *xit* with respect to **u***<sup>t</sup>* (or structural shocks obtained by a linear transformation of **u***t*) requires the estimation of a VAR for the singular *I*(1) vector **F***t*. On the other hand, the latter is necessarily cointegrated with cointegration rank *c* at least equal to *r* − *q* (the rank of the spectral density of (1 − *L*)**F***<sup>t</sup>* does not exceed *q* at all frequencies and, therefore, at frequency zero).

Singular vectors of factors in an *I*(1) DFM and *I*(1) singular vectors in DSGE models provide strong motivation for studying singular *I*(1) vectors in a general time-series context. The main contributions of the paper are:

(I) A generalization of Johansen's proof of the Granger Representation Theorem (from MA to AR), this is Proposition 2. Consider an *I*(1) singular vector **y***t*, with dimension *r*, rank *q* < *r*, and cointegrating rank *c* ≥ *r* − *q*. Assuming that (1 − *L*)**y***<sup>t</sup>* has an ARMA structure, **S**(*L*)(1 − *L*)**y***<sup>t</sup>* = **B**(*L*)**u***<sup>t</sup>* and that some simple additional conditions hold, **y***<sup>t</sup>* has a representation as a vector error correction mechanism (VECM) with *c* error correction terms:

$$\mathbf{A}(L)\mathbf{y}\_{l} = \mathbf{A}^\*(L)(1-L)\mathbf{y}\_{l} + \mathfrak{a}(\mathcal{J}^\prime \mathbf{y}\_{t-1} - \mathbf{w}) = \mathbf{B}(0)\mathbf{u}\_{l},\tag{4}$$

where *α* and *β* are both *r* × *c* and full rank, *β* **y***<sup>t</sup>* − **w** is *I*(0), **A**(*L*) and **A**∗(*L*) are *r* × *r* rational matrices in *L*. Under the additional assumption that unity is the only zero of **B**(*L*), i.e., if *z* = 1 then **B**(*z*) is full rank, **A**(*L*) and **A**∗(*L*) are finite-degree matrix polynomials.

(II) Assuming that the parameters of **S**(*L*) and **B**(*L*) may vary in an open subset of R*λ*, see Section 3.2 for the definition of *λ*, in Proposition 3 we show that all the assumptions used to obtain (4), and also the assumption that unity is the only possible zero of **B**(*L*), hold for generic values of

<sup>1</sup> Usually orthonormality is assumed. This is convenient but not necessary in the present paper.

the parameters. This implies that the matrices **A**(*L*) and **A**∗(*L*) are generically of finite degree, which is obviously not the case for nonsingular vectors.2

The paper is organized as follows. Section 2 is preliminary. We firstly recall recent results for stationary singular stochastic vectors with rational spectral density, see (Anderson and Deistler 2008a, 2008b). Secondly, we discuss cointegration and the cointegrating rank for *I*(1) singular stochastic vectors.

In Section 3 we prove our main results. We also obtain the permanent-transitory shock representation in the singular case: **y***<sup>t</sup>* is driven by *r* − *c* permanent shocks, i.e., *r* minus the cointegrating rank, the usual result. However, the number of transitory shocks is *c* − (*r* − *q*), not *c* as in the nonsingular case.

Section 3 also contains an exercise carried out with simulated singular *I*(1) vectors. We compare the results obtained by estimating an unrestricted VAR in the levels and a VECM. Though limited to a simple example, the results confirm what has been found for nonsingular vectors, that under cointegration the long-run features of impulse-response functions are better estimated using a VECM rather than an unrestricted VAR in the levels (Phillips 1998).

In Section 4 we analyse cointegration of the observable variables *xit* in a DFM. Our results on cointegration of the singular vector **F***<sup>t</sup>* have the implication that *p*-dimensional subvectors of the *n*-dimensional common-component vector *χt*, with *p* > *r* − *c*, are cointegrated. As a consequence, stationarity of the idiosyncratic components would imply that all *p*-dimensional subvectors of the *n*-dimensional dataset **x***<sup>t</sup>* are cointegrated if *p* > *r* − *c*. For example, if *q* = 3 and *d* = 1, then all 3-dimensional subvectors in the dataset are cointegrated, a kind of regularity that we do not observe in actual large macroeconomic datasets. This suggests that an estimation strategy robust to the assumption that the idiosyncratic components can be *I*(1) has to be preferred (for this aspect we refer to Barigozzi et al. 2019). Section 5 concludes. Some proofs, a discussion of some non-uniqueness problems arising with singularity and details on the simulations are collected in the Appendix.

#### **2. Stationary and I**(**1**) **Singular Vectors**

#### *2.1. Stationary Singular Vectors*

As in this paper we only consider representation issues it is convenient to assume that all stochastic processes are defined for *<sup>t</sup>* <sup>∈</sup> <sup>Z</sup>. Accordingly, the lag operator *<sup>L</sup>* is defined as *<sup>L</sup>***y***<sup>t</sup>* <sup>=</sup> **<sup>y</sup>***t*−<sup>1</sup> for *<sup>t</sup>* <sup>∈</sup> <sup>Z</sup> (Bauer and Wagner (2012) also study *<sup>I</sup>*(1) and cointegrated processes for *<sup>t</sup>* <sup>∈</sup> <sup>Z</sup>).

We start by introducing results on singular vectors with an ARMA structure from (Anderson and Deistler 2008a, 2008b). Some preliminary definitions are needed.

#### **Definition 1. (Zeros and Poles)**

*(A) When considering matrices* **<sup>V</sup>**(*z*) *whose entries are rational functions of <sup>z</sup>* <sup>∈</sup> <sup>C</sup> *we always assume that numerator and denominator of each entry have no common roots. If* **V**(*z*) *is an r* × *q matrix of rational functions, we say that z*∗ *is a pole of* **V**(*z*) *if it is a pole of some entry of* **V**(*z*)*.*

*(B) Suppose that* **<sup>V</sup>**(*z*) *is an <sup>r</sup>* <sup>×</sup> *<sup>q</sup> matrix whose entries are polynomial functions of <sup>z</sup>* <sup>∈</sup> <sup>C</sup>*, with <sup>q</sup>* <sup>≤</sup> *r. We say that <sup>z</sup>*<sup>∗</sup> <sup>∈</sup> <sup>C</sup> *is a zero of* **<sup>V</sup>**(*z*) *if rank*(**V**(*z*∗)) <sup>&</sup>lt; *q, and that* **<sup>V</sup>**(*z*) *is zeroless if it has no zeros, i.e., rank*(**V**(*z*)) = *q for all z* <sup>∈</sup> <sup>C</sup>*.*

<sup>2</sup> To our knowledge, the present paper is the first to study cointegration and error correction representations for *I*(1) singular vectors, the factors of *I*(1) dynamic factor models in particular. An error correction model in the DFM framework is studied in (Banerjee et al. 2014, 2017). However, their focus is on the relationship between the observable variables and the factors. Their error correction term is a linear combination of the variables *xit* and the factors **F***t*, which is stationary if the idiosyncratic components are stationary (so that the *x*'s and the factors are cointegrated). Because of this and other differences their results are not directly comparable to those in the present paper.

With a minor abuse of language, we may speak of zeros and poles of the corresponding matrix **V**(*L*). When a *r* × *r* polynomial matrix **S**(*L*) has all its zeros outside the unit circle we say that **S**(*L*) is stable.

All the stationary vector processes considered have an ARMA structure. Precisely, the *r*-dimensional process **y***<sup>t</sup>* has an ARMA structure with rank *q*, *q* ≤ *r*, if there exist


$$\mathbf{y}\_t = \mathbf{V}(L)\mathbf{u}\_{t\prime} \tag{5}$$

where **V**(*L*) = **S**(*L*)−1**B**(*L*).

Suppose that **y***<sup>t</sup>* has also the representation **y***<sup>t</sup>* = **S˜**(*L*)−<sup>1</sup>**B˜**(*L*)**u˜***t*, where **u˜***<sup>t</sup>* is a *q*˜-dimensional nonsingular white noise. Denoting by **Σ***y*(*θ*) the spectral density of **y***t*,

$$
\Sigma\_{\mathcal{Y}}(\theta) = (2\pi)^{-1} \mathbf{V}(\varepsilon^{-i\theta}) \Sigma\_u \mathbf{V}'(\varepsilon^{i\theta}),
$$

so that the rank of **Σ***y*(*θ*) is *q* for all *θ*, with the exception of a finite subset of [−*π*, *π*]. As the spectral density is independent of the ARMA representation, *q* = *q*˜ and **B˜**(*z*) has rank *q* except for a finite subset of C.

**Remark 1.** *Let us recall that the equation*

$$\mathbf{S}(L)\mathbf{\tilde{\xi}}\_{\ell} = \mathbf{B}(L)\mathbf{u}\_{\ell\prime}$$

*in the unknown vector process ζt, where* **S**(*L*) *is stable, has only one stationary solution, and this is* **y***<sup>t</sup>* = **S**(*L*)−1**B**(*L*)**u***t*. *Thus the ARMA process* **y***<sup>t</sup> can also be defined as the stationary solution of* **S**(*L*)*ζ<sup>t</sup>* = **B**(*L*)**u***t.*

**Definition 2. (Genericity)** *Suppose that a statement <sup>Q</sup> depends on* **<sup>p</sup>** ∈ A*, where* <sup>A</sup> *is an open subset of* <sup>R</sup>*λ. We say that Q holds generically in* A*, or that Q holds for generic values of* **p** ∈ A*, if the subset* N *of* A *where it does not hold is nowhere dense in* A*, i.e., the closure of* N *in* A *has no internal points.*

For example, assuming that **<sup>p</sup>** ∈ A <sup>=</sup> <sup>R</sup>, the statement "The roots of the polynomial *<sup>x</sup>*<sup>2</sup> <sup>+</sup> **<sup>p</sup>***<sup>x</sup>* <sup>+</sup> <sup>1</sup> are distinct" holds generically in A.

**Definition 3. (Rational reduced-rank family of filters)** *Assume that r* > *q and let* G *be a set of ordered couples* (**S**(*L*),**B**(*L*))*, where:*


*We say that* G *is a rational reduced-rank family of filters with parameter set* Π*.*

The notation **Sp**(*L*), **Bp**(*L*), though more rigorous, would be heavy and not really necessary. We use it only in Appendix A.1.

#### **Proposition 1.** *Assume that r* > *q.*

*(I) Suppose that* **V**(*L*) *is an r* × *q matrix polynomial in L. If* **V**(*z*) *is zeroless then* **V**(*L*) *has an r* × *r finite-degree stable left inverse, i.e., there exists a finite-degree polynomial r* × *r matrix* **W**(*L*) *such that:* (a) **W**(0) = **I***r,* (b) det(**W**(*z*)) = 0 *implies* |*z*| > 1*,* (c) **W**(*L*)**V**(*L*) = **V**(0)*. Let* **y***<sup>t</sup> be the stationary solution of* **S**(*L*)*ζ<sup>t</sup>* = **B**(*L*)**u***<sup>t</sup> and suppose that* **B**(*L*) *is zeroless. Then* **y***<sup>t</sup> has a finite vector autoregressive representation (VAR)* **A**(*L*)**y***<sup>t</sup>* = **B**(0)**u***t, where* **A**(*L*) = **N**(*L*)**S**(*L*) *and* **N**(*L*) *is a finite-degree left inverse of* **B**(*L*)*.*

*(II) Assume that* **y***<sup>t</sup> is the stationary solution of* **S**(*L*)*ζ<sup>t</sup>* = **B**(*L*)**u***t, where* (**S**(*L*), **B**(*L*)) *belongs to a rational reduced-rank family of filters with parameter set* Π*. For generic values of the parameters in* Π*,* **B**(*L*) *is zeroless so that* **y***t has a finite VAR representation.*

For statement (I) see Anderson and Deistler (2008a), Theorem 3. Statement (II) is a modified version of their Theorem 2, see for a proof Forni et al. (2009), p. 1327.

#### *2.2. Fundamentalness*

Assume that the *r*-dimensional vector **y***t* has an ARMA structure, rank *q* and the moving average representation (5). If rank(**B**(*z*)) = *<sup>q</sup>* for |*z*| < 1, then **<sup>u</sup>***<sup>t</sup>* belongs to the space spanned by **<sup>y</sup>***t*−*k*, with *k* ≥ 0, and representation (5), as well as **u***t*, is called fundamental (for these definitions and results see e.g., Rozanov (1967), pp. 43–7). Note that if (5) is fundamental rank(**B**(0)) = *q*. Note also that when *q* = *r*, the condition that rank(**B**(*z*)) = *q* for |*z*| < 1 becomes det(**B**(*z*)) = 0 for |*z*| < 1.

**Remark 2.** *Note that in Proposition 1, part (II), we do not assume that* **u***<sup>t</sup> is fundamental for* **y***t. However, Proposition 1, (II), states that for generic values of* **p** ∈ Π *the matrix* **B**(*L*) *is zeroless and therefore* **u***<sup>t</sup> is fundamental for* **y***t.*

#### *2.3. I*(1) *Singular Vectors*

To analyze cointegration and the autoregressive representations of singular non-stationary vectors let us first recall the definitions of *I*(0), *I*(1) and cointegrated vectors. This requires some preliminary definitions and results.

We denote by *L*2(Ω, F, *P*) the space of the square-integrable functions on the probability space (Ω, F, *P*). Let **z***<sup>t</sup>* = (*z*1*<sup>t</sup> z*2*<sup>t</sup>* ··· *zrt*) , *zht* ∈ *L*2(Ω, F, *P*), be an *r*-dimensional stochastic process and consider the difference equation

$$(1 - L)\mathbb{J}\_t = \mathbf{z}\_t,\tag{6}$$

in the unknown *r*-dimensional process *ζ <sup>t</sup>*. A solution of (6) is

$$\tilde{\boldsymbol{\Psi}}\_{t} = \begin{cases} \mathbf{z}\_{1} + \mathbf{z}\_{2} + \cdots + \mathbf{z}\_{t}, \text{for } t > 0 \\ \mathbf{0}, \text{for } t = 0 \\ -(\mathbf{z}\_{0} + \mathbf{z}\_{-1} \cdot \cdots + \mathbf{z}\_{t+1}), \text{for } t < 0, \end{cases}$$

see e.g., Gregoir (1999), p. 439, Franchi and Paruolo (2019). All the solutions of (6) are *ψ<sup>t</sup>* = *ψ***˜***<sup>t</sup>* + *φt*, where *φ<sup>t</sup>* = (*φ*1*<sup>t</sup> φ*2*<sup>t</sup>* ··· *φrt*) , *φht* ∈ *L*2(Ω, F, *P*), is a solution of the homogeneous equation (1 − *L*)*ζ<sup>t</sup>* = **<sup>0</sup>**, so that *<sup>φ</sup><sup>t</sup>* <sup>=</sup> **<sup>K</sup>**, for some *<sup>r</sup>*-dimensional stochastic vector **<sup>K</sup>**, for all *<sup>t</sup>* <sup>∈</sup> <sup>Z</sup>. We say that the process *φ<sup>t</sup>* = **K** is a constant stochastic process. Obviously a constant stochastic process *φ<sup>t</sup>* = **K** is weakly stationary. Its spectral measure has the jump **Σ***<sup>K</sup>* at frequency zero. Thus *φ<sup>t</sup>* has a spectral density (has an absolutely continuous spectral measure) if and only if **Σ***<sup>K</sup>* = **0**, i.e., if and only if *φt*(*ω*) = **k**, where **<sup>k</sup>** <sup>∈</sup> <sup>R</sup>*<sup>r</sup>* , for *ω* almost everywhere in Ω.

#### **Definition 4. (I(0), I(1) and Cointegrated vectors)**

**I(0).** *An r-dimensional ARMA* **y***<sup>t</sup> with spectral density* **Σ***y*(*θ*) *is I*(0) *if* **Σ***y*(0) = **0***.*

**I(1).** *The r-dimensional vector stochastic process* **y***<sup>t</sup> is I*(1) *if it is a solution* (1 − *L*)*ζ<sup>t</sup>* = **z***<sup>t</sup> where* **z***<sup>t</sup> is an r-dimensional I*(0) *process. The rank of* **y***<sup>t</sup> is defined as the rank of* **z***t.*

**Cointegration.** *Assume that the r-dimensional stochastic vector* **y***<sup>t</sup> is I*(1) *and denote by* **Σ**Δ*y*(*θ*) *the spectral density of* (1 − *L*)**y***t. The vector* **y***<sup>t</sup> is cointegrated with cointegrating rank c, with* 0 < *c* < *r, if* rank(ΣΔ*y*(0)) = *r* − *c.*

If *q* is the rank of **y***<sup>t</sup>* and *r* ≥ *q*, then *c* = *r* − *q* + *d*, where *q* > *d* > 0. Thus in the singular case, *r* > *q*, **y***<sup>t</sup>* is necessarily cointegrated with cointegrating rank at least equal to *r* − *q*.

If **y***<sup>t</sup>* is *I*(1) and cointegrated with cointegrating rank *c*, there exist *c* linearly independent *r* × 1 vectors **c***j*, *j* = 1, ... , *c*, such that the spectral density of **c** *j* (1 − *L*)**y***<sup>t</sup>* vanishes at frequency zero. The vectors **c***<sup>j</sup>* are called cointegrating vectors and the set **c***j*, *j* = 1, ... , *c*, a complete set of cointegrating vectors. Of course a complete set of cointegrating vectors **c***j*, *j* = 1, ... , *c*, can be replaced by the set **d***j*, *j* = 1, . . . , *c*, where the vectors **d***<sup>j</sup>* are *c* independent linear combinations of the vectors **c***j*.

**Lemma 1.** (I) *Assume that* **y***<sup>t</sup> has an ARMA structure and has the rational representation* (5)*:* **y***<sup>t</sup>* = **V**(*L*)**u***t. Then* **y***<sup>t</sup> is I*(0) *if and only if* **V**(1) = **0***.*

(II) *Assume* (1 − *L*)**y***<sup>t</sup> has an ARMA structure and has the rational representation*

$$((1 - L)\mathbf{y}\_t = \mathbf{V}(L)\mathbf{u}\_t. \tag{7}$$

*The process* **y***<sup>t</sup> is I*(1) *if and only if* **V**(1) = **0***.*

(III) *If* **y***<sup>t</sup> is I*(1)*, cointegrated and has representation* (7)*, the cointegrating rank of* **y***<sup>t</sup> is c if and only if the rank of* **V**(1) *is r* − *c. Moreover* **c** *is a cointegrating vector for* **y***<sup>t</sup> if and only if* **c V**(1) = **0***.*

(IV) *Assume that* **y***<sup>t</sup> is I*(1)*.* **c** *is a cointegrating vector for* **y***<sup>t</sup> if and only if a scalar stochastic variable w* ∈ *L*2(Ω, F, *P*) *can be determined such that* **c y***<sup>t</sup>* − *w is stationary with an ARMA structure.*

**Proof.** (I) is an immediate consequence of **Σ***y*(0)=(2*π*)−1**V**(1)**Γ***u***V**(1) , where **Γ***<sup>u</sup>* is the nonsingular covariance matrix of **u***t*. (II) and (III) are obtained in the same way from **Σ**Δ*y*(0)=(2*π*)−1**V**(1)**Γ***u***V**(1) . (IV) The process **y***<sup>t</sup>* solves (6) with **z***<sup>t</sup>* = **V**(*L*)**u***t*, so that, defining

$$\mu\_t = \begin{cases} \mathbf{u}\_1 + \mathbf{u}\_2 + \dots + \mathbf{u}\_{l\prime} \text{ for } t > 0 \\ \mathbf{0}, \text{for } t = 0 \\ -(\mathbf{u}\_0 + \mathbf{u}\_{-1} \cdot \dots + \mathbf{u}\_{l+1}), \text{ for } t < 0, \end{cases} \tag{8}$$

we have

$$\mathbf{y}\_t = \mathbf{V}(L)\boldsymbol{\mu}\_t + \mathbf{K} = \left[\mathbf{V}(1) + (1 - L)\frac{\mathbf{V}(L) - \mathbf{V}(1)}{1 - L}\right]\boldsymbol{\mu}\_t + \mathbf{K} = \mathbf{V}(1)\boldsymbol{\mu}\_t + \mathbf{V}^\*(L)\mathbf{u}\_t + \mathbf{K}\_t$$

where (i) the entries of **V**∗(*L*)=(**V**(*L*) − **V**(1))/(1 − *L*) are rational functions of *L* with no poles of modulus less or equal to unity, (ii) **K** is a constant *r*-dimensional stochastic process. We have:

$$\mathbf{c}'\mathbf{y}\_t = \mathbf{c}'\mathbf{V}(1)\boldsymbol{\mu}\_t + \mathbf{c}'\mathbf{V}^\*(L)\mathbf{u}\_t + \mathbf{c}'\mathbf{K}.\tag{9}$$

If **c** is a cointegrating vector of **y***t* we have **c V**(1) = **0**, so that

$$\mathbf{c}'\mathbf{y}\_t = \mathbf{c}'\mathbf{V}^\*(L)\mathbf{u}\_t + \mathbf{c}'\mathbf{K}.$$

Setting *w* = **c K**, the process **c y***<sup>t</sup>* − *w* = **c V**∗(*L*)**u***<sup>t</sup>* has the desired properties. Note that *w* has the equivalent definition *w* = **c y**<sup>0</sup> − **c V**∗(*L*)**u**0. Conversely, suppose that *w* is such that **c y***<sup>t</sup>* − *w* has an ARMA structure. By (9),

$$\mathbf{c}'\mathbf{y}\_t - w = \mathbf{c}'\mathbf{V}(1)\mu\_t + \mathbf{c}'\mathbf{V}^\*(L)\mathbf{u}\_t + \mathbf{c}'\mathbf{K} - w\_{\prime\prime}$$

so that

$$
\sqrt{\mathbb{E}(\mathbf{c}'\mathbf{y}\_t - w)^2} + \sqrt{\mathbb{E}(\mathbf{c}'\mathbf{V}^\*(L)\mathbf{u}\_t)^2} + \sqrt{\mathbb{E}(\mathbf{c}'\mathbf{K} - w)^2} \ge \sqrt{\mathbf{c}'\mathbf{V}(1)\Sigma\_\mu\mathbf{V}'(1)\mathbf{c}}.
$$

The three terms on the left-hand side are finite and independent of *t*. As **Σ***μ<sup>t</sup>* = |*t*|**Σ***<sup>u</sup>* and **Σ***<sup>u</sup>* is positive definite, the right-hand side diverges for |*t*| → ∞ unless **c V**(1) = **0**.

Lemma 1 shows that our definitions of *I*(0) and *I*(1) processes are equivalent to Definitions 3.2, and 3.3 in Johansen (1995), p. 35, with two minor differences: (i) our assumption of rational spectral density, (ii) the time span of the stochastic processes is *t* = 0, 1, ... in Johansen's book, *<sup>t</sup>* <sup>∈</sup> <sup>Z</sup> in the present paper. Also, under the assumption that (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)**y***<sup>t</sup>* has an ARMA structure, our definition of cointegration is equivalent to that in Johansen (1995), p. 37.

#### **3. Representation Theory for Singular I**(**1**) **Vectors**

In Section 3.1 we prove our generalization to singular vectors of the Granger representation theorem (from MA to AR). We closely follow the proof in Johansen (1995), Theorem 4.5, p. 55–57. In Section 3.2 we show that, under a suitable parameterization, the matrix of the autoregressive representation is generically of finite degree.

#### *3.1. The Granger Representation Theorem (MA to AR)*

Suppose that *r* ≥ *q*, *c* > 0 and *r* > *c* ≥ *r* − *q*. Let **B**(*L*) be an *r* × *q* polynomial matrix of degree *s*<sup>1</sup> ≥ 0 and **S**(*L*) an *r* × *r* polynomial matrix of degree *s*<sup>2</sup> ≥ 0 with **S**(0) = **I***r*.

**Assumption 1. S**(*L*) *is stable.*

**Assumption 2.** *If z*<sup>∗</sup> *is a zero of* **B**(*z*) *(i.e.* rank(**B**(*z*∗)) < *q) then either z*<sup>∗</sup> = 1 *or* |*z*∗| > 1*.*

Assumption 2 implies that the rank of **B**(0) is *q*. The next is a stronger version of Assumption 2:

**Assumption 3.** *If z*∗ *is a zero of* **B**(*z*) *then z*∗ = 1*.*

**Assumption 4.** rank(**B**(1)) = *r* − *c.*

Under Assumption 1, let **y***t* be a solution of the equation

$$\mathbf{S}(1-L)\mathbf{J}\_{\mathfrak{t}} = \mathbf{S}(L)^{-1}\mathbf{B}(L)\mathbf{u}\_{\mathfrak{t}}.\tag{10}$$

We have

$$\mathbf{y}\_t = \mathbf{S}(L)^{-1}\mathbf{B}(L)\boldsymbol{\mu}\_l + \mathbf{K}\_\prime \tag{11}$$

where *<sup>μ</sup><sup>t</sup>* is defined in (8) and **<sup>K</sup>** is a constant stochastic process. By Assumption 4, **<sup>S</sup>**(1)−1**B**(1) <sup>=</sup> **<sup>0</sup>**, so that **y***<sup>t</sup>* is *I*(1) with cointegrating rank *c*, see Lemma 1, (II) and (III).

Consider the finite Taylor expansion of **B**(*z*) around *z* = 1:

$$\mathbf{B}(z) = \mathbf{B}(1) - (1 - z)\mathbf{B}'(1) + (1 - z)^2 \mathbf{B}''(1) + \dotsb \dots$$

Assumption 4 implies that

$$\mathbb{B}(1) = \xi \eta',$$

where *ξ* is *r* × (*r* − *c*) of rank *r* − *c*, *η* is *q* × (*r* − *c*) of rank *r* − *c*, see Lancaster and Tismenetsky (1985, p. 97, Proposition 3). The Taylor expansion above can be rewritten as

$$\mathbf{B}(z) = \mathbf{\dot{y}}\mathbf{\dot{y}}^{\prime} + (1 - z)\mathbf{B}^{\*} + (1 - z)^{2}\mathbf{E}(z),\tag{12}$$

where **B**<sup>∗</sup> = −**B** (1) and **E**(*z*) is a polynomial matrix.

Let *ξ*<sup>⊥</sup> be an *r* × *c* matrix whose columns are orthogonal to all columns of *ξ*: (i) the columns of *<sup>ξ</sup>*<sup>⊥</sup> are a complete set of cointegrating vectors for **<sup>B</sup>**(*L*)**u***t*, (ii) the columns of the matrix **<sup>S</sup>** (1)*ξ*<sup>⊥</sup> are a complete set of cointegrating vectors for **y***t*. Regarding (i), using (11) and (12), we have

$$\mathbf{g}\_{\perp}^{\prime}\mathbf{S}(L)\mathbf{y}\_{l} = \mathbf{g}\_{\perp}^{\prime}\mathbf{B}(L)\boldsymbol{\mu}\_{l} + \mathbf{g}\_{\perp}^{\prime}\mathbf{S}(1)\mathbf{K} = (\mathbf{g}\_{\perp}^{\prime}\mathbf{B}^{\*} + (1 - L)\mathbf{g}\_{\perp}^{\prime}\mathbf{E}(L))\mathbf{u}\_{l} + \mathbf{g}\_{\perp}^{\prime}\mathbf{S}(1)\mathbf{K},\tag{13}$$

so that *ξ* <sup>⊥</sup>**S**(*L*)**y***<sup>t</sup>* <sup>−</sup> *<sup>ξ</sup>* ⊥**S**(1)**<sup>K</sup>** has an ARMA structure. Regarding (ii), see the proof of Proposition 2.

$$\text{Assumption 5. rank} \left[ \begin{pmatrix} \mathbf{g}'\_{\perp} \mathbf{B}^\* \\ \mathbf{\eta}' \end{pmatrix} \right] = \text{rank} \left[ \begin{pmatrix} \mathbf{g}'\_{\perp} \mathbf{B}^\* \\ \mathbf{g}'\_{\mathbf{g}} \mathbf{\eta}' \end{pmatrix} \right] = q.c.$$

$$\text{Define } \mathbf{S}^\*(L) = \frac{\mathbf{S}(L) - \mathbf{S}(1)}{1 - L}.$$

**Assumption 6.** *ξ* ⊥(**B**<sup>∗</sup> <sup>−</sup> **<sup>S</sup>**∗(1)**S**(1)−1*ξη* ) = **0***.*

**Remark 3.** *Let* **y***<sup>t</sup> be a solution of* (10) *so that* (1 − *L*)**y***<sup>t</sup> is stationary and* **S**(*L*)[(1 − *L*)**y***t*] = **B**(*L*)**u***t. Assumption 2, and therefore 3, implies that* **u***<sup>t</sup> is fundamental for* (1 − *L*)**y***t, see Section 2.2.*

We are now ready for our main representation result.

**Proposition 2.** *(I) Weak form. Suppose that Assumptions 1, 2, 4, 5 and 6 hold and let* **y***<sup>t</sup> be a solution of the difference Equation* (10)*, so that* **y***<sup>t</sup>* = **S**(*L*)−1**B**(*L*)*μ<sup>t</sup>* + **K***, with μ<sup>t</sup> defined in* (8) *and* **K** *a constant stochastic process. Set β* = **S**(1) *ξ*⊥*. Then a c-dimensional stochastic vector* **w** *can be determined such that* (i) *β* **y***<sup>t</sup>* − **w** *is I*(0)*,* (ii) **y***<sup>t</sup> has the error correction representation*

$$\mathbf{A}(L)\mathbf{y}\_t = \mathbf{A}^\*(L)(1-L)\mathbf{y}\_t + \mathfrak{a}(\mathcal{J}^\prime \mathbf{y}\_{t-1} - \mathbf{w}) = \mathbf{B}(0)\mathbf{u}\_{t\prime} \tag{14}$$

*where* **A**(*L*) *is a rational r* × *r matrix with no poles in or on the unit circle,* **A**(1) = **I***r,* **A**∗(*L*)=(**A**(*L*) − **<sup>A</sup>**(1)*L*)(<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)−1, *<sup>α</sup> is r* <sup>×</sup> *c and full rank, αβ* <sup>=</sup> **<sup>A</sup>**(1)*.*

*(II) Strong form. Under Assumptions 1, 3, 4, 5 and 6, statement (I) holds with an r* × *r stable, finite-degree matrix polynomial* **A**(*L*)*.*

**Proof.** Multiply both sides of (1 − *L*)**S**(*L*)**y***<sup>t</sup>* = **B**(*L*)**u***<sup>t</sup>* by the *r* × *r* invertible matrix **Ξ** = *ξ* ⊥ *ξ* .

We obtain

$$\begin{split} &(1-L)\mathbf{E}\mathbf{S}(L)\mathbf{y}\_{t} = \mathbf{E}\mathbf{B}(L)\mathbf{u}\_{t} \\ &= \left\{ \begin{pmatrix} \mathbf{0}\_{\mathbb{C}\times\eta} \\ \mathbf{\xi}^{\mathsf{T}}\xi\eta \end{pmatrix} + (1-L) \begin{pmatrix} \mathbf{\xi}^{\mathsf{T}}\_{\bot}\mathbf{B}^{\star} \\ \mathbf{\xi}^{\mathsf{T}}\mathbf{B}^{\star} \end{pmatrix} + (1-L)^{2} \begin{pmatrix} \mathbf{\xi}^{\mathsf{T}}\_{\bot}\mathbf{E}(L) \\ \mathbf{\xi}^{\mathsf{T}}\mathbf{E}(L) \end{pmatrix} \right\} \mathbf{u}\_{t} \\ &= \begin{pmatrix} (1-L)\mathbf{I}\_{\mathsf{L}} & \mathbf{0} \\ \mathbf{0} & \mathbf{I}\_{\mathsf{T}-\mathsf{c}} \end{pmatrix} \left\{ \begin{pmatrix} \mathbf{\xi}^{\mathsf{T}}\_{\bot}\mathbf{B}^{\star} \\ \mathbf{\xi}^{\mathsf{T}}\mathbf{E}\mathbf{y}^{\mathsf{T}} \end{pmatrix} + (1-L) \begin{pmatrix} \mathbf{\xi}^{\mathsf{T}}\_{\bot}\mathbf{E}(L) \\ \mathbf{\xi}^{\mathsf{T}}\mathbf{B}^{\star} \end{pmatrix} + (1-L)^{2} \begin{pmatrix} \mathbf{0}\_{\mathsf{C}\times\mathsf{q}} \\ \mathbf{\xi}^{\mathsf{T}}\mathbf{E}(L) \end{pmatrix} \right\} \mathbf{u}\_{t}. \end{split} \tag{15}$$

Taking the first *c* rows in (15),

$$((1-L)\mathfrak{F}'\_{\perp}\mathbf{S}(L)\mathbf{y}\_{\ell} = (1-L)\left(\mathfrak{F}'\_{\perp}\mathbf{B}^\* + (1-L)\mathfrak{F}'\_{\perp}\mathbf{E}(L)\right)\mathbf{u}\_{\ell}.$$

This implies that

$$\mathbf{g}'\_{\perp}\mathbf{S}(L)\mathbf{y}\_{l} = \left(\mathbf{g}'\_{\perp}\mathbf{B}^\* + (1 - L)\mathbf{g}'\_{\perp}\mathbf{E}(L)\right)\mathbf{u}\_{l} + \mathbf{w}\_{l} \tag{16}$$

where **w** is a *c*-dimensional constant stochastic vector. Comparing with (13), **w** = *ξ* ⊥**S**(1)**K**. On the other hand,

$$\begin{split} \mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}(1) \mathbf{y}\_{t} - \mathbf{w} &= (\mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}(L) \mathbf{y}\_{t} - \mathbf{w}) - (\mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}(L) \mathbf{y}\_{t} - \mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}(1) \mathbf{y}\_{t}) \\ &= (\mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}(L) \mathbf{y}\_{t} - \mathbf{w}) - \mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}^{\*}(L) (1 - L) \mathbf{y}\_{t} \\ &= (\mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}(L) \mathbf{y}\_{t} - \mathbf{w}) - \mathbf{\tilde{g}}\_{\perp}^{\prime} \mathbf{S}^{\*}(L) \mathbf{S}(L)^{-1} \mathbf{B}(L) \mathbf{u}\_{t} \\ &= \left\{ \mathbf{\tilde{g}}\_{\perp}^{\prime} (\mathbf{B}^{\*} - \mathbf{S}^{\*}(1) \mathbf{S}(1)^{-1} \mathbf{f} \mathbf{y}^{\prime}) + (1 - L) \mathbf{H}(L) \right\} \mathbf{u}\_{t}, \end{split} \tag{17}$$

where the last equality has been obtained using (16) and H(*L*) is a suitable polynomial matrix. Thus *β* **y***<sup>t</sup>* − **w** = *ξ* <sup>⊥</sup>**S**(1)**y***<sup>t</sup>* <sup>−</sup> **<sup>w</sup>** has an ARMA structure. Moreover, by Assumption 6, *<sup>β</sup>* **y***<sup>t</sup>* − **w** is *I*(0).

Joining (16) with the last *r* − *c* rows of (15),

$$
\begin{pmatrix}
\mathbf{I}\_{\varepsilon} & \mathbf{0} \\
\mathbf{0} & (1 - L)\mathbf{I}\_{r - \varepsilon}
\end{pmatrix} \boldsymbol{\Xi} \mathbf{S}(L) \mathbf{y}\_{t} - \begin{pmatrix}
\mathbf{I}\_{\varepsilon} \\
\mathbf{0}\_{(r - \varepsilon) \times \varepsilon}
\end{pmatrix} \mathbf{w} = \mathbf{M}(L) \mathbf{u}\_{t} \tag{18}
$$

where

$$\mathbf{M}(L) = \begin{pmatrix} \mathfrak{z}'\_{\perp} \mathbf{B}^\* \\ \mathfrak{z}' \mathfrak{z} \mathfrak{y}' \end{pmatrix} + (1 - L) \begin{pmatrix} \mathfrak{z}'\_{\perp} \mathbf{E}(L) \\ \mathfrak{z}' \mathbf{B}^\* \end{pmatrix} + (1 - L)^2 \begin{pmatrix} \mathbf{0}\_{c \times q} \\ \mathfrak{z}' \mathbf{E}(L) \end{pmatrix}. \tag{19}$$

By (15) and (19),

$$\mathbf{B}(L) = \boldsymbol{\Xi}^{-1} \begin{pmatrix} (1 - L)\mathbf{I}\_{\boldsymbol{\varepsilon}} & \mathbf{0} \\ \mathbf{0} & \mathbf{I}\_{r - \boldsymbol{\varepsilon}} \end{pmatrix} \mathbf{M}(L). \tag{20}$$

By Assumption 5, **M**(*z*) has no zero at *z* = 1, see (19). On the other hand, (i) if *z*∗ is a zero of **M**(*z*) then *z*<sup>∗</sup> is a zero of **B**(*z*), (ii) if *z*<sup>∗</sup> is a zero of **B**(*z*), *z*<sup>∗</sup> = 1, then *z*<sup>∗</sup> is a zero of **M**(*z*). Therefore, Assumption 3 implies that **M**(*z*) is zeroless and viceversa. Under Assumption 2, the zeros of **M**(*z*) lie outside the unit circle. In order to conclude the proof we need inverting **M**(*L*) in (18).

(I) Under Assumption 3, Proposition 1, part (I), states that there exists an *r* × *r* stable, finite-degree polynomial matrix **<sup>N</sup>**(*L*) = **<sup>I</sup>***<sup>r</sup>* <sup>+</sup> **<sup>N</sup>**1*<sup>L</sup>* <sup>+</sup> ··· <sup>+</sup> **<sup>N</sup>***pLp*, for some *<sup>p</sup>*, such that: (i) **<sup>N</sup>**(0) = **<sup>I</sup>***r*, (ii) **N**(*L*)**M**(*L*) = **M**(0).

(II) Under Assumption 2, by a standard procedure we remove all the zeros of **M**(*z*) which lie outside the unit circle4, then use Proposition 1, part (I), to left-invert the residual zeroless polynomial, thus obtaining an *r* × *r* rational matrix **N**(*L*) such that (i) **N**(*L*) has no poles in or on the unit circle (possible poles of **N**(*L*) are the zeros of **M**(*L*), which lie outside the unit circle), (ii) **N**(0) = **I***r*, (iii) **N**(*L*)**M**(*L*) = **M**(0). See also Deistler et al. (2010).

Defining

$$\mathbf{A}(L) = \boldsymbol{\Xi}^{-1} \mathbf{N}(L) \begin{pmatrix} \mathbf{I}\_{\varepsilon} & \mathbf{0} \\ \mathbf{0} & (1 - L)\mathbf{I}\_{\mathrm{r} - \varepsilon} \end{pmatrix} \boldsymbol{\Xi} \mathbf{S}(L) = \boldsymbol{\Xi}^{-1} \mathbf{N}(L) \begin{pmatrix} \boldsymbol{\mathfrak{E}}'\_{\perp} \\ (1 - L)\mathbf{\mathfrak{E}}' \end{pmatrix} \mathbf{S}(L)$$

and using **M**(0) = **ΞB**(0), we have

$$\mathbf{A}(L)\mathbf{y}\_t - \boldsymbol{\Xi}^{-1}\mathbf{N}(1)\begin{pmatrix} \mathbf{I}\_{\boldsymbol{\varepsilon}} \\ \mathbf{0}\_{(r-\boldsymbol{\varepsilon})\times\boldsymbol{\varepsilon}} \end{pmatrix}\mathbf{w} = \mathbf{B}(0)\mathbf{u}\_{t\times t}$$

<sup>3</sup> In the square case, *r* = *q*, Assumption 3 holds if and only if **M**(*z*) is unimodular.

<sup>4</sup> If *<sup>z</sup>*<sup>∗</sup> is a zero of **<sup>M</sup>**(*z*), multiply **<sup>M</sup>**(*z*) by an invertible *<sup>r</sup>* <sup>×</sup> *<sup>r</sup>* matrix **<sup>Q</sup>***z*<sup>∗</sup> such that *<sup>z</sup>*<sup>∗</sup> is a zero of, say, the first row of **<sup>Q</sup>***z*∗**M**(*z*). Then multiply by the *<sup>r</sup>* <sup>×</sup> *<sup>r</sup>* diagonal matrix with (*<sup>z</sup>* <sup>−</sup> *<sup>z</sup>*∗)−<sup>1</sup> in position (1, 1) and unity elsewhere on the main diagonal. Iterating, all the zeros of **M**(*z*) are removed.

with **<sup>A</sup>**(0) = **<sup>I</sup>***r*. Defining **<sup>A</sup>**∗(*L*)=(**A**(*L*) <sup>−</sup> **<sup>A</sup>**(1)*L*)(<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)−1,

$$\mathbf{A}^\*(L)(1-L)\mathbf{y}\_t + \mathbf{A}(1)\mathbf{y}\_{t-1} - \boldsymbol{\Xi}^{-1}\mathbf{N}(1)\begin{pmatrix} \mathbf{I}\_{\boldsymbol{\varepsilon}} \\ \mathbf{0}\_{(r-\boldsymbol{\varepsilon})\times\boldsymbol{\varepsilon}} \end{pmatrix}\mathbf{w} = \mathbf{B}(0)\mathbf{u}\_t.$$

Defining

$$\mathfrak{a} = \Xi^{-1} \mathbf{N}(1) \begin{pmatrix} \mathbf{I}\_c \\ \mathbf{0}\_{(r-c)\times c} \end{pmatrix}.$$

we see that **A**(1) = *αβ* and

$$\mathbf{A}^\*(L)(1-L)\mathbf{y}\_t + \mathfrak{a}(\mathcal{J}^\prime \mathbf{y}\_{t-1} - \mathbf{w}) = \mathbf{B}(0)\mathbf{u}\_t.$$

Some remarks are in order.

**Remark 4.** (I) *Under our assumption of an ARMA structure, Assumption 1 corresponds to Definition 3.1 in Johansen's book, see p. 34. Assumption 2 is Johansen's Assumption 1 (see p. 14), adapted for singularity. Assumption 3 has no counterpart in Johansen's nonsingular framework. In Section 3.2 we show that under the parameterization adopted in Definition 5, Assumption 3 holds generically.*

(II) *Simplifying the model by taking* **S**(*L*) = **I***r, Assumption 5 generalizes to the singular case Johansen's assumption that ξ* <sup>⊥</sup>**C**∗*η*<sup>⊥</sup> *is full rank (see Theorem 4.5, p. 55;* **<sup>C</sup>**<sup>∗</sup> *corresponds to our* **<sup>B</sup>**∗*). For, assuming that <sup>r</sup>* <sup>=</sup> *q, multiplying the matrix in Assumption <sup>5</sup> by the nonsingular matrix η*<sup>⊥</sup> *η , we obtain that Assumption 5 holds if and only if ξ* <sup>⊥</sup>**B**∗*η*<sup>⊥</sup> *is full rank. Assumption <sup>5</sup> is used in the proof of Proposition <sup>2</sup> to invert the matrix* **M**(*L*)*, which remains on the right-hand side after the removal of the unit roots, see Equation (18), which is the same rôle played by Johansen's assumption in his proof.*

(III) *Under* **S**(*L*) = **I***r, assumption 6 simplifies to ξ* ⊥**B**<sup>∗</sup> <sup>=</sup> **<sup>0</sup>***. If <sup>d</sup>* <sup>&</sup>gt; <sup>0</sup> *Assumption <sup>6</sup> is a consequence of Assumption 5. For, if d* > 0 *then r* − *c* = *q* − *d* < *q. On the other hand, r* − *c is the number of rows of η , so that Assumption 5 holds only if Assumption 6 holds. In particular, if r* = *q and c* = *d* > 0*, Assumption 6 is redundant. However if r* > *q and d* = 0*, so that the rank of η is q, then Assumption 5 holds even if ξ* ⊥**B**<sup>∗</sup> <sup>=</sup> **<sup>0</sup>***. Assumption 6 is necessary in Proposition 2 to prove that the error correction term is I*(0)*, not only stationary.*

**Remark 5.** *Uniqueness issues arise with autoregressive representations of singular vectors. For example, suppose that c* = *r* − *q* > 0*, so that d* = 0*. Representation* (14) *has an* (*r* − *q*)*-dimensional error correction term β* **y***<sup>t</sup>* − **w***. On the other hand, in this case* **B**(1) *has full rank q, so that Proposition 1* (I) *applies and, in spite of cointegration,* **y***t has an autoregressive representation in differences*

$$\mathbf{D}(L)\mathbf{S}(L)(1-L)\mathbf{y}\_{\mathbf{f}} = \mathbf{B}(0)\mathbf{u}\_{\mathbf{f}}.$$

*In Appendix B.1 we sketch a proof of the statement that in general,* **y***<sup>t</sup> has VECM representations with a number of error correction terms ranging from d to c. However, as we show in Appendix B.2, different autoregressive representations of* **y***t produce the same impulse-response functions. Both in this and the companion paper Barigozzi et al. (2019) the number of error correction terms in the error correction representation for reduced-rank I(1) vectors is always the maximum c. It is worth reporting that, in our experiments with simulated data, the best results in estimation of singular VECMs are obtained using c as the number of error correction terms.*

**Remark 6.** *Assume for simplicity that* **S**(*L*) = **I***r. From equation* (17)*:*

$$\mathbf{e}\_t = \boldsymbol{\mathfrak{z}}'\mathbf{y}\_t - \mathbf{w} = \boldsymbol{\mathfrak{z}}'\_{\perp}\mathbf{y}\_t - \mathbf{w} = \left\{ \boldsymbol{\mathfrak{z}}'\_{\perp}\mathbf{B}^\* + (1 - L)\mathbf{\mathcal{H}}(L) \right\} \mathbf{u}\_t.$$

*If r* = *q, Assumption 5 implies that ξ* ⊥**B**<sup>∗</sup> *has rank c, so that no c-dimensional vector* **<sup>d</sup>** <sup>=</sup> **<sup>0</sup>** *can be determined such that some of the coordinates of* **de***<sup>t</sup> is stationary but not I*(0)*. Thus, according to the definition introduced in Franchi and Paruolo (2019), p. 1181, the error term* **e***<sup>t</sup> is a "non-cointegrated I*(0) *process." When r* > *q and c* ≤ *q, i.e., r* ≤ 2*q* − *d, elementary examples can be produced in which* **e***<sup>t</sup> is an I*(0) *but not a non-cointegrated I*(0) *process (one is given in Appendix A.2). Thus Assumption 6 only implies that* **e***<sup>t</sup> is I*(0)*. Of course, under c* ≤ *q, the assumption that ξ* ⊥(**B**<sup>∗</sup> <sup>−</sup> **<sup>S</sup>**∗(1)**S**(1)−1*ξη* ) *has rank c, an enhancement of Assumption 6, implies that* **e***<sup>t</sup> is a non-cointegrated I*(0) *process. On the other hand, if c* > *q, i.e., r* > 2*q* − *d,* **e***<sup>t</sup> cannot be a non-cointegrated I*(0) *process.*

#### *3.2. Generically,* **A**(*L*) *Is a Finite-Degree Polynomial*

Suppose that the couple (**S**(*L*),**B**(*L*)) is parameterized as in Definition 3. It easy to see that **B**(1) has generically rank *q*, so that generically the cointegrating rank of **y***<sup>t</sup>* is *r* − *q*. In particular, if *r* = *q* cointegration is non generic.

It is quite easy to see that this paradoxical result only depends on the choice of a parameter set that is unfit to study cointegration. Our starting point here is that a specific value of *c* between *r* − *q* and *r* − 1 has a motivation in economic theory or in statistical inference, and must be therefore built in the parameter set. Thus in Definition 5 below the family of filters is redefined so that generically the cointegrating rank is equal to a given *c* between *r* − *q* and *r* − 1.

**Definition 5. (Rational reduced-rank family of filters with cointegrating rank c)** *Assume that r* > *q, c* > 0 *and r* > *c* ≥ *r* − *q. Let* G *be a set of couples* (**S**(*L*),**B**(*L*))*, where:*

*(i) The matrix* **B**(*L*) *has the parameterization*

$$\mathbf{B}(L) = \xi \eta' + (1 - L)\mathbf{B}^\* + (1 - L)^2 \mathbf{E}(L)\_{\prime\prime}$$

*where ξ and η are r* × (*r* − *c*) *and q* × (*r* − *c*) *respectively,* **B**<sup>∗</sup> *is an r* × *q matrix and* **E**(*L*) *is an r* × *q matrix polynomial of degree s*<sup>1</sup> ≥ 0*.*


*We say that* G *is a rational reduced-rank family of filters with cointegrating rank c.*

**Proposition 3.** *Assume that r* > *q. Let* **y***<sup>t</sup> be a I*(1) *solution of Equation* (10)*, where* (**S**(*L*),**B**(*L*)) *belongs to a rational reduced-rank family of filters with cointegrating rank c. For generic values of the parameters in* Π*, Assumptions 1, 3, 4, 5 and 6 hold. Thus the Strong Form of Proposition 2 holds and* **y***<sup>t</sup> has an error correction representation*

$$\mathbf{A}(L)\mathbf{y}\_{l} = \mathbf{A}^\*(L)(1-L)\mathbf{y}\_{l} + \mathfrak{a}(\boldsymbol{\mathcal{S}}^\prime\mathbf{y}\_{t-1} - \mathbf{w}) = \mathbf{B}(0)\mathbf{u}\_{l}.$$

*where* **A**(*L*) *is a finite-degree polynomial matrix.*

**Proof.** Part (iii) of Definition 5 implies that Assumptions 1 and 4 hold for all **p** ∈ Π. The sets where Assumptions 5 and 6 do not hold are the intersections of the open set Π with the algebraic varieties

$$\mathbf{f}(\mathbf{a}) \text{ rank} \left[ \begin{pmatrix} \mathbf{g}'\_{\perp} \mathbf{B}^\* \\ \mathbf{\eta}' \end{pmatrix} \right] < q, \quad \text{(b)} \ \mathbf{g}'\_{\perp} \left( \mathbf{B}^\* - \mathbf{S}^\*(1) \mathbf{S}(1)^{-1} \mathbf{f} \mathbf{\eta}' \right) = \mathbf{0}$$

(the variety described by (a) is obtained by equating to zero the determinant of all the *q* × *q* submatrices of the *r* × *q* matrix between brackets). It is easy to see that the varieties (a) and (b) are not trivial, i.e., that their dimension is lower than *λ*. Thus Assumptions 5 and 6 hold generically. The same result holds for Assumption 3. The points of Π where it is not fulfilled belong to a lower-dimensional algebraic variety. This is proved in A.1, see in particular Lemma A4.

**Remark 7.** *It is easy to see that, assuming that c* ≤ *q,* rank(*ξ* ⊥(**B**<sup>∗</sup> <sup>−</sup> **<sup>S</sup>**∗(1)**S**(1)−1*ξη* ) = *c holds generically in* Π*. Thus, in that case, the error term β***y***<sup>t</sup>* − **w** *is generically a non-cointegrated I*(0) *process, see Remark 6.*

**Remark 8.** *A general comment on genericity results is in order. Theorems like Proposition 3 or Proposition 1, part* (II)*, show that the subset where some statement does not hold belong to some algebraic variety of lower dimension (see the proof of Proposition 3 in particular), and is therefore negligible from a topological point of view. This suggests the working hypothesis that such subset is negligible from an economic or statistical point of view as well. If, for example, economic theory produces a singular vector* **y***t with cointegrationg rank c, we may find reasonable to conclude that* **y***<sup>t</sup> has representation* (14) *with a finite autoregressive polynomial. However, a greater degree of certainty is obtained by checking that the parameters of* (**S**(*L*), **B**(*L*))*, that are implicit in the theory, do not necessarily lie in one of the three algebraic varieties described in the proof of Proposition 3.*

Definition 5 does not assume that **B**(*L*) has no zeros inside the unit circle. Thus we have not assumed that **u***<sup>t</sup>* is fundamental for (1 − *L*)**F***t*, see Section 2.2. However, Proposition 3 shows that for generic values of the parameters in Π, the assumptions of Proposition 2, strong form, hold, Assumption 3 in particular, so that **B**(*L*) has no zeros of non-unit modulus and therefore inside the unit circle. Thus:

**Proposition 4.** *Assume that r* > *q. Let* **y***<sup>t</sup> be a solution of Equation* (10)*, where* (**S**(*L*),**B**(*L*)) *belongs to a rational reduced-rank family of filters with cointegrating rank c. For generic values of the parameters in* Π*,* **u***<sup>t</sup> is fundamental for* (1 − *L*)**y***t.*

**Remark 9.** *Note that Propositions 3 and 4 do not hold in the nonsingular case, where no genericity argument can be used to rule out non-unit zeros of* **B**(*L*)*, either inside or outside the unit circle. In particular, fundamentalness of* **u***<sup>t</sup> for* (1 − *L*)**y***<sup>t</sup> is not generic if r* = *q.*

#### *3.3. Permanent and Transitory Shocks*

Let *η*<sup>⊥</sup> be a *q* × *d* matrix whose columns are independent and orthogonal to the columns of *η*, and let

$$
\overline{\eta} = \eta(\eta'\eta)^{-1}, \quad \overline{\eta}\_{\perp} = \eta\_{\perp}(\eta'\_{\perp}\eta\_{\perp})^{-1}.
$$

Defining **v**1*<sup>t</sup>* = *η* <sup>⊥</sup>**u***t*, and **<sup>v</sup>**2*<sup>t</sup>* <sup>=</sup> *<sup>η</sup>* **u***t*, we have

$$\mathbf{u}\_{\mathbf{f}} = \overline{\eta}\_{\perp} \mathbf{v}\_{1t} + \overline{\eta} \mathbf{v}\_{2t} = \begin{pmatrix} \overline{\eta}\_{\perp} & \overline{\eta} \end{pmatrix} \begin{pmatrix} \mathbf{v}\_{1t} \\ \mathbf{v}\_{2t} \end{pmatrix}.$$

We have

$$\mathbf{B}(L)\mathbf{u}\_{t} = \left[\mathbf{B}(L)\left(\overline{\eta}\_{\perp}\,\overline{\eta}\right)\right] \begin{pmatrix} \mathbf{v}\_{1t} \\ \mathbf{v}\_{2t} \end{pmatrix} = \left(1 - L\right)\mathbf{G}\_{1}(L)\mathbf{v}\_{1t} + \left(\xi + (1 - L)\mathbf{G}\_{2}(L)\right)\mathbf{v}\_{2t}.$$

where **<sup>G</sup>**1(*L*) = (**B**<sup>∗</sup> + (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)**E**(*L*)) *<sup>η</sup>*⊥, and **<sup>G</sup>**2(*L*) = (**B**<sup>∗</sup> + (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)**E**(*L*)) *<sup>η</sup>*. All the solutions of the difference equation (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)**y***<sup>t</sup>* <sup>=</sup> **<sup>S</sup>**(*L*)−1**C**(*L*)**u***<sup>t</sup>* are

$$\mathbf{y}\_t = \mathbf{S}(L)^{-1} \left[ \mathbf{G}\_1(L)\mathbf{v}\_{1t} + \mathbf{G}\_2(L)\mathbf{v}\_{2t} + \mathbf{T}\_l \right] + \mathbf{K},\tag{21}$$

where **K** is a constant stochastic process, and

$$\mathbf{T}\_t = \begin{cases} \mathfrak{F}(\mathbf{v}\_{21} + \mathbf{v}\_{22} + \dots + \mathbf{v}\_{2t}), \text{ for } t > 0 \\ \mathbf{0}, \text{ for } t = 0 \\ -\mathfrak{F}(\mathbf{v}\_{20} + \mathbf{v}\_{2,-1} + \dots + \mathbf{v}\_{2,t+1}), \text{ for } t < 0. \end{cases}$$

As *ξ* is full rank, we see that **y***<sup>t</sup>* is driven by the *q* − *d* = *r* − *c* permanent shocks **v**2*t*, and by the *d* temporary shocks **v**1*t*. In representation (21), the component **T***<sup>t</sup>* is the common-trend of Stock and Watson (1988). Note that the number of permanent shocks is obtained as *r* minus the cointegrating rank, as usual. However, the number of transitory shocks is only *d* = *c* − (*r* − *q*), as though *r* − *q* transitory shocks had a zero coefficient.

#### *3.4. VECMs and Unrestricted VARs in The Levels*

Several papers have addressed the issue of whether and when an error correction model or an unrestricted VAR in the levels should be used for estimation in the case of nonsingular cointegrated vectors: Sims et al. (1990) have shown that the parameters of a cointegrated VAR are consistently estimated using an unrestricted VAR in the levels; on the other hand, Phillips (1998) shows that if the variables are cointegrated, the long-run features of the impulse-response functions are consistently estimated only if the unit roots are explicitly taken into account, that is within a VECM specification. The simulation exercise described below provides evidence in favour of the VECM specification in the singular case.


The results are shown in Table 1. We see that the RMSE of both the VECM and the LVAR decreases as *T* increases. However, for all values of *T*, the RMSE of the VECM stabilizes as the lag increases, whereas it deteriorates for the LVAR, in line with the claim that the long-run response of the variables are better estimated with the VECM. The performance of the misspecified DVAR is uniformly poor with the exception of lag zero.


**Table 1.** Monte Carlo Simulations. VECM: *r* = 4, *q* = 3, *c* = 3.

Root mean squared errors at different lags, when estimating the impulse-response functions of the simulated variables **y***<sup>t</sup>* to the shocks **u***t*. Estimation is carried out using three different autoregressive representations: a VAR for (1 − *L*)**y***<sup>t</sup>* (DVAR), a VAR for **y***<sup>t</sup>* (LVAR), and a VECM with *c* = *r* − *q* + *d* error terms (VECM). The results are based on 1000 replications. For the data generating process see Appendix C. The RMSEs are obtained averaging over all replications and all 4 × 3 responses.

#### **4. Cointegration of the Observable Variables in a DFM**

Consider again the factor model *xit* = *χit* + *it*, rewritten here as

$$\mathbf{x}\_{l} = \mathbf{x}\_{l} + \mathbf{e}\_{l}, \quad \mathbf{x}\_{l} = \boldsymbol{\Lambda} \mathbf{F}\_{l} \tag{22}$$

where **Λ** is *n* × *r*, with *n* > *r*. The relationship between cointegration of the factors **F***<sup>t</sup>* and cointegration of the variables *xit* is now considered.

Let us recall that the the common factors **F***jt* are assumed to be orthogonal to the idiosyncratic components *ks* for all *i*, *j*, *t*,*s*, i.e., E*χt <sup>s</sup>* = **0***n*×*n*. for all *t*,*s*, see the Introduction. The other assumptions on model (22) are asymptotic, see e.g., Forni et al. (2000); Forni and Lippi (2001); (Stock and Watson 2002a, 2002b), and put no restriction on the matrix **Λ** and the vector *<sup>t</sup>* for a given finite *n*. In particular, the first *r* eigenvalues of the matrix **ΛΛ** must diverge as *n* → ∞, but this has no implications on the rank of the matrix **Λ** corresponding to, say, *n* = 10. However, as we see in Proposition 5 (iii), if the idiosyncratic components are *I*(0), then, independently of **Λ**, all *p*-dimensional subvectors of **x***<sup>t</sup>* are cointegrated for *p* > *q* − *d*, which is at odds with what is observed in the macroeconomic datasets analyzed in the empirical Dynamic Factor Model literature. This motivates assuming that *<sup>t</sup>* is *I*(1). In that case, see Proposition 5 (i), cointegration of **x***<sup>t</sup>* requires that both the common and the idiosyncratic components are cointegrated. Some results are collected in the statement below.

**Proposition 5.** *Let* **x** (*p*) *<sup>t</sup>* <sup>=</sup> *<sup>χ</sup>*(*p*) *<sup>t</sup>* + (*p*) *<sup>t</sup>* = **<sup>Λ</sup>**(*p*)**F***<sup>t</sup>* + (*p*) *<sup>t</sup> be a p-dimensional subvector of* **x***t, p* ≤ *n. Denote by c p <sup>χ</sup> and c p the cointegrating rank of <sup>χ</sup>*(*p*) *<sup>t</sup> and* (*p*) *<sup>t</sup> respectively. Both range from p, stationarity, to* 0*, no cointegration.*


**Proof.** Because *χit* and *js* are orthogonal for all *i*, *j*, *t*,*s*, the spectral densities of (1 − *L*)**x** (*p*) *<sup>t</sup>* , (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)*χ*(*p*) *<sup>t</sup>* , (1 − *L*) (*p*) *<sup>t</sup>* fulfill:

$$
\Delta\_{\Lambda\mathfrak{x}}^{(p)}(\theta) = \Delta\_{\Lambda\chi}^{(p)}(\theta) + \Delta\_{\Lambda\mathfrak{c}}^{(p)}(\theta) \quad \theta \in [-\pi, \pi]. \tag{23}
$$

Now, (23) implies that

$$
\lambda\_p \left( \Sigma\_{\Lambda x}^{(p)}(0) \right) \ge \lambda\_p \left( \Sigma\_{\Lambda \chi}^{(p)}(0) \right) + \lambda^{(p)} \left( \Sigma\_{\Lambda x}^{(p)}(0) \right), \tag{24}
$$

where *λp*(*A*) denotes the smallest eigenvalue of the hermitian matrix *A*; this is one of the Weyl's inequalities, see Franklin (2000), p. 157, Theorem 1. Because the spectral density matrices are non-negative definite, the right hand side in (24) vanishes if and only if both terms on the right hand side vanish, i.e., the spectral density of Δ**x** (*p*) *<sup>t</sup>* is singular at zero if and only if the spectral densities of <sup>Δ</sup>*χ*(*p*) *<sup>t</sup>* and Δ (*p*) *<sup>t</sup>* are singular at zero. By definition 4, (i) is proved.

Without loss of generality we can assume that **S**(*L*) = **I***r*. By substituting (21) in (22), we obtain

$$\mathbf{x}\_{l} = \mathbf{A}\left[ (\mathbf{G}\_{1}(L)\mathbf{v}\_{1t} + \mathbf{G}\_{2}(L)\mathbf{v}\_{2t} + \mathbf{T}\_{l}) + \mathbf{K} \right] + \mathbf{e}\_{l},\tag{25}$$

where on the right hand side the only non-stationary terms are **T***<sup>t</sup>* and possibly *t*. By recalling that **T***<sup>t</sup>* = *ξ* ∑*<sup>t</sup> <sup>s</sup>*=<sup>1</sup> **v**2*<sup>s</sup>* where *ξ* is of dimension *r* × (*q* − *d*) and rank *q* − *d*, and by defining *G<sup>t</sup>* = **<sup>Λ</sup>**[**G**1(*L*)**v**1*<sup>t</sup>* + **<sup>G</sup>**2(*L*)**v**2*<sup>t</sup>* + **<sup>K</sup>**] and *<sup>T</sup><sup>t</sup>* = <sup>∑</sup>*<sup>t</sup> <sup>s</sup>*=<sup>1</sup> **v**2*s*, we can rewrite (25) as

$$\mathbf{x}\_t = \Lambda \mathbf{\mathcal{J}} \mathcal{T}\_t + \mathbf{\mathcal{G}}\_t + \boldsymbol{\epsilon}\_t.$$

For **x** (*p*) *<sup>t</sup>* :

$$\mathbf{x}\_t^{(p)} = \mathbf{x}\_t^{(p)} + \mathbf{e}\_t^{(p)} = \boldsymbol{\Lambda}^{(p)} \boldsymbol{\xi} \mathcal{T}\_t + \mathbf{g}\_t^{(p)} + \boldsymbol{\varepsilon}\_t^{(p)},$$

where **<sup>Λ</sup>**(*p*) and *<sup>G</sup>*(*p*) *<sup>t</sup>* have an obvious definition. Of course cointegration of the common components *χ*(*p*) *<sup>t</sup>* is equivalent to cointegration of **<sup>Λ</sup>**(*p*)*ξTt*, which in turn is equivalent to rank(**Λ**(*p*)*ξ*) <sup>&</sup>lt; *<sup>p</sup>*. Statement (ii) follows from

$$\text{rank}\left(\boldsymbol{\Lambda}^{(p)}\boldsymbol{\xi}\right) \le \min\left(\text{rank}(\boldsymbol{\Lambda}^{(p)}), \text{rank}(\boldsymbol{\xi})\right).$$

The first part of (iii) is obvious. Assume now that *p* > *q* − *d*. If *c p <sup>χ</sup>* + *c p* = dim(*Vχ*) + dim(*V*) = *p* − (*q* − *d*) + *c p* > *p*, i.e., if *c p* <sup>&</sup>gt; *<sup>q</sup>* <sup>−</sup> *<sup>d</sup>*, then the intersection between *<sup>V</sup><sup>χ</sup>* and *<sup>V</sup>* is non-trivial, so that **x** (*p*) *<sup>t</sup>* is cointegrated.

#### **5. Summary and Conclusions**

The paper studies representation theory for singular *I*(1) stochastic vectors, the factors of an *I*(1) Dynamic Factor Model in particular. Singular *I*(1) vectors are cointegrated, with a cointegrating rank *c* equal to *r* − *q*, the dimension of **y***<sup>t</sup>* minus its rank, plus *d*, with 0 ≤ *d* < *q*.

If (1 − *L*)**y***<sup>t</sup>* has rational spectral density, under assumptions that generalize to the singular case those in Johansen (1995), we show that **y***<sup>t</sup>* has an error correction representation with *c* error terms, thus generalizing the Granger representation theorem (from MA to AR) to the singular case. Important consequences of singularity are that generically: (i) the autoregressive matrix polynomial of the error correction representation is of finite degree, (ii) the white noise vector driving (1 − *L*)**y***<sup>t</sup>* is fundamental.

We find that **y***<sup>t</sup>* is driven by *r* − *c* permanent shocks and *d* = *c* − (*r* − *q*) transitory shocks, not *c* as in the nonsingular case.

Using simulated data generated by a simple singular VECM, confirms previous results, obtained for nonsingular vectors, showing that under cointegration the long-run features of impulse-response functions are better estimated using a VECM rather than a VAR in the levels.

In Section 4 we argue that stationarity of the idiosyncratic components in a DFM produce an amount of cointegration for the observable variables *xit* that is not observed in the datasets that are standard in applied Dynamic Factor Model literature. Thus the idiosyncratic vector in those datasets is likely to be *I*(1), so that an estimation strategy robust to the assumption that some of the idiosyncratic variables *it* are *I*(1) should be preferred.

The results in this paper are the basis for estimation of *I*(1) Dynamic Factor Models with cointegrated factors, which is developed in the companion paper (Barigozzi et al. 2019).

**Author Contributions:** All authors contributed equally to the paper. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** Dietmar Bauer, Manfred Deistler, Massimo Franchi, Martin Wagner, three anonymous referees and the Editors of this Special Issue gave important suggestions for improvements. We also thank the participants to the Workshop on Estimation and Inference Theory for Cointegrated Processes in the State Space Representation, Technische Universität Dortmund, January 2016. Part of this paper was written while Matteo Luciani was chargé de recherches F.R.S.- F.N.R.S., and he gratefully acknowledges their financial support. Of course we are responsible for any remaining errors.

**Disclaimer**: The views expressed in this paper are those of the authors and do not necessarily reflect those of the Board of Governors or the Federal Reserve System.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A. Proofs**

*Appendix A.1. Assumption 3 Holds Generically*

Proving that Assumption 3 holds generically is equivalent to proving that **M**(*z*) is generically zeroless, see the argument below Equation (20).

We need some preliminary results. Lemma A1, though quite easy, is not completely standard and is therefore carefully stated and proved below. Regarding notation, to avoid possible misunderstandings, let us recall that vectors and matrices are always denoted by boldface symbols, while light symbols denote scalars, see Lemmas A1 and A2 in particular.

**Lemma A1.** *Let Aj, <sup>j</sup>* <sup>=</sup> 1, ... ,*s, be scalar polynomials defined on* <sup>R</sup>*λ, let* **<sup>p</sup>** <sup>∈</sup> <sup>R</sup>*<sup>λ</sup> and <sup>Q</sup>*(**p**) *be the statement*

$$A\_j(\mathbf{p}) = 0, \quad \text{for } j = 1, \dots, s\_{\prime}$$

*for example the statement that all the q* × *q minors of* **M**(1) *vanish, i.e., that* rank(**M**(1)) < *q. Let* Π *be an open subset of* <sup>R</sup>*λ. If Q is false for one point* **<sup>p</sup>**<sup>∗</sup> <sup>∈</sup> <sup>R</sup>*λ, then Q is generically false in* <sup>Π</sup>*.*

**Proof.** Let N be the closure in Π (in the topology of Π) of the subset of Π where *Q* is true. Suppose that *Q* is not generically false in Π. Then the interior of N in Π, call it N ◦, is not empty. As Π is open, <sup>N</sup> ◦ is open both in the topology of <sup>Π</sup> and of <sup>R</sup>*λ*. On the other hand a polynomial function defined on R*<sup>λ</sup>* vanishes on an open set if and only if it vanishes on the whole R*λ*, which contradicts the existence of a point in R*<sup>λ</sup>* where *Q* is false.

**Lemma A2.** *Consider the scalar polynomials*

$$A(z) = a\_0 z^n + a\_1 z^{n-1} + \dots + a\_m \quad B(z) = b\_0 z^m + b\_1 z^{m-1} + \dots + a\_{m+1}$$

*with a*<sup>0</sup> = 0 *and b*<sup>0</sup> = 0*, and let αi, i* = 1, ... , *n and βj, j* = 1, ... , *m, be the roots of A and B, respectively. Then: (i)*

$$a\_0^m b\_0^n \prod\_{i,j} (a\_i - \beta\_j) = \mathcal{R}(a\_0, a\_1, \dots, a\_n; b\_0, b\_1, \dots, b\_m),$$

*where R is a polynomial function which is called the resultant of A and B. (ii) The resultant vanishes if and only if A and B have a common root. (iii) Suppose that the coefficients ai and bj are polynomial functions of* **<sup>p</sup>** <sup>∈</sup> <sup>Π</sup>*, where* <sup>Π</sup> *is an open subset of* <sup>R</sup>*λ. If there exists a point* **<sup>p</sup>**<sup>∗</sup> <sup>∈</sup> <sup>R</sup>*<sup>λ</sup> such that <sup>a</sup>*0(**p**∗) <sup>=</sup> <sup>0</sup>*, <sup>b</sup>*0(**p**∗) <sup>=</sup> <sup>0</sup>*, and R*(**p**∗) = 0*, then generically in* Π *the polynomials A and B have no common roots.*

**Proof.** For (i) and (ii) see van der Waerden (1953, pp. 83-8). Statement (iii) is an obvious consequence of (ii) and Lemma A1.

**Lemma A3.** *Recall that a zero of* **M**(*z*) *is a complex number z*∗ *such that* rank(**M**(*z*∗)) < *q. If* **M**(*z*) *has two q* × *q submatrices whose determinants have no common roots, then* **M**(*z*) *is zeroless.*

**Proof.** If *z*<sup>∗</sup> is a zero of **M**(*z*), then *z*<sup>∗</sup> is a zero of all the *q* × *q* submatrices of **M**(*z*).

For the statement and proof of our last result it is convenient to make explicit the dependence of the matrix **M**(*z*) and its submatrices on the vector **p**. Thus we use **Mp**(*z*), etc. The parameters of the matrix **S**(*L*) play no role here. Hence, with no loss of generality, we assume *s*<sup>2</sup> = 0, so that *λ* = (*r* − *c*)(*r* + *q*) + *rq*(*s*<sup>1</sup> + 2). Lemmas A2–A4 below imply that Assumption 3 holds generically in Π.

**Lemma A4.** *Let* **M<sup>p</sup>** <sup>1</sup> (*z*), **<sup>M</sup><sup>p</sup>** <sup>2</sup> (*z*), ... *be all the <sup>q</sup>* <sup>×</sup> *<sup>q</sup> submatrices of* **<sup>M</sup>p**(*z*) *and let* <sup>L</sup>**<sup>p</sup>** *<sup>i</sup> be the leading coefficient of* det**M<sup>p</sup>** *<sup>i</sup>* (*z*) *and R***<sup>p</sup>** *ij is the resultant of* det**Mp** *<sup>i</sup>* (*z*) *and* det **Mp** *<sup>j</sup>* (*z*)*. There exist i, j,* **<sup>p</sup>**<sup>∗</sup> <sup>∈</sup> <sup>R</sup>*<sup>λ</sup> such that*

$$
\mathcal{L}\_i^{\mathbb{P}^\*} \mathcal{L}\_j^{\mathbb{P}^\*} \neq 0
$$

*and*

$$\mathcal{R}\_{ij}^{\mathbb{P}^\*} \neq 0.$$

**Proof.** Assume that *r* = *q* + 1. To each **p** ∈ Π there corresponds the matrix

$$\mathbf{M}^{\mathbf{P}}(z) = \begin{pmatrix} \mathfrak{E}'\_{\perp} \mathbf{B}^\* \\ \mathfrak{E}' \mathfrak{E} \mathfrak{p}' \end{pmatrix} + (1 - z) \begin{pmatrix} \mathfrak{E}'\_{\perp} \mathbf{E}(z) \\ \mathfrak{E}' \mathbf{B}^\* \end{pmatrix} + (1 - z)^2 \begin{pmatrix} \mathbf{0}\_{\mathbb{C} \times q} \\ \mathfrak{E}' \mathbf{E}(z) \end{pmatrix}.$$

Of course, the definition of **<sup>M</sup>p**(*z*) makes sense for all **<sup>p</sup>** <sup>∈</sup> <sup>R</sup>*λ*, see Equation (19). Let **<sup>M</sup><sup>p</sup>** <sup>1</sup> (*z*) and **M<sup>p</sup>** <sup>2</sup> (*z*) be the matrices obtained from **<sup>M</sup>p**(*z*) by removing the first and the last row respectively. We have:

$$\begin{aligned} \text{degree}[\det(\mathbf{M}\_1^\mathbf{P}(z))] &\le (q-d)(s\_1+2) + d(s\_1+1) = d\_1, \\ \text{degree}[\det(\mathbf{M}\_2^\mathbf{P}(z))] &\le (q-d-1)(s\_1+2) + (d+1)(s\_1+1) = d\_2. \end{aligned}$$

We will construct a point **<sup>p</sup>**<sup>∗</sup> <sup>∈</sup> <sup>R</sup>*<sup>λ</sup>* such that: (A) the coefficient of *<sup>z</sup>d*<sup>1</sup> in det(**Mp**<sup>∗</sup> <sup>1</sup> (*z*)) and the coefficient of *zd*<sup>2</sup> in det(**Mp**<sup>∗</sup> <sup>2</sup> (*z*)) (the leading coefficients) do not vanish, (B) the resultant of det(**Mp**<sup>∗</sup> <sup>1</sup> (*z*)) and det(**Mp**<sup>∗</sup> <sup>2</sup> (*z*)) does not vanish.

Let us firstly define a family of matrices, denoted by **M**(*z*), obtained by specifying *η*, *ξ*, *ξ* ⊥, **<sup>B</sup>**<sup>∗</sup> and **E**(*z*) in the following way:

$$\begin{aligned} \underline{\mathbf{u}}' &= \begin{pmatrix} \mathbf{0}\_{(q-d)\times d} & \mathbf{I}\_{q-d} \end{pmatrix}, \underline{\mathbf{g}} = \begin{pmatrix} \mathbf{I}\_{q-d} \\ \mathbf{0}\_{c\times(q-d)} \end{pmatrix}, \underline{\mathbf{g}}'\_{\perp} = \begin{pmatrix} \mathbf{K} \\ \mathbf{H} \end{pmatrix}, \\\underline{\mathbf{B}}" &= \begin{pmatrix} \mathbf{H}' & \mathbf{0}\_{(q+1)\times(q-d)} \end{pmatrix}, \underline{\mathbf{E}}(z) = \begin{pmatrix} \mathbf{E}\_{1}(z) \\ \mathbf{E}\_{2}(z) \\ \mathbf{E}\_{3}(z) \end{pmatrix}, \end{aligned}$$

where:

$$\begin{aligned} \mathbf{K} &= \begin{pmatrix} \mathbf{0}\_{1 \times (q-d)} & 1 & \mathbf{0}\_{1 \times d} \end{pmatrix}, \quad \mathbf{H} = \begin{pmatrix} \mathbf{0}\_{d \times (q+1-d)} & \mathbf{I}\_d \end{pmatrix}, \\\\ \mathbf{E}\_1(z) &= \begin{pmatrix} k\_1(z) & h\_1(z) & \cdots & 0 \\\\ \mathbf{0}\_{(q-d) \times d} & \ddots & \ddots & \\\\ & & \ddots & \ddots & \\\\ & & 0 & \cdots & k\_{q-d}(z) \end{pmatrix}, \quad \mathbf{E}\_2(z) = \begin{pmatrix} e(z) & \mathbf{0}\_{1 \times (q-1)} \end{pmatrix}, \\\\ &\mathbf{0} & \cdots & k\_{q-d}(z) \end{pmatrix}, \quad \mathbf{E}\_3(z) = \begin{pmatrix} e(z) & \mathbf{0}\_{1 \times (q-1)} \end{pmatrix}, \quad \mathbf{E}\_4(z) = \begin{pmatrix} \mathbf{0}\_{d \times (q-1)} & \mathbf{0}\_{d \times (q-1)} \end{pmatrix}, \quad \mathbf{E}\_5(z) = \begin{pmatrix} \mathbf{0}\_{d \times (q-1)} & \mathbf{0}\_{d \times (q-1)} \end{pmatrix}, \quad \mathbf{E}\_6(z) = \begin{pmatrix} \mathbf{0}\_{d \times (q-1)} & \mathbf{0}\_{d \times (q-1)} \end{pmatrix}, \quad \mathbf{E}\_7(z) = \begin{pmatrix} \mathbf{0}\_{d \times (q-1)} & \mathbf{0}\_{d \times (q-1)} \end{pmatrix}. \end{aligned}$$

the entries *e*, *ki*, *hi*, *fi* and *gi* being scalar polynomials of degree *s*1.

0 ··· *fd*(*z*) *gd*(*z*)

We denote by **q1** the vector including the coefficients of the polynomials *fi*, *i* = 1, ... , *d* and *ki*, *i* = 1, ... ,(*q* − *d*), a total of *q*(*s*<sup>1</sup> + 1) coefficients, by **q**<sup>2</sup> the vector including the coefficients of the polynomials *e*, *gi*, *i* = 1, ... , *d* and *hi*, *i* = 1, ... ,(*q* − *d* − 1), a total of *q*(*s*<sup>1</sup> + 1) coefficients, by **q**<sup>0</sup> the vector including the zeros and the ones in the definition of *ξ*, *η*, **B**∗, **E**, and define **q** = (**q**<sup>0</sup> **q**<sup>1</sup> **q**2), which is a *λ*-dimensional parameter vector. We put no restriction on **q**<sup>1</sup> and **q**2, so that both can take any value in R*ν*, with *ν* = *q*(*s*<sup>1</sup> + 1). Note that **q** does not necessarily belong to Π. We have:

$$\underline{\mathbf{M}}^{\mathbf{q}}(z) = \begin{pmatrix} \mathbf{0}\_{1 \times d} & \mathbf{0}\_{1 \times (q-d)} \\ \mathbf{I}\_d & \mathbf{0}\_{d \times (q-d)} \\ \mathbf{0}\_{(q-d) \times d} & \mathbf{I}\_{q-d} \end{pmatrix} + (1-z) \begin{pmatrix} \mathbf{E}\_2(z) \\ \mathbf{E}\_3(z) \\ \mathbf{0}\_{(q-d) \times q} \end{pmatrix} + (1-z)^2 \begin{pmatrix} \mathbf{0}\_{1 \times q} \\ \mathbf{0}\_{d \times q} \\ \mathbf{E}\_1(z) \end{pmatrix}.\tag{A1}$$

The matrix **Mq**(*z*) has zero entries except for the diagonal joining the positions (1, 1) and (*q*, *q*), and the diagonal joining (2, 1) and (*q* + 1, *q*). The matrices **M<sup>q</sup>** <sup>1</sup> (*z*) and **<sup>M</sup><sup>q</sup>** <sup>2</sup> (*z*) are upper- and lower-triangular, respectively, and

$$\begin{split} \det(\underline{\mathbf{M}}\_1^{\mathbf{q}}(z)) &= \left[1 + (1 - z)f\_1(z)\right] \cdots \left[(1 + (1 - z)f\_d(z))\right] \\ &\quad \times \left[1 + (1 - z)^2 k\_1(z)\right] \cdots \left[1 + (1 - z)^2 k\_{q - d}(z)\right] = \mathcal{L}\_{1,d\_1}^{\mathbf{q}} z^{d\_1} + \cdots + \mathcal{L}\_{1,0}^{\mathbf{q}} \\ \det(\underline{\mathbf{M}}\_2^{\mathbf{q}}(z)) &= (1 - z)^{2q - d - 1} c(z) [\underline{g}\_1(z) \cdots \underline{g}\_d(z)] [h\_1(z) \cdots h\_{q - d - 1}(z)] = \mathcal{L}\_{2,d\_2}^{\mathbf{q}} z^{d\_1} + \cdots + \mathcal{L}\_{2,0}^{\mathbf{q}}. \end{split}$$

Note that det(**Mq** <sup>1</sup> (*z*)) does not depend on **<sup>q</sup>**2, while det(**M<sup>q</sup>** <sup>2</sup> (*z*)) does not depend on **q**1. Thus we use the notation *δ* **q**1 <sup>1</sup> (*z*) = det(**M<sup>q</sup>** <sup>1</sup> (*z*)), *δ* **q**2 <sup>2</sup> (*z*) = det(**M<sup>q</sup>** <sup>2</sup> (*z*)), <sup>M</sup>**q**<sup>1</sup> 1,*d*<sup>1</sup> <sup>=</sup> <sup>L</sup>**<sup>q</sup>** 1,*d*<sup>1</sup> , <sup>M</sup>**q**<sup>2</sup> 2,*d*<sup>2</sup> <sup>=</sup> <sup>L</sup>**<sup>q</sup>** 2,*d*<sup>2</sup> . Now:


$$\det(\mathbf{M}\_1^{\mathbf{P}^\*}(z)) = \delta\_1^{\mathbf{q}\_1^\*}(z)), \quad \det(\mathbf{M}\_2^{\mathbf{P}^\*}(z)) = \delta\_2^{\mathbf{q}\_2^\*}(z).$$

Using (i) and (ii), (A) the leading coefficients of det(**Mp**<sup>∗</sup> <sup>1</sup> (*z*)) and det(**Mp**<sup>∗</sup> <sup>2</sup> (*z*)) do not vanish, (B) det(**Mp**<sup>∗</sup> <sup>1</sup> (*z*)) and det(**Mp**<sup>∗</sup> <sup>2</sup> (*z*)) have no root in common so that their resultant does not vanish. This proves the proposition for *r* = *q* + 1.

Generalizing this result to *r* > *q* + 1 is easy. Let us define the family **N**(*z*) in the following way: (a) specify *η* , *ξ*, **E**1(*z*) and **E**3(*z*) as in the definition of **M**(*z*), (b) then let

,

$$\begin{aligned} \mathbf{K} &= \begin{pmatrix} \mathbf{0}\_{(r-q)\times(r-d-1)} & \begin{pmatrix} \mathbf{0}\_{1\times(r-q-1)} \mathbf{1} \end{pmatrix}' & \mathbf{0}\_{(r-q)\times d} \end{pmatrix}, \ \mathbf{H} = \begin{pmatrix} \mathbf{0}\_{d\times(r-d)} & \mathbf{I}\_d \end{pmatrix}, \\\ \mathbf{f}'\_{\perp} &= \begin{pmatrix} \mathbf{K} \\ \mathbf{H} \end{pmatrix}, \ \mathbf{\underline{D}} = \begin{pmatrix} \mathbf{H}' & \mathbf{I}\_{r\times(q-d)} \end{pmatrix}, \ \mathbf{E}\_2(z) = \begin{pmatrix} \mathbf{0}\_{(r-q)\times q} \\ \begin{pmatrix} c(z) & \mathbf{0}\_{1\times(q-1)} \end{pmatrix} \end{pmatrix}. \end{aligned}$$

We have:

$$\mathbf{M}(z) = \begin{pmatrix} \mathbf{0}\_{(r-q)\times d} & \mathbf{0}\_{(r-q)\times(q-d)} \\ \mathbf{I}\_d & \mathbf{0}\_{d\times(q-d)} \\ \mathbf{0}\_{(q-d)\times d} & \mathbf{I}\_{q-d} \end{pmatrix} + (1-z) \begin{pmatrix} \mathbf{E}\_2(z) \\ \mathbf{E}\_3(z) \\ \mathbf{0}\_{(q-d)\times q} \end{pmatrix} + (1-z)^2 \begin{pmatrix} \mathbf{0}\_{(r-q)\times q} \\ \mathbf{0}\_{d\times q} \\ \mathbf{E}\_1(z) \end{pmatrix}.$$

It is easy to see that the (*<sup>q</sup>* <sup>+</sup> <sup>1</sup>) <sup>×</sup> *<sup>q</sup>* lower submatrix of **<sup>N</sup>**(*z*) is identical to the matrix **Mq**(*z*) in (A1).

*Appendix A.2. if R* > *Q and C* ≤ *Q, Assumptions 5 and 6 Do Not Imply That* **e***<sup>t</sup> Is a Non-Cointegrated I*(0) *Process.*

$$\text{Let } r = 3, q = 2, \mathbf{S}(L) = \mathbf{l}\_{3'} $$

$$\mathfrak{F} = \begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix}, \quad \eta = \begin{pmatrix} 0 \\ 1 \end{pmatrix}, \quad \mathfrak{F}\_{\perp} = \begin{pmatrix} 0 & 0 \\ 1 & 0 \\ 0 & 1 \end{pmatrix}, \quad \mathbb{B}^\* = \begin{pmatrix} a & b \\ 1 & 0 \\ 1 & 0 \end{pmatrix}.$$

In this case *c* = 2 and *d* = 1, so that *c* = *q* (see Remark 6). We have

$$
\begin{pmatrix} \mathfrak{F}'\_{\perp} \mathbf{B}^\* \\ \mathfrak{h}' \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 1 & 0 \\ \hline 0 & 1 \end{pmatrix} \dots
$$

We see that Assumptions 5 and 6 hold. However, rank(*ξ* ⊥**B**∗) = 1, so that **<sup>e</sup>***t*, though being *<sup>I</sup>*(0), is not a non-cointegrated *I*(0) process. On the other hand, if the (3, 2) entry of **B**∗ is 1 instead of 0, **e***t* is non-cointegrated.

#### **Appendix B. Non Uniqueness**

In Proposition 3 we prove that a singular *I*(1) vector with cointegrating rank *c* has a finite error correction representation with *c* error terms. On the other hand, as we have seen in Remark 5, when *c* = *r* − *q* the singular vector **y***<sup>t</sup>* has also an autoregressive representation in the differences, i.e., a representation with zero error terms. In Appendix B.1 we give an example hinting that **y***t* has error correction representations with any number of error terms between *d* and *c*. However, in Appendix B.2 we show that all such representations produce the same impulse-response functions.

*Appendix B.1. Alternative Representations with Different Numbers of Error Terms*

Let **S**(*L*) = **I***<sup>r</sup>* and consider the following example, with *r* = 3, *q* = 2, *c* = 2, so that *d* = 1:

$$\begin{aligned} \mathfrak{F}' &= \begin{pmatrix} 1 & 1 & 1 \end{pmatrix} \\ \mathfrak{p}' &= \begin{pmatrix} 1 & 2 \end{pmatrix} \\ \mathfrak{F}'\_{\perp} &= \begin{pmatrix} 1 & -1 & 0 \\ 0 & 1 & -1 \end{pmatrix} \end{aligned}$$

We have,

$$(1 - L) \begin{pmatrix} \mathfrak{F}'\_{\perp} \\ \mathfrak{F}'\_{\perp} \end{pmatrix} \mathbf{y}\_{\ell} = \begin{pmatrix} 1 - L & 0 & 0 \\ 0 & 1 - L & 0 \\ 0 & 0 & 1 \end{pmatrix} \left\{ \begin{pmatrix} b\_{11}^{\*} - b\_{21}^{\*} & b\_{12}^{\*} - b\_{22}^{\*} \\ b\_{21}^{\*} - b\_{31}^{\*} & b\_{22}^{\*} - b\_{32}^{\*} \\ 3 & 6 \end{pmatrix} + (1 - L)\mathbf{\hat{E}}(L) \right\} \mathbf{u}\_{\ell}.$$

where (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)**Eˆ**(*L*) gathers the second and third terms in **<sup>M</sup>**(*L*). If the assumptions of Proposition <sup>2</sup> hold, we obtain an error correction representation with error terms

$$
\xi'\_{\perp} \mathbf{y}\_t = \begin{pmatrix} y\_{1t} - y\_{2t} \\ y\_{2t} - y\_{3t} \end{pmatrix}.
$$

However, we also have

$$\begin{aligned} &(1-L)\begin{pmatrix} \mathfrak{E}'\_{\perp} \\ \mathfrak{E}' \end{pmatrix} \mathbf{y}\_t = \begin{pmatrix} 1-L & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} \\ &\times \left\{ \begin{pmatrix} b\_{11}^\* - b\_{21}^\* & b\_{12}^\* - b\_{22}^\* \\ (1-L)(b\_{21}^\* - b\_{31}^\*) & (1-L)(b\_{22}^\* - b\_{32}^\*) \\ 3 & 6 \end{pmatrix} + (1-L)\mathbb{E}(L) \right\} \mathbf{u}\_t = \begin{pmatrix} 1-L & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} \mathbf{\breve{N}}(L)\mathbf{u}\_t. \end{aligned}$$

Under suitable assumptions on the coefficients *b*∗ *ij* and **Eˇ**(*L*), assuming in particular that the matrix

$$
\begin{pmatrix} b\_{11}^\* - b\_{21}^\* & b\_{12}^\* - b\_{22}^\* \\ 3 & 6 \end{pmatrix}
$$

is nonsingular, the matrix **Mˇ** (*L*) is zeroless and has therefore a finite-degree left inverse. Proceeding as in Proposition 2, we obtain an alternative error correction representation with just one error term, namely *y*1*<sup>t</sup>* − *y*2*t*.

This example should be sufficient to convey the idea that **y***t* admits error correction representations with a minimum *d* and a maximum *c* = *r* − *q* + *d* of error terms.

The problem of error correction representations, with different numbers of error terms, has been recently addressed in Deistler and Wagner (2017). An implication of their main result (see Theorem 1, p. 41) is that if **y***t* has the error correction representation

$$\mathbf{\tilde{A}}(L)\mathbf{y}\_t = \mathbf{\tilde{A}}^\*(L)(1-L)\mathbf{y}\_t + \mathbf{\tilde{A}}(1)\mathbf{y}\_{t-1} = \mathbf{\tilde{B}}\mathbf{\tilde{u}}\_{t-1}$$

and rank(**A˜** (1)) < *c* (the number of error terms is not the maximum), then **A˜** (*L*) and **B˜** are not left coprime.

The consequences of Deistler and Wagner's paper have not yet been developed. In Propositions 2 and 3 we have only considered representations with *c* error terms. On non-uniqueness of autoregressive representations for singular vectors with rational spectral density see also Chen et al. (2011); Anderson et al. (2012); Forni et al. (2015).

#### *Appendix B.2. Uniqueness of Impulse-Response Functions*

Suppose that the assumptions of Proposition 2, weak form, hold. Let **y***<sup>t</sup>* be a solution of Equation (10), so that

$$(1 - L)\mathbf{y}\_t = \mathbf{S}(L)^{-1}\mathbf{B}(L)\mathbf{u}\_{t\prime} \tag{A2}$$

and suppose that **y***t* has the autoregressive representation

$$
\tilde{\mathbf{A}}(L)\mathbf{y}\_l = \tilde{\mathbf{B}}\tilde{\mathbf{u}}\_{l'} \tag{A3}
$$

where **A˜** (*L*) is a rational matrix with poles outside the unit circle, **A˜** (0) = **I***r*, **u˜***<sup>t</sup>* is a nonsingular *<sup>q</sup>*-dimensional white noise, **B˜** is a full rank *<sup>r</sup>* <sup>×</sup> *<sup>q</sup>* matrix5. We have

$$
\tilde{\mathbf{A}}(L)[(1-L)\mathbf{y}\_l] = (1-L)\tilde{\mathbf{B}}\tilde{\mathbf{u}}\_l.\tag{A4}
$$

The assumption that **B˜** is full rank and the argument used e.g., in Brockwell and Davis (1991), p. 111, Problem 3.8, imply that **u˜***<sup>t</sup>* is fundamental for (1 − *L*)**y***t*. Thus **u˜***<sup>t</sup>* = **Qu***t*, where **Q** is a nonsingular *<sup>q</sup>* <sup>×</sup> *<sup>q</sup>* matrix (see Rozanov (1967), p. 57), and **B ˜ ˜ <sup>u</sup>***<sup>t</sup>* = [**BQ˜** ]**u***t*.

On the other hand, from (A2) and (A4):

$$\mathbf{\tilde{A}}(L)\mathbf{S}(L)^{-1}\mathbf{B}(L)\mathbf{u}\_{l} = (1-L)[\mathbf{\tilde{B}}\mathbf{Q}]\mathbf{u}\_{l}.\tag{A5}$$

As **<sup>u</sup>***<sup>t</sup>* is nonsingular, **A˜** (*L*)**S**(*L*)−1**B**(*L*)=(<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)[**BQ˜** ]. Setting *<sup>L</sup>* <sup>=</sup> 0 we have **BQ˜** <sup>=</sup> **<sup>B</sup>**(0), so that (A3) becomes

$$
\tilde{\mathbf{A}}(L)\mathbf{y}\_t = \mathbf{B}(0)\mathbf{u}\_t \tag{A6}
$$

while (A5) becomes

$$\mathbf{\tilde{A}}(L)\mathbf{S}(L)^{-1}\mathbf{B}(L)\mathbf{u}\_{l} = (1-L)\mathbf{B}(0)\mathbf{u}\_{l}.\tag{A7}$$

The impulse-response function of **y***<sup>t</sup>* to **u***<sup>t</sup>* resulting from (A6) is **H**(*L*)**B**(0), where **H**(*L*)**A˜** (*L*) = **I***r*. Multiplying both sides of (A7) by **H**(*L*) we obtain

$$\mathbf{S}(L)^{-1}\mathbf{B}(L) = (1 - L)\mathbf{H}(L)\mathbf{B}(0),$$

so that **H**(*L*)**B**(0) is obtained by cumulating **S**(*L*)−1**B**(*L*) and is therefore independent of **A˜** (*L*).

#### **Appendix C. Data Generating Process for the Simulations**

The simulation results of Section 3.4 are obtained using the following specification of (14):

$$\mathbf{A}(L)\mathbf{y}\_t = \mathbf{A}^\*(L)(1-L)\mathbf{y}\_t + \mathfrak{a}\beta^\prime\mathbf{y}\_{t-1} = \mathbf{C}(0)\mathbf{u}\_t = \mathbf{G}\mathbf{H}\mathbf{u}\_{t\wedge\tau}$$

where *r* = 4, *q* = 3, *c* = 3, the degree of **A**(*L*) is 2, so that the degree of **A**∗(*L*) is 1. **A**(*L*) is generated using the factorization

$$\mathbf{A}(L) = \mathcal{U}(L)\mathcal{M}(L)\mathcal{V}(L)\_{\prime}$$

where *U*(*L*) and *V*(*L*) are *r* × *r* matrix polynomials with all their roots outside the unit circle, and

$$\mathcal{M}(L) = \begin{pmatrix} (1-L)\mathbf{I}\_{r-c} & \mathbf{0} \\ \mathbf{0} & \mathbf{I}\_c \end{pmatrix}$$

<sup>5</sup> Multiplying both sides of (A3) by (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*) and using (A2), we obtain **A˜** (*L*)**S**(*L*)−1**B**(*L*)**u***<sup>t</sup>* = (<sup>1</sup> <sup>−</sup> *<sup>L</sup>*)**B ˜ ˜ <sup>u</sup>***t*. Comparing the spectral densities of right- and left-hand terms, it is easy to prove that **u˜***<sup>t</sup>* must be a *q*-dimensional, nonsingular white noise and the rank of **B˜** must be *q*.

(see Watson 1994). To get a VAR(2) we set *U*(*L*) = **I***<sup>r</sup>* − *U*1*L*, and *V*(*L*) = **I***r*, and then, by rewriting *M*(*L*) = **I***<sup>r</sup>* −*M*1*L*, we get **A**<sup>1</sup> = *M*<sup>1</sup> + *U*1, and **A**<sup>2</sup> = −*M*1*U*1.

Regarding the generation of the data, the diagonal entries of the matrix *U*<sup>1</sup> are drawn from a uniform distribution between 0.5 and 0.8, while the extra–diagonal entries are drawn from a uniform distribution between 0 and 0.3. *U*<sup>1</sup> is then multiplied by a scalar so that its largest eigenvalue is 0.6. The matrix **<sup>G</sup>** is generated as in Bai and Ng (2007): (1) **G˜** is an *<sup>r</sup>* <sup>×</sup> *<sup>r</sup>* diagonal matrix of rank *<sup>q</sup>* where *g*˜*ii* is drawn from the uniform distribution between 0.8 and 1.2, (2) **Gˇ** is obtained by orthogonalizing an *<sup>r</sup>* <sup>×</sup> *<sup>r</sup>* uniform random matrix, (3) **<sup>G</sup>** is equal to the first *<sup>q</sup>* columns of the matrix **Gˇ G˜** 1/2. Lastly, the orthogonal matrix **H** is such that the upper 3 × 3 submatrix of **GH** is lower triangular. The results are based on 1000 replications. The matrices *U*1, **G** and **H** are generated only once (the numerical values are available on request) so that the set of impulse responses to be estimated is the same for all replications, whereas the vector **u***<sup>t</sup>* is redrawn from N (**0**,**I**4) at each replication.

#### **References**


Canova, Fabio. 2007. *Methods for Applied Macroeconomics*. Princeton: Princeton University Press.


Van der Waerden, Bartel Leendert. 1953.*Modern Algebra*, 2nd ed. New York: Frederick Ungar, vol. I.

Watson, Mark W. 1994. Vector autoregressions and cointegration. In *Handbook of Econometrics*. Edited by Robert F. Engle and Daniel L. McFadden. Amsterdam: North Holland, Elsevier, vol. 4, chp. 47, pp. 2843–915.

© 2020 by Matteo Barigozzi and Marco Lippi. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
