Skip to main content

ORIGINAL RESEARCH article

Front. Appl. Math. Stat., 17 September 2020
Sec. Mathematical Finance
This article is part of the Research Topic Long-Memory Models in Mathematical Finance View all 5 articles

Modeling and Estimation of Multivariate Discrete and Continuous Time Stationary Processes

  • Department of Mathematics and Systems Analysis, Aalto University School of Science, Espoo, Finland

In this paper, we give an autoregressive model of order 1 type of characterization covering all multivariate strictly stationary processes indexed by the set of integers. Consequently, under square integrability, we derive continuous time algebraic Riccati equations for the parameter matrix of the characterization. This provides us with a natural way to define the corresponding estimator. In addition, we show that the estimator inherits consistency from autocovariances of the stationary process. Furthermore, the limiting distribution is given by a linear function of the limiting distribution of the autocovariances. We also present the corresponding existing results of the continuous time setting paralleling them to the discrete case treated in this paper.

1. Introduction

Stationary stochastic processes provide a significant instrument for modeling numerous temporal phenomena related to different fields of science. In particular, due to the evidence of long dependence structures in the real financial data, stationary processes possessing long-memory have been widely applied in mathematical finance.

When discrete time is considered, stationary data is typically modeled by applying autoregressive–moving-average (ARMA) processes or their extensions. One focal reason for popularity of ARMA processes is that for every stationary process with a vanishing autocovariance γ(·) and for every n ∈ ℕ there exists an ARMA process X such that γX(k) = γ(k) for |k| ≤ n. For a comprehensive overview of ARMA processes we mention [13]. The immense ARMA family include for example SARIMA (seasonal autoregressive integrated moving-average) models, where a seasonal ARMA process is obtained by differencing the original data. We mention also different GARCH (generalized autoregressive conditional heteroskedasticity) models originating from [4] and [5] that are commonly used in financial modeling taking account of the time-dependent volatility. ARMA processes, their extensions and estimation in these models have been concerned e.g., in [616], to name but a few. Moreover, in [17] we showed that all univariate strictly stationary processes indexed by the integers are characterized by the AR(1) (autoregressive model of order 1) equation

Xt-ϕXt-1=Zt,     t.

However, contrary to the classical AR(1), the noise Z belonging to a certain class of stationary processes is not necessarily white. Established on the characterization, we proposed an estimation method for ϕ in the case of a square integrable stationary process. This method has several advantages over conventional ones such as maximum likelihood and least squares fitting of ARMA models. Furthermore, in [18] we applied our method in estimation of a generalization of the ARCH model involving a covariate process that can be interpreted as the liquidity of an asset. We would like to emphasize that the proposed method is applicable in estimation of essentially any square integrable one-dimensional stationary process. Hence, it also covers stationary solutions, which are often of central interest, of models of the vast ARMA family.

In the case of continuous time, the Ornstein-Uhlenbeck process X given by the Langevin equation

dXt=-θXtdt+dBt,     t,    (1)

where θ > 0 and B is a two-sided Brownian motion, can be seen as the analog of the discrete time AR(1) process. By posing a suitable initial condition, (1) yields a stationary solution. The foregoing can be generalized, for example, by replacing Brownian motion with other stationary increment processes satisfying certain integrability conditions. One popular option as the replacement is fractional Brownian motion recovering the fractional Ornstein-Uhlenbeck process introduced in [19]. This kind of generalized Ornstein-Uhlenbeck processes are applied e.g., in mathematical finance to describe mean-reverting systems under the influence of shocks, and they are a highly active topic of research. Equations of type (1) with varying driving forces, and estimation in such models have been concerned e.g., in [2032], to mention but a few. Furthermore, in [33] we showed that a generalized version of (1) characterizes all multivariate strictly stationary processes with continuous paths. Consequently, we proposed an estimation method for the parameter matrix of (1). The method is based on continuous time algebraic matrix Riccati equations (CAREs) and it is applicable in estimation of essentially any square integrable multivariate stationary process.

Algebraic Riccati equations, occurring naturally e.g., in optimal control and filtering theory, is an intensively studied topic in the literature on its own. In many applications, real-valued CAREs often take the symmetric form

BA+AB-ACA+D=0,    (2)

where C and D are symmetric, and symmetric solutions A are to be found. For a general approach to algebraic Riccati equations the reader may consult for example [34]. The existence and uniqueness of a solution to (2) is a well-studied topic, especially when C and D are positive semidefinite [see e.g., [3537]].

The main contribution of this paper is to extend the characterization and the consequent estimation method of discrete stationary processes of [17] to multivariate settings. Serving as an instrument for the characterization, we also define a multivariate discrete Lamperti transform giving a one-to-one correspondence between stationary and H-self-similar processes with a positive definite matrix H. The characterization leads to covariance based symmetric CAREs for the related parameter matrix providing us with a novel estimation method of multivariate discrete stationary processes. In outline, this results in the following correspondence between the continuous time case. When the concepts of noises are extended from the conventional ones, the reputed analogous AR(1) and Langevin equations characterize discrete and continuous time stationary processes, respectively. The characterizations provide us with models of stationary processes yielding symmetric CAREs for the corresponding parameter matrices. Furthermore, these equations can be applied in a similar manner in estimation in both cases. For readers' interest and to highlight the obtained connection between discrete and continuous time, we also give the key results of [33].

The rest of the paper is organized as follows. In section 2.1, we first give an AR(1) type of characterization covering all multivariate strictly stationary processes indexed by the integers. Consequently, under the assumption of square integrability, we obtain a set of symmetric CAREs for the model parameter matrix serving as a basis for estimation. Finally, we state theorems for consistency and asymptotic distribution of the parameter (matrix) estimator. In section 2.2, we present the main results of [33], while at the same time comparing them to the results obtained in discrete time. For the reader's convenience, all technical proofs are postponed to section 3.

2. Main Results

We begin with some preliminaries and a short notational introduction. The processes we consider in this paper are n-dimensional, real-valued and indexed by I ∈ {ℤ, ℝ}. For such a process Y we write Y = (Yt)tI, where the ith component of the random vector Yt is denoted by Yt(i). Equality of the distributions of two random vectors Yt and Zt is denoted by Yt=lawZt. Similarly, equality of two processes Y and Z in the sense of finite dimensional distributions is denoted by Y=(Yt)tI=law(Zt)tI=Z. Throughout the paper, we investigate strictly stationary processes meaning that (Xt+s)tI=law(Xt)tI for every sI. Consequently, we omit the word “strictly” and simply say that X is stationary. By writing A ≥ 0 or A > 0 we mean that the matrix A is positive semidefinite or positive definite, respectively. We denote an eigendecomposition of a symmetric matrix by A = QΛQ, where Λ = diag(λi). Furthermore, the L2 vector norm and the corresponding induced matrix norm is denoted by ‖·‖.

In sections 2.1 and 2.2, we introduce models for discrete and continuous time stationary processes, respectively. Consequently, this leads to symmetric CAREs for the model parameter matrix of the form

BA+AB-ACA+D=0,    (3)

where C, D ≥ 0, and we are solving the equation for a positive definite A. There exists a vast amount of literature on existence and uniqueness of a solution [see e.g., [35] or [36]] in the described setting. In particular, if C, D > 0, then there exists a unique positive semidefinite solution to (3). Furthermore, there exists several numerical methods for finding the positive semidefinite solution of (3) [see e.g., [38, 39] or monograph [40]]. Hence, provided that the solution to (3) is unique, a prospective estimation method of A can be based on the equation.

2.1. Discrete Time

In this subsection, we extend the characterization of discrete (I = ℤ) stationary processes of [17] to multivariate settings. Consequently, we derive quadratic equations for the corresponding model parameter matrix providing us with a natural way to define an estimator for the parameter. Finally, we pose theorems for consistency and asymptotic distribution. A strong analog with the continuous time case I = ℝ covered in [33] is obtained. We start by providing some definitions.

DEFINITION 2.1. Let G = (Gt)t∈ℤ be an n-dimensional stationary increment process. We define a stationary process ΔG = (ΔtG)t∈ℤ by

ΔtG=Gt-Gt-1.

Next, we define a class of stationary increment processes having sub-exponentially deviating sample paths. These processes serve as the noise in the subsequent AR(1) type of characterization of stationary processes.

DEFINITION 2.2. Let H > 0 be a positive definite n × n matrix, and let G = (Gt)t∈ℤ be an n-dimensional stochastic process with stationary increments and G0 = 0. If

liml-k=l0ekHΔkG    (4)

exists in probability and defines an almost surely finite random variable, we denote GGH.

REMARK 2.3. Lemma 3.1 shows that existence of a logarithmic moment is sufficient for GGH for all H > 0. Particularly, this is the case if G is square integrable. On the other hand, an example of an one-dimensional stationary increment process G with G0 = 0, but GGH for any H > 0 was provided in [41]. Moreover, Definition 2.2 could also be stated without the assumption of positive definiteness. However, in the one-dimensional case with H ≤ 0, the convergence of (4) would imply G ≡ 0. As expected, we encounter a similar kind of dimensional degeneracy when considering e.g., symmetric matrices with non-positive eigenvalues. Hence, the assumption H > 0 may be regarded as natural. See also Remark 3.3.

The next theorem characterizes all multivariate stationary processes, including processes possessing long-memory.

THEOREM 2.4. Let H > 0 be a positive definite n × n matrix, and let X = (Xt)t∈ℤ be an n-dimensional stochastic process. Then X is stationary if and only if limt-etHXt=0 and

ΔtX=(e-H-I)Xt-1+ΔtG    (5)

for GGH and t ∈ ℤ. Moreover, the process GGH is unique.

COROLLARY 2.5. Let H > 0 be a positive definite n × n matrix, and let X be stationary. Then X admits an AR(1) type of representation

Xt-ΦXt-1=ΔtG,    (6)

where Φ = eH and GGH.

By using (6) and the expression (17) from the proof of Theorem 2.4, it is straightforward to show that ΔG is centered and square integrable if and only if X is centered and square integrable, respectively. In what follows, we assume these two attributes and write γ(t)=𝔼(XtX0) and r(t)=𝔼[(ΔtG)(Δ0G)]. Furthermore, since

Gt=k=1tΔkG,     t1,

in this case also G is centered and square integrable, and we denote v(t)=cov(Gt)=𝔼(GtGt). We would like to point out that centeredness can be assumed without loss of generality (see Remark 2.11).

Under the discussed assumptions, we obtain an expression for γ(t) in terms of the noise process.

REMARK 2.6. The autocovariance function γ(t) is given by

γ(t)=e-tHk=-tj=-0ekHr(k-j)ejH.

Furthermore, if G has independent components, we obtain

γ(t)=etH2k=tj=0ekH(v(kj+1)+v(kj1)                 2v(kj))ejH.

The following lemma writes the quadratic equations for the model parameter Φ = eH presented in [17] in our multivariate setting.

LEMMA 2.7. Let H > 0 be a positive definite n × n matrix, and let X be stationary of the form (6). Then

r(t)=Φγ(t)Φ-γ(t+1)Φ-Φγ(t-1)+γ(t)    (7)

for every t ∈ ℤ.

REMARK 2.8. In the proof of the lemma, we utilize the increment process ΔG similarly as in [17] yielding quadratic equations for Φ in terms of r and γ. One could consider an estimation method for Φ based on the above equations. However, for a general stationary X, (7) is a symmetric CARE only if t = 0, and even in this case, existence of a unique positive semidefinite solution is not guaranteed.

By adapting the approach of [33], we obtain a set of symmetric CAREs on which we construct an estimator for the model parameter Φ = eH. For this, we define the following matrix coefficients.

DEFINITION 2.9. We set

Bt=k=1tγ(k-1)-γ(k)Ct=k=1tj=1tγ(k-j)Dt=v(t)-2γ(0)+γ(t)+γ(t)

for every t ∈ ℕ.

THEOREM 2.10. Let H > 0 be a positive definite n × n matrix, and set Θ = IeH. Let X = (Xt)t∈ℤ be stationary of the form (5). Then the CARE

BtΘ+ΘBt-ΘCtΘ+Dt=0    (8)

is satisfied for every t ∈ ℕ.

REMARK 2.11. Equations (7) and (8) are covariance based. Consequently, they hold also when X and G in Theorem 2.4 are not centered.

REMARK 2.12. It is worth to emphasize the interplay between the parameter matrix H and the noise G in the context of a fixed stationary X. The noise in Theorem 2.4 and Corollary 2.5 is unique only after H > 0 is fixed. On the other hand, if we fix some covariances of the noise process such that the equations of Theorem 2.10 (or Lemma 2.7) yield a positive definite solution, then we may set H according to this solution. Furthermore, the corresponding noise is now given by Corollary 2.5 and its covariance satisfies the set prerequisites. In practice, the parameter matrix H is estimated by assuming some information on the noise, and by estimating the autocovariance of the observed stationary process X. The estimation will discussed in more detail at the end of this subsection. See also the following examples illustrating how to obtain familiar noise terms in the case of two stationary ARMA type of processes.

We give a couple of examples on how some basic multivariate processes of ARMA type can be presented in the form (6), and how to derive the corresponding noise G together with its covariance function v.

EXAMPLE 2.13. Let X be an n-dimensional stationary AR(1) type of process given by

Xt-ϕXt-1=ϵt,

with 0 < ϕ = QΛQ, ‖ϕ‖ < 1 and ϵ ~ iid(0, Σ). Then, we may set H=-Qdiag(logλi)Q giving Φ = ϕ. Now ΔG = ϵ and Gt=k=1tϵk. Furthermore, v(t)=k=1tcov(ϵk)=tΣ for t ≥ 1.

EXAMPLE 2.14. Let X be an n-dimensional stationary ARMA(1, q) type of process given by

Xt-ϕXt-1=ϵt+θ1ϵt-1++θqϵt-q,

with 0 < ϕ, ‖ϕ‖ < 1 and ϵ ~ iid(0, Σ). Similarly as above, we may set Φ = ϕ and now ΔG equals to the MA(q) process on the right. Consequently, for t ≥ 1,

Gt=k=1tϵk+θ1ϵk-1++θqϵk-q

and

v(t)=i,j=0qmax(0,t-|i-j|)θiΣθj,

where θ0 = I.

In [17], we proposed an estimation method of one-dimensional stationary processes based on equations (7). In particular, we showed that the method is applicable except in some special class of stationary processes. In [42], we provided a comprehensive analysis of the class, and proved that it consists of highly degenerate processes. On the other hand, due to the strong dependence structure, the failure of different estimation methods is expected. Fundamentally, a stationary process X belongs to the class if there exists two values H and H~ such that the corresponding processes ΔG and ΔG~ in (6) have identical autocovariance functions. Next, we state a lemma showing that these degenerate processes have a special characteristic also under the new set of Equations (8).

LEMMA 2.15. Let X be a one-dimensional stationary process and let H > 0 be fixed. Set Φ = eH and

Xt-ΦXt-1=ΔtG,     t.

If the equation

γ(t)Φ2-(γ(t+1)-γ(t-1))Φ+γ(t)-r(t)=0    (9)

yields the same two solutions Φ,Φ~>0 for every t ∈ ℤ, then also the equation

CtΘ2-2BtΘ-Dt=0    (10)

yields the same two solutions Θ = 1 − Φ and Θ~=1-Φ~ for every t ∈ ℕ.

The following remark illuminates the connection between the coefficient matrices of Definition 2.9 and uniqueness of the solution to (8), which we discussed in the beginning of section 2.

REMARK 2.16. By Lemma 3.8, the matrix Θ is positive definite. Since

Ct=𝔼[k=1tXk-1(k=1tXk-1)]=cov(k=1tXk-1),

the matrix Ct is positive semidefinite. Furthermore, if the smallest eigenvalue of v(t) grows enough in time, Dt becomes positive definite [see [33]]. This is the case e.g., when the noise G has independent components with growing variances.

For estimation, it is desirable that (8) admits a unique positive semidefinite solution. In this case, the solution is the correct parameter matrix Θ by the construction. Moreover, it guarantees that a convergent numerical scheme has the desired limit Θ. In the sequal, we simply assume that t is chosen in such a way that Ct, Dt > 0 ensuring the uniqueness of the solution, and we omit the subindex t from Equation (8). We have justified the assumption of positive definiteness in detail in continuous time [see section 2.1 and Remark 2.11 in [33]]. Furthermore, we assume that v(t) is known and the stationary process X is observed up to the time T > t, and the coefficient matrices B, C, and D are estimated from these observations by replacing the autocovariances γ(·) with some estimators γ^T(·). The coefficient estimators are denoted by B^T,C^T, and D^T, and we set

ΔTB=B^T-B,     ΔTC=C^T-C,     ΔTD=D^T-D.

Next, we define an estimator Θ^T for the matrix Θ = I − Φ = IeH. The proofs of the related asymptotic results allow a certain amount of flexibility in the definition. Thus, we give a definition that probably is the most convenient from the practical point of view. Consistency and the rate of convergence of Θ^T are inherited from autocovariance estimators γ^T(·) of the observed stationary process. In addition, the limiting distribution is obtained as a linear function of the limiting distribution of the autocovariance estimators.

DEFINITION 2.17. The estimator Θ^T is defined as the unique positive semidefinite solution to the perturbed CARE

B^TΘ^T+Θ^TB^T-Θ^TC^TΘ^T+D^T=0

whenever C^T,D^T>0. Otherwise, we set Θ^T=0.

THEOREM 2.18. Let C, D > 0. Assume that

maxs{0,1,,t}γ^T(s)-γ(s)0.

Then

Θ^T-Θ0,

where Θ^T is given by Definition 2.17.

THEOREM 2.19. Let l(T) be a rate function. If

l(T)[vec(γ^T(0)-γ(0))vec(γ^T(1)-γ(1))vec(γ^T(t)-γ(t)) ]lawZ,

where Z is a (t + 1)n2-dimensional random vector, then:

(1) Let Z~ be the permutation of elements of Z corresponding to the order of elements of

[vec((γ^T(0)-γ(0)))vec((γ^T(1)-γ(1)))vec((γ^T(t)-γ(t))).]

Define a linear mapping L1:(t+1)n23n2 by

L1(Z)=[k=0t-1(t-k)[Z(kn2+1)Z((k+1)n2)]+k=1t-1(t-k)[Z~(kn2+1)Z~((k+1)n2)]k=1t[Z((k-1)n2+1)-Z~(kn2+1)Z(kn2)-Z~((k+1)n2)][Z(tn2+1)+Z~(tn2+1)Z((t+1)n2)+Z~((t+1)n2)]-2[Z(1)Z(n2)]],

where 10 is an empty sum. Then

l(T)vec(ΔTC,ΔTB,ΔTD)lawL1(Z).

(2) If D, C > 0 and Θ^T is given by Definition 2.17, then

l(T)vec(Θ^T-Θ)lawL2(L1(Z)),

where L2:3n2n2 is a linear mapping expressible in terms of Θ, t and r.

2.2. Continuous Time

We have collected the main results (Theorems 2.21, 2.28, 2.32 and 2.33) of [33] concerning continuous time multivariate stationary processes into this subsection. In addition, in order to complete the analog between discrete and continuous time, we derive quadratic equations for the model parameter (Proposition 2.24) that correspond to the equations of Lemma 2.7. Throughout the subsection, we assume that the considered processes have continuous paths almost surely and hence, the related stochastic integrals can be interpreted as pathwise Riemann-Stieltjes integrals. Again, we start by defining the class GH of stationary increment processes for H > 0.

DEFINITION 2.20. Let H > 0 be a positive definite n × n matrix, and let G = (Gt)t∈ℝ be an n-dimensional stochastic process with stationary increments and G0 = 0. If

lims-s0eHudGu

exists in probability and defines an almost surely finite random variable, we denote GGH.

As in discrete time, it can be shown that existence of some logarithmic moments ensure that GGH for all H > 0. In particular, square integrability of G suffices, which is the case in our second moment based estimation method.

The next theorem is the continuous time counterpart of Theorem 2.4 showing that all stationary processes are characterized by the Langevin equation, whereas in discrete time, the characterization was given by an AR(1) type of equation.

THEOREM 2.21. Let H > 0 be a positive definite n × n matrix, and let X = (Xt)t∈ℝ be an n-dimensional stochastic process. Then X is stationary if and only if

X0=-0eHudGu

and

dXt=-HXtdt+dGt,         (11)

for GGH and t ∈ ℝ. Moreover, the process GGH is unique.

COROLLARY 2.22. From Theorem 2.21 it follows that X is the unique stationary solution

Xt=e-Ht-teHudGu    (12)

to (11).

In order to apply Theorem 2.21 in estimation, we pose the assumption

sups[0,1]𝔼(Gs2)<.

This guarantees that GGH for all H > 0, and square integrability of X and G. On the other hand, if X is square integrable, then G is also. In addition and without loss of generality, we assume that the processes are centered. Again, we write γ(t)=𝔼(XtX0) and v(t)=𝔼(GtGt). Now, the autocovariance function of the following stationary process is well-defined.

DEFINITION 2.23. Let G = (Gt)t∈ℝ be a centered square integrable stationary increment process and let δ > 0. We define a stationary process ΔδG=(ΔtδG)t by

ΔtδG=Gt-Gt-δ

and the corresponding autocovariance function rδ by

rδ(t)=𝔼[(ΔtδG)(Δ0δG)].

As in discrete time (Lemma 2.7), we obtain quadratic equations for the model parameter H in terms of rδ and γ. The equations could potentially be used to construct an estimator for H and they might be also of independent interest.

PROPOSITION 2.24. Let H > 0 be a positive definite n × n matrix, and let X be of the form (12). Then

rδ(t)=2γ(t)γ(t+δ)γ(tδ)+(tt+δγ(s)dstδtγ(s)ds)H+H(tδtγ(s)dstt+δγ(s)ds)+H(tδt(st+δ)γ(s)ds                                                                    +tt+δ(ts+δ)γ(s)ds)H

for every t ∈ ℝ.

REMARK 2.25. The advantage of the equations above is that we have to consider γ(s) only for s ∈ [t − δ, t + δ], but as in the discrete case, for a general stationary X we obtain a symmetric CARE only when t = 0. In addition, similarly as above, we could set in discrete time ΔtkG: =Gt-Gt-k,k. However, this would lead to more complicated equations in Lemma 2.7.

A significant difference compared to the discrete time Equations (7) occurs in the univariate case. Namely, the first order term with respect to H vanishes.

COROLLARY 2.26. The univariate case yields

rδ(t)=2γ(t)γ(t+δ)γ(tδ)+H2(tδt(st+δ)γ(s)ds                  +tt+δ(ts+δ)γ(s)ds)

for every t ∈ ℝ.

One could base a univariate estimation method on the above equations without the concern of existence of a unique positive solution. However, since we wish to treat also multivariate settings, we present the most central results of [33] that are obtained from Theorem 2.21 by considering the noise G directly. First, we define matrix coefficients corresponding to Definition 2.9. Consequently, we write symmetric CAREs for the parameter H that are similar to the CAREs (8) for the discrete time parameter Θ.

DEFINITION 2.27. We set

Bt=0tγ(s)-γ(s)dsCt=0t0tγ(s-u)dudsDt=v(t)-2γ(0)+γ(t)+γ(t)

for every t ≥ 0.

THEOREM 2.28. Let H > 0 be a positive definite n × n matrix, and let X = (Xt)t∈ℝ be stationary of the form (12). Then the CARE

BtH+HBt-HCtH+Dt=0    (13)

is satisfied for every t ≥ 0.

REMARK 2.29. As in discrete time, the Equations (13) and (14) are covariance based and hence, they hold also for a non-centered stationary X.

REMARK 2.30. Contrary to the discrete time Equations (8), the first order term vanishes in the univariate setting as in (13).

Again, we assume that t is chosen in such a way that Ct, Dt > 0 ensuring the existence of a unique positive semidefinite solution. We have discussed this assumption in detail in [33]. We define an estimator H^T for the model parameter matrix H identically to the discrete time by replacing the autocovariances γ(·) in the matrix coefficients with their estimators γ^T(·). The below given definition differs slightly from the definition in [33], but the same asymptotic results still apply.

DEFINITION 2.31. The estimator H^T is defined as the unique positive semidefinite solution to the perturbed CARE

B^TH^T+H^TB^T-H^TC^TH^T+D^T=0

whenever C^T,D^T>0. Otherwise, we set H^T = 0.

As in discrete time, asymptotic properties of H^T are inherited from the autocovariance estimators. However, due to the continuous time setting, instead of pointwise convergence, we have to consider functional form of convergence of γ^T(·). In [33], we have provided sufficient conditions in the case of Gaussian noise G with independent components, under which the assumptions of the following theorems are satisfied. In particular, the results are valid for fractional Brownian motion that is widely applied in the field of mathematical finance.

THEOREM 2.32. Let C, D > 0. Assume that

sups[0,t]γ^T(s)-γ(s)0.

Then

HT^-H0,

where H^T is given by Definition 2.31.

THEOREM 2.33. Let Y = (Ys)s∈[0, t] be an n2-dimensional stochastic process with continuous paths almost surely and let l(T) be a rate function. If

l(T)vec(γ^T(s)-γ(s))lawYs

in the uniform topology of continuous functions, then:

(1) Let Y~s be the permutation of elements of Ys that corresponds to the order of elements of vec((γ^(s)γ(s))). Then

l(T)vec(ΔTC,ΔTB,ΔTD)law[0t(t-s)(Ys+Y~s)ds0t(Ys-Y~s)dsYt+Y~t-2Y0]                                                                                             =:L1(Y).

(2) If C, D > 0 and H^T is given by Definition 2.31, then

l(T)vec(H^T-H)lawL2(L1(Y)),

where L2:3n2n2 is a linear mapping expressible in terms of H, t and the covariance function of G.

3. Proofs

In the following, we denote the smallest eigenvalue of H > 0 by λmin. Consequently ekH=Qdiag(eλik)Q=eλmink for a negative k.

3.1. Discrete Time

The proof of the next lemma follows the lines of the proof of Theorem 2.2. in [41] that concerns the one-dimensional continuous time case. However, in our setting, we obtain a weaker sufficient condition for GGH for all H > 0.

LEMMA 3.1. Let G = (Gt)t∈ℤ be an n-dimensional stationary increment process with G0 = 0. Assume that

𝔼(|logG1𝟙{G1>1}|1+δ)<

for some δ > 0. Then GGH for all positive definite n × n matrices H.

PROOF. Let H > 0. We apply the Borel-Cantelli lemma together with Markov's inequality to show that ekHΔkG0 almost surely as k → −∞. Let ϵ > 0 be fixed below.

(ekHΔkG>ϵ)(eλminkΔkG>ϵ)=(eλminkG1>ϵ)=(G1>ϵeλmink)                                  =(logG1>log ϵ-λmink),

since ΔG is stationary and G0 = 0. Furthermore, k(log ϵk-λmin)-Ck for some C > 0 and kkϵ. Thus, for kkϵ,

(ekHΔkG>ϵ)(logG1-Ck)                  =(logG1𝟙{G1>1}-Ck)𝔼|logG1𝟙{G1>1}|1+δ(-Ck)1+δc1(-k)1+δ

giving the wanted result. We conclude the proof by noting that

k=-0ekHΔkGk=-0e12kHe12kHΔkGsupke12kHΔkGk=-0e12λmink<

almost surely.

Next, we extend the concept of self-similarity to discrete time multivariate processes. While the following two definitions are natural, to the best of our knowledge they are not widely acknowledged in the literature.

DEFINITION 3.2. Let H > 0 be a positive definite n × n matrix, and let Y=(Yet)t be an n-dimensional stochastic process. Then Y is H-self-similar if

(Yet+s)t=law(esHYet)t

for every s ∈ ℤ.

REMARK 3.3. Similarly as in Definition 2.2, positive definite matrices serve as natural counterparts of the conventional positive scalar valued exponents of self-similarity. In addition, liml-Yel does not converge (see the proof of the auxiliary Lemma 3.7) in the case of a general H-self-similar Y with non-positive definite H. It seems intuitive to expect the term to be convergent. However, posing such requirement leads to degeneracy of Y e.g., when H is symmetric with negative eigenvalues [cf. Remark 3.2 in [33]]. Furthermore, this indicates a possibility of dimensional reduction. For details on continuous time self-similar processes, we refer to [43]. See also the discussion on non-positive exponents of self-similarity in [44].

The following transform and the corresponding theorem giving one-to-one correspondence between self-similar and stationary processes were originally introduced by Lamperti in the univariate continuous time setting [45].

DEFINITION 3.4. Let H > 0 be a positive definite n × n matrix, and let X = (Xt)t∈ℤ and Y=(Yet)t be n-dimensional stochastic processes. We define

(LHX)et=etHXt

and

(LH-1Y)t=e-tHYet.

THEOREM 3.5. The operator LH together with its inverse LH-1 define a bijection between n-dimensional stationary processes and n-dimensional H-self-similar processes.

PROOF. First, let X be stationary and set Zet=(LHX)et. Then

[Zet1+sZet2+sZetm+s]=[e(t1+s)HXt1+se(t2+s)HXt2+se(tm+s)HXtm+s]=law[e(t1+s)HXt1e(t2+s)HXt2e(tm+s)HXtm]=[esHZet1esHZet1esHZetm]

for every m ∈ ℕ, t ∈ ℤm and s ∈ ℤ. Hence, Z is H-self-similar.

Now, let Y be H-self-similar and set Zt=(LH-1Y)t. Then

[Zt1+sZt2+sZtm+s]=[e-(t1+s)HYet1+se-(t2+s)HYet2+se-(tm+s)HYetm+s]=law[e-t1HYet1e-t2HYet2e-tmHYetm]=[Zt1Zt2Ztm]

for every m ∈ ℕ, t ∈ ℤm and s ∈ ℤ. Hence, Z is stationary completing the proof.

REMARK 3.6. The form of the conventional continuous time Lamperti transform is not directly applicable in the discrete time setting due to the scaling of time. Hence, in order to stay within the given discrete parameter sets, we use exponential clocks with self-similar processes in Definitions 3.2 and 3.4.

Before the proof of Theorem 2.4 we state an auxiliary lemma.

LEMMA 3.7. Let H > 0 be a positive definite n × n matrix, and let (Yet)t be an n-dimensional H-self-similar process. We define a process G = (Gt)t∈ℤ by

Gt={k=1te-kHΔkYek,     t10,     t=0-k=t+10e-kHΔkYek,     t-1.

Then GGH.

PROOF. It is straightforward to verify that ΔtG=e-tHΔtYet for every t ∈ ℤ. In addition

liml-k=l0ekHΔkG=liml-k=l0ΔkYek=Y1-liml-Yel-1,

where by self-similarity of Y

(Yel-1ϵ)=(eH(l-1)Y1ϵ)(eλmin(l-1)Y1ϵ)                                      0.

Hence, we set

liml-k=l0ekHΔkG=Y1.

PROOF OF THEOREM 2.4. Assume that limt-etHXt=0 and (5) holds for GGH. Then, by using (5) repeatedly

Xt=e-HXt-1+ΔtG=e-(n+1)HXt-n-1+j=0ne-jHΔt-jG                   =e-(n+1)HXt-n-1+e-tHk=t-ntekHΔkG=e-tH                                            (e(t-n-1)HXt-n-1+k=t-ntekHΔkG)

for every n ∈ ℕ. Since, as n → ∞, the limit of the sum above is well-defined, and limne(t-n-1)HXt-n-1=0, we obtain that

Xt=e-tHk=-tekHΔkG.    (14)

Let m ∈ ℕ, t ∈ ℤm, and s ∈ ℤ. Then, by stationary increments of G, we have

[e-(t1+s)Hj=-Mt1e(j+s)HΔj+sGe-(tm+s)Hj=-Mtme(j+s)HΔj+sG]=law[e-t1Hj=-Mt1ejHΔjGe-tmHj=-MtmejHΔjG]

for every −M < min{ti}. Since the random vectors above converge in probability as M → ∞, we obtain that

[Xt1+sXtm+s]=[e-(t1+s)Hj=-t1e(j+s)HΔj+sGe-(tm+s)Hj=-tme(j+s)HΔj+sG]              =law[e-t1Hj=-t1ejHΔjGe-tmHj=-tmejHΔjG]=[Xt1Xtm]

and hence, X is stationary.

Next, assume that X is stationary. Then, by Theorem 3.5 there exists a H-self-similar Y such that

ΔtX=e-tHYet-e-(t-1)HYet-1=(e-H-I)Xt-1+e-tHΔtYet.

Defining G as in Lemma 3.7 completes the proof of the other direction.

To prove uniqueness, we use (15). Assume that, for G,G~GH,

etHXt=k=-tekHΔkG=k=-tekHΔkG~

for every t ∈ ℤ. Then

etHXt-e(t-1)HXt-1=etHΔtG=etHΔtG~.

Since etH is invertible and both processes start from zero, we conclude that G=G~.

PROOF OF LEMMA 2.7. We have that

ΔtG(Δ0G)=(Xt-ΦXt-1)(X0-X-1Φ).

Taking expectations yields

r(t)=Φγ(t)Φ-γ(t+1)Φ-Φγ(t-1)+γ(t).

PROOF OF THEOREM 2.10. Let

ΔtX=-ΘXt-1+ΔtG.

Then for t ∈ ℕ we have

Gt=k=1tΔkG=k=1tΔkX+Θk=1tXk-1=Xt-X0+Θk=1tXk-1.

Hence

        cov(Gt)=cov(Xt-X0)+𝔼[(Xt-X0)k=1tXk-1]Θ+Θ𝔼[k=1tXk-1(Xt-X0)]+Θ𝔼[k=1tXk-1(k=1tXk-1)]Θ

giving (8) since

𝔼[(Xt-X0)k=1tXk-1]=k=1tγ(t-k+1)-γ(-k+1)                                                   =k=1tγ(k)-γ(k-1)

and

cov(Xt-X0)=𝔼[(Xt-X0)(Xt-X0)]=2γ(0)-γ(t)                            -γ(-t)=2γ(0)-γ(t)-γ(t)

PROOF OF LEMMA 2.15. Assume that Φ~=e-H~ satisfies (9) for every t ∈ ℤ and set

Xt-Φ~Xt-1=ΔtG~,     t.

Consequently

r~(t)=γ(t)Φ~2-(γ(t+1)-γ(t-1))Φ~+γ(t)=r(t),     t,

where r~(t) is the autocovariance function of (ΔtG~)t. Now, since G0=G~0=0, we obtain that

var(Gt)=var(k=1tΔkG)=k,j=1tcov(ΔkG,ΔjG)=k,j=1tr(k-j)               =var(k=1tΔkG~)=var(G~t)

for all t ∈ ℕ. Hence, both Θ and Θ~ are solutions to (10).

LEMMA 3.8. The matrix Θ = IeH is positive definite.

PROOF. Let a be a real vector of length n, and let H = QΛQ be an eigendecomposition of H. Then

a(I-e-H)a=a2-ae-Ha,

where

|ae-Ha|a2e-λmin<a2

completing the proof.

In order to show that Θ^T is consistent, we simply need to find suitable bounds for ΔTB, ΔTC and ΔTD in terms of the autocovariance estimators. After that, the same strategy as in [33] can be applied.

LEMMA 3.9. Set

Mt,T=maxs{0,1,,t}γ^T(s)-γ(s).

Then the coefficients of the perturbed CARE satisfy

ΔTD4Mt,TΔTCt2Mt,TΔTB2tMt,T,

PROOF. First, we recall first that

γ^T(-s)-γ(-s)=γ^T(s)-γ(s)=γ^T(s)-γ(s).

Now, since v(t) is known,

ΔTD2γ^T(0)-γ(0)+γ^T(t)-γ(t)+γ^T(t)                                                                                     -γ(t)4Mt,T.

Moreover

ΔTCk=1tj=1tγ^T(k-j)-γ(k-j)t2Mt,T.

Finally

ΔTBk=1tγ^T(k-1)-γ(k-1)+γ(k)-γ^T(k)                                                                                                         2tMt,T.

PROOF OF THEOREM 2.18. The result follows by replacing sups[0,t]γ^T(s)-γ(s) with Mt,T in Corollary 3.14 and in the proof of Theorem 2.9 of [33]. The details are left to the reader.

PROOF OF THEOREM 2.19. For the first part of the theorem, we notice that

Ct=k=1tj=1tγ(k-j)=k=1tl=k-tk-1γ(l)=l=1-t-1k=1t+lγ(l)      +        l=0t-1k=l+1tγ(l)=1-t-1(t+l)γ(l)+l=0t-1(t-l)γ(l)      =                                           l=0t-1(t-l)γ(l)+l=1t-1(t-l)γ(l),

where 0-1 and 10 are interpreted as empty sums. Now we have that

ΔTC=k=0t-1(t-k)(γ^T(k)-γ(k))           +k=1t-1(t-k)(γ^T(k)-γ(k))ΔTB=k=1tγ^T(k-1)-γ(k-1)-γ^T(k)+γ(k)ΔTD=2(γ(0)-γ^T(0))+γ^T(t)-γ(t)+γ^T(t)-γ(t)

and furthermore

                                                                                         l(T)vec(ΔTC,ΔTB,ΔTD)=l(T)[     k=0t1(tk)vec(γ^(k)γ(k))+k=1t1(tk)vec((γ^(k)γ(k)))k=1tvec(γ^(k1)γ(k1))vec((γ^(k)γ(k)))2vec(γ^(0)γ(0))+vec(γ^(t)γ(t))+vec((γ^(t)γ(t)))]                                                             =L1[l(T)[vec(γ^T(0)γ(0))vec(γ^T(1)γ(1))vec(γ^T(t)γ(t))]]lawL1(Z)

by the continuous mapping theorem. For the second part of the theorem, the proof of the continuous time case of [33] can be applied just by replacing sups[0,t]γ^T(s)-γ(s) with Mt,T in the definition of the set AT.

3.2. Continuous Time

We only provide the proof of Proposition 2.24, while the other proofs can be found from [33]

PROOF OF PROPOSITION 2.24. Integrating (11) from 0 to t gives

Gt=Xt-X0+H0tXsds.

Hence

ΔtδG(Δ0δG)=(Xt-Xt-δ+Ht-δtXsds)                                  (X0-X-δ+-δ0XsdsH).

Taking expectations yields

rδ(t)=2γ(t)-γ(t+δ)-γ(t-δ)+-δ0γ(t-s)-γ(t-δ-s)dsH+Ht-δtγ(s)-γ(s+δ)ds+Ht-δt-δ0γ(s-u)dudsH,

where the first order terms can be treated with a simple change of variables. For the second order term we obtain that

                 t-δt-δ0γ(s-u)duds=t-δtss+δγ(x)dxds                 =t-δtt-δxγ(x)dsdx+tt+δx-δtγ(x)dsdx=t-δt(x-t+δ)γ(x)dx+tt+δ(t-x+δ)γ(x)dx.

Author Contributions

The author confirms being the sole contributor of this work and has approved it for publication.

Conflict of Interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

I would like to thank Lauri Viitasaari for his comments and suggestions.

References

1. Brockwell PJ, Davis RA. Time Series: Theory and Methods. New York, NY: Springer Science & Business Media (1991). doi: 10.1007/978-1-4419-0320-4

CrossRef Full Text | Google Scholar

2. Hamilton JD. Time Series Analysis. Vol. 2. Princeton: Princeton University Press (1994).

3. Neusser K. Time Series Econometrics. Springer (2016). doi: 10.1007/978-3-319-32862-1

CrossRef Full Text | Google Scholar

4. Engle RF. Autoregressive conditional heteroskedasticity with estimates of the variance of United Kingdom inflation. Econometr J Econometr Soc. (1982) 50:987–1007. doi: 10.2307/1912773

CrossRef Full Text | Google Scholar

5. Bollerslev T. Generalized autoregressive conditional heteroskedasticity. J. Econometr. (1986) 31:307–27. doi: 10.1016/0304-4076(86)90063-1

CrossRef Full Text | Google Scholar

6. Hannan EJ. The asymptotic theory of linear time-series models. J Appl Probabil. (1973) 10:130–45. doi: 10.1017/S0021900200042145

CrossRef Full Text | Google Scholar

7. Tiao GC, Tsay RS. Consistency properties of least squares estimates of autoregressive parameters in ARMA models. Ann Stat. (1983) 11:856–71. doi: 10.1214/aos/1176346252

CrossRef Full Text | Google Scholar

8. Mikosch T, Gadrich T, Kluppelberg C, Adler RJ. Parameter estimation for ARMA models with infinite variance innovations. Ann Stat. (1995) 23:305–26. doi: 10.1214/aos/1176324469

CrossRef Full Text | Google Scholar

9. Mauricio JA. Exact maximum likelihood estimation of stationary vector ARMA models. J Am Stat Assoc. (1995) 90:282–91. doi: 10.1080/01621459.1995.10476511

CrossRef Full Text | Google Scholar

10. Ling S, McAleer M. Asymptotic theory for a vector ARMA-GARCH model. Econometr Theory. (2003) 19:280–310. doi: 10.1017/S0266466603192092

CrossRef Full Text | Google Scholar

11. Francq C, Zakoian JM. Maximum likelihood estimation of pure GARCH and ARMA-GARCH processes. Bernoulli. (2004) 10:605–37. doi: 10.3150/bj/1093265632

CrossRef Full Text | Google Scholar

12. Bollerslev T. Glossary to ARCH (GARCH). CREATES Res Pap. (2008) 49:1–44. doi: 10.2139/ssrn.1263250

CrossRef Full Text | Google Scholar

13. Han H, Kristensen D. Asymptotic theory for the QMLE in GARCH-X models with stationary and nonstationary covariates. J Bus Econ Stat. (2014) 32:416–29. doi: 10.1080/07350015.2014.897954

CrossRef Full Text | Google Scholar

14. Baillie RT, Bollerslev T, Mikkelsen HO. Fractionally integrated generalized autoregressive conditional heteroskedasticity. J Econometr. (1996) 74:3–10. doi: 10.1016/S0304-4076(95)01749-6

CrossRef Full Text | Google Scholar

15. Ling S, Li WK. On fractionally integrated autoregressive moving-average time series models with conditional heteroscedasticity. J Am Stat Assoc. (1997) 92:1184–94. doi: 10.1080/01621459.1997.10474076

CrossRef Full Text | Google Scholar

16. Zhang MY, Russell JR, Tsay RS. A nonlinear autoregressive conditional duration model with applications to financial transaction data. J Econometr. (2001) 104:179–207. doi: 10.1016/S0304-4076(01)00063-X

CrossRef Full Text | Google Scholar

17. Voutilainen M, Viitasaari L, Ilmonen P. On model fitting and estimation of strictly stationary processes. Modern Stochast Theory Appl. (2017) 4:381–406. doi: 10.15559/17-VMSTA91

CrossRef Full Text | Google Scholar

18. Voutilainen M, Viitasaari L, Ilmonen P, Torres S, Tudor C. On the ARCH model with stationary liquidity. (2020). Available online at: https://link.springer.com/content/pdf/10.1007/s00184-020-00779-x.pdf

Google Scholar

19. Cheridito P, Kawaguchi H, Maejima M. Fractional Ornstein-Uhlenbeck processes. Electron J Probabil. (2003) 8:1–44. doi: 10.1214/EJP.v8-125

CrossRef Full Text | Google Scholar

20. Hu Y, Nualart D. Parameter estimation for fractional Ornstein-Uhlenbeck processes. Stat Probabil Lett. (2010) 80:1030–8. doi: 10.1016/j.spl.2010.02.018

CrossRef Full Text | Google Scholar

21. Kleptsyna ML, Le Breton A. Statistical analysis of the fractional Ornstein-Uhlenbeck type process. Stat Inference Stochast Process. (2002) 5:229–48. doi: 10.1023/A:1021220818545

CrossRef Full Text | Google Scholar

22. Azmoodeh E, Viitasaari L. Parameter estimation based on discrete observations of fractional Ornstein-Uhlenbeck process of the second kind. Stat Inference Stochast Process. (2015) 18:205–27. doi: 10.1007/s11203-014-9111-8

CrossRef Full Text | Google Scholar

23. Bajja S, Es-Sebaiy K, Viitasaari L. Least squares estimator of fractional Ornstein-Uhlenbeck processes with periodic mean. J Korean Stat Soc. (2017) 46:608–22. doi: 10.1016/j.jkss.2017.06.002

CrossRef Full Text | Google Scholar

24. Balde MF, Es-Sebaiy K, Tudor CA. Ergodicity and drift parameter estimation for infinite-dimensional fractional Ornstein-Uhlenbeck process of the second kind. Appl Math Opt. (2020) 81:785–814. doi: 10.1007/s00245-018-9519-4

CrossRef Full Text | Google Scholar

25. Brouste A, Iacus SM. Parameter estimation for the discretely observed fractional Ornstein-Uhlenbeck process and the Yuima R package. Comput Stat. (2013) 28:1529–47. doi: 10.1007/s00180-012-0365-6

CrossRef Full Text | Google Scholar

26. Douissi S, Es-Sebaiy K, Tudor CA. Hermite Ornstein-Uhlenbeck processes mixed with a Gamma distribution. Publicationes Mathematicae Debrecen. (2020) 96:1–22. doi: 10.5486/PMD.2020.8443

CrossRef Full Text | Google Scholar

27. Hu Y, Nualart D, Zhou H. Parameter estimation for fractional Ornstein-Uhlenbeck processes of general Hurst parameter. Stat Inference Stochast Process. (2019) 22:111–42. doi: 10.1007/s11203-017-9168-2

CrossRef Full Text | Google Scholar

28. Nourdin I, Tran TD. Statistical inference for Vasicek-type model driven by Hermite processes. Stochast Process Appl. (2019) 129:3774–91. doi: 10.1016/j.spa.2018.10.005

CrossRef Full Text | Google Scholar

29. Sottinen T, Viitasaari L. Parameter estimation for the Langevin equation with stationary-increment Gaussian noise. Stat Inference Stochast Process. (2018) 21:569–601. doi: 10.1007/s11203-017-9156-6

CrossRef Full Text | Google Scholar

30. Tanaka K. Maximum likelihood estimation for the non-ergodic fractional Ornstein-Uhlenbeck process. Stat Inference Stochast Process. (2015) 18:315–32. doi: 10.1007/s11203-014-9110-9

CrossRef Full Text | Google Scholar

31. Applebaum D, et al. Infinite dimensional Ornstein-Uhlenbeck processes driven by Lévy processes. Probabil Surv. (2015) 12:33–54. doi: 10.1214/14-PS249

CrossRef Full Text | Google Scholar

32. Magdziarz M. Fractional Ornstein-Uhlenbeck processes. Joseph effect in models with infinite variance. Phys A Stat Mech Appl. (2008) 387:123–33. doi: 10.1016/j.physa.2007.08.016

CrossRef Full Text | Google Scholar

33. Voutilainen M, Viitasaari L, Ilmonen P, Torres S, Tudor C. Vector-valued generalised Ornstein-Uhlenbeck processes. arXiv[Preprint].arXiv:1909.02376. (2019).

Google Scholar

34. Lancaster P, Rodman L. Algebraic Riccati Equations. New York, NY: Oxford University Press (1995).

Google Scholar

35. Kucera V. A contribution to matrix quadratic equations. IEEE Trans Automat Control. (1972) 17:344–7. doi: 10.1109/TAC.1972.1099983

CrossRef Full Text | Google Scholar

36. Wonham WM. On a matrix Riccati equation of stochastic control. SIAM J Control. (1968) 6:681–97. doi: 10.1137/0306044

CrossRef Full Text | Google Scholar

37. Sun Jg. Perturbation theory for algebraic Riccati equations. SIAM J Matrix Anal Appl. (1998) 19:39–65. doi: 10.1137/S0895479895291303

CrossRef Full Text | Google Scholar

38. Byers R. Solving the algebraic Riccati equation with the matrix sign function. Linear Algebra Appl. (1987) 85:267–79. doi: 10.1016/0024-3795(87)90222-9

CrossRef Full Text | Google Scholar

39. Laub A. A Schur method for solving algebraic Riccati equations. IEEE Trans Automat Control. (1979) 24:913–21. doi: 10.1109/TAC.1979.1102178

CrossRef Full Text | Google Scholar

40. Bini DA, Iannazzo B, Meini B. Numerical Solution of Algebraic Riccati Equations. Vol. 9. Philadelphia: SIAM (2012). doi: 10.1137/1.9781611972092

CrossRef Full Text | Google Scholar

41. Viitasaari L. Representation of stationary and stationary increment processes via Langevin equation and self-similar processes. Stat Probabil Lett. (2016) 115:45–53. doi: 10.1016/j.spl.2016.03.020

CrossRef Full Text | Google Scholar

42. Voutilainen M, Viitasaari L, Ilmonen P. Note on AR(1)-characterisation of stationary processes and model fitting. Modern Stochast Theory Appl. (2019) 6:195–207. doi: 10.15559/19-VMSTA132

CrossRef Full Text | Google Scholar

43. Embrechts P, Maejima M. Selfsimilar Processes. Princeton: Princeton University Press (2002). doi: 10.1515/9781400825103

CrossRef Full Text | Google Scholar

44. Pipiras V, Taqqu MS. Long-Range Dependence and Self-Similarity. Vol. 45. Cambridge: Cambridge University Press (2017). doi: 10.1017/CBO9781139600347

CrossRef Full Text | Google Scholar

45. Lamperti J. Semi-stable stochastic processes. Trans Am Math Soc. (1962) 104:62–78. doi: 10.1090/S0002-9947-1962-0138128-7

CrossRef Full Text | Google Scholar

Keywords: time-series analysis, stationary processes, characterization, multivariate Ornstein-Uhlenbeck processes, generalized Langevin equation, algebraic Riccati equations, estimation, consistency

Citation: Voutilainen M (2020) Modeling and Estimation of Multivariate Discrete and Continuous Time Stationary Processes. Front. Appl. Math. Stat. 6:43. doi: 10.3389/fams.2020.00043

Received: 04 November 2019; Accepted: 03 August 2020;
Published: 17 September 2020.

Edited by:

Tommi Sottinen, University of Vaasa, Finland

Reviewed by:

Foad Shokrollahi, University of Vaasa, Finland
Ehsan Azmoodeh, Ruhr University Bochum, Germany

Copyright © 2020 Voutilainen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Marko Voutilainen, marko.voutilainen@aalto.fi

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.