RESEARCH PAPER
FRACTIONAL SKELLAM PROCESSES WITH APPLICATIONS TO FINANCE
Alexander Kerss 1, Nikolai N. Leonenko 2, Alla Sikorskii 3
Abstract The recent literature on high frequency financial data includes models that use the difference of two Poisson processes, and incorporate a Skellam distribution for forward prices. The exponential distribution of inter-arrival times in these models is not always supported by data. Fractional gener- alization of Poisson process, or fractional Poisson process, overcomes this limitation and has Mittag-Leffler distribution of inter-arrival times. This paper defines fractional Skellam processes via the time changes in Poisson and Skellam processes by an inverse of a standard stable subordinator. An application to high frequency financial data set is provided to illustrate the advantages of models based on fractional Skellam processes. MSC 2010 : Primary 60E05, 60G22; Secondary 60G51, 26A33 Key Words and Phrases: fractional Poisson process, fractional Skellam process, Mittag-Leffler distribution, high frequency financial data
1. Introduction The advent of high frequency financial data has spurred new model- ing techniques to describe characteristics of trade by trade data. Recent literature on the subject includes [1, 2, 3, 10] where models based on the difference of two point processes are proposed. Difference of Poisson pro- cesses is considered in [3, 10], and Hawkes processes are discussed in [1, 2]. This paper extends the models where the forward price of a risky asset c 2014 Diogenes Co., Sofia pp. 532–551 , DOI: 10.2478/s13540-014-0184-2 FRACTIONAL SKELLAM PROCESSES . . . 533 is modeled via the difference of two independent Poisson processes, also known as Skellam processes. A drawback of the existing models, which may be at odds with empirical facts, is exponential inter-arrival time, or time between trades. Mainardi et al. [21, 22] studied the fractional Poisson process, where the exponential waiting time distribution is replaced by a Mittag-Leffler distribution, see also [4, 17, 28, 32]. Meerschaert et al. [23] showed that the same fractional Poisson process can also be obtained via an inverse stable time change. Using a time change in models for financial data has been popularized in the last decade based on the idea that it is not the calendar time that drives the changes in price, but rather the information flow or activity time is what matters for modeling of the prices [15, 14, 18]. In this paper we define fractional Skellam processes via time changes in ordinary Skellam processes. The resulting fractional Skellam models incor- porate Mittag-Leffler distribution of inter-arrival times and may provide better fit to high frequency financial data.
2. Preliminaries This section collects definitions and some results on the Skellam process, subordinators and the fractional Poisson process. These results will be used in the next section for the construction of the fractional Skellam processes. 2.1. Skellam processes Definition 2.1. A Skellam process is defined as S(t)=N (1)(t) − N (2)(t),t≥ 0, where N (1)(t), t ≥ 0andN (2)(t), t ≥ 0 are two independent homogeneous Poisson processes with intensities λ(1) > 0andλ(2) > 0, respectively.
The Skellam distribution has been introduced in [30] and [16], and the Skellam processes are considered in [3]. The probability mass function of S(t)isoftheform (1) k/2 (1) (2) λ √ −t(λ +λ ) (1) (2) sk(t)=P(S(t)=k)=e I|k| 2t λ λ , λ(2) (2.1) k ∈ Z = {0, ±1, ±2,...}, where Ik is the modified Bessel function of the first kind [31, p. 114] ∞ z/ 2n+k I z ( 2) . k( )= n n k n=0 !( + )! The Skellam process is a L´evy process with
E[exp{−θS(t)}]=exp{−tψS(1)(θ)}, 534 A. Kerss, N. N. Leonenko, A. Sikorskii and the L´evy exponent is given by ∞ −θy ψS(1)(θ)= (1 − e )ν(dy), −∞ where L´evy measure ν is the linear combination of two Dirac measures: (1) (2) ν = λ δ{1} + λ δ{−1}. It is easy to see that the moment generating function of the Skellam process is
(1) (2) (1) θ (2) −θ M(θ,t):=E[eθS(t)]=e−t(λ +λ −λ e −λ e ),θ∈ R. (2.2) Themeanandthevarianceare ES(t)=(λ(1) − λ(2))t, VarS(t)=(λ(1) + λ(2))t, (2.3) and the covariance function Cov(S(t),S(s)) = (λ(1) + λ(2))min(t, s),t,s>0. The next result on the Skellam processes is straightforward, but to the best of our knowledge, it has not appeared in the literature.
Lemma 2.1. The Skellam process is a stochastic solution of the fol- lowing system of differential equations: d s t λ(1) s t − s t − λ(2) s t − s t ,k∈ Z dt k( )= ( k−1( ) k( )) ( k( ) k+1( )) (2.4) with the initial conditions s0(0) = 1 and sk(0) = 0 for k =0 . The moment generating function of the Skellam process solves the dif- ferential equation dM θ,t ( ) M θ,t λ(1) eθ − λ(2) e−θ − ,θ∈ R dt = ( )( ( 1) + ( 1)) (2.5) with the initial condition M(θ,0) = 1.
P r o o f. Using the properties of the modified Bessel function [31, p. 115] Ik(x)=I−k(x) for any integer n and all x,and
dIν (z) 1 = (Iν−1(z)+Iν+1(z)) , dz 2 differentiate equation (2.1) to get d s t λ(1) s t − s t − λ(2) s t − s t ,k∈ Z. dt k( )= ( k−1( ) k( )) ( k( ) k+1( )) FRACTIONAL SKELLAM PROCESSES . . . 535
Thus (2.4) holds. Now multiply both sides of (2.4) by eθk and sum over k to get dM θ,t ( ) M θ,t λ(1) eθ − λ(2) e−θ − dt = ( )( ( 1) + ( 1)) with the initial condition M(θ,0) = 1, thus equation (2.5) holds. This equation clearly has solution (1) (2) (1) s (2) −s M(s, t)=e−t(λ −λ −λ e +λ e ),θ∈ R, which agrees with equation (2.2) above. 2 2.2. Subordinators and inverse subordinators AL´evy subordinator L(t), t ≥ 0 is a homogeneous positive nondecreas- ing L´evy process with Laplace transform E[e−θL(t)]=e−tφL(1)(θ),θ≥ 0, where the Laplace exponent ψL(1)(·)isgivenby −θx φL(1)(θ)=bθ + (1 − e )ν(dx) (0, ∞) with drift b ≥ 0andtheL´evy measure ν satisfying (0, ∞)(1∧x)ν(dx) < ∞. A standard α-stable subordinator D(t)isaL´evy subordinator with α Laplace transform E[e−θD(t)]=e−tθ , θ>0, t ≥ 0withα ∈ (0, 1). The inverse α-stable subordinator E(t) is defined as the inverse or first passage time of a stable subordinator D(t), that is E(t)=inf{u ≥ 0:D(u) >t},t≥ 0, see, for example, [7], [25, p. 101]. Note that E(t), t ≥ 0 is non-Markovian with non-stationary and non-independent increments. From [8], the mo- ments are tαkk! E[Ek(t)] = . (2.6) Γ(αk +1) The covariance function of this process is computed in [33], see also [19, equation (9)]: Cov[E(t),E(s)] = 1 min(t,s) (st)α ((t − τ)α +(s − τ)α) τ α−1dτ − ,t,s≥ 0. Γ(1 + α)Γ(α) 0 Γ(1 + α) (2.7) The Laplace transform of the inverse stable subordinator is −θE(t) α E[e ]=Eα(−θt ),θ>0,t≥ 0,α∈ (0, 1), (2.8) where 536 A. Kerss, N. N. Leonenko, A. Sikorskii