Confidence interval for correlation estimator
between latent processes

arXiv:1710.06683v1 [math.ST] 18 Oct 2017

Akitoshi Kimura

∗

March 14, 2018

Abstract
Kimura and Yoshida [26] treated a model in which the finite variation
part of a two-dimensional semimartingale is expressed by time-integration
of latent processes. They proposed a correlation estimator between the
latent processes and proved its consistency and asymptotic mixed normality. In this paper, we discuss the confidence interval of the correlation
estimator to detect the correlation. We propose two types of estimators
for asymptotic variance of the correlation estimator and prove their consistency in a high frequency setting. Our model includes doubly stochastic
Poisson processes whose intensity processes are correlated Itô processes.
We compare our estimators based on the simulation of the doubly stochastic Poisson processes.

1

Introduction

Integrated correlation is an important index in high frequency financial data
analysis. Epps [16] pointed out that the sample correlation between the returns
of two different stocks decreases as the sampling frequency of data increases.
It is considered that non-synchronicity and market microstructure of trading
cause this phenomenon.
Non-synchronous covariance estimation schemes have been developed: Fourier
analytic approach (Malliavin and Mancino [32], Malliavin et al. [34]) and the
cumulative covariance estimator (Hayashi and Yoshida [23, 21, 22], Mykland
[40]).
The market microstructure is modeled as the noise added to the latent price
process. This modeling was successful and denoising techniques have also been
developed: sub-sampling (Zhang et al. [49], Zhang [48]), pre-averaging (Podolskij and Vetter [43], Jacod et al. [25]), and others (Zhou [50]). There are many
studies that treat both non-synchoronicity and market microstructure noise:
Malliavin and Mancino [33], Mancino and Sanfelici [35], Park and Linton [42],
Voev and Lunde [47], Griffin and Oomen [19], Christensen et al. [12, 13], Koike
[27, 29, 28], Aı̈t-Sahalia et al. [4], Barndorff-Nielsen et al. [6], Bibinger [8, 9].
The followings are proposed as the cause of market microstructure: bid-ask
spread (Roll [44]), discretization error (Gottlieb and Kalay [18]), and asymmetric information (Glosten and Milgrom [17]). On the other hand, the relationship
∗ Graduate School of Mathematical Sciences, University of Tokyo, 3-8-1 Komaba, Meguroku, Tokyo 153-8914, Japan, Japan Science and Technology Agency, CREST

1

between the additive noise modeling and the market microstructure is not so
clearly explained.
In ultra high frequency sampling, recently, the market microstructure is
modeled as dynamics of the limit order book (LOB) rather than the noise. Many
studies adopt approaches modeling LOB with Poisson processes (Cont et al. [15],
Abergel and Jedidi [2], Muni Toke [46, 37], Smith et al. [45], Muni Toke and
Yoshida [38]), Hawkes processes (Hewlett [24], Large [31], Bowsher [10], Bacry
et al. [5], Abergel and Jedidi [3], Muni Toke and Pomponio [39], Muni Toke
[36], Clinet and Yoshida [14], Ogihara and Yoshida [41]), and doubly stochastic
Poisson processes (Abergel et al. [1], Guilbaud and Pham [20], Chertok et al.
[11], Korolev et al. [30]). These approaches can also treat the non-synchronicity
of trading.
In this stream, Kimura and Yoshida [26] focused on the integrated correlation between intensity processes of doubly stochastic Poisson processes. They
introduced the following model to treat it in a generalized form. This paper is
also based on the model.
Now, we consider a stochastic basis B = (Ω, F , F, P ), F = (Ft )t∈[0,T ] . On B,
let X = (X 1 , X 2 ) be an R2 -valued Itô process given by
Xt = X0 +

Z

t

0

X0s ds +

Z

t

0

X1s dws (t ∈ [0, T ]),

(1.1)

where w is an r-dimensional F -Wiener process, X0 is an F0 -measurable random
variable, X0 is a two-dimensional F-adapted process, and X1 is an R2 ⊗ Rr valued F-adapted process satisfies condition [A] mentioned later. Let an be a
positive number depending on n. On B, consider a two-dimensional measurable
process Yn = (Y n,1 , Y n,2 ) having a decomposition
Ynt = Yn0 +

Z

0

t

an Xs ds + Mnt (t ∈ [0, T ]),

n

where M is a two-dimensional measurable process with Mn0 = 0 satisfies con♯
dition [B], [B’], or [B♯ ] mentioned later.
In this model, Y can be a counting process with intensity process an X where
X is R2+ -valued, i.e. Y is a two-dimensional doubly stochastic Poisson process.
This counting process models the high frequency counting data of the orders or
transactions in the active market, for example.
Let Ij = [tj−1 , tj ) for a sampling design Π = (tj )j=0,...,bn with 0 = t0 < t1 <
· · · < tbn = T and hj = tj − tj−1 .
Kimura and Yoshida [26] introduced the estimator of the covariance hX 1 , X 2 i
from the sampled data of Y as follows. For α, β = 1, 2,
Snαβ =

bn  n,α
X
Ytj − Ytn,α
j−1
j=2

an hj

−

 Y n,β − Y n,β
− Ytn,α
Ytn,α
tj−1
tj
j−2
j−1
an hj−1

an hj

−


− Ytn,β
Ytn,β
j−2
j−1
an hj−1

.

Here, Sn = (Sn12 , Sn11 , Sn22 )⋆ is the (co)variance estimator. Sn depends on the
scaling parameter an which is not derived
p from the data of Y. Therefore, treat
the correlation estimator Cn12 = Sn12 / Sn11 Sn22 which does not depend on the
scaling parameter an .

2

Let
U αβ =

2 α,c β,c
2
hX , X iT =
3
3

T

Z

0

Xtα1 · Xtβ1 dt

(α, β = 1, 2)

and U = (U 12 , U 11 , U 22 )⋆ , where X α,c is the continuous part of X α and Xtα1 is
the α-th row of X1t . Let
γ

(α1 ,β1 ),(α2 ,β2 )

=

Z

T

0

α2 1 β2 1
α2 1 β2 1
α1 1 β1 1
α1 1 β1 1
r
X
Xj,s + Xj,s
Xi,s
Xi,s
Xj,s + Xj,s
Xi,s Xi,s
ds
2
2
i,j=1

α1
and Γ = (γ pq )p,q=(1,2),(1,1),(2,2) , where Xi,s
is the (α, i)-element of X1s .
♯

Under [A] and any one of [B], [B’], and [B♯ ], Kimura and Yoshida [26]
proved the the asymptotic mixed normality of the covariance estimator and the
correlation estimator.
 T −1/2
(Sn − U ) →ds Γ1/2 ζ
bn
as n → ∞, where ζ is an R3 -valued standard normal√variable independent of F .
and ds means F -stable convergence. Let R = U 12 / U 11 U 22 and suppose that
U 11 U 22 6= 0 a.s.
 T −1/2
bn

(Cn12 − R) →ds Ξ1/2 ζ

as n → ∞, where


−U 12
1
−U 12
, p
Ξ1/2 := √
Γ1/2 .
, p
U 11 U 22 2 (U 11 )3 U 22 2 U 11 (U 22 )3

Ξ = (Ξ1/2 )(Ξ1/2 )⋆ is the asymptotic variance of the correlation estimator.
In this paper, we introduce two types of estimators for Γ and prove their
consistency in order to obtain the consistent estimator for Ξ. Therefore, we can
detect the correlation between latent processes based on hypothesis testing.
This paper is organized as follows. In Section 2, we list the assumptions. In
Section 3, we introduce the estimators for Γ and Ξ, and state the main result:
their consistency. In Section 4, the results in Section 3 are proved. In Section
5, an example and simulation studies are given.

2

Assumptions

For simplicity, we assume tj = j(T /bn ), hj = T /bn =: δn . We write ∆j V =
Vtj − Vtj−1 for a process V . Consider the following conditions.
[A] Process X admits the representation (1.1) for an R2 -valued F0 -measurable
random variable X0 and coefficients Xκ (κ = 0, 1) such that X0 is a cádlág
F-adapted process and that X1 has a representation
X1t

=

X10

+

Z

0

t

X10
s

ds +

Z

t

′

r
X

0 κ′ =1

3

′

dw̃sκ
X1κ
s

′

(t ∈ [0, T ]),

where X10 is an R2 ⊗ Rr -valued F0 -measurable random variable, w̃ =
′
(w̃1 , ..., w̃r ) is an r′ -dimensional F-Wiener process (not necessary indepen′
dent of w), and X1κ
(κ′ = 0, 1, ..., r′ ) are R2 ⊗ Rr -valued cádlág F-adapted
s
processes.
[B] Mn = (M n,α )α=1,2 is a two-dimensional F-local martingale with Mn0 = 0
and such that
5/2

(i) limn→∞ bn /an = 0.
Pbn
1/2
n 2
n
(ii)
for a
j=1 |∆j M | = Op (an ) as n → ∞, supt∈[0,1] |∆M | ≤ can
constant c independent of n.
(iii) The absolutely continuous (w.r.t. the Lebesgue measure a.s.) mapping [0, T ] ∋ t 7→ hMn , wit ∈ R2 ⊗Rr satisfies supt∈[0,T ] |dhMn , wit /dt| =
Op (bn ) as n → ∞.
Here, hMn , wit is the 2 × r matrix of angle brackets hM n,α , wk it for Mn =
(M n,α )α=1,2 and w = (wk )k=1,...,r .
[B′ ] Mn is a two-dimensional F-local martingale with Mn0 = 0, satisfies [B] (i),
(iii), and
Pbn
(ii′ ) E[ j=1
|∆j Mn |2 ] = O(an ) as n → ∞.
♯

[B♯ ] Mn is a two-dimensional F-local martingale with Mn0 = 0, satisfies [B] (ii),
(iii), and
♯

(i♯ ) limn→∞ b3n /an = 0.
Pbn
[C] E[ j=1
|∆j Mn |4 ] = O(a2n b−1
n ) as n → ∞

♯

Remark 2.1 (A). , [B], and [B′ ] are the same condition in [26]. [B♯ ] is a bit
stronger condition than [B♯ ] in [26].

3

Results

We need the estimator of Γ to obtain the asymptotic variance estimator. Here,
we give two types of the gamma estimator.
Before that, we introduce some notations to simplify the description. Let
A⊗(i,j) is an (i, j)-element of A⊗ , A⊗ = A ⊗ A = AA⋆ , and ⋆ denotes transpose.
˜ = ((xi yj + xj yi )/2) ∈ Rr ⊗ Rr for x = (xi ), y = (yi ) ∈ Rr , and let
Let x⊗y
P
˜
⊗(α,β)
α
˜ β
x
= xα
) ∈ R2 ⊗Rr . We write x·y = ri=1 xi yi for x = (xi ),
· ⊗x· for x = (xiP
r
y = (yi ) ∈ Rr , and x · y = i,j=1 xi,j yi,j for x = (xi,j ), y = (yi,j ) ∈ Rr ⊗ Rr .
Write
Ṽk :=

∆k V
∆k V
=
.
an hk
an δn

for a stochastic process V . For example, we can write
Snαβ

=

bn
X
j=2

(Ỹnj

−

Ỹnj−1 )⊗(α,β) ,

γ

pq

=

Z

0

4

1

˜

˜

(X1s )⊗p · (X1s )⊗q ds.

Γ estimator

3.1

In the similar way in Barndorff-Nielsen and Shephard [7], we define
b

pq
γ̂n,1

=

n
9n X
(Ỹnk − Ỹnk−1 )⊗p (Ỹnk − Ỹnk−1 )⊗q
8

k=2

−

bn −2 
o
1 X
(Ỹnk − Ỹnk−1 )⊗p (Ỹnk+2 − Ỹnk+1 )⊗q + (Ỹnk+2 − Ỹnk+1 )⊗p (Ỹnk − Ỹnk−1 )⊗q
2
k=2

By the similar idea, we define
pq
=
γ̂n,2

bn −2 n
o
on
9 X
1
(Ỹnk+2 − Ỹnk+1 )⊗p − (Ỹnk − Ỹnk−1 )⊗p (Ỹnk+2 − Ỹnk+1 )⊗q − (Ỹnk − Ỹnk−1 )⊗q .
8
2
k=2

3.2

Kernel based Γ estimator

In the similar way in Hayashi and Yoshida [22], we define the kernel based
estimator.
α

β

∂h {X , X }k :=

k
X

(∆l X α ∆l X β )h−1 .

l=(k−n(h)+1)+
−1/2 −1
hn

where h = hn is a parameter satisfying hn → 0 and bn
maxm {tm ≤ h}.

→ 0, and n(h) :=

b

pq
γ̂n,h
=

n n
o T 2
9X
∂h {Ỹ n,α1 , Ỹ n,α2 }k ∂h {Ỹ n,β1 , Ỹ n,β2 }k + ∂h {Ỹ n,α1 , Ỹ n,β2 }k ∂h {Ỹ n,β1 , Ỹ n,α2 }k
8
bn

k=2

3.3

Ξ estimator

Here, we define


−Sn12
−Sn12
1
1/2
1/2
Ξ̂n,∗ = p
, p
, p
Γ̂n,∗
Sn11 Sn22 2 (Sn11 )3 Sn22 2 Sn11 (Sn22 )3
1/2

and Ξ̂n,∗ = (Ξ̂n,∗ )⊗ , for ∗ = 1, 2, h.

3.4

consistency of the estimators

pq
Theorem 3.1. (a) (T /bn )−1 γ̂n,i
= γ pq + op (1) under [A] and any one of [B],
♯

[B′ ] and [B♯ ], (i = 1, 2).
♯
−1/2
pq
(b) (T /bn )−1 γ̂n,i
= γ pq + Op (bn ) under [A] and [B♯ ], (i = 1, 2).
pq
Theorem 3.2. (T /bn )−1 γ̂n,h
→p γ pq under [A] and any one of [B], [B′ ] and
♯

[B♯ ].

corollary 3.3. (1a) (T /bn )−1 Ξ̂n,i = Ξ + op (1) under [A] and any one of [B],
♯
[B′ ] and [B♯ ], (i = 1, 2).
♯
−1/2
(1b) (T /bn )−1 Ξ̂n,i = Ξ + Op (bn ) under [A] and [B♯ ], (i = 1, 2).
♯
(2) (T /bn )−1 Ξ̂n,h = Ξ + op (1) under [A] and any one of [B], [B′ ] and [B♯ ].
5

4

Proof

Proof of Theorem 3.1
pq
First, we approximate γ̂n,i
(i = 1, 2) by the following quantities. Let
b

Gpq
n,1

n
9n X
(χ̃nk − χ̃nk−1 )⊗p (χ̃nk − χ̃nk−1 )⊗q
8

=

k=2

−

bn −2 
o
1 X
(χ̃nk − χ̃nk−1 )⊗p (χ̃nk+2 − χ̃nk+1 )⊗q + (χ̃nk+2 − χ̃nk+1 )⊗p (χ̃nk − χ̃nk−1 )⊗q
2
k=2

and
Gpq
n,2 =

bn −2 n
o
on
9 X
1
(χ̃nk+2 − χ̃nk+1 )⊗p − (χ̃nk − χ̃nk−1 )⊗p (χ̃nk+2 − χ̃nk+1 )⊗q − (χ̃nk − χ̃nk−1 )⊗q ,
8
2
k=2

where χn = (χn,1 , χn,2 ) and χn,α
=
t

Rt
0

an Xsα ds.

pq
Lemma 4.1. (a) (T /bn )−1 γ̂n,i
= (T /bn )−1 Gpq
n,i + op (1) under any one of [B],
♯

[B′ ] and [B♯ ], (i=1,2).
♯
−1/2
pq
(b) (T /bn )−1 γ̂n,i
= (T /bn )−1 Gpq
) under [B♯ ], (i=1,2).
n,i + op (bn
Proof. Write (V 1 ⊗V 2 ⊗· · ·⊗V n )i1 ,i2 ,...,in = Vi11 Vi22 · · · Vinn for vectors V 1 , V 2 , . . . , V n
and confirm |V 1 ⊗V 2 ⊗· · ·⊗V n | = |V 1 ||V 2 | · · · |V n |. Let sk = (Ỹnk −Ỹnk−1 )⊗ , tk =
α1 ,β1 ,α2 ,β2
α1 ,β1 ,α2 ,β2
)α1 ,β1 ,α2 ,β2 =1,2
)α1 ,β1 ,α2 ,β2 =1,2 , and Gn,i = (Gn,i
(χ̃nk −χ̃nk−1 )⊗ , γ̂n,i = (γ̂n,i
(i = 1, 2). It holds that
 9 −1  T −1
|γ̂n,1 − Gn,1 |
8
bn
bn
bn −2 
 T −1 X

1 X
=
(sk ⊗ sk − tk ⊗ tk ) −
(sk ⊗ sk+2 − tk ⊗ tk+2 ) + (sk+2 ⊗ sk − tk+2 ⊗ tk )
bn
2
k=2

≤

bn
 T −1 n X

bn

k=2

k=2

|sk ⊗ sk − tk ⊗ tk | +

bn −2
bn −2
o
1 X
1 X
|sk ⊗ sk+2 − tk ⊗ tk+2 | +
|sk+2 ⊗ sk − tk+2 ⊗ tk |
2
2
k=2

k=2

and
 9 −1  T −1
=

|γ̂n,2 − Gn,2 |
8
bn
bn −2
1  T −1 X
(sk+2 ⊗ sk+2 − tk+2 ⊗ tk+2 ) − (sk+2 ⊗ sk − tk+2 ⊗ tk )
2 bn
k=2

−(sk ⊗ sk+2 − tk ⊗ tk+2 ) + (sk ⊗ sk − tk ⊗ tk )
≤

bn −2
bX
n −2
1  T −1 n X
|sk+2 ⊗ sk+2 − tk+2 ⊗ tk+2 | +
|sk+2 ⊗ sk − tk+2 ⊗ tk |
2 bn
k=2

+

bX
n −2
k=2

|sk ⊗ sk+2 − tk ⊗ tk+2 | +

k=2

bX
n −2
k=2

6

o
|sk ⊗ sk − tk ⊗ tk | .

Now, let ρ = 0, 1/2 and see the typical term.
b1+ρ
n

bX
n −2
k=2

=

b1+ρ
n

|sk+2 ⊗ sk − tk+2 ⊗ tk |

bX
n −2

(sk+2 − tk+2 ) ⊗ (sk − tk ) + (sk+2 − tk+2 ) ⊗ tk + tk+2 ⊗ (sk − tk )

bX
n −2

(sk+2 − tk+2 ) ⊗ (sk − tk ) + (sk+2 − tk+2 ) ⊗ tk + tk+2 ⊗ (sk − tk )

k=2

≤

b1+ρ
n

=

b1+ρ
n

k=2

bX
n −2
k=2

≤

|sk+2 − tk+2 ||sk − tk | + |sk+2 − tk+2 ||tk | + |tk+2 ||sk − tk |

bX
bX
n −2
n −2



b1/2+ρ/2
|sk+2 − tk+2 | b1/2+ρ/2
|sk − tk |
n
n
k=2

k=2

n −2
n −2
 bX

 
 bX

+ bρn
|sk+2 − tk+2 | bn max |tk | + bn max |tk+2 | bρn
|sk − tk |

k

k=2

k

k=2

♯

Under [B♯ ], the same argument of the proof of Lemma 1 (b) in [26] under the
assumption [B] yields
b3/4
n

bn
X

k=2

|sk − tk | →p 0.
♯

Under any one of [B], [B′ ] and [B♯ ], it holds that
b1/2
n

bn
X

k=2

|sk − tk | →p 0

by the proof of Lemma 1 (b) in [26]. By inequality (4) in [26], and chebychev’s
inequality or Markov’s inequality, we have
bn max |tk | = Op (1).
k

Therefore, the typical term converges to 0 in probability. The other terms are
evaluated in the same way. We obtain the conclusion.
pq
Proposition 4.2. (T /bn )−1 γ̂n,i
= (T /bn )−1 Gpq
n,i + op (bn
and [C] (i=1,2).

−1/2

) under [B] or [B’],

Proof. We just confirm the convergence of the easiest term.
b1+ρ
n

bX
n −2
k=2

=

b1+ρ
n

bX
n −2
k=2

|sk ⊗ sk − tk ⊗ tk |
(sk − tk ) ⊗ (sk − tk ) + (sk − tk ) ⊗ tk + tk ⊗ (sk − tk )

7

By definition,
sk − tk = (M̃nk − M̃nk−1 )⊗ + (M̃nk − M̃nk−1 ) ⊗ (χ̃nk − χ̃nk−1 ) + (χ̃nk − χ̃nk−1 ) ⊗ (M̃nk − M̃nk−1 ).
Therefore, it is enough to prove the following equations.
b1+ρ
n

bX
n −2

=

op (1)

(M̃nk − M̃nk−1 )⊗⊗ ⊗ (χ̃nk − χ̃nk−1 ) =

op (1)

(M̃nk − M̃nk−1 )⊗ ⊗ (χ̃nk − χ̃nk−1 )⊗

=

op (1)

=

op (1),

k=2

bX
n −2

b1+ρ
n

k=2

b1+ρ
n

(M̃nk − M̃nk−1 )⊗⊗⊗

bX
n −2
k=2

b1+ρ
n

bX
n −2
k=2

(M̃nk − M̃nk−1 ) ⊗ (χ̃nk − χ̃nk−1 )⊗⊗

where for vector V , V ⊗ = V ⊗ V , V ⊗⊗ = V ⊗ V ⊗ V , V ⊗⊗⊗ = V ⊗ V ⊗ V ⊗ V .
By assumption [C], it holds that




bX
bX
n −2
n −2
n
n
⊗⊗⊗
1+ρ
n
n
4
E b1+ρ
(
M̃
−
M̃
)
≤
E
b
|
M̃
−
M̃
|
n
k
k−1
n
k
k−1
k=2

≤

k=2



bX
n −2
n 4
−2
16E b1+ρ
|
M̃
|
= O(b4+ρ
n
k
n an ) = o(1).
k=2

Lemma 1 (b) and inequality (4) in [26] yields
b1+ρ
n

bX
n −2
k=2

≤



bρn

bX
n −2
k=2

(M̃nk − M̃nk−1 )⊗ ⊗ (χ̃nk − χ̃nk−1 )⊗

|M̃nk − M̃nk−1 |2



bn |χ̃nk − χ̃nk−1 |2 = op (1)Op (1) = op (1).

and
b1+ρ
n

bX
n −2
k=2

bρn

≤
=

(M̃nk − M̃nk−1 ) ⊗ (χ̃nk − χ̃nk−1 )⊗⊗

bX
n −2
k=2

(M̃nk − M̃nk−1 ) ⊗ (χ̃nk − χ̃nk−1 ) bn (χ̃nk − χ̃nk−1 )⊗

op (1)Op (1) = op (1).

From the above inequalities, it is derived that
b1+ρ
n

bX
n −2
k=2

≤
=



b1+ρ
n

(M̃nk − M̃nk−1 )⊗⊗ ⊗ (χ̃nk − χ̃nk−1 )

bX
n −2

(M̃nk

k=2

−

M̃nk−1 )⊗⊗⊗

1/2 
1/2
bX
n −2
1+ρ
n
n
⊗
n
n
⊗
bn
(M̃k − M̃k−1 ) ⊗ (χ̃k − χ̃k−1 )
k=2

op (1)op (1) = op (1).
8

Now, we let
Aα
k

Xtα1
k−2

=

Bkα

=

1
hk

Z

Ckα

=

1
hk

Z

Dkα

=

1
hk

Z

Z

Hk (s) dws +

Ik

Ik−1
Z t

Ik

Z

Ik

tk−2
t

tk−2
Z t
tk−2

Xtα1
k−2

Z

Kk (s) dws
Ik

1

(Xsα1 − Xtα1
) dws dt −
k−2
1

Xtα0
ds dt −
k−2

hk−1

Z

Ik−1

(Xsα0 − Xtα0
)ds dt −
k−2

hk−1
Z t

tk−2

1
hk−1

Z

Ik−1

Z

t

tk−2

(Xsα1 − Xtα1
) dws dt
k−2

Xtα0
ds dt
k−2

Z

Ik−1

Z

t

tk−2

(Xsα0 − Xtα0
)ds dt.
k−2

By Lemma 2 and inequality (4) in [26], we have
 T −1
bn

Gpq
n,1

bn
 T −1 9 n X

=

bn

−

8

β1
α2 β2
1
(Aα
k Ak )(Ak Ak )

k=2

bX
n −2 

1
2

o
β1
β2
β1
α2
α1
α2 β2
1
(Aα
)(A
)
+
(A
A
)(A
A
)
+ op (bn−1/2 )
A
A
k
k
k+2 k+2
k
k+2 k+2
k

k=2

=

 T −1

=

n −2 n
 T −1 9 1 bX

bn

−1/2
R1,pq
)
n,1 + op (bn

and
 T −1
bn

Gpq
n,2

bn

82

k=2

β1
β2
β1
α2
α1
α2 β2
1
(Aα
k+2 Ak+2 )(Ak+2 Ak+2 ) − (Ak+2 Ak+2 )(Ak Ak )

β1
β2
α2
1
−(Aα
k Ak )(Ak+2 Ak+2 )

 T −1

=

bn

o
β1
α2 β2
−1/2
1
+ (Aα
)
k Ak )(Ak Ak ) + op (bn

−1/2
R1,pq
).
n,2 + op (bn

Here, we write
qk2,αβ

=

Xtα1
· Xtβ1
k−2
k−2

hk
hk−1
+ Xtα1
· Xtβ1
k−2
k−2
3
3

and
qk3,αβ

˜

= 2(X1tk−2 )⊗(α,β) ·
+

Z

Ik

Z

Z

Ik−1

Z

t

tk−2

Hk (s) dws ⊗ Hk (t) dwt

t

tk−1

Kk (s) dws ⊗ Kk (t) dwt +

Z

Ik−1

Hk (s) dws ⊗

Z

Ik


Kk (s) dws .

Then, we have
 T −1
bn

R1,pq
n,1

=

bn
 T −1 9 n X

bn

−

1
2

8

bX
n −2 
k=2

(qk2,p + qk3,p )(qk2,q + qk3,q )

k=2

o
2,q
3,q
2,p
3,p
(qk2,p + qk3,p )(qk+2
+ qk+2
) + (qk+2
+ qk+2
)(qk2,q + qk3,q )
9

and
 T −1
bn

R1,pq
n,2

n −2 n
 T −1 9 1 bX

=

bn

82

−(qk2,p

k=2

2,p
3,p
2,p
3,p
2,p
3,p
(qk+2
+ qk+2
)(qk+2
+ qk+2
) − (qk+2
+ qk+2
)(qk2,p + qk3,p )

2,p
qk3,p )(qk+2

+

o
3,p
+ qk+2
) + (qk2,p + qk3,p )(qk2,p + qk3,p )

by Itô’s formula and simple calculus. We only see the typical terms. We have
bn
 T −1 X

bn

b

k=2

k=2

bn
 T −1 X

bn

n
T
4X
11
21
11
21
(Xtαk−2
· Xtβk−2
)(Xtαk−2
· Xtβk−2
) ,
9
bn

qk2,p qk2,q =

bn
X
2 α1 1
11
(X
· Xtβk−2
)qk3,q = Op (bn−1/2 ),
3 tk−2

qk2,p qk3,q =

k=2

k=2

and
n −2
 T −1 bX

bn

3,p 3,q
qk+2
qk = Op (bn−1/2 ).

k=2

By Lemma 4 in [26] and simple calculus, it holds that
bn
 T −1 X

bn

qk3,p qk3,q

=

k=2

bn n X
r
 T −1 X

bn
Z

k=2
t

Ik−1

+

Z

tk−2

Z

Z

I

+

Z k−1

Z

t

tk−2

Hk (s) dwsi Hk (t) dwtj +

dwsi

bn n X
r
 T −1 X

bn

k=2

Z

Ik

Kk (s) dwsj

Z

Ik

Kk (s)

dwsj

Ik tk−1
r
X

+

T

0

=

Z

Kk (s) dwsi Kk (t) dwtj

β2 1
β2 1
α2 1
α2 1
(Xi,t
Xj,t
+ Xj,t
Xi,t
)
k−2
k−2
k−2
k−2

i,j=1

Z Z
Ik

t
tk−1

Kk (s) dwsi Kk (t) dwtj

o

β1 1
β1 1
α1 1
α1 1
(Xi,t
Xj,t
+ Xj,t
Xi,t
)
k−2
k−2
k−2
k−2

i,j=1

on Z
β2 1
β2 1
α2 1
α2 1
(Xi,t
X
+
X
X
)
j,tk−2
j,tk−2 i,tk−2
k−2
Z

t

Z Z

on

Hk (s) dwsi Hk (t) dwtj +

Hk (s)

Ik−1

=

i,j=1

Hk (s) dwsi

Ik−1

β1 1
β1 1
α1 1
α1 1
(Xi,t
Xj,t
+ Xj,t
Xi,t
)
k−2
k−2
k−2
k−2

t

2

2

(K̃k (s)) ds(K̃k (t)) dt +
0

Z

T

0
T

0

Z

t

(H̃k (s))2 ds(H̃k (t))2 dt

0
2

(H̃k (s)) ds

Z

0

T

o
(H̃k (s))2 ds + Op (bn−1/2 )

bn n
oT
4X
11
21
21
11
11
21
11
21
(Xtαk−2
· Xtαk−2
)(Xtβk−2
· Xtβk−2
) + (Xtαk−2
· Xtβk−2
)(Xtβk−2
· Xtαk−2
)
9
bn
k=2

+Op (bn−1/2 ),

10

where H̃k (s) = 1Ik−1 (s)Hk (s) and K̃k (s) = 1Ik (s)Kk (s). Therefore, we have
 bn
 T −1
9 4X
T
1,pq
11
21
11
21
Rn,1 =
(Xtαk−2
· Xtβk−2
)(Xtαk−2
· Xtβk−2
)
bn
8 9
bn
k=2

+

−

=
and
 T −1
bn

R1,pq
n,2

=

bn n
4X

9

k=2

oT
11
21
21
11
11
21
11
21
(Xtαk−2
· Xtαk−2
)(Xtβk−2
· Xtβk−2
) + (Xtαk−2
· Xtβk−2
)(Xtβk−2
· Xtαk−2
)
bn


bn −2
bn −2
1n4 X
4 X
T
To
11
21
11
21
(Xtαk−2
· Xtβk−2
)(Xtαk2 1 · Xtβk2 1 ) +
(Xtαk1 1 · Xtβk1 1 )(Xtαk−2
· Xtβk−2
)
2 9
bn
9
bn

k=2
−1/2
+Op (bn )
γ pq + op (bn−1/2 )

k=2

+ Op (bn−1/2 )

 bn −2
91 4 X
T
(Xtαk1 1 · Xtβk1 1 )(Xtαk2 1 · Xtβk2 1 )
82 9
bn
k=2

+

−
+

4
9

bX
n −2 n

4
9

bX
n −2

4
9

bX
n −2

4
+
9

=
−1/2

The last op (bn
calculus.

k=2

k=2

k=2

oT
11
21
21
11
11
21
11
21
(Xtαk−2
· Xtαk−2
)(Xtβk−2
· Xtβk−2
) + (Xtαk−2
· Xtβk−2
)(Xtβk−2
· Xtαk−2
)
bn

21
21
(Xtαk1 1 · Xtβk1 1 )(Xtαk−2
· Xtβk−2
)

bn −2
4 X
T
T
11
11
−
(Xtαk−2
· Xtβk−2
)(Xtαk2 1 · Xtβk2 1 )
bn
9
bn

11
21
11
21
(Xtαk−2
· Xtβk−2
)(Xtαk−2
· Xtβk−2
)

bX
n −2 n

11
(Xtαk−2

k=2
+Op (bn−1/2 ).
γ pq + op (bn−1/2 )

·

11
21
Xtαk−2
)(Xtβk−2

·

k=2

T
bn

21
Xtβk−2
)

+

11
(Xtαk−2

·

21
11
Xtβk−2
)(Xtβk−2

·

oT 

21
Xtαk−2
)

bn

+ Op (bn−1/2 ).

)s are derived from the proof of Lemma 3 in [26] and simple

Proof of Theorem 3.2
In the similar way to the proof of Therem 3.1, let
b

Gpq
n,h

n n
o T 2
9X
.
=
∂h {χ̃n,α1 , χ̃n,α2 }k ∂h {χ̃n,β1 , χ̃n,β2 }k + ∂h {χ̃n,α1 , χ̃n,β2 }k ∂h {χ̃n,β1 , χ̃n,α2 }k
8
bn

k=2

pq
Lemma 4.3. (a) (T /bn )−1 γ̂n,h
= (T /bn )−1 Gpq
n,h + op (1) under any one of [B],
♯

[B′ ] and [B♯ ].
♯
−1/2
pq
(b) (T /bn )−1 γ̂n,h
= (T /bn )−1 Gpq
) under [B♯ ].
n,h + op (bn
Proof. We write
šk,h :=

k
X

sl h−1 ,

ťk,h :=

k
X

l=(k−n(h)+1)+

l=(k−n(h)+1)+

11

tl h−1 ,

pq
γ̂n,h = (γ̂n,h
)p,q=(1,2),(2,1),(1,1),(2,2) , and Gn,h = (Gpq
n,h )p,q=(1,2),(2,1),(1,1),(2,2) . It
holds that
 9 −1  T −1
|γ̂n,h − Gn,h |
bρn
8
bn
bn
 T −1 X
 T 2
= bρn 2
(šk,h ⊗ šk,h − ťk,h ⊗ ťk,h )
bn
bn
k=2

=

bρn 2

≤

bρn 2

=

bρn 2

≤

bn
 T −1 X

bn

(šk,h − ťk,h ) ⊗ (šk,h − ťk,h ) + (šk,h − ťk,h ) ⊗ ťk,h + ťk,h ⊗ (šk,h − ťk,h )

k=2

bn n
 T −1 X

bn

k=2

bn
 T −1 X

bn

k=2

 T 2

(šk,h − ťk,h ) ⊗ (šk,h − ťk,h ) + (šk,h − ťk,h ) ⊗ ťk,h + ťk,h ⊗ (šk,h − ťk,h )

|šk,h − ťk,h |

2

 T 2
bn

+

bρn 4

bn
 T −1 X

bn

k=2

|ťk,h ||šk,h − ťk,h |

 T 2
bn


bn
bn
 T −1/2 X
o X
n
 T 2
 T 
ρ
.
2 bρ/2
+
4
max
|
ť
|
b
|š
−
ť
|
|š
−
ť
|
k,h
k,h
k,h
k,h
k,h
n
n
k
bn
bn
bn
k=2

k=2

We obtain the conclusion from the inequalities
bn
X

k=2

|šk,h − ťk,h |

T 
bn

≤

bn
X

k
X

k=2 l=(k−n(h)+1)+

|sl − tl | h

−1

T 
bn

≤

bn
X

k=2

|sk − tk |,

max |ťk,h | ≤ max |tk |,
k

k

and the proof of Lemma 4.1.
pq
Proposition 4.4. (T /bn )−1 γ̂n,h
= (T /bn )−1 Gpq
n,h + op (bn
and [C].

−1/2

) under [B] or [B’],

Proof. The conclusion is derived from the proof of Lemma 4.3, that of Proposition 4.2, and the inequality
bρn

bn
 T −1 X

bn

k=2

|šk,h − ťk,h |2

 T 2
bn

≤ bρn

bn
 T −1 X

bn

k=2

|sk − tk |2 .

.
Lemma 4.5. bρn

Pk

bn

l=(k−n(h)+1)+
β
1, 2}2 \{(Aα
l , Al ); α, β = 1, 2}.

α
α
α
Fl Gl h−1 →p 0 for all pairs of (Fl, Gl ) ∈ {Aα
k , Bk , Ck , Dk ; α =

The proof of this lemma is very similar to proof of Lemma 2 in [26].
Pk
β
Proof. Let Qn,h = bρn l=(k−n(h)+1)+ Fl Gl h−1 . For Fl = Aα
l and Gl = Cl ,
−2 −3
E[(|Qn,h |)2 ] = O(b2ρ
bn ) = O(h−1 bn2ρ−2 ) = o(1). For Fl = Aα
n n(h)h
l and
ρ−1/2
β
ρ
−1 −3/2
Gl = Dl , E[|Qn,h |] = O(bn n(h)h bn )o(1) = o(bn
) = o(1). For Fl = Aα
l
β
2
2ρ
−2 −3
−1 2ρ−2
and Gl = Bl , E[(|Qn,h |) ] = O(bn n(h)h bn ) = O(h bn ) = o(1).
12

o T 2
bn

Now, we let
k
X

Q1,αβ
k,h =

β −1
Aα
.
l Al h

l=(k−n(h)+1)+

By Lemma 4.3. and (4) (esp. Q1,αβ
k,h = Op (1)), it holds that
 T −1
bn

Gpq
n,h

=

bn
 T −1 X
9

=:

 T −1

bn

k=2

bn

 T 2
1 α2
1 β2
1 β2
1 α2
)
Q1,β
+ Q1,α
Q1,β
(Q1,α
+ op (bn−1/2 )
k,h
k,h
k,h
k,h
8
bn

−1/2
R1,pq
).
n,h + op (bn

2,αβ
3,αβ
Simple calculus yields Q1,αβ
k,h = Qk,h + Qk,h , where
k
X

Q2,αβ
k,h =

ql2,αβ h−1 ,

k
X

Q3,αβ
k,h =

ql3,αβ h−1 .

l=(k−n(h)+1)+

l=(k−n(h)+1)+

2
−2
−1 −1
By E[|Q3,αβ
n(h)b−2
bn ) and Q2,αβ
n ) = O(h
k,h | ] = O(h
k,h = Op (1), we have

 T −1
bn

R1,pq
n,h

=

bn
 T −1 X
9

 T 2
1 α2
1 β2
1 β2
1 α2
)
Q2,β
+ Q2,α
Q2,β
(Q2,α
+ op (bn−1/2 h−1/2 )
k,h
k,h
k,h
k,h
8
bn

bn
k=2
 T −1
−1/2 −1/2
R2,pq
h
).
=:
n,h + op (bn
bn

Now, we confirm Q2,αβ
k,h =
derived that
Q2,αβ
k,h

2 α1
3 Xtk−2

· Xtβ1
+ op (1). From simple calculus, it is
k−2



αβ
−1/2
−
U
+
o(b
)
h−1 Utαβ
t
n
k−2
(k−2−n(h))+
Z tk−2
2
X α1 · Xsβ1 ds + op (bn−1/2 h−1 )
h−1
3 tk−2 −h s
2 α1
X
· Xtβ1
+ Op (h1/2 ) + op (bn−1/2 h−1 ).
k−2
3 tk−2

=
=
=

Therefore, it holds that
 T −1
bn

R2,pq
n,h

=

bn 
 2
 2

 2
2 α1 1
9T X
11
21
11
21
11
21
21
Xtk−2 · Xtαk−2
Xtβk−2
· Xtβk−2
+
Xtαk−2
Xtβk−2
· Xtαk−2
· Xtβk−2
8 bn
3
3
3
3
k=2

+op (1)
=

bn  X
r
r
 2 X

2
9T X
β1 1
β2 1
α1 1
α2 1
Xi,t
X
X
X
j,t
i,t
j,t
k−2
k−2
k−2
k−2
8 bn
3 i=1
3 j=1
k=2

+

=

r
2 X

3

β2 1
α1 1
Xi,t
Xi,t
k−2
k−2

i=1

9  2 2 T
4 3 bn

bn
X

r
X

r
 2 X
j=1

β1 1
α1 1
Xi,t
Xj,t
k−2
k−2

β1 1
β2 1
β2 1
α1 1
α2 1
α2 1
+ Xj,t
Xi,t
Xi,t
Xj,t
+ Xj,t
Xi,t
k−2
k−2
k−2
k−2
k−2
k−2

2

k=2 i,j=1

+op (1)
=

3


β1 1
α2 1
+ op (1)
Xj,t
Xj,t
k−2
k−2

γ pq + op (bn−1/2 ) + op (1).
13

2

5

Simulation

We consider the following example. X is an R2+ -valued Itô process satisfying
Xt1
Xt2

= X01 +
=

X02

+

Z

t

Z

t

0

0

Xs1 µ1 ds +

Z

Xs2 µ2

Z

ds +

t
0
t
0

Xs1 σ1 dws1
Xs2 ρσ2 dws1 +

Z

0

t

Xs2

p
1 − ρ2 σ2 dws2 ,

where w1 and w2 are independent Wiener processes, and t ∈ [0, T ]. Y is a
counting process with intensity process an X. Now, we fix the parameters µ1 =
0.2, µ2 = 0.3, σ1 = 0.2, σ2 = 0.3, ρ = 0.7, X01 = 1, X02 = 2, T = 1. We examine
the asymptotic variance estimators in various an , bn settings where h = b−0.45
.
n
Calculate these with 1000 paths and .compare these by mean squared error
(MSE).

1
2
h

Table 1: MSE of asymptotic variance
a n=10^5
a n=10^6
a n=10^7
0.61604509 0.15684688 0.01680340
0.58219026 0.14576924 0.01574452
1.18562040 0.56408060 0.27658070

estimator under bn = 102
a n=10^8
a n=10^9
0.01497601 0.01496691
0.01435753 0.01463455
0.25139000 0.24816300

a n=10^10
0.01391545
0.01330052
0.24773520

1
2
h

Table 2: MSE of asymptotic variance
a n=10^5
a n=10^6
a n=10^7
0.68418095 0.69442275 0.60087796
0.68064701 0.69036001 0.59723320
1.60337550 1.60723040 1.46665960

estimator under bn = 103
a n=10^8
a n=10^9
0.11642533 0.00371270
0.11544003 0.00363454
0.61081140 0.26640340

a n=10^10
0.00156057
0.00155093
0.23161730

1
2
h

Table 3: MSE of asymptotic variance
a n=10^7
a n=10^8
a n=10^9
0.68759049 0.68810632 0.60240947
0.68718650 0.68770141 0.60207140
2.28807526 2.29027915 2.12823010

estimator under bn = 104
a n=10^10 a n=10^11
0.11224316 0.00230202
0.11214653 0.00229482
0.98841491 0.50342850

a n=10^12
0.00016181
0.00016117
0.45111698

pq
pq
Tables 1, 2, and 3 show the MSE of the Ξ estimators derived from γ̂n,1
,γ̂n,2
,
pq
and γ̂n,h respectively. We can see that the kernel based version (indexed by h)
converges more slowly than others. This result is imagined from the proof of
Theorem 3.2 and the choice of parameter h.

Acknowledgment
The author would like to thank his supervisor Prof. Nakahiro Yoshida who
provided valuable suggestions and comments.

14

References
[1] F. Abergel, C. Huré, and H. Pham. Algorithmic trading in a microstructural limit order book model. ArXiv e-prints, May 2017.
[2] Frédéric Abergel and Aymen Jedidi. A mathematical approach to order
book modeling. International Journal of Theoretical and Applied Finance,
16(05):1350025, 2013.
[3] Frédéric Abergel and Aymen Jedidi. Long-time behavior of a hawkes
process–based limit order book. SIAM Journal on Financial Mathematics, 6(1):1026–1043, 2015.
[4] Yacine Aı̈t-Sahalia, Jianqing Fan, and Dacheng Xiu. High-frequency covariance estimates with noisy and asynchronous financial data. Journal of
the American Statistical Association, 105(492):1504–1517, 2010.
[5] Emmanuel Bacry, Sylvain Delattre, Marc Hoffmann, and Jean-François
Muzy. Modelling microstructure noise with mutually exciting point processes. Quantitative Finance, 13(1):65–77, 2013.
[6] Ole E Barndorff-Nielsen, Peter Reinhard Hansen, Asger Lunde, and
Neil Shephard. Multivariate realised kernels: consistent positive semidefinite estimators of the covariation of equity prices with noise and nonsynchronous trading. Journal of Econometrics, 162(2):149–169, 2011.
[7] Ole E. Barndorff-Nielsen and Neil Shephard. Econometric analysis of realized covariation: High frequency based covariance, regression, and correlation in financial economics. Econometrica, 72(3):885–925, 2004.
[8] Markus Bibinger. Efficient covariance estimation for asynchronous noisy
high-frequency data. Scandinavian Journal of Statistics, 38(1):23–45, 2011.
[9] Markus Bibinger. An estimator for the quadratic covariation of asynchronously observed itô processes with noise: asymptotic distribution theory. Stochastic Processes and their Applications, 122(6):2411–2453, 2012.
[10] Clive G Bowsher. Modelling security market events in continuous time: Intensity based, multivariate point process models. Journal of Econometrics,
141(2):876–912, 2007.
[11] A. V. Chertok, V. Yu. Korolev, and A. Yu. Korchagin. Modeling highfrequency non-homogeneous order flows by compound cox processes*. Journal of Mathematical Sciences, 214(1):44–68, Apr 2016.
[12] Kim Christensen, Silja Kinnebrock, and Mark Podolskij. Pre-averaging
estimators of the ex-post covariance matrix in noisy diffusion models with
non-synchronous data. Journal of Econometrics, 159(1):116–133, 2010.
[13] Kim Christensen, Mark Podolskij, and Mathias Vetter. On covariation
estimation for multivariate continuous itô semimartingales with noise in
non-synchronous observation schemes. Journal of Multivariate Analysis,
120:59–84, 2013.

15

[14] Simon Clinet and Nakahiro Yoshida. Statistical inference for ergodic point
processes and application to limit order book. Stochastic Processes and
their Applications, 127(6):1800 – 1839, 2017.
[15] Rama Cont, Sasha Stoikov, and Rishi Talreja. A stochastic model for order
book dynamics. Operations research, 58(3):549–563, 2010.
[16] Thomas W Epps. Comovements in stock prices in the very short run.
Journal of the American Statistical Association, 74(366a):291–298, 1979.
[17] Lawrence R. Glosten and Paul R. Milgrom. Bid, ask and transaction prices
in a specialist market with heterogeneously informed traders. Journal of
Financial Economics, 14(1):71 – 100, 1985.
[18] Gary Gottlieb and Avner Kalay. Implications of the discreteness of observed
stock prices. The Journal of Finance, 40(1):135–153, 1985.
[19] Jim E Griffin and Roel CA Oomen. Covariance measurement in the presence of non-synchronous trading and market microstructure noise. Journal
of Econometrics, 160(1):58–68, 2011.
[20] F. Guilbaud and H. Pham. Optimal High Frequency Trading with limit
and market orders. ArXiv e-prints, June 2011.
[21] Takaki Hayashi and Nakahiro Yoshida. Asymptotic normality of a covariance estimator for nonsynchronously observed diffusion processes. Annals
of the Institute of Statistical Mathematics, 60(2):367–406, 2008.
[22] Takaki Hayashi and Nakahiro Yoshida. Nonsynchronous covariation process and limit theorems. Stochastic Processes and their Applications,
121(10):2416–2454, 2011.
[23] Takaki Hayashi, Nakahiro Yoshida, et al. On covariance estimation of nonsynchronously observed diffusion processes. Bernoulli, 11(2):359–379, 2005.
[24] Patrick Hewlett. Clustering of order arrivals, price impact and trade path
optimisation. In Workshop on Financial Modeling with Jump processes,
Ecole Polytechnique, pages 6–8, 2006.
[25] Jean Jacod, Yingying Li, Per A Mykland, Mark Podolskij, and Mathias
Vetter. Microstructure noise in the continuous case: the pre-averaging
approach. Stochastic processes and their applications, 119(7):2249–2276,
2009.
[26] Akitoshi Kimura and Nakahiro Yoshida. Estimation of Correlation Between
Latent Processes, pages 131–146. Springer International Publishing, Cham,
2016.
[27] Yuta Koike. Estimation of integrated covariances in the simultaneous presence of nonsynchronicity, microstructure noise and jumps. Econometric
Theory, pages 1–79, 2013.
[28] Yuta Koike. An estimator for the cumulative co-volatility of asynchronously
observed semimartingales with jumps. Scandinavian Journal of Statistics,
41(2):460–481, 2014.
16

[29] Yuta Koike. Limit theorems for the pre-averaged hayashi–yoshida estimator with random sampling. Stochastic Processes and their Applications,
124(8):2699–2753, 2014.
[30] V. Y. Korolev, A. V. Chertok, A. Y. Korchagin, and A. I. Zeifman. Modeling high-frequency order flow imbalance by functional limit theorems for
two-sided risk processes. ArXiv e-prints, October 2014.
[31] Jeremy Large. Measuring the resiliency of an electronic limit order book.
Journal of Financial Markets, 10(1):1–25, 2007.
[32] Paul Malliavin and Maria Elvira Mancino. Fourier series method for measurement of multivariate volatilities. Finance and Stochastics, 6(1):49–61,
2002.
[33] Paul Malliavin, Maria Elvira Mancino, et al. A fourier transform method for
nonparametric estimation of multivariate volatility. The Annals of Statistics, 37(4):1983–2010, 2009.
[34] Paul Malliavin, Maria Elvira Mancino, and Maria Cristina Recchioni. A
non-parametric calibration of the hjm geometry: an application of itô calculus to financial statistics. Japanese Journal of Mathematics, 2(1):55–77,
Mar 2007.
[35] Maria Elvira Mancino and Simona Sanfelici. Estimating covariance via
fourier method in the presence of asynchronous trading and microstructure
noise. Journal of financial econometrics, 9(2):367–408, 2011.
[36] I. Muni Toke. “Market making” behaviour in an order book model and its
impact on the bid-ask spread. ArXiv e-prints, March 2010.
[37] I. Muni Toke. Stationary distribution of the volume at the best quote in a
Poisson order book model. ArXiv e-prints, February 2015.
[38] I. Muni Toke and N. Yoshida. Modelling intensities of order flows in a limit
order book. ArXiv e-prints, February 2016.
[39] Ioane Muni Toke and Fabrizio Pomponio. Modelling trades-through in a
limited order book using hawkes processes. Economics discussion paper,
(2011-32), 2011.
[40] Per A Mykland. A gaussian calculus for inference from high frequency data.
Annals of finance, 8(2-3):235–258, 2012.
[41] T Ogihara and N Yoshida. Quasi likelihood analysis for point process
regression models. preprint, 2015.
[42] Sujin Park, Oliver B Linton, et al. Estimating the quadratic covariation
matrix for an asynchronously observed continuous time signal masked by
additive noise. 2012.
[43] Mark Podolskij, Mathias Vetter, and Margit Sommer. Estimation of volatility functionals in the simultaneous presence of microstructure noise and
jumps. Available at SSRN 950344, 2007.

17

[44] Richard Roll. A simple implicit measure of the effective bid-ask spread in
an efficient market. The Journal of Finance, 39(4):1127–1139, 1984.
[45] Eric Smith, J Doyne Farmer, L szl Gillemot, Supriya Krishnamurthy, et al.
Statistical theory of the continuous double auction. Quantitative finance,
3(6):481–514, 2003.
[46] I. Muni Toke. The order book as a queueing system: average depth and
influence of the size of limit orders. Quantitative Finance, 15(5):795–808,
2015.
[47] Valeri Voev and Asger Lunde. Integrated covariance estimation using highfrequency data in the presence of noise. Journal of Financial Econometrics,
5(1):68–104, 2007.
[48] Lan Zhang et al. Efficient estimation of stochastic volatility using noisy
observations: A multi-scale approach. Bernoulli, 12(6):1019–1043, 2006.
[49] Lan Zhang, Per A Mykland, and Yacine Aı̈t-Sahalia. A tale of two time
scales. Journal of the American Statistical Association, 2012.
[50] Bin Zhou. High-frequency data and volatility in foreign-exchange rates.
Journal of Business & Economic Statistics, 14(1):45–52, 1996.

18

