<?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.2016.61016</article-id><article-id pub-id-type="publisher-id">JMF-63958</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>
 
 
  Statistical Arbitrage in S&amp;P500
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>tefanos</surname><given-names>Drakos</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>International Centre for Computational Engineering, Rhodes, Greece</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>stefanos.drakos@gmail.com</email></corresp></author-notes><pub-date pub-type="epub"><day>05</day><month>02</month><year>2016</year></pub-date><volume>06</volume><issue>01</issue><fpage>166</fpage><lpage>177</lpage><history><date date-type="received"><day>8</day>	<month>January</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>24</month>	<year>February</year>	</date><date date-type="accepted"><day>29</day>	<month>February</month>	<year>2016</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>
 
 
   A methodology to create statistical arbitrage in stock Index S&amp;P500 is presented. A synthetic asset based on the cointegration relationship of the stocks with Index was constructed. In order to capture the dynamic of the market time adaptive algorithms have been developed and discussed. The pair trading strategy was applied in different periods between S&amp;P500 and synthetic asset and the results were evaluated. Different metrics have shown that the Multvariate Kalman Algorithm creates statistical arbitrage in index with much lower Maximum Drawdown and higher profit. The algorithm is neutral as the beta is close to zero and the Sharp Ratio remains high in all cases. 
 
</p></abstract><kwd-group><kwd>Statistical Arbitrage</kwd><kwd> Mean Reverting</kwd><kwd> Pair Trading</kwd><kwd> Kalman Filter</kwd><kwd> Trading Algorithms</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Financial markets are based on the general trading rule: buy with low price and sell with high price. The aim is the development of strategies with low risk and succeeds this general rule. Pure Arbitrage is a category of strategies with zero risk. As an example we can refer the case of buying and selling a stock at the same time with a different value in two different exchanges. The profit results from the difference in prices, breaking the law of one price.</p><p>Another category is the statistical arbitrage which is not risk free at all. Strategies of this type are aimed at the expected gain which is greater than the risk. The profit results from the mispricing of the stocks. To achieve this, one needs to assess whether the price of a stock is overvalued or undervalued relative to the actual value which is really hard to determine. The fundamentals of the stock, the demand of each period and the general economic environment are some of the factors that make the fair value evaluation difficult. Part of this category is the Pair Trading and is presented in several references as in [<xref ref-type="bibr" rid="scirp.63958-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.63958-ref10">10</xref>] . The Pairs trading is a strategy which is based on the relative pricing of stocks without actually interested in the true value of them. The relative pricing is based on the idea that two assets with the same features can be priced about the same price as the basis of the law of one price. When at a given time the stock prices are different then one is overvalued and the other undervalued relative to the actual price. The classic trading Pairs strategy derives from the characteristics of this incorrect pricing (mispricing) between the two stocks.</p><p>In the current work instead of the two stocks the mispricing between the index S&amp;P500 and a sub-set of stocks belonging to this is considered. This forms a portfolio consisting of one unit of the Index in long (or short) and the corresponding number γ (Hedge ratio) of subset’s stocks in the opposite position. Three methods adopted to determine the number γ (Hedge ratio): a) the ordinary, b) the rolling ordinary least squares (OLS) regression and c) the Kalman filter process as it is presented below. According to the strategy, the spread of the index S&amp;P500 and the stocks prices’ combination are computed and when this deviates from its historical average value then the investor bets on the return to the historical with selling and buying respectively the stocks of the portfolio. In practice to construct the synthetic asset is necessary to investigate the appropriate stock exhibiting long relationship between them with the index. The technique used for this purpose is the Cointegration as this presented in [<xref ref-type="bibr" rid="scirp.63958-ref1">1</xref>] according to which two time series <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x6.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x7.png" xlink:type="simple"/></inline-formula> are cointegrated if <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x8.png" xlink:type="simple"/></inline-formula> for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x9.png" xlink:type="simple"/></inline-formula> and notation I(d) means integrated order d.</p><p>In [<xref ref-type="bibr" rid="scirp.63958-ref1">1</xref>] concluded that pairs who cointegrate in sample period behave better in the out-of-sample period than those not cointegrate in the sample period. Based on the previous, the spread of the linear combination of the cointegrated stocks and the S&amp;P500 index have to be stationary process. To analyze this relationship the augmented Dickey-Fuller test (ADF) was used as presented in [<xref ref-type="bibr" rid="scirp.63958-ref2">2</xref>] . According to the strategy the stocks are not restricted to cointegrate in the out-of-sample period but it is an indication that they will present mean-reverting behavior.</p><p>Additionally the log of prices used is also instead the prices as in [<xref ref-type="bibr" rid="scirp.63958-ref3">3</xref>] . The main reason is that log-returns are time additive. So, in order to calculate the return over n periods using real returns we need to calculate the product of n numbers:<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x10.png" xlink:type="simple"/></inline-formula>.</p><p>If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x11.png" xlink:type="simple"/></inline-formula> defined as:</p><disp-formula id="scirp.63958-formula438"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x12.png"  xlink:type="simple"/></disp-formula><p>And</p><disp-formula id="scirp.63958-formula439"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x13.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.63958-formula440"><graphic  xlink:href="http://html.scirp.org/file/16-1490401x14.png"  xlink:type="simple"/></disp-formula><p>Consequently the profit of spread over a period is equal to:</p><disp-formula id="scirp.63958-formula441"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x15.png"  xlink:type="simple"/></disp-formula><p>The classic application of the trading pair strategy has published in many works in the past and has been applied in practice. In [<xref ref-type="bibr" rid="scirp.63958-ref4">4</xref>] the author went a further step and presented an algorithm for the pair trading of a stock basket cointegrated with the S&amp;P500 using constant cointegrating coefficient with a single value for all stocks belonging to the basket. To capture the real behavior of the markets this work presents new Time Adaptive algorithms using rolling ordinary least squares (OLS) regression and Multivariate Kalman Filter process where the time dependent hedge ratio is computed separately for each of the stocks forming the synthetic asset creating thereby statistical arbitrage conditions in index S&amp;P500 and increasing the strategy performance.</p></sec><sec id="s2"><title>2. Period of Methodology</title><p>The proposed methodology and the trading algorithm designed based on that divided in two different spaces. The first refers to the in sample period which is used to make all the appropriate test and construct the synthetic asset and the other to the out of sample period (Trading Period) where the synthetic asset trading based on the specific rules (<xref ref-type="fig" rid="fig1">Figure 1</xref>).</p><sec id="s2_1"><title>2.1. In Sample Period</title><p>The data of the in sample period used for the synthetic asset construction. Working on a daily data domain a year of closing prices was chosen as the in sample period for determine the set of cointegrated stocks with S&amp;P500 and create the synthetic asset.</p>Synthetic Asset Construction.<p>The paper presents different algorithms in order to create statistical arbitrage in S&amp;P500. The S&amp;P500, based on themarket capitalizations of 500 large companies equity indices, and many consider it one of the best representations of the U.S. stock market. The initial choice of stocks for cointegration test with S&amp;P500 is made by the S&amp;P100 which is a sub-set of the S&amp;P500, and measures the performance of large cap companies in the United States. Constituents of the S&amp;P100 are selected for sector balance and represent about 57% of the market capitalization of the S&amp;P500 and almost 45% of the market capitalization of the U.S. equity markets. The stocks in the S&amp;P100 tend to be the largest and most established companies in the S&amp;P500 (Wikipedia).</p><p>Using the data of a selected In Sample Period for each stock <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x16.png" xlink:type="simple"/></inline-formula> the second step of Engle and Granger’ approach adopted. Using the logarithmic price of stocks and S&amp;P500 the OLS regression shown below is performed:</p><disp-formula id="scirp.63958-formula442"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x17.png"  xlink:type="simple"/></disp-formula><p>Τhe Augmented Dickey-Fuller unit root test applied for the stationarity of the OLS residuals. According to the successful stationary result a subset of the S&amp;P100 is created and the components of this set are the candidate for the synthetic asset construction. Still working in the sample period a new OLS regression was carried out:</p><disp-formula id="scirp.63958-formula443"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x18.png"  xlink:type="simple"/></disp-formula><p>Or</p><disp-formula id="scirp.63958-formula444"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x19.png"  xlink:type="simple"/></disp-formula><p>Using again the Augmented Dickey-Fuller the stationarity of the new OLS residuals is examined.</p><p>If the stationarity exist then this an evidence of mean reverting long term behavior of the spread</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x20.png" xlink:type="simple"/></inline-formula>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x21.png" xlink:type="simple"/></inline-formula> and S is the set of stock where individually and as logarithmic</p><p>sum cointegrated with S&amp;P500. The dimension of S is dim(S).</p></sec><sec id="s2_2"><title>2.2. Out of Sample Period-Trading Period</title><p>According to the previous period different algorithms for synthetic asset trading are developed. The first was</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Period of study</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/16-1490401x22.png"/></fig><p>designed assumed that the cointegration coefficient is constant during the trading period.</p><p>In that case and using the Equation (3) the profit of the strategy during the period t &#247; t + h arises from the following equation:</p><disp-formula id="scirp.63958-formula445"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x23.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s3"><title>3. Time Adaptive Coefficient γ</title><p>In reality the system of trading is dynamic and updated as new information get to in and the cointegration coefecient (or the hedge ratio) cannot stay constant during the trading period. For that reason time adaptive algorithms are developed to capture the real conditions of the markets.</p><sec id="s3_1"><title>3.1. Rolling Ordinary Least Squares (OLS) Regression</title><p>The first one considering a rolling ordinary least squares (OLS) regression. The frequency of regression calculations raised by an optimization procedure and the cointegration coefficient calculated at each step by the regres-</p><p>sion of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x24.png" xlink:type="simple"/></inline-formula> against the<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x25.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3_2"><title>3.2. Kalman Filter Process</title><p>The Kalman filter process can be described by three different steps: the prediction the observation and the correction. A new approach were developed using a Multivariate Kalman filter process. Based on that the hedge ratio calculated separately for each stock owned in the synthetic asset and the computed vector of the calculated parameters at each time step has dimensions (N + 1) &#215; 1 where N = dim(S) while the dimensions of the covariance matrix is (N + 1) &#215; (N + 1). The aim of these algorithms is to calculate at each time step the updated hedge ratio of the synthetic asset. Assuming that the hedge ratio and the premium follow a random walk we have:</p><disp-formula id="scirp.63958-formula446"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x26.png"  xlink:type="simple"/></disp-formula><p>where:</p><p>y<sub>t</sub>: is the current state of the of the parameters.</p><p>y<sub>t</sub><sub>-1</sub>: is the previous state of the of the parameters.</p><disp-formula id="scirp.63958-formula447"><graphic  xlink:href="http://html.scirp.org/file/16-1490401x27.png"  xlink:type="simple"/></disp-formula><p>where for the multivariate Kalman filter process</p><disp-formula id="scirp.63958-formula448"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x28.png"  xlink:type="simple"/></disp-formula><p>With:</p><disp-formula id="scirp.63958-formula449"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x29.png"  xlink:type="simple"/></disp-formula><p>h: is the cointegration coefficient from in sample period and μ<sub>0</sub> coming from the same period.</p><p>The vector of logarithmic price of stocks:</p><disp-formula id="scirp.63958-formula450"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x30.png"  xlink:type="simple"/></disp-formula><p>And <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x31.png" xlink:type="simple"/></inline-formula></p><p>The process following the steps as below:</p><p>Prediction state where the next system state is predicted based on the knowledge of the previous state</p><disp-formula id="scirp.63958-formula451"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x32.png"  xlink:type="simple"/></disp-formula><p>The covariance of prediction state is given by:</p><disp-formula id="scirp.63958-formula452"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x33.png"  xlink:type="simple"/></disp-formula><p>The next step concerns the measurement prediction. Given the price of the synthetic asset and the predicted hedge ratio the measurement prediction are given as:</p><disp-formula id="scirp.63958-formula453"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x34.png"  xlink:type="simple"/></disp-formula><p>The residual of measurement and real value at each step calculated as:</p><disp-formula id="scirp.63958-formula454"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x35.png"  xlink:type="simple"/></disp-formula><p>The variance of measurement error is equal to:</p><disp-formula id="scirp.63958-formula455"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x36.png"  xlink:type="simple"/></disp-formula><p>The Kalman Gain is the filter, which tells how much the predictions should be corrected on time step is given as:</p><disp-formula id="scirp.63958-formula456"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x37.png"  xlink:type="simple"/></disp-formula><p>The last step of process is the update step where:</p><p>The updated state is estimated as following:</p><disp-formula id="scirp.63958-formula457"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x38.png"  xlink:type="simple"/></disp-formula><p>And the Updated state covariance is equal to</p><disp-formula id="scirp.63958-formula458"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x39.png"  xlink:type="simple"/></disp-formula><p>All the process repeated at every time step of out of sample period. The estimation of V<sub>w</sub> and V<sub>e</sub> has been discussed in [<xref ref-type="bibr" rid="scirp.63958-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.63958-ref12">12</xref>] .</p></sec><sec id="s3_3"><title>3.3. Profit of Strategies</title><p>In the case of Time Adaptive coefficient γ of the linear regression case the profit of the pair trading strategy raised by the following equations:</p><disp-formula id="scirp.63958-formula459"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x40.png"  xlink:type="simple"/></disp-formula><p>In the case of Multivariate Kalman Filter where the hedge ratio is different for each stock and for each time step in the synthetic asset it can been shown that:</p><disp-formula id="scirp.63958-formula460"><graphic  xlink:href="http://html.scirp.org/file/16-1490401x41.png"  xlink:type="simple"/></disp-formula><p>Similarly</p><disp-formula id="scirp.63958-formula461"><graphic  xlink:href="http://html.scirp.org/file/16-1490401x42.png"  xlink:type="simple"/></disp-formula><p>And finally the profit during a period <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x43.png" xlink:type="simple"/></inline-formula> is equal to:</p><disp-formula id="scirp.63958-formula462"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x44.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s4"><title>4. The Pair Trading Strategy</title><p>The algorithm of the pair trading strategy is based on the distance of the spread from its historical mean value and its mean-reverting behavior. To measure this distance a normalized variable called z-score introduced as:</p><disp-formula id="scirp.63958-formula463"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/16-1490401x45.png"  xlink:type="simple"/></disp-formula><p>where:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x46.png" xlink:type="simple"/></inline-formula>: is the mean value of the spread over a lookback period.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x47.png" xlink:type="simple"/></inline-formula>: is the standard deviation of the spread over the same period.</p><p>The trading it takes place when this variable exceeds some limits based on the spread mean reverting behavior. Thus:</p><p>Open-long position if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x48.png" xlink:type="simple"/></inline-formula>.</p><p>Open-short position if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x49.png" xlink:type="simple"/></inline-formula>.</p><p>Exit-long position if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x50.png" xlink:type="simple"/></inline-formula>.</p><p>Exit-short position if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x51.png" xlink:type="simple"/></inline-formula>.</p><p>When a long (short) position is opened we buy (sell) one unit of S&amp;P500 and sell (buy) the following amount</p><p>of stocks from synthetic asset as: <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x52.png" xlink:type="simple"/></inline-formula>in the case of constant hedge ratio, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x53.png" xlink:type="simple"/></inline-formula>when the rolling OLS regression is applied, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/16-1490401x54.png" xlink:type="simple"/></inline-formula> if the Kalman Filter algorithm applied. The performance</p><p>of strategy evaluated using the following metrics: 1) Cumulative return, 2) Annualized Return, 3) Sharpe Ratio, 4) Maximum drawdown, 5) beta.</p></sec><sec id="s5"><title>5. Back Testing-Performance Evaluation</title><p>Five different time periods was studied. In each case one year data collection was used in order to make the entire test and construct the synthetic asset. After that the algorithm starts trading with ending day for all cases the 30/12/2015. All the sample periods started at the first day of the year and ending at the last of the same year. The trading started on the first day the next year. In <xref ref-type="table" rid="table1">Table 1</xref> the back testing periods are presented. In <xref ref-type="fig" rid="fig2">Figure 2</xref> and <xref ref-type="fig" rid="fig3">Figure 3</xref> the results of cumulative return of Multivariate Kalman Filter algorithm against S&amp;P500 and its maximum drawdown are shown. In <xref ref-type="fig" rid="fig4">Figure 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref> the cumulative return of rolling OLS regression algorithm against S&amp;P500 and the cumulative return of rolling OLS regression algorithm against constant hedge ratio are illustrated. In t <xref ref-type="table" rid="table2">Table 2</xref> the sets of synthetic assets are presented while in <xref ref-type="table" rid="table3">Table 3</xref> the name of symbols is given. In Tables 4-8 all the metrics of the algorithms are displayed.</p><p>It can be shown that the dimension of the set and the constituents are different from period to period. This is an evidence of the non constant cointegration behavior of the stocks with the index but as we can see from the graphs the synthetic asset of each period continues to trade with profit and good metrics results.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Back testing periods</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >2006</th><th align="center" valign="middle" >2007</th><th align="center" valign="middle" >2008</th><th align="center" valign="middle" >2009</th><th align="center" valign="middle" >2010</th><th align="center" valign="middle" >2011</th><th align="center" valign="middle" >2012</th><th align="center" valign="middle" >2013</th><th align="center" valign="middle" >2014</th><th align="center" valign="middle" >2015</th></tr></thead><tr><td align="center" valign="middle" >In Sample</td><td align="center" valign="middle"  colspan="9"  >Trading</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2007</td><td align="center" valign="middle" >2008</td><td align="center" valign="middle" >2009</td><td align="center" valign="middle" >2010</td><td align="center" valign="middle" >2011</td><td align="center" valign="middle" >2012</td><td align="center" valign="middle" >2013</td><td align="center" valign="middle" >2014</td><td align="center" valign="middle" >2015</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >In Sample</td><td align="center" valign="middle"  colspan="8"  >Trading</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >2008</td><td align="center" valign="middle" >2009</td><td align="center" valign="middle" >2010</td><td align="center" valign="middle" >2011</td><td align="center" valign="middle" >2012</td><td align="center" valign="middle" >2013</td><td align="center" valign="middle" >2014</td><td align="center" valign="middle" >2015</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >In Sample</td><td align="center" valign="middle"  colspan="7"  >Trading</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >2009</td><td align="center" valign="middle" >2010</td><td align="center" valign="middle" >2011</td><td align="center" valign="middle" >2012</td><td align="center" valign="middle" >2013</td><td align="center" valign="middle" >2014</td><td align="center" valign="middle" >2015</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >In Sample</td><td align="center" valign="middle"  colspan="6"  >Trading</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >2010</td><td align="center" valign="middle" >2011</td><td align="center" valign="middle" >2012</td><td align="center" valign="middle" >2013</td><td align="center" valign="middle" >2014</td><td align="center" valign="middle" >2015</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >In Sample</td><td align="center" valign="middle"  colspan="5"  >Trading</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Sets of stock in synthetic asset on each period</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Trading Period</th><th align="center" valign="middle" >Stock of Synthetic Asset (Set S)</th><th align="center" valign="middle" >Dim (S)</th></tr></thead><tr><td align="center" valign="middle" >2007-2015</td><td align="center" valign="middle" >{BRK.B, DVN, IBM, SPG}</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >2008-2015</td><td align="center" valign="middle" >{CAN, BAX, BMY, DELL, EMC, GE, HAL, MDLZ, NSC, T, TXN, VZ}</td><td align="center" valign="middle" >12</td></tr><tr><td align="center" valign="middle" >2009-2015</td><td align="center" valign="middle" >{DVN, FCX, IBM, MON, SPG}</td><td align="center" valign="middle" >5</td></tr><tr><td align="center" valign="middle" >20010-2015</td><td align="center" valign="middle" >{APA, AXP, C, CAT, COF, FCX, IBM, MMM}</td><td align="center" valign="middle" >8</td></tr><tr><td align="center" valign="middle" >2011-2015</td><td align="center" valign="middle" >{COST, FDX, GS, QCOM}</td><td align="center" valign="middle" >4</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Union set of companies raised from cointegration test from all the period of study</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >APA</th><th align="center" valign="middle" >APACHE CORP</th><th align="center" valign="middle" >FDX</th><th align="center" valign="middle" >FEDEX CORPORATION</th></tr></thead><tr><td align="center" valign="middle" >AXP</td><td align="center" valign="middle" >AMERICAN EXPRESS COMPANY</td><td align="center" valign="middle" >GE</td><td align="center" valign="middle" >GENERAL ELECTRIC CO</td></tr><tr><td align="center" valign="middle" >BAX</td><td align="center" valign="middle" >BAXTER INTERNATIONAL INC</td><td align="center" valign="middle" >GS</td><td align="center" valign="middle" >GOLDMAN SACHS GROUP INC</td></tr><tr><td align="center" valign="middle" >BMY</td><td align="center" valign="middle" >RISTOL MYERS SQUIBB COMPANY</td><td align="center" valign="middle" >HAL</td><td align="center" valign="middle" >HALLIBURTON CO (HOLDING CO)</td></tr><tr><td align="center" valign="middle" >BRK.B</td><td align="center" valign="middle" >BERKSHIRE HATHWY INC (HLDG CO) B</td><td align="center" valign="middle" >IBM</td><td align="center" valign="middle" >INTL BUSINESS MACHINES CORP</td></tr><tr><td align="center" valign="middle" >C</td><td align="center" valign="middle" >CITIGROUP</td><td align="center" valign="middle" >MDLZ</td><td align="center" valign="middle" >MONDELEZ INTERNATIONAL INC</td></tr><tr><td align="center" valign="middle" >CAT</td><td align="center" valign="middle" >CATERPILLAR INC</td><td align="center" valign="middle" >MMM</td><td align="center" valign="middle" >3M COMPANY</td></tr><tr><td align="center" valign="middle" >COF</td><td align="center" valign="middle" >CAPITAL ONE FINANCIAL CORP</td><td align="center" valign="middle" >MON</td><td align="center" valign="middle" >MONSANTO COMPANY</td></tr><tr><td align="center" valign="middle" >COST</td><td align="center" valign="middle" >COSTCO WHOLESALE CORP</td><td align="center" valign="middle" >NSC</td><td align="center" valign="middle" >NORFOLK SOUTHERN CORP</td></tr><tr><td align="center" valign="middle" >DELL</td><td align="center" valign="middle" >DELL INC</td><td align="center" valign="middle" >QCOM</td><td align="center" valign="middle" >QUALCOMM INC</td></tr><tr><td align="center" valign="middle" >DVN</td><td align="center" valign="middle" >DEVON ENERGY CORP (NEW)</td><td align="center" valign="middle" >SPG</td><td align="center" valign="middle" >SIMON PROPERTIES GROUP INC</td></tr><tr><td align="center" valign="middle" >EMC</td><td align="center" valign="middle" >EMC CORPORATION</td><td align="center" valign="middle" >T</td><td align="center" valign="middle" >AT&amp;T INC. COM</td></tr><tr><td align="center" valign="middle" >FCX</td><td align="center" valign="middle" >FREEPORT-MCMORAN COPPER&amp;GOLD B</td><td align="center" valign="middle" >TXN</td><td align="center" valign="middle" >EXAS INSTRUMENTS INC</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >VZ</td><td align="center" valign="middle" >VERIZON COMMUNICATIONS</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Statistical metrics performance for the trading period 01/01/2007-30/12/2015</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >In Sample Period</th><th align="center" valign="middle"  colspan="3"  >01/01/2006-31/12/2006</th><th align="center" valign="middle"  rowspan="4"  >S&amp;P500</th></tr></thead><tr><td align="center" valign="middle" >Out of Sample Period</td><td align="center" valign="middle"  colspan="3"  >01/01/2007-30/12/2015</td></tr><tr><td align="center" valign="middle"  rowspan="2"  ></td><td align="center" valign="middle"  colspan="3"  >Algorithm</td></tr><tr><td align="center" valign="middle" >MultiVariate Kalman Filter</td><td align="center" valign="middle" >Rolling (OLS) regression</td><td align="center" valign="middle" >Constant Hedge Ratio</td></tr><tr><td align="center" valign="middle" >Cumulative Return %</td><td align="center" valign="middle" >387.069</td><td align="center" valign="middle" >24.637</td><td align="center" valign="middle" >12.678</td><td align="center" valign="middle" >45.592</td></tr><tr><td align="center" valign="middle" >Annual Return %</td><td align="center" valign="middle" >19.261</td><td align="center" valign="middle" >2.481</td><td align="center" valign="middle" >1.337</td><td align="center" valign="middle" >4.268</td></tr><tr><td align="center" valign="middle" >Sharpe Ratio</td><td align="center" valign="middle" >3.298</td><td align="center" valign="middle" >0.892</td><td align="center" valign="middle" >0.483</td><td align="center" valign="middle" >0.303</td></tr><tr><td align="center" valign="middle" >Beta</td><td align="center" valign="middle" >0.019</td><td align="center" valign="middle" >0.017</td><td align="center" valign="middle" >0.027</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Maximum Drawdown %</td><td align="center" valign="middle" >12.642</td><td align="center" valign="middle" >2.925</td><td align="center" valign="middle" >5.666</td><td align="center" valign="middle" >62.453</td></tr><tr><td align="center" valign="middle" >Duration of Maximum Drawdown</td><td align="center" valign="middle" >73.0</td><td align="center" valign="middle" >147.0</td><td align="center" valign="middle" >390.0</td><td align="center" valign="middle" >1365.0</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Statistical metrics performance for the trading period 01/01/2008-30/12/201</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >In Sample Period</th><th align="center" valign="middle"  colspan="3"  >01/01/2007-31/12/2007</th><th align="center" valign="middle" ></th></tr></thead><tr><td align="center" valign="middle" >Out of Sample Period</td><td align="center" valign="middle"  colspan="3"  >01/01/2008-30/12/2015</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle"  rowspan="2"  ></td><td align="center" valign="middle"  colspan="3"  >Algorithm</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Multi Variate Kalman Filter</td><td align="center" valign="middle" >Rolling (OLS) regression</td><td align="center" valign="middle" >Constant Hedge Ratio</td><td align="center" valign="middle" >S&amp;P500</td></tr><tr><td align="center" valign="middle" >Cumulative Return %</td><td align="center" valign="middle" >171.671</td><td align="center" valign="middle" >22.347</td><td align="center" valign="middle" >6.855</td><td align="center" valign="middle" >42.470</td></tr><tr><td align="center" valign="middle" >Annual Return %</td><td align="center" valign="middle" >13.321</td><td align="center" valign="middle" >2.556</td><td align="center" valign="middle" >0.833</td><td align="center" valign="middle" >4.528</td></tr><tr><td align="center" valign="middle" >Sharpe Ratio</td><td align="center" valign="middle" >1.135</td><td align="center" valign="middle" >0.738</td><td align="center" valign="middle" >0.228</td><td align="center" valign="middle" >0.312</td></tr><tr><td align="center" valign="middle" >Beta</td><td align="center" valign="middle" >−0.010</td><td align="center" valign="middle" >−0.005</td><td align="center" valign="middle" >−0.013</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Maximum Drawdown %</td><td align="center" valign="middle" >28.636</td><td align="center" valign="middle" >6.730</td><td align="center" valign="middle" >5.301</td><td align="center" valign="middle" >52.909</td></tr><tr><td align="center" valign="middle" >Duration of Maximum Drawdown</td><td align="center" valign="middle" >175.0</td><td align="center" valign="middle" >404.0</td><td align="center" valign="middle" >628.0</td><td align="center" valign="middle" >154.0</td></tr></tbody></table></table-wrap>
<table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Statistical metrics performance for the trading period 01/01/2010-30/12/201</title></caption></table-wrap></sec></body>
<back><ref-list><title>References</title><ref id="scirp.63958-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Engle, R. and Granger, C. (1987) Co-Integration and Error Correction: Representation, Estimation, and Testing. Econometrica, 55, 251-276. &lt;/br&gt;http://dx.doi.org/10.2307/1913236</mixed-citation></ref><ref id="scirp.63958-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Vidyamurthy, G. (2004) Pairs Trading, Quantitative Methods and Analysis. John Wiley &amp; Sons, Hoboken.</mixed-citation></ref><ref id="scirp.63958-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Infantino, L. and Itzhaki, S. (2010) Developing High-Frequency Equities Trading Models. Master of Business Administration, Massachusetts Institute of Technology, Cambridge.</mixed-citation></ref><ref id="scirp.63958-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Ernie, C. (2013) Algorithmic Trading: Winning Strategies and Their Rationale. John Wiley &amp; Sons, Hoboken.</mixed-citation></ref><ref id="scirp.63958-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Caldeira, J. and Moura, G.V. (2013) Selection of a Portfolio of Pairs Based on Cointegration: A Statistical Arbitrage Strategy. &lt;/br&gt;http://dx.doi.org/10.2139/ssrn.2196391</mixed-citation></ref><ref id="scirp.63958-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Dunis, C.L. and Ho, R. (2005) Cointegration Portfolios of European Equities for Index Tracking and Market Neutral Strategies. Journal of Asset Management, 6, 33-52.&lt;/br&gt;http://dx.doi.org/10.1057/palgrave.jam.2240164</mixed-citation></ref><ref id="scirp.63958-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Dunis, C.L, Giorgioni, G., Laws, J. and Rudy, J. (2010) Statistical Arbitrage and High-Frequency Data with an Application to Eurostoxx 50 Equities. CIBEF Working Papers, CIBEF.</mixed-citation></ref><ref id="scirp.63958-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Elliott, R., van der Hoek, J. and Malcolm, W. (2005) Pairs Trading. Quantitative Finance, 5, 271-276. &lt;/br&gt;http://dx.doi.org/10.1080/14697680500149370</mixed-citation></ref><ref id="scirp.63958-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Gatev, E., Goetzmann, W.N. and Rouwenhorst, K.G. (2006) Pairs Trading: Performance of a Relative-Value Arbitrage Rule. Review of Financial Studies, 19, 797-827. &lt;/br&gt;http://dx.doi.org/10.1093/rfs/hhj020</mixed-citation></ref><ref id="scirp.63958-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Khandani, A. and Lo, A.W. (2007) What Happened to the Quants in August 2007. &lt;/br&gt;http://web.mit.edu/Alo/www/Papers/august07.pdf</mixed-citation></ref><ref id="scirp.63958-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Rajamani, M. (2007) Data-Based Techniques to Improve State Estimation in Model Predictive Control. PhD Thesis, University of Wisconsin-Madison, Madison.</mixed-citation></ref><ref id="scirp.63958-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Rajamani, M.R. and Rawlings, J.B. (2009) Estimation of the Disturbance Structure from Data Using Semidefinite Programming and Optimal Weighting. Automatica, 45, 142-148. &lt;/br&gt;http://dx.doi.org/10.1016/j.automatica.2008.05.032</mixed-citation></ref></ref-list></back></article>