Skip to main content

ORIGINAL RESEARCH article

Front. Appl. Math. Stat., 11 May 2015
Sec. Mathematical Finance

Multivariate tempered stable model with long-range dependence and time-varying volatility

  • College of Business, The State University of New York at Stony Brook, Stony Brook, NY, USA

High-frequency financial return time series data have stylized facts such as the long-range dependence, fat-tails, asymmetric dependence, and volatility clustering. In this paper, a multivariate model which describes those stylized facts is presented. To construct the model, a multivariate ARMA-GARCH model is considered along with fractional Lévy process. The fractional Lévy process in this paper is defined by the stochastic integral with a tempered stable driving process. Parameters of the new model are fit to high-frequency returns for five U.S stocks. Approximated form of portfolio value-at-risk and average value-at-risk are provided and portfolio optimization is discussed under the model.

1. Introduction

The long-range dependence, fat-tail property, and volatility clustering effect are important issues for modeling high-frequency return time series in finance. The fractional Brownian motion introduced by Mandelbrot and Ness [1] can explain short-range or long-range dependence but it cannot explain the volatility clustering effect. The volatility clustering effect can be captured by the autoregressive conditional heteroskedastic (ARCH) and the generalized ARCH (GARCH) models formulated by Engle [2] and Bollerslev [3], respectively. However, GARCH models based on the normal distribution have not performed well in explaining high-frequency data analysis (see [4] and [5]), since the normal distribution does not capture the fat-tail property in the empirical innovation. To capture the fat-tail property in high-frequency return data, [5] suggested the univariate ARMA-GARCH model with tempered stable innovations without considering long-range dependence. The long-range dependence in high-frequency data is empirically reported in Sun et al. [6]. Sun et al. [4] provide a univariate model having long-range dependence, fat-tail property, and volatility clustering by taking the ARMA-GARCH model with fractional stable noise residuals exhibits, and show that the model has superior performance in high-frequency returns.

In this paper, a new multivariate market model that describes the long-range dependence, fat-tail property, and volatility clustering effect is developed. The new market model is constructed by taking the fractional tempered stable innovations on the multivariate ARMA-GARCH model. Univariate fractional tempered stable process was defined by the stochastic integral for the Volterra kernel in Houdre and Kawai [7] based on subclasses of Rosinski's tempered stable processes [8]. In order to construct multivariate model, we would better use normal tempered stable (NTS) process rather than Rosinski's tempered stable processes, since the NTS process is defined by the time-changed Brownian motion, and hence we can easily obtain a multivariate model which allows fatter tails than the multivariate Gaussian distribution and an asymmetric dependence structure. The NTS process has been discussed in many literatures including [912], and [13]. Although NTS is not included in Rosinski's tempered stable processes (see [14]), the fractional NTS process is redefined in this paper. The fractional NTS process is a multivariate process having the long-range dependence in time and asymmetric dependence between elements1.

To verify the performance of the new model, empirical illustration is provided using high-frequency stock return data. Useful simulation and parameter estimation methods are provided, and the goodness-of-fit tests are performed for the estimated parameters. The long-range dependence and fat-tail property of the high frequency stock return data are observed by this investigation.

To apply the new model for the financial risk management and portfolio management, portfolio value-at-risk (VaR), average VaR (AVaR), and portfolio optimization based on the new market model are discussed. VaR has a number of well-known limitations as a risk measure, nevertheless the VaR measure has been popularly used as a standard risk measure in the financial industry. The AVaR is the average of VaRs exceeding the VaR for a given confidence level2 AVaR is a superior alternative to VaR because it is a coherent risk measure3. and it is consistent with preference relations of risk-averse investors (see [20]).

A major contribution to the portfolio theory is the mean-variance model presented by Markowitz [21]. The importance of the model cannot be overstated, but some of assumptions underlying the model have been challenged since its introduction. One of the assumptions is that asset returns follow a Gaussian (normal) distribution and another is that the variance is a measure of risk ignoring higher-order moments. In this paper, the portfolio optimization is discussed based on the Markowitz's theory, but the Gaussian assumption is replaced by the ARMA-GARCH model with fractional NTS innovations and the variance risk measure is superseded by VaR and AVaR risk measures.

The remainder of this paper is organized as follows. The NTS process is reviewed in Section 2. The definition of the multivariate fractional NTS process and its simulation method are presented in Sections 3, 4, respectively. Multivariate ARMA-GARCH model having long-range dependence is defined in Section 5 along with the empirical illustration. In Section 6, portfolio risk will be assessed and the optimal portfolio will be found based on the new ARMA-GARCH model using high-frequency market data. In Section 7, the principal findings are summarized. Volterra kernel is briefly reviewed in the appendix.

2. Normal Tempered Stable Process

In the remainder of this paper, we assume that N is a positive integer standing for the dimension and α ∈ (0, 2), θ > 0, β = (β1, β2, …, βN)T, γ = (γ1, γ2, …, γN)T such that γn > 0 for n ∈ {1, 2, …, N}, and Σ = [σm,n]m,n∈{1, 2, …, N} is a given correlation matrix such that σn,n = 1 for n ∈ {1, 2, …, N}. The pure jump Lévy process T = (T(t))t ≥ 0 whose characteristic function ϕT(t) is equal to

ϕT(t)(u)=exp(θ1α2tΓ(1α2)0(eiux1)eθxxα/2+1dx).

is a subordinator and referred to as the tempered stable subordinator with parameters (α, θ). Solving the integration in the last Equation, we can obtain the following formula,

ϕT(t)(u)=exp(2θ1α2tα((θiu)α2θα2)).(1)

Consider a N-dimensional Brownian motion B = (B(t))t ≥ 0 such that B(t) = (B1(t), B2(t), …, BN(t))T, and suppose that

cov(Bm(t),Bn(t))=σm,nt

for all m, n ∈ {1, 2, …, N}.

Suppose T is independent of B. Consider a N-dimensional process Z = (Z(t))t ≥ 0 such that Z(t) = (Z1(t), Z2(t), …, ZN(t))T. For n ∈ {1, 2, …, N}, define (Zn(t))t ≥ 0 by the time-changed Brownian motion as

Zn(t)=βn(T(t)t)+γnBn(T(t)).(2)

Then the process Z is referred to as the N-dimensional NTS process with parameters (α, θ, β, γ, Σ) and we denote by Z~ NTSN (α, θ, β, γ, Σ).

By composing characteristic functions of B(t) and T(t), we obtain the characteristic function of Zn(t) as follows:

ϕZn(t)(u)=       exp(βnuti2θ1α2αt((θiβnu+γn2u22)α2θα2)).(3)

The mean of Zn(t) is equal to zero for n ∈ {1, 2, …, N}. Covariance between Zm(t) and Zn(t) is given by

cov(Zm(t),Zn(t))=tγmγnσm,n+tβmβn(2α2θ)(4)

for m, n ∈ {1, 2, …, N}. Moreover the variance of Zn(t) is equal to

Var(Zn(t))=t(γn2+βn2(2α2θ))(5)

The linear combination of elements of Z is also NTS as follows.

Proposition 1. Let w = (w1, w2, …, wN)T ∈ ℝN, Z ~ NTSN (α, θ, β, γ, Σ), and Y(t) = Y(t)=n=1NwnZn(t). Then (Y(t))t ≥ 0~ NTS1(α, θ, β, γ, 1), where

β¯=n=1Nwnβn,and γ¯=m=1Nn=1Nwmwnγmγnσm,n.

Proof. See [12].

If γn=1βn2(2α2θ) and |βn|<2θ2α for n ∈ {1, 2, …, N} then var(Zn(t)) = t. In this case the process Z is referred to as the N-dimensional standard NTS process with parameters (α, θ, β, Σ) and denoted by Z ~ stdNTSN (α, θ, β, Σ).

3. Fractional Normal Tempered Stable Process

Let KH(t, s) be the Voltera kernel and Z ~ NTSN (α, θ, β, γ, Σ). The N-dimensional fractional NTS (fNTS) process generated by Z is defined by the process of vector X = (X(t))t ≥ 0 with X(t) = (X1(t), X2(t), …, XN(t))T such that

Xn(t)=lim||P||0j=1MKH(t,tj1)(Zn(tj)Zn(tj1))

in distribution sense for n ∈ {1, 2, …, N}, where

P:0=t0<t1<<tM=t

is a partition of the interval [0, t] and

||P||=max{tjtj1|j=1,2,,M}.

In this case we denote that Xn(t)=0tKH(t,s)dZn(s) and X~fNTSN(H,α,θ,β,γ,Σ). Since we have

E[(iuj=1MKH(t,tj1)(Zn(tj)Zn(tj1)))]=j=1ME[exp(iuKH(t,tj1)(Zn(tj)Zn(tj1)))]=j=1Mexp(βnuiKH(t,tj1)(tjtj1)2θ1α2α((θiβnuKH(t,tj1)+γn2u2(KH(t,tj1))22)α2θα2)(tjtj1))=exp(βnuij=1MKH(t,tj1)(tjtj1)2θ1α2αj=1M((θiβnuKH(t,tj1)+γn2u2(KH(t,tj1))22)α2θα2)(tjtj1)),

the characteristic function of Xn(t) is given by

ϕXn(t)(u)=lim||P||0exp(βnuij=1MKH(t,tj1)(tjtj1)2θ1α2αj=1M((θiβnuKH(t,tj1)+γn2u2(KH(t,tj1))22)α2θα2)(tjtj1))=exp(βnui0tKH(t,s)ds2θ1α2α0t((θiβnuKH(t,s)+γn2u2(KH(t,s))22)α2θα2)ds).

Proposition 2. For n ∈ {1, 2, …, N}, the covariance between Xn(s) and Xn(t) is equal to

cov(Xn(s),Xn(t))=12(γn2+βn2(2α2θ))(t2H+s2H|ts|2H),   s,t>0.(6)

Proof. Let P be a partition such that

P:  0=t0<t1<<tM1<tM=st<tM+1<<tM*=st.

Then we have

cov(Xn(s),Xn(t))=E[Xn(s)Xn(t)]  =lim||P||0E[j=1MKH(s,tj1)(Zn(tj)2ptZn(tj1))     k=1M*KH(t,tk1)(Zn(tk)2ptZn(tk1))]  =lim||P||0j=1Mk=1M*KH(s,tj1)KH(t,tk1)    E[(Zn(tj)Zn(tj1))(Zn(tk)Zn(tk1))]

By the property of the NTS process Zn, we have

E[(Zn(tj)Zn(tj1))(Zn(tk)Zn(tk1))]={(tjtj1)Var(Zn(1)) if j=k0 if jk

Hence, we obtain

cov(Xn(s),Xn(t))=lim||P||0j=1MKH(s,tj1)KH(t,tj1)(tjtj1)Var(Zn(1))=Var(Zn(1))0stKH(s,u)KH(t,u)du.

Hence we obtain Equation (6) by Equation (5) and Equation (21) in Appendix.

Proposition 3. For m, n ∈ {1, 2, …, N}, the covariance between Xm(t) and Xn(t) is equal to

cov(Xm(t),Xn(t))=t2H(σm,nγmγn+βmβn(2α2θ)),t>0.(7)

Proof. Let P be a partition such that

P:0=t0<t1<<tM1<tM=t.

We have

cov(Xm(t),Xn(t))=E[Xm(t)Xn(t)]=lim||P||0E[j=1MKH(t,tj1)(Zm(tj)Zm(tj1)) k=1MKH(t,tk1)(Zn(tk)Zn(tk1))]=lim||P||0j=1Mk=1MKH(t,tj1)KH(t,tk1)E[(Zm(tj)Zm(tj1))(Zn(tk)Zn(tk1))]

By the property of the NTS process Zn, we have

E[(Zm(tj)Zm(tj1))(Zn(tk)Zn(tk1))]={(tjtj1)cov(Zm(1),Zn(1)) if j=k0 if jk

Hence, we obtain

cov(Xm(t),Xn(t))=cov(Zm(1),Zn(1))0t(KH(t,u))2du.

Hence we obtain Equation (7) by Equation (4) and Equation (22) in Appendix.

For a given stochastic process Y = (Y(t))t ≥ 0, the summation

j=1E[(Y(1)Y(0))(Y(j+1)Y(j))]

diverges, then we say that Y exhibits long-range dependence (See [22]). By Proposition 2 and L'Hopital's rule, we have

E[Xn(1)(Xn(j+1)Xn(j))]=v2((j+1)2H2j2H+(j1)2H)=v2j2H2(j2((1+1j)2H2+(11j)2H))vH(2H1)j2H2 as j,

where v=(γn2+βn2(2α2θ)). Hence, j=1E[Xn(1)(Xn(j+1)Xn(j))] diverges, i.e., the process (Xn(t))t ≥ 0 has long-range dependence, when ½ < H < 1.

Since the NTS has an asymmetric dependence structure, X has also asymmetric dependence. By Proposition 1, we can prove the following proposition.

Proposition 4. Let w = (w1, w2, …, wN)T ∈ ℝN, X ~ fNTSN (H, α, θ, β, γ, Σ), and Y(t)=n=1NwnXn(t). Then (Y(t))t ≥ 0 ~ fNTS1(H, α, θ, β, γ, 1), where

β¯=n=1Nwnβn,and γ¯=m=1Nn=1Nwmwnγmγnσm,n.

When Z ~ stdNTSN (α, θ, β, Σ), the multivariate fractional NTS process X generated by Z is referred to as the fractional standard NTS process, and we denote that X ~ fstdNTSN (H, α, θ, β, Σ). In this case, we have

cov(Xn(s),Xn(t))=12(t2H+s2H|ts|2H),   s,t>0.(8)

and

cov(Xm(t),Xn(t))=t2H(σm,n1βm2(2α2θ)1βn2(2α2θ)  +βmβn(2α2θ)),   t>0.(9)

By Proposition 4, we can prove the following corollary.

Corollary 5. Let w = (w1, w2, …, wN)T ∈ ℝN, X ~ fstdNTSN (H, α, θ, β, Σ), and Y(t)=n=1NwnXn(t). Then (Y(t))t ≥ 0 ~ fNTS1(H, α, θ, β, γ, 1), where

β¯=n=1Nwnβn,

and

γ¯=m=1Nn=1Nwmwnσm,n1βm2(2α2θ)1βn2(2α2θ).

4. Simulation

In this section, a numerical method is provided to generate the sample path for the multivariate fractional normal tempered stable process. By Theorem 5.3 (i) in Rosiński [8]4, we obtain series representation for the tempered stable subordinator T as follows:

T(t)=limMj=1M1(0,t)(τj)((αξjΓ(α2)2θ1α2T)2αejuj2αθ),t[0,T],(10)

where

• {uj} is an iid sequence of random variables on (0, 1),

• {ej} and {ej} are an iid sequence of exponential random variables with parameter 1,

• ξj = e1 + e2 + … + ej,

• {τj} be an independent and identically distributed uniform random variable in [0, yes], where yes > 0 is fixed.

• and assume that {uj}, {ej}, {ej} and {τi} are independent.

Let LΣ be the lower triangular matrix obtained by the Cholesky decomposition for Σ with Σ = LΣLTΣ, where Σ is the correlation matrix in Equation (2). Then we have B(t) = LΣ B(t) where B(t) = (B1(t), B2(t), …, BN(t))T is a mutually independent vector of Brownian motions.

Sample path of Z ~ NTSN (α, θ, β, γ, Σ) is generated as follows. For a given partition {t0, t1, …, tM} of the interval [0, yes] with t0 = 0, tM = yes and tj < tk for j < k, we have

B¯n(T(tj))B¯n(T(tj1))=T(tj)T(tj1)ϵj,n, n{1,2,,N}

where ϵj,n ~ N(0, 1). Therefore, we have

Z(tk)=(T(tk)tk)β+diag(γ)j=1kT(tj)T(tj1)LΣϵj,(11)

where ϵj = (ϵj, 1, ϵj, 2 …, ϵj, N)T, ϵj,n ~ N(0,1), and ϵj,m is independent of ϵj,n for all m, n ∈ {1, 2, …, N}, and j ∈ {1, 2, …, M}.

Finally, sample paths of X ~ fNTSN (H, α, θ, β, γ, Σ) is generated as follows:

X(tk)=j=0k1KH(tk,tj)(Z(tj+1)Z(tj)),     k1,2,,M.(12)

Figure 1 presents one example sample path of the tempered stable subordinator T with parameters α = 1.2 and θ = 0.8. Figure 2 exhibits one pair of simulated sample path of the 2 dimensional standard fractional NTS process X = (X1(t), X2(t))t ≥ 0 with parameters α = 1.2, θ = 0.8, β = (0.1,−0.3), σ1, 2 = σ2, 1 = 0.7, and H = 0.6525. In the process Figure 2A (X1(t))t ≥ 0 and (X2(t))t ≥ 0 are drawn on the time line. Figure 2B presents 2 dimensional movements of X.

FIGURE 1
www.frontiersin.org

Figure 1. Simulated sample path of the tempered stable subordinator T with parameters α = 1.2 and θ = 0.8.

FIGURE 2
www.frontiersin.org

Figure 2. Simulated sample path of the 2 dimensional standard fractional NTS process X = (X1(t), X2(t))t ≥ 0 with parameters α = 1.2, θ = 0.8, β = (0.1, −0.3), σ1, 2 = σ2, 1 = 0.7, and H = 0.6525.

5. ARMA-GARCH Model with fNTS Innovations and Empirical Illustration

Let X ~ fstdNTSN (H, α, θ, β, Σ) generated by Z ~ stdNTSN (α, θ, β, Σ). A N-dimensional discrete time process Y = (Y(k))k ∈ {0, 1, 2, … } with Y(k) =(Y1(k), Y2(k), …, YN(k)) is referred to as the N-dimensional ARMA-GARCH model with fNTS innovations when it is given by the ARMA(1,1) - GARCH(1,1) model as follows: Yn(0) = 0, εn(0) = 0, and

{ Yn(k+1)=cn+anYn(k)+bnσn(k)εn(k)+σn(k+1)εn(k+1)(σn(k+1))2=κn+ξn(σn(k)εn(k))2+ζn(σn(k))2 ,(13)

where εn(k + 1) = Xn(k + 1) − Xn(k) and n ∈ {1, 2, …, N}. This model describes volatility clustering effect by GARCH(1,1) model, the fat-tails and the asymmetric dependence between elements by the standard NTS process Z, and the long-range dependence by fractional NTS process X.

Since we have

Xn(k)j=0k1KH(k,j)(Zn(j+1)Zn(j))

The increment Xn(k + 1)−Xn(k) can be approximated as follows:

Xn(k+1)Xn(k)KH(k+1,k)(Zn(k+1)Zn(k))+j=0k1(KH(k+1,j)KH(k,j))(Zn(j+1)Zn(j))(14)

Let M be the number of time steps in the sample and N be the number of assets in the portfolio. ARMA(1,1)-GARCH(1,1) model is fit to the data and extract εn(tk), for n = 1, 2, …, N and k = 1, …, M. Since Equation (8) is the same as the covariance of the fractional Brownian motion, set Wn(tk)=j=1kεn(tj) and estimate the Hurst index Hn of the process (Wn(tk))tk ≥ 0 using the wavelet details regression estimator method by Flandrin [23] and Abry et al. [24]. The parameter H is obtained finally as the mean of Hn for n = 1, 2, …, N.

Suppose we have estimated Hurst index H. We estimate the parameters of the model as follows in this investigation.

1. Estimate ARMA(1,1)-GARCH(1,1) parameters an, bn, cn, κn, ξn, and ζn with standard normal innovations by maximum likelihood estimation (MLE) with assumption (σn(0))2 = κn/(1 − ξn − ζn) for n = 1, 2, …, N.

2. Extract residuals using the estimated parameters.

3. Put Xn(tk)=j=1kεn(tj) and extract {Zn(tk)| k = 1, 2, …, M} for n = 1, 2, …, N as follows.

Zn(1)=Xn(1)/KH(1,0) andZn(k)=Zn(k1)+Xn(k)Xn(k1)KH(k,k1)j=0k2KH(k,j)KH(k1,j)KH(k,k1)(Zn(j+1)Zn(j))(15)

for k = 2, 3, …, M.

4. Estimate parameters αn, θn, and βn of the standard NTS process using {Zn(k)| k = 1, 2, …, M} extracted in the step 4 by curve-fitting in least-squares sense. Set α=n=1NαnN and θ=n=1NθnN Estimate parameters βn again using {Zn(k)| k = 1, 2, …, M} by means of MLE for n = 1, 2, …, N.

5. Calculate the covariance between (Zm(1)) and (Zn(1)) for m,n ∈ {1, 2, …, N} using data {(Zm(k), Zn(k))| k =1, 2, …, M} extracted in the step 4. Estimate

Σ=[σm,n]m,n{1,2,,N}

by Equation (4) and cov(Zm(1), Zn(1)).

The parameter H is estimated using 2,158 observed 1 min returns for five stocks (Apple Inc., Google Inc., IBM Co., A&T, Wal-Mart Stores Inc.) from November 21 to November 29, 2011. Set dt = 1/390, since 1 day has 6 h and 30 min (from 9:30 to 16:00) trading time, that is 390 min, on New York Stock Exchange. The parameters Hn are presented in Table 1 and parameter H = 0.5806 is obtained finally as the mean of Hn.

TABLE 1
www.frontiersin.org

Table 1. Estimated Hurst index Hn.

The ARMA(1,1)-GARCH(1,1) parameters and estimated fNTS parameters α, θ and β are reported in Table 2 for each stocks. The estimated matrix Σ is presented in Table 3.

TABLE 2
www.frontiersin.org

Table 2. Parameters the ARMA-GARCH model with (non-fractional) NTS innovations and fNTS innovations.

TABLE 3
www.frontiersin.org

Table 3. Estimated Σ = [σm,n]m,n ∈ {1, 2, …, N}.

The Kolmogorov-Smirnov (KS) and Anderson-Darling (AD) tests are used for goodness-of-fit tests. The KS and AD2 statistic are given by

KS=x|F^(x)F(x)|,   AD2=s(F^(x)F(x))2F(x)(1F(x))dF(x),

where F^(x) is the empirical sample distribution, F(x) is the estimated theoretical distribution, and s is the number of observed samples. According to the p-values of KS statistic and AD2 statistic, estimated parameters are not rejected. Calculating p-values for KS and AD2 statistic are explained in Marsaglia et al. [25] and Marsaglia and Marsaglia [26]. The KS and AD2 statistics are calculated between the standard NTS distribution with estimated parameters (α, θ, βn), and the empirical cumulative distribution of {Zn(tk) − Zn(tk − 1)| k = 1, 2, …, M} where Zn(tk) is the extracted process by Equation (15). According to the p-values of KS and AD2 statistic in Table 4, all estimated parameters are not rejected at the 1% significance level for the five stock returns investigated.

TABLE 4
www.frontiersin.org

Table 4. Goodness-of-fit test for the ARMA-GARCH model with normal innovations.

If H = 0.5 then the ARMA-GARCH model with fNTS innovations becomes the ARMA-GARCH model with non-fractional NTS innovations. The estimated parameters of the non-fractional standard NTS process are presented also in Table 2 for each stocks. To compare the performance of parameter estimation for ARMA-GARCH model with fNTS innovations with non-fractional NTS innovation, we provide KS and AD2 statistic with p-values for the ARMA-GARCH model with non-fractional NTS innovations in Table 4 as well. Unfortunately, the estimated parameters for IBM Co. and AT&T are rejected by KS and AD tests at the 1% significance level.

If we assume that (εn(tk))k = 1, 2, …, M is independent and identically distributed and εn(tk) ~ N(0, 1), then the Yn follows the ARMA-GARCH model with normal innovations. The ARMA-GARCH model with normal innovations is rejected by both KS test and and AD test for all considered stocks at the 1% significance level.

Therefore, we can conclude that the ARMA-GARCH model with fNTS innovations describes behavior of the high frequency return time series investigated in this section. While, the ARMA-GARCH models with non-fractional NTS innovations and normal innovations do not perfectly explain the behavior of the high frequency return data.

6. Assessment Risk on the ARMA-GARCH Model with fNTS Innovations

In this section, portfolio VaR and AVaR on the ARMA-GARCH model with fNTS innovations are discussed, and they are applied to the portfolio optimization.

Let yes > 0 be a time horizon. Assume that we have a portfolio with N stocks. Suppose the stock return vector process Y = (Y(tk))k ∈ {0, 1, 2, …, M} with Y(tk) = (Y1(tk), Y2(tk), …, YN(tk)) follows the ARMA-GARCH model with fNTS innovations as defined in Section 5 and X ~ fstdNTSN (H, α, θ, β, Σ) generated by Z~ stdNTSN (α, θ, β, Σ), where {t0, t1,…, tM} is a given discrete time such that tk = k · Δt with Δt =yes/M for k ∈ {0, 1, 2, …, M}. Then a portfolio return process R = (R(tk))k ∈ {0, 1, 2, …, M} with allocation weight vector w = (w1, w2, …, wN)T with n=1Nwn=1 is given by R(tk)=n=1NwnYn(tk).

We have

R(k+1)=n=1NwnYn(k+1)=n=1Nwn(cn+anYn(k)+bnσn(k)εn(k))+n=1Nwnσn(k+1)εn(k+1)(16)

and

n=1Nwnσn(k+1)εn(k+1)=n=1Nwnσn(k+1)Xn(k+1)Xn(k)(tk+1tk)H=(Δt)Hn=1Nwnσn(k+1)(Xn(k+1)Xn(k)).

By Equation (14), we obtain

n=1Nwnσn(k+1)εn(k+1)(Δt)HKH(tk+1,tk)n=1Nwnσn(k+1)(Zn(k+1)Zn(k))+(Δt)Hn=1Nj=0k1wnσn(k+1)(KH(tk+1,tj)KH(tk,tj))(Zn(tj+1)Zn(tj))(17)

Let (yes(tk))k ∈ {1, 2, …, M} be the natural filtration generated by Y. Then σn(k + 1) and

n=1Nj=0k1wnσn(k+1)(KH(tk+1,tj)KH(tk,tj))(Zn(tj+1)Zn(tj))

are yes(tk)-measurable. Moreover, since Zn has stationary increments, we have

yes

Hence, we have

yes

where (Ξ(t))t ≥ 0 ~ NTS1(α, θ, β, γ, 1) with

β¯=n=1Nwnσn(k+1)βn, γ¯=m=1Nn=1Nwmwnσm(k+1)σn(k+1)σm,nγmγn,

and

γn=1βn2(1α2θ),    n=1,2,,N,

by Proposition 1. Therefore, we can discuss VaR and AVaR as follows. The VaR and AVaR for R(k + 1) with the significance level η under information until time tk are defined by

yes

and

yes

respectively. By Equation (16), we have

yes

and by Equation (17), we obtain

yes

Thus, by Equation (18) we obtain

yes

By the same argument, we obtain AVaR as follows

yes

The closed-form solutions of VaRη (Ξ(Δt)) and AVaRη (Ξ(Δt)) for NTS process (Ξ(t))t ≥ 0 is presented in Kim et al. [27].

The Figures 3A,B exhibit the forecasted 1 min ahead VaR and AVaR values, respectively, for an equally weighted portfolio for the five stocks in this study based on the parameters in Table 2 and information we discussed in the Section 5. The equally weighted portfolio is the portfolio having allocation weight as w = (1/N, 1/N, …, 1/N), where N is the number of stocks in the portfolio. The forecasted 1-min VaR and AVaR values are calculated by Equation (19) and Equation (20), respectively, at confidence levels from 0.5 to 10%. The VaR and AVaR values for the portfolio based on the ARMA-GARCH model with non-fractional NTS innovations are presented in the figure. The VaR (AVaR) values of the ARMA-GARCH model with fNTS innovations are larger than VaR (AVaR) values of the model with non-fractional NTS innovations.

FIGURE 3
www.frontiersin.org

Figure 3. Forecasted 1 min (A) VaR and (B) AVaR under the ARMA-GARCH model with non-fractional NTS and fNTS innovations. VaR and AVaR values for the ARMA-GARCH model with fNTS innovations are given by solid curves, and those values for the ARMA-GARCH model with non-fractional NTS innovations are given by dash curves.

7. Portfolio Optimization and Risk Budgeting on the ARMA-GARCH Model with fNTS Innovations

Using VaR and AVaR values by Equation (19) and Equation (20), we can find the optimal portfolio. The VaR minimizing portfolio is obtained by solving the following optimization problem:

yes

where μ0 is the benchmark expected return. Similarly, the AVaR minimizing portfolio is obtained by solving the following optimization problem:

yes

Table 5 presents the VaR and AVaR minimizing portfolios with the benchmark return μ0 = 2.5924 · 10−4. To measure performance of the optimal portfolio, we use the VaR ratio and stable tail-adjusted return ratio (STARR), which are defined by

yes

and

yes

respectively5. The following is concluded from the results reported in Table 5:

• The VaR minimizing portfolio has the best expected return among the three portfolios.

• The AVaR minimizing portfolio has smaller VaR than the case of the VaR minimizing portfolio, but that is not surprising since the AVaR minimizing portfolio has less expected return than the case of the VaR minimizing portfolio.

• The AVaR minimizing portfolio has the largest VaR Ratio among the three VaR ratios.

• The AVaR minimizing portfolio has the largest STARR among the three STARR.

TABLE 5
www.frontiersin.org

Table 5. Optimal allocation weight for the portfolio with the five stocks and performance measures.

8. Conclusion

The multivariate ARMA-GARCH model with fNTS innovations exhibits fat-tails, asymmetric dependence, volatility clustering, and long-range dependence. Comparing with two ARMA-GARCH models with non-fractional NTS innovation and normal innovations, the ARMA-GARCH model with fNTS innovations has better performance in parameter estimation for 1-min stock return data investigated in this paper. That means the fNTS process describes the behavior of the residual process of 1-min returns better than the non-fractional NTS process or Brownian motion. The portfolio VaR and AVaR are calculated by the approximation method under the model, and those risk measures are used for portfolio optimization. In this investigation, we obtain the fact that the risk measures of the ARMA-GARCH model with fNTS innovations are more conservative than those of the model with non-fractional NTS innovations. The AVaR minimizing portfolio performs better than the VaR minimizing portfolio for the case considered in this paper.

Conflict of Interest Statement

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.

Footnotes

1. ^Kim [15] defines another different multivariate fractional NTS process using the time-changed fractional Brownian motion.

2. ^AVaR is also known as conditional value-at-risk (CVaR). See Pflug [16] and Rockafellar and Uryasev [17, 18].

3. ^VaR is not is not a coherent risk measure. The notion of a coherent risk measure was introduced by Artzner et al. [19].

4. ^The tempered stable subordinator T with parameter (α, θ) is included in the class of tempered stable processes provided by Rosinski [8]. The parameter α/2 of the tempered stable subordinator corresponds to the parameter α in Rosinski's tempered stable class.

5. ^Many literatures define the VaR ratio and STARR as

yes

and

yes

respectively, where Rf is a market index return. In this paper, we assume Rf = 0 for considering absolute performance.

References

Mandelbrot BB, Ness JWV. Fractional brownian motions, fractional noises applications. SIAM Rev. (1968) 10, 422–37. doi: 10.1137/1010093

CrossRef Full Text | Google Scholar

Engle R. Autoregressive conditional heteroskedasticity with estimates of the variance of United Kingdom inflation. Econometrica (1982) 50, 987–1007. doi: 10.2307/1912773

CrossRef Full Text | Google Scholar

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

Sun W, Rachev ST, Fabozzi FJ. Long-range dependence, fractal processes, and intra-daily data. In: Seese D, Weinhardt C, Schlottmann F. editors. Handbook on Information Technology in Finance (2008). (Springer), p. 543–85.

Beck A, Kim YS, Rachev ST, Feindt M, Fabozzi FJ. Empirical analysis of ARMA-GARCH models in market risk estimation on high-frequency U.S. data. Stud Nonlin Dyn Econometri. (2012). 17:167–77. doi: 10.1515/snde-2012-0033

CrossRef Full Text | Google Scholar

Sun W, Rachev ST, Fabozzi FJ. Fractals or I.I.D.: evidence of long-range dependence and heavy tailedness from modeling German equity market returns. J Econ Bus. (2007) 59:575–95. doi: 10.1016/j.jeconbus.2007.02.001

CrossRef Full Text | Google Scholar

Houdre C, Kawai R. On fractional tempered stable motion. Stochast Process Appl. (2006) 116:1161–84. doi: 10.1016/j.spa.2006.01.008

CrossRef Full Text | Google Scholar

Rosiński J. Tempering stable processes. Stochast Proce Appl. (2007) 117:677–707. doi: 10.1016/j.spa.2006.10.003

CrossRef Full Text

Barndorff-Nielsen OE, Shephard N. Normal Modified Stable Processes. Economics Series Working Papers from University of Oxford, Department of Economics. (2001) 72.

Google Scholar

Barndorff-Nielsen OE, Levendorskii S. Feller processes of normal inverse gaussian type. Quantitat Finan. (2001) 1:318–31. doi: 10.1088/1469-7688/1/3/303

CrossRef Full Text | Google Scholar

Eberlein E, Madan DB. On correlating lévy processes. J Risk. (2010) 13:3–16.

Kim YS, Giacometti R, Rachev ST, Fabozzi FJ, Mignacca D. Measuring financial risk and portfolio optimization with a non-Gaussian multivariate model. Anna Operat Res. (2012) 201:325–43. doi: 10.1007/s10479-012-1229-8

CrossRef Full Text | Google Scholar

Kim YS, Volkmann D. Normal tempered stable copula. Appl Math Lett. (2013). 26:676–80. doi: 10.1016/j.aml.2013.01.009

CrossRef Full Text | Google Scholar

Kim YS, Rachev ST, Chung DM, Bianchi ML. The modified tempered stable distribution, GARCH-models and option pricing. Probabil Math Stat. (2009) 29:91–117.

Google Scholar

Kim YS. Fractional multivariate normal tempered stable process. Appl Math Lett. (2012) 25:2396–401. doi: 10.1016/j.aml.2012.07.011

CrossRef Full Text | Google Scholar

Pflug G. Some Remarks on the value-at-risk and the Conditional value-at-risk. In: Uryasev S, editor. Probabilistic Constrained Optimization: Methodology and Applications. Kluwer Academic Publishers (2000). p. 272–281.

Rockafellar RT, Uryasev S. Optimization of conditional value-at-risk. J Risk (2000) 2:21–41.

Google Scholar

Rockafellar RT, Uryasev S. Conditional value-at-risk for general loss distributions. J Bank Finance. (2002) 26:1443–71. doi: 10.1016/S0378-4266(02)00271-6

CrossRef Full Text | Google Scholar

Artzner P, Delbaen F, Eber JM, Heath D. Coherent measures of risk. Math Finance. (1999) 9:203–28. doi: 10.1111/1467-9965.00068

CrossRef Full Text | Google Scholar

Rachev ST, Stoyanov S, Fabozzi FJ. Advanced Stochastic Models, Risk Assessment, and Portfolio Optimization: The Ideal Risk, Uncertainty, and Performance Measures. Hoboken, NJ: John Wiley & Sons (2007).

Markowitz HM. Portfolio selection. J Finance. (1952) 7:77–91.

Google Scholar

Samorodnitsky G, Taqqu MS. Stable Non-Gaussian Random Processes. Chapman & Hall/CRC (1994). (Baca Raton, London, New York, Washinton, DC).

Google Scholar

Flandrin P. Wavelet analysis and synthesis of fractional Brownian motion. IEEE Trans Inf Theory. (1992) 38:910–17. doi: 10.1109/18.119751

CrossRef Full Text | Google Scholar

Abry P, Flandrin P, Taqqu MS, Veitch D. Self-similarity and long-range dependence through the wavelet lens. In: Theory and Applications of Long-Range Dependence. Birkhäuser (2003). eds D. Paul, O. George and T. Murad p. 527–556.

Marsaglia G, Tsang WW, Wang G. Evaluating kolmogorov's distribution. J Stat Softw. (2003) 8:1–4.

Marsaglia G, Marsaglia J. Evaluating the anderson-darling distribution. J Stat Softw. (2004) 9:1–5.

Google Scholar

Kim YS, Rachev ST, Bianchi ML, Fabozzi FJ. Computing VaR and AVaR in infinitely divisible distributions. Probab Math Stat. (2010) 30:223–45.

Google Scholar

Nualart D. Stochastic integration with respect to fractional Brownian motion and applications. Contemp Math. (2003) 336:3–39. doi: 10.1090/conm/336/06025

CrossRef Full Text | Google Scholar

Appendix

To define a process with long-range dependence, we use the Volterra kernel KH:[0, ∞) × [0, ∞) → [0, ∞), given by

cH=(H(12H)Γ(12H)Γ(22H)Γ(H+12))12.

with

0tsKH(t,u)KH(s,u)du=12(t2H+s2H|ts|2H),

and H ∈ (0, 1). According to Nualart [28] and Houdre and Kawai [7], we have the following facts:

• We have

0tKH(t,s)2ds=t2(H12)+1.

and

0tKH(t,s)2ds=t2(H12)+1.(22)

• If H ∈ (½, 1) then

KH(t,s)=cH(H12)s12Hst(us)H32uH12du1[0,t](s).

Moreover, K½(t, s) = 1[0, t](s).

• Let t > 0 and let p ≥ 2. KH(t, ·) ∈ Lp([0, t]) if and only if !!i!!. When KH(t, ·) ∈ Lp([0, t]), we have

0tKH(t,s)pds=CH,ptp(H12)+1,

where

CH,p=cHp01vp(12H)×[(1v)H12(H12)v1wH32(wv)H12dw]pdv.(23)

Keywords: multivariate fractional normal tempered stable process, long-range dependence, fractional Brownian motion, fractional, Lévy processes, high-frequency market, intraday trading, volatility clustering, asymmetric dependence

JEL: C13, C32, C58, C61, G11, G32

Citation: Kim YS (2015) Multivariate tempered stable model with long-range dependence and time-varying volatility. Front. Appl. Math. Stat. 1:1. doi: 10.3389/fams.2015.00001

Received: 16 January 2015; Accepted: 05 February 2015;
Published online: 11 May 2015.

Edited by:

Yong Hyun Shin, Sookmyung Women's University, South Korea

Reviewed by:

Kum-Hwan Roh, Hannam University, South Korea
Younhee Lee, Chungnam National University, South Korea

Copyright © 2015 Kim. 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) or licensor 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: Young Shin Kim, College of Business, Stony Brook University, Harriman Hall 315, Stony Brook, 11794-3775 NY, USA aaron.kim@stonybrook.edu

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.