<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMF</journal-id><journal-title-group><journal-title>Journal of Mathematical Finance</journal-title></journal-title-group><issn pub-type="epub">2162-2434</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmf.2015.51002</article-id><article-id pub-id-type="publisher-id">JMF-53602</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Business&amp;Economics</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Modeling Returns and Unconditional Variance in Risk Neutral World for Liquid and Illiquid Market
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>vivi</surname><given-names>Joseph Mwaniki</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>School of Mathematics, Statistics, Actuarial and Finance Division, University of Nairobi, Nairobi, Kenya</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>jimwaniki@uonbi.ac.ke</email></corresp></author-notes><pub-date pub-type="epub"><day>20</day><month>01</month><year>2015</year></pub-date><volume>05</volume><issue>01</issue><fpage>15</fpage><lpage>25</lpage><history><date date-type="received"><day>6</day>	<month>January</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>25</month>	<year>January</year>	</date><date date-type="accepted"><day>28</day>	<month>January</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  This article seeks to model daily asset returns using log-ARCH-L&#233;vy type model which is expected to reproduce most of the stylized features of financial time series data (such as volatility clustering, leptokurtic nature of log returns, joint covariance structure and aggregational Gaussianity) that are empirically found in different types of market. In addition, unconditional variance of daily log returns in risk neutral world of different conditional heteroscedastic models is derived. A key observation is that liquid markets and illiquid market may not have the same underlying dynamics. For instance empirical analysis based on S&amp;P500 index log returns as a liquid market do not have autoregressive part in their first moments while in Nairobi Securities Exchange NSE20 index there is strong presence of autoregressive dynamics of order three,
  <em> i.e.</em> AR(3). Higher moments of both markets are serially correlated.
 
</p></abstract><kwd-group><kwd>AR-APARCH</kwd><kwd> L&amp;eacute;vy Increments</kwd><kwd> Generalized Hyperbolic Distribution</kwd><kwd> Normal Inverse Gaussian</kwd><kwd>  Illiquid Market</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>It is well known that the stock price changes are neither independent nor identically distributed. There are linear and nonlinear dependencies between successive price changes. Distributional assumptions concerning risky asset log returns play a key role in option pricing. According to research finding of Mandelbrot [<xref ref-type="bibr" rid="scirp.53602-ref1">1</xref>] , evidence indicates that the empirical distributions of daily stock returns differ significantly from the traditional Gaussian model. In search of satisfactory descriptive models for financial data, large number of distributions have been tried (see for example, [<xref ref-type="bibr" rid="scirp.53602-ref2">2</xref>] - [<xref ref-type="bibr" rid="scirp.53602-ref6">6</xref>] ).</p><p>The deviations from normality become more severe when more frequent data are used to calculate stock returns. Various studies have shown that the normal distribution does not accurately describe observed stock return data. Over the past several decades, some stylized facts have emerged about the statistical behavior of speculative market returns such as aggregational Gaussianity, volatility clustering, etc see [<xref ref-type="bibr" rid="scirp.53602-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.53602-ref8">8</xref>] . On the same note, most of the literature for example [<xref ref-type="bibr" rid="scirp.53602-ref9">9</xref>] - [<xref ref-type="bibr" rid="scirp.53602-ref12">12</xref>] and references therein, assume that daily log returns, can be modeled by exponential L&#233;vy processes and geometric L&#233;vy process.</p><p>There are two important directions in the literature regarding these type of stochastic volatility models. Continuous-time stochastic volatility process represented in general by a bivariate diffusion process, and the discrete time autoregressive conditionally heteroscedastic (ARCH) model of [<xref ref-type="bibr" rid="scirp.53602-ref13">13</xref>] or its generalization (GARCH) as first defined by [<xref ref-type="bibr" rid="scirp.53602-ref14">14</xref>] . Option pricing in GARCH models has been typically done using the local risk neutral valuation relationship (LRNVR) pioneered by [<xref ref-type="bibr" rid="scirp.53602-ref15">15</xref>] . The crucial assumptions in his construction are the con- ditional, normal distribution of the asset returns under the underlying probability space and the invariance of the conditional volatility to the change of measure. The empirical performance of these normal option pricing models has been studied extensively, for example in [<xref ref-type="bibr" rid="scirp.53602-ref16">16</xref>] , [<xref ref-type="bibr" rid="scirp.53602-ref17">17</xref>] .</p><p>The main focus of this paper is to develop a ARCH type L&#233;vy model which attempts to capture some of the stylized features observed in demeaned log returns from any market data. More so we derive unconditional variance of daily log returns in risk neutral world of different ARCH type models, and an in-depth empirical study in liquid and illiquid market. All parameters are estimated from historical data, i.e. for S&amp;P500 index from January 3, 1990 to January 18, 2008 and NSE20 index from March 2, 1998 to July 11, 2007.</p><p>The article is organized as follows. Section 2 provides a brief overview of ARCH type models and L&#233;vy increments resulting to parameter estimation of observed salient features. In Section 3 which is our major con- tribution, unconditional variance of different ARCH type models is presented. Filtered Leptokurtic residuals of L&#233;vy increments are calibrated. Conclusions are drawn in Section 4. Appendix is in the last section.</p></sec><sec id="s2"><title>2. ARCH Type Models</title><p>ARCH-type models are in general, discrete models used to estimate volatility of financial time series data such stock returns, interest rates and foreign exchange rates. Let</p><disp-formula id="scirp.53602-formula69"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x5.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x6.png" xlink:type="simple"/></inline-formula> denotes the price of stock at time<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x7.png" xlink:type="simple"/></inline-formula>. Define the following equation</p><disp-formula id="scirp.53602-formula70"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x8.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.53602-formula71"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x9.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x10.png" xlink:type="simple"/></inline-formula> is the GARCH(p, q) volatility process. If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x11.png" xlink:type="simple"/></inline-formula> then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x12.png" xlink:type="simple"/></inline-formula> is ARCH(p). [<xref ref-type="bibr" rid="scirp.53602-ref18">18</xref>] and [<xref ref-type="bibr" rid="scirp.53602-ref19">19</xref>] provide a general specifications of volatility dynamic that nest most ARCH type models. In this connection volatility dynamics can be written as</p><disp-formula id="scirp.53602-formula72"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x13.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x14.png" xlink:type="simple"/></inline-formula> is the innovation function. Different GARCH models are mainly characterized by the following specifications of the innovation function<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x15.png" xlink:type="simple"/></inline-formula>.</p><disp-formula id="scirp.53602-formula73"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x16.png"  xlink:type="simple"/></disp-formula><p>The innovation function is used to model asymmetry and news impact to say the least. These GARCH models can be generalized to allow non-linearity of volatility dynamics by using Box-Cox transformation as follows</p><disp-formula id="scirp.53602-formula74"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x17.png"  xlink:type="simple"/></disp-formula><p>which implies modeling news and power, will nest most of the proposed GARCH models in Literature. Note that the leverage parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x18.png" xlink:type="simple"/></inline-formula> shifts the innovation function, the news parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x19.png" xlink:type="simple"/></inline-formula> tilts the innovation, and the power parameters <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x20.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x21.png" xlink:type="simple"/></inline-formula> flatten or steepen the innovation function. Such a model (4) is the Asym- metric Power Autoregressive Conditional Heteroscedastic model i.e. APARCH model defined in (5).</p><p>The APARCH(m, n) model of can be written as follows</p><disp-formula id="scirp.53602-formula75"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x22.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula76"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x23.png"  xlink:type="simple"/></disp-formula><p>subject to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x24.png" xlink:type="simple"/></inline-formula> for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x25.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x26.png" xlink:type="simple"/></inline-formula>, for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x27.png" xlink:type="simple"/></inline-formula>. and</p><disp-formula id="scirp.53602-formula77"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x28.png"  xlink:type="simple"/></disp-formula><p>The model introduces a Box-Cox power transformation on the conditional standard deviation process and on</p><p>the asymmetric innovations, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x29.png" xlink:type="simple"/></inline-formula>, adds flexibility of a varying exponent with an asymmetry co-</p><p>efficient to take the leverage effect into account. The properties of APARCH model have been studied, see [<xref ref-type="bibr" rid="scirp.53602-ref20">20</xref>] . The model nests seven other ARCH extensions as special cases.</p><p>・ ARCH model of [<xref ref-type="bibr" rid="scirp.53602-ref13">13</xref>] when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x30.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x31.png" xlink:type="simple"/></inline-formula></p><p>・ GARCH model of [<xref ref-type="bibr" rid="scirp.53602-ref14">14</xref>] when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x32.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x33.png" xlink:type="simple"/></inline-formula></p><p>・ GJR-GARCH Model of [<xref ref-type="bibr" rid="scirp.53602-ref21">21</xref>] when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x34.png" xlink:type="simple"/></inline-formula></p><p>・ TARCH Model of [<xref ref-type="bibr" rid="scirp.53602-ref22">22</xref>] when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x35.png" xlink:type="simple"/></inline-formula></p><p>Note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x36.png" xlink:type="simple"/></inline-formula> denote the conditional mean given the information set <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x37.png" xlink:type="simple"/></inline-formula> available at time t − 1. The innovation process for the conditional mean is then given by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x38.png" xlink:type="simple"/></inline-formula> with corresponding unconditional variance <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x39.png" xlink:type="simple"/></inline-formula> and zero unconditional mean. The conditional variance is defined as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x40.png" xlink:type="simple"/></inline-formula></p><sec id="s2_1"><title>2.1. Empirical Data</title><p>For simplicity, we focus on daily closing indices <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x41.png" xlink:type="simple"/></inline-formula> as reported in Nairobi Securities Exchange for NSE20 share index and S&amp;P500 index in New-York Stock Exchange. Daily log-returns <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x42.png" xlink:type="simple"/></inline-formula> of S&amp;P500 index are computed from January 3, 1990 to January 18, 2008 for a total of 4550 daily observations. While for NSE20, share indexes are computed from March 2, 1998 to July 11, 2007 for a total of 2317 daily observations.</p><p>All return series exhibit strong conditional heteroscedasticity. The Ljung and Box test rejects the hypothesis of homoscedasticity at all common levels both for returns in S&amp;P500 index and AR(3) residuals of linear re- gression in NSE20 share index. We estimate GARCH type models assuming conditional normality. With re- spect to the absolute value of parameter estimates, we find that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x43.png" xlink:type="simple"/></inline-formula> but different for both indices (NSE20<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x44.png" xlink:type="simple"/></inline-formula>, S&amp;P500<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x45.png" xlink:type="simple"/></inline-formula>), indicating the typical higher per- sistence of shocks in volatility in New York Stock exchange compared to Nairobi Securities Exchange. Model (5) is estimated using Pseudo Maximum Likelihood estimator based on the assumption of conditional normal in- novations. The parameter estimates of (8) are reported in <xref ref-type="table" rid="table1">Table 1</xref> and AR-ARCH residual calibrations of GH distribution (9) are presented in <xref ref-type="table" rid="table2">Table 2</xref>. Empirical and kernel densities of fitted distributions for both indices are compared in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><disp-formula id="scirp.53602-formula78"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x46.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula79"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x47.png"  xlink:type="simple"/></disp-formula></sec><sec id="s2_2"><title>2.2. L&#233;vy Increments</title><p>Suppose <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x48.png" xlink:type="simple"/></inline-formula> is the characteristic function of a distribution. If for every positive integer n, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x49.png" xlink:type="simple"/></inline-formula>is the</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> GARCH and GJR model estimates for the indices</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >NSE20</th><th align="center" valign="middle" ></th><th align="center" valign="middle" >S&amp;P500</th><th align="center" valign="middle" ></th></tr></thead><tr><td align="center" valign="middle" >Parameter</td><td align="center" valign="middle" >GARCH</td><td align="center" valign="middle" >GJR<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x50.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >GARCH</td><td align="center" valign="middle" >GJR<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x51.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x52.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.18915 (0.024496)</td><td align="center" valign="middle" >0.18136 (0.02424)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x53.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.16451 (0.023785)</td><td align="center" valign="middle" >0.16245 (0.02352)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x54.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.11388 (0.023413)</td><td align="center" valign="middle" >0.11516 (0.02308)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x55.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.03549 (0.006902)</td><td align="center" valign="middle" >0.03458 (0.00647)</td><td align="center" valign="middle" >0.006577 (0.001645)</td><td align="center" valign="middle" >0.01088 (0.00204)</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x56.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.15023 (0.017978)</td><td align="center" valign="middle" >0.18578 (0.02528)</td><td align="center" valign="middle" >0.056461 (0.0067528)</td><td align="center" valign="middle" >0.00322 (0.00512)</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x57.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.78763 (0.024753)</td><td align="center" valign="middle" >0.79045 (0.02373)</td><td align="center" valign="middle" >0.937566 (0.0074845)</td><td align="center" valign="middle" >0.93202 (0.0079)</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x58.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >−0.07332 (0.02592)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.10558 (0.0123)</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x59.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >9.3468 (0.2287)</td><td align="center" valign="middle" >8.8337 (0.2648)</td><td align="center" valign="middle" >16.5309 (0.08541)</td><td align="center" valign="middle" >15.2862 (0.1220)</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x60.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >7.1689 (0.5739)</td><td align="center" valign="middle" >8.46159 (0.38973)</td><td align="center" valign="middle" >6.8918 (0.54835)</td><td align="center" valign="middle" >5.9298 (0.6551)</td></tr><tr><td align="center" valign="middle" >lgl</td><td align="center" valign="middle" >−8363.5</td><td align="center" valign="middle" >−8367.7</td><td align="center" valign="middle" >−15090.9</td><td align="center" valign="middle" >−15090.9</td></tr><tr><td align="center" valign="middle" >n</td><td align="center" valign="middle" >2316</td><td align="center" valign="middle" >2316</td><td align="center" valign="middle" >4549</td><td align="center" valign="middle" >4549</td></tr></tbody></table></table-wrap><p>Notes: standard errors are in parenthesis. lgl is the log likelihood.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x61.png" xlink:type="simple"/></inline-formula>power of a characteristic function, we say that the distribution is infinitely divisible. One can define for every such infinitely divisible distribution a stochastic process <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x62.png" xlink:type="simple"/></inline-formula> called a L&#233;vy process, which starts at zero, has independent and stationary increments and such that the distribution of an increment over <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x63.png" xlink:type="simple"/></inline-formula> has <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x64.png" xlink:type="simple"/></inline-formula> is the characteristic function. For more detailed treatment of L&#233;vy process, see [<xref ref-type="bibr" rid="scirp.53602-ref23">23</xref>] .</p><p>Definition 2.1 The probability density function of the one-dimensional Generalized Hyperbolic distribution is given by the following:</p><disp-formula id="scirp.53602-formula80"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x65.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x66.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x67.png" xlink:type="simple"/></inline-formula> is the modified Bessel function of third kind, with the index <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x68.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.53602-formula81"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x69.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x70.png" xlink:type="simple"/></inline-formula>is the location parameter and can take any real value, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x71.png" xlink:type="simple"/></inline-formula>is a scale parameter; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x72.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x73.png" xlink:type="simple"/></inline-formula> determine the distribution shape and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x74.png" xlink:type="simple"/></inline-formula> defines the subclasses of GH and is related to the tail flatness.</p><p>The mean and variance of GH distribution are given respectively by the followings</p><disp-formula id="scirp.53602-formula82"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x75.png"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.53602-formula83"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x76.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x77.png" xlink:type="simple"/></inline-formula>. Note that, if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x77.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x78.png" xlink:type="simple"/></inline-formula>, then</p><disp-formula id="scirp.53602-formula84"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x79.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula85"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x80.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula86"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x81.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula87"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x82.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula88"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x83.png"  xlink:type="simple"/></disp-formula><p>For more information about GH distribution, see [<xref ref-type="bibr" rid="scirp.53602-ref24">24</xref>] .</p></sec></sec><sec id="s3"><title>3. Modeling the Underlying</title><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x84.png" xlink:type="simple"/></inline-formula> be a stochastic basis describing the uncertainty of the economy. We refer to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x85.png" xlink:type="simple"/></inline-formula> as the</p><p>physical probability measure and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x86.png" xlink:type="simple"/></inline-formula> represent the information flow driven by Brownian motion <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x87.png" xlink:type="simple"/></inline-formula></p><p>and L&#233;vy proces<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x88.png" xlink:type="simple"/></inline-formula>. Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x89.png" xlink:type="simple"/></inline-formula> be the price of a stock at time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x90.png" xlink:type="simple"/></inline-formula> adapted to the natural filtration<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x91.png" xlink:type="simple"/></inline-formula>.</p><p>Define daily log return as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x92.png" xlink:type="simple"/></inline-formula> It is well known from our empirical studies that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x93.png" xlink:type="simple"/></inline-formula> can he represented as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x94.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x95.png" xlink:type="simple"/></inline-formula> is a mean function and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x96.png" xlink:type="simple"/></inline-formula> are the two components of the error term. Moreover, define a <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x97.png" xlink:type="simple"/></inline-formula> order autoregressive process <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x98.png" xlink:type="simple"/></inline-formula> with APARCH(m,n) error as</p><disp-formula id="scirp.53602-formula89"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x99.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x100.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x100.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x101.png" xlink:type="simple"/></inline-formula> are identically and independently distributed random variables. A general time series model for log returns would be</p><disp-formula id="scirp.53602-formula90"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x102.png"  xlink:type="simple"/></disp-formula><sec id="s3_1"><title>3.1. Risk Neutralization</title><p>In this section, we construct risk neutral probability measure in the context of [<xref ref-type="bibr" rid="scirp.53602-ref15">15</xref>] and [<xref ref-type="bibr" rid="scirp.53602-ref19">19</xref>] . Duan [<xref ref-type="bibr" rid="scirp.53602-ref15">15</xref>] intro- duced the GARCH option pricing model by generalizing the traditional risk neutral valuation methodology to the case of conditional heteroscedasticity, the so called Local Risk Neutral Valuation Relationship (LRNVR).</p><p>Definition 3.1 A pricing measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x103.png" xlink:type="simple"/></inline-formula> is said to satisfy the locally risk-neutral valuation relationship (LRNVR) if measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x104.png" xlink:type="simple"/></inline-formula> is equivalent to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x105.png" xlink:type="simple"/></inline-formula>, and</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Calibration of AR-GARCH(1,1) residuals to a class of infinitely divisible distributions</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >NSE20</th><th align="center" valign="middle" >GH</th><th align="center" valign="middle" >HY</th><th align="center" valign="middle" >NIG</th><th align="center" valign="middle" >S&amp;P500</th><th align="center" valign="middle" >GH</th><th align="center" valign="middle" >HY</th><th align="center" valign="middle" >NIG</th></tr></thead><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x106.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >−1.79233</td><td align="center" valign="middle" >1.0000</td><td align="center" valign="middle" >−0.5000</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x107.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >2.38336</td><td align="center" valign="middle" >1.0000</td><td align="center" valign="middle" >−0.500</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x108.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.98225</td><td align="center" valign="middle" >1.15813</td><td align="center" valign="middle" >0.66862</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x109.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.14671</td><td align="center" valign="middle" >1.68640</td><td align="center" valign="middle" >1.33977</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x110.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >−0.05226</td><td align="center" valign="middle" >−0.06604</td><td align="center" valign="middle" >−0.05864</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x111.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >−0.14279</td><td align="center" valign="middle" >−0.14976</td><td align="center" valign="middle" >−0.15755</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x112.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >1.79373</td><td align="center" valign="middle" >0.45207</td><td align="center" valign="middle" >1.18530</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x113.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.04052</td><td align="center" valign="middle" >1.04004</td><td align="center" valign="middle" >1.59588</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x114.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.12296</td><td align="center" valign="middle" >0.13923</td><td align="center" valign="middle" >0.13014</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x115.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.14292</td><td align="center" valign="middle" >0.15130</td><td align="center" valign="middle" >0.16032</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Empirical and kernel densities of standardized GARCH filtered L&#233;vy increments of NSE20 index (left) S&amp;P500 index (right) calibrated vs. density of fitted infinitely divisible distributions and normal distributions</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-1490293x116.png"/></fig><disp-formula id="scirp.53602-formula91"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x117.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula92"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x118.png"  xlink:type="simple"/></disp-formula><p>almost surely with respect to measure<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x119.png" xlink:type="simple"/></inline-formula>.</p><p>For some commonly used assumptions concerning utility functions and distributions of change of con- sumption, [<xref ref-type="bibr" rid="scirp.53602-ref15">15</xref>] shows that a representative agent maximizes his expected utility using the LRNVR measure<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x120.png" xlink:type="simple"/></inline-formula>. Risk neutralization should leave the variance unchanged and should transform the conditional expectation so that the discounted expected price of the underlying asset becomes a martingale. It is worth noting that in the case of homoscedasticity process, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x120.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x121.png" xlink:type="simple"/></inline-formula>, the conditional variances become the same constant and the LRNVR reduces to conventional risk neutral valuation relationship.</p><p>Consider the general model of daily log returns under the data generating probability measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x122.png" xlink:type="simple"/></inline-formula> as</p><disp-formula id="scirp.53602-formula93"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x123.png"  xlink:type="simple"/></disp-formula><p>where the parameters <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x125.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x126.png" xlink:type="simple"/></inline-formula> and given<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x127.png" xlink:type="simple"/></inline-formula>. The sequence <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x128.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x129.png" xlink:type="simple"/></inline-formula> are conditionally independent, while <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x129.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x130.png" xlink:type="simple"/></inline-formula> is the past information set. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x129.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x131.png" xlink:type="simple"/></inline-formula>represents the conditional expectation of returns.</p><p>The pricing measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x132.png" xlink:type="simple"/></inline-formula> shifts the error term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x133.png" xlink:type="simple"/></inline-formula> by some measurable function<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x134.png" xlink:type="simple"/></inline-formula>, so that the conditional expectation of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x135.png" xlink:type="simple"/></inline-formula> becomes equal to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x135.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x136.png" xlink:type="simple"/></inline-formula>. In the case of AR(1)APARCH(1,1)-L&#233;vy filter, we follow the [<xref ref-type="bibr" rid="scirp.53602-ref25">25</xref>] argument. Therefore under the equivalent martingale measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x135.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x136.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x137.png" xlink:type="simple"/></inline-formula> the model (16) translates to</p><disp-formula id="scirp.53602-formula94"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x138.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula95"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x139.png"  xlink:type="simple"/></disp-formula><p>The LRNVR implies that under the risk neutral measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x140.png" xlink:type="simple"/></inline-formula> the return process evolves as</p><disp-formula id="scirp.53602-formula96"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x141.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula97"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x142.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula98"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x143.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula99"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x144.png"  xlink:type="simple"/></disp-formula><p>It follows quite easily that</p><disp-formula id="scirp.53602-formula100"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x145.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3_2"><title>3.2. Unconditional Variance</title><p>The following propositions provide the unconditional variance for the process <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x146.png" xlink:type="simple"/></inline-formula> under <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x146.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x147.png" xlink:type="simple"/></inline-formula></p><p>Proposition 3.1 Consider AR(3) APARCH(1,1) L&#233;vy filter, with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x148.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x148.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x149.png" xlink:type="simple"/></inline-formula> which implies AR(3)- GARCH(1,1) L&#233;vy model, the unconditional variance of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x148.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x149.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x150.png" xlink:type="simple"/></inline-formula> under the LRNVR equivalent measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x148.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x149.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x150.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x151.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.53602-formula101"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x152.png"  xlink:type="simple"/></disp-formula><p>Proof: See Appendix. W</p><p>Proposition 3.2 A special case of AR(1)GARCH(1,1)L&#233;vy filter the unconditional variance under the LRNVR equivalent measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x153.png" xlink:type="simple"/></inline-formula> is given by</p><disp-formula id="scirp.53602-formula102"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x154.png"  xlink:type="simple"/></disp-formula><p>Proof: See Appendix. W</p><p>Example 3.1 In case of Hyperbolic distribution we substitute mean and variance respectively into (25). Where the parameters used maximize the likelihood function of Hyperbolic distribution. i.e. Let</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x155.png" xlink:type="simple"/></inline-formula>then,</p><disp-formula id="scirp.53602-formula103"><label>(25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x156.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula104"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x157.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula105"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x158.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula106"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x159.png"  xlink:type="simple"/></disp-formula><p>Consider a discrete time economy, where interest rates and returns are paid after each time interval of equal spaced length. Suppose there is a price for risk, measured in terms of a risk premium that is added to the risk free interest rate r to build the expected next period return. As in Duan [<xref ref-type="bibr" rid="scirp.53602-ref15">15</xref>] , we adopt and extend the ARCH-M model of [<xref ref-type="bibr" rid="scirp.53602-ref26">26</xref>] with the risk premium being linear functional of the conditional standard deviation, hence the following model under<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x160.png" xlink:type="simple"/></inline-formula>,</p><disp-formula id="scirp.53602-formula107"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x161.png"  xlink:type="simple"/></disp-formula><p>The parameters <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x162.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x162.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x163.png" xlink:type="simple"/></inline-formula> are constant parameters satisfying stationarity and positivity conditions, while the constant parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x162.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x164.png" xlink:type="simple"/></inline-formula> may be interpreted as the unit price for risk. If we change the function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x162.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x164.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x165.png" xlink:type="simple"/></inline-formula> in (29) to model news impact, we get threshold GARCH model of [<xref ref-type="bibr" rid="scirp.53602-ref21">21</xref>] where</p><disp-formula id="scirp.53602-formula108"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x166.png"  xlink:type="simple"/></disp-formula><p>hence the resulting TGARCH L&#233;vy filter model</p><disp-formula id="scirp.53602-formula109"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x167.png"  xlink:type="simple"/></disp-formula><p>Proposition 3.3 The unconditional variance of the GARCH-M L&#233;vy filter model under the LRNVR equivalent martingale measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x168.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.53602-formula110"><label>(32)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x169.png"  xlink:type="simple"/></disp-formula><p>Proof: See Appendix. W</p><p>Proposition 3.4 The unconditional variance of the TGARCH-M L&#233;vy filter model under equivalent martingale measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x170.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.53602-formula111"><label>(33)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x171.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.53602-formula112"><label>(34)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-1490293x172.png"  xlink:type="simple"/></disp-formula><p>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x173.png" xlink:type="simple"/></inline-formula> denoting the cumulative standard normal distribution function.</p><p>Proof: See Appendix. W</p></sec></sec><sec id="s4"><title>4. Concluding Remarks</title><p>This article develops an log-ARCH-L&#233;vy type risk neutral model. The proposed method delivers predictive dis- tribution of the payoff function for a given econometric model. As a result, the probability distribution could be useful to market participants who wish to compare the model predictions to the potential prices in liquid and illiquid markets.</p><p>Any effective option pricing model is expected to be consistent with distributional and time series properties of the underlying asset. The proposed model accommodates most of the observed stylistic fact about financial time series data i.e. skewness and leptokurtic nature of demeaned GARCH filtered log returns and perhaps aggregational Gaussianity. In summary,</p><p>・ developed markets and emerging markets may not have the same underlying dynamics. It would be incorrect to assume that a universal model for the underlying process for all markets.</p><p>・ The presence of linear autoregressive dynamics AR(3)-GARCH(1,1) effects in NSE20 index affects the un- conditional variance in risk neutral world. S&amp;P500 index was found to follow GARCH(1,1) plus leptokurtic residual which was calibrated in one class of generalized hyperbolic distributions,say for example, Normal inverse Gaussian (NIG).</p><p>・ The presence of autoregressive dynamics, i.e. AR(3)-GARCH(1,1) model of NSE20 index as an example of illiquid market would have an impact in pricing options, if the index were to be used as an underlying process.</p><p>The log-ARCH-L&#233;vy model is very tractable compared to other jump-diffusion or stochastic volatility models. It attempts to addresses the drawbacks of local volatilities. Further refinements and extensions are left for future research.</p></sec><sec id="s5"><title>Acknowledgements</title><p>Comments from the Editor and the anonymous referee are acknowledged. Financial support from International Science Progam (Sweden)/EAUMP is greatly appreciated.</p></sec><sec id="s6"><title>Appendix</title><p>Proof of proposition 3.1</p><p>Given <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x174.png" xlink:type="simple"/></inline-formula> We note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x174.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x175.png" xlink:type="simple"/></inline-formula> and</p><disp-formula id="scirp.53602-formula113"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x176.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula114"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x177.png"  xlink:type="simple"/></disp-formula><p>after rearranging and simple algebra</p><disp-formula id="scirp.53602-formula115"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x178.png"  xlink:type="simple"/></disp-formula><p>Thus under stationarity, the unconditional expectations are independent of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x179.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.53602-formula116"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x180.png"  xlink:type="simple"/></disp-formula><p>Therefore, the unconditional variance of AR(3)GARCH(1,1)Levy filter model under LRNVR equivalent mar- tingale measure is</p><disp-formula id="scirp.53602-formula117"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x181.png"  xlink:type="simple"/></disp-formula><p>Proof of proposition 3.2</p><p>This is a special case of (3.1) with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x182.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x182.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x183.png" xlink:type="simple"/></inline-formula>.</p><p>Proof of proposition 3.3</p><p>It is a special case of proposition 3.4 when we take <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x184.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x184.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x185.png" xlink:type="simple"/></inline-formula></p><p>Proof of proposition 3.4</p><p>Under measure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x186.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.53602-formula118"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x187.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x188.png" xlink:type="simple"/></inline-formula> is the risk premium and</p><disp-formula id="scirp.53602-formula119"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x189.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula120"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x190.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula121"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x191.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula122"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x192.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53602-formula123"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x193.png"  xlink:type="simple"/></disp-formula><p>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x194.png" xlink:type="simple"/></inline-formula> denoting the cumulative standard normal distribution. Note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x194.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-1490293x195.png" xlink:type="simple"/></inline-formula> and</p><disp-formula id="scirp.53602-formula124"><graphic  xlink:href="http://html.scirp.org/file/2-1490293x196.png"  xlink:type="simple"/></disp-formula><p>Therefore, for positive support</p></sec></body><back><ref-list><title>References</title><ref id="scirp.53602-ref1"><label>1</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Mandelbrot</surname><given-names> B. </given-names></name>,<etal>et al</etal>. (<year>1963</year>)<article-title>The Variation of Certain Speculative Prices</article-title><source> International Statistical Review</source><volume> 36</volume>,<fpage> 394</fpage>-<lpage>419</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.53602-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Clark, P. (1973) A Surbodinated Stochastic Process Model with Finite Variance for Speculative Prices. Econometrica, 41, 135-155. http://dx.doi.org/10.2307/1913889</mixed-citation></ref><ref id="scirp.53602-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Madan, D. and Seneta, E. (1990) The Variance Gamma (V.G.) Model for Share Markets. Journal of Business, 63, 511-524. http://dx.doi.org/10.1086/296519</mixed-citation></ref><ref id="scirp.53602-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Eberlein, E. and Keller, U. (1995) Hyperbolic Distributions in Finance. Bernolli, 1, 281-299. 
http://dx.doi.org/10.2307/3318481</mixed-citation></ref><ref id="scirp.53602-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Berndorff-Nielson, O. (1998) Process of Normal Inverse Gaussian Type. Finance and Stochastics, 2, 41-68. 
http://dx.doi.org/10.1007/s007800050032</mixed-citation></ref><ref id="scirp.53602-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Mwaniki, I.J. (2010) On APARCH L&amp;eacute;vy Filter Option Pricing Formula for Developed and Emerging Markets. PhD Thesis, University of Nairobi, Nairobi.</mixed-citation></ref><ref id="scirp.53602-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Rydberg, T. (2000) Realistic Statistical Modeling of Financial Data. International Statistical Review, 68, 233-258. 
http://dx.doi.org/10.1111/j.1751-5823.2000.tb00329.x</mixed-citation></ref><ref id="scirp.53602-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Cont, R. (2001) Empirical Properties of Asset Returns: Stylized Facts and Statistical Issues. Quantitative Finance, 1, 223-236. http://dx.doi.org/10.1080/713665670</mixed-citation></ref><ref id="scirp.53602-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Carr, P. and Madan, D. (1998) Option Valuation Using the Fast Fourier Transform. Journal of Computational Finance, 2, 61-73.</mixed-citation></ref><ref id="scirp.53602-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Barndorff-Nielson, O. (1998) Process of Normal Inverse Gaussian Type. Finance Stochastics, 2, 41-68. 
http://dx.doi.org/10.1007/s007800050032</mixed-citation></ref><ref id="scirp.53602-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Chan, T. (1999) Pricing Contingent Claims on Stock Driven by L&amp;eacute;vy Processes. The Annals of Applied Probability, 9, 504-528. http://dx.doi.org/10.1214/aoap/1029962753</mixed-citation></ref><ref id="scirp.53602-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Carr, P., German, H., Madan, D. and Yor, M. (2002) The Fine Structure of Asset Returns: An Empirical Investigation. Journal of Business, 75, 305-332. http://dx.doi.org/10.1086/338705</mixed-citation></ref><ref id="scirp.53602-ref13"><label>13</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Engle</surname><given-names> R. </given-names></name>,<etal>et al</etal>. (<year>1982</year>)<article-title>Autoregressive Conditional Heteroscedasticity with Estimates of Variance of United Kingdom Inflation</article-title><source> Journal of Business and Economic Statistics</source><volume> 9</volume>,<fpage> 987</fpage>-<lpage>1008</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.53602-ref14"><label>14</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bollerslev</surname><given-names> T. </given-names></name>,<etal>et al</etal>. (<year>1986</year>)<article-title>Generalized Autoregressive Conditional Heteroskedasticity</article-title><source> Journal of Econometrics</source><volume> 31</volume>,<fpage> 307</fpage>-<lpage>327</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.53602-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Duan, J. (1995) The GARCH Option Pricing Model. Mathematical Finance, 5, 13-32.  
http://dx.doi.org/10.1111/j.1467-9965.1995.tb00099.x</mixed-citation></ref><ref id="scirp.53602-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">H&amp;auml;rdle, W. and Hafner, C. (2000) Discrete Time Option Pricing with Flexible Volatility Estimation. Finance and Stochastics, 4, 189-207. http://dx.doi.org/10.1007/s007800050011</mixed-citation></ref><ref id="scirp.53602-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Christoffersen, P. and Jacobs, K. (2004) Which GARCH Model for Option Valuation? Management Science, 50, 1204-1221. http://dx.doi.org/10.1287/mnsc.1040.0276</mixed-citation></ref><ref id="scirp.53602-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Ding, Z., Granger, W. and Engle, R. (1993) A Long Memory Property of Stock Markets Returns and a New Model. Journal of Empirical Finance, 1, 83-106. http://dx.doi.org/10.1016/0927-5398(93)90006-D</mixed-citation></ref><ref id="scirp.53602-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Hentschel, L. (1995) All in the Family Nesting Symmetric and Asymmetric GARCH Models. Journal of Financial Economics, 39, 71-104. http://dx.doi.org/10.1016/0304-405X(94)00821-H</mixed-citation></ref><ref id="scirp.53602-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Laurent, S. (2004) Analytical Derivatives of the APARCH Model. Computational Economics, 24, 51-57.  
http://dx.doi.org/10.1023/B:CSEM.0000038851.72226.76</mixed-citation></ref><ref id="scirp.53602-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Glosten, L., Jagannathan, R. and Runkle, D. (1993) The Relationship between Expected Value and the Volatility of the Nominal Excess Returns on Stocks. Journal of Finance, 48, 1779-1801.  
http://dx.doi.org/10.1111/j.1540-6261.1993.tb05128.x</mixed-citation></ref><ref id="scirp.53602-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Zakoian, J. (1994) Threshold Heteroskedastic Models. Journal of Economic Dynamics and Control, 18, 931-955.  
http://dx.doi.org/10.1016/0165-1889(94)90039-6</mixed-citation></ref><ref id="scirp.53602-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Sato, K. (1999) L&amp;eacute;vy Process and Infinitely Divisible Distributions. Cambridge University Press, Cambridge.</mixed-citation></ref><ref id="scirp.53602-ref24"><label>24</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Barndorff-Nielsen</surname><given-names> O. </given-names></name>,<etal>et al</etal>. (<year>1977</year>)<article-title>Exponentially Decreasing Distributions for Logarithms of Particle Size</article-title><source> Proceedings of the Royal Society London Series A</source><volume> 353</volume>,<fpage> 401</fpage>-<lpage>419</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.53602-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Hafner, C. and Herwartz, H. (2001) Option Pricing under Linear Autoregressive Dynamics, Heteroskedasticity, and Conditional Leptokurtosis. Journal of Empirical Finance, 8, 1-34. http://dx.doi.org/10.1016/S0927-5398(00)00024-4</mixed-citation></ref><ref id="scirp.53602-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Engle, R., Lilian, D. and Robins, R. (1987) Estimating Time Varying Premia in Term Structure: The ARCH-M Model. Econometrica, 55, 391-407. http://dx.doi.org/10.2307/1913242</mixed-citation></ref></ref-list></back></article>