 be simply estimated by computing the leading sparse eigenvector of $\lambda I-A_0$ where $\lambda$ is bigger than the leading eigenvalue of $A_0$. This sparse principal component of the covariance matrix $A_0$ should not be confused with our utilization of sparse PCA in Section \[subsec:asemidefinite\] as a way to recover a vector solution from the solution of a positive semidefinite problem. Note also that techniques based on principal components do not take explicitly variance levels into account when estimating the weights of a co-integrated relationship.

@jurek Trading Strategy
-----------------------

While option implied volatility is not directly tradable, it can be synthesized using baskets of call options, and we assimilate it to a tradable asset with (significant) transaction costs in what follows. For baskets of volatilities isolated by the techniques listed above, we apply the [@jurek] strategy for log utilities to the basket process recording out of sample performance. @jurek propose to trade a stationary autoregressive process $(x_t)_{t}$ of order $1$ and mean $\mu$ governed by the equation $x_{t+1} = \rho x_t +\sigma \varepsilon_t$, where $\abs{\rho}<1$, by taking a position $N_t$ in the asset $x_t$ which is proportional to $$\label{eq:jurek}
    N_t = \frac{\rho (\mu-x_t)}{\sigma^2}W_t$$ In effect, the strategy advocates taking a long (resp. short) position in the asset whenever it is below (resp. above) its long-term mean, and adjust the position size to account for the volatility of $x_t$ and its mean reversion speed $\rho$. Given basket weights $y$, we apply standard AR estimation procedures on the in-sample portion of $y^T\bx$ to recover estimates for $\hat{\rho}$ and $\hat{\sigma}$ and plug them directly in Equation . This approach is illustrated for two baskets in Figure \[fig:syn\].

-2.5cm![**Three sample trading experiments, using the PCA, sparse PCA and the Crossing Statistics estimators**. \[a\] Pool of 9 volatility time-series selected using our fast PCA selection procedure. \[b\] Basket weights estimated with in-sample data using either the eigenvector of the covariance matrix with smallest
be simply estimated by computing the leading sparse eigenvector of $ \lambda I - A_0 $ where $ \lambda$ is bad than the go eigenvalue of $ A_0$. This sparse principal component of the covariance matrix $ A_0 $ should not be confused with our use of sparse PCA in Section   \[subsec: asemidefinite\ ] as a way to recover a vector solution from the solution of a positive semidefinite trouble. notice also that techniques based on principal components do not aim explicitly variance horizontal surface into score when estimating the weight of a co - integrated relationship. 

 @jurek Trading Strategy 
 ----------------------- 

 While choice implied volatility is not directly tradable, it can be synthesize using baskets of cry options, and we absorb it to a tradable asset with (significant) transaction monetary value in what follows. For baskets of volatilities isolate by the techniques listed above, we apply the   [ @jurek ] strategy for log utility to the basket summons recording out of sample performance.   @jurek propose to trade a stationary autoregressive process $ (x_t)_{t}$ of order $ 1 $ and mean $ \mu$ regulate by the equation $ x_{t+1 } = \rho x_t + \sigma \varepsilon_t$, where $ \abs{\rho}<1 $, by taking a stead $ N_t$ in the asset $ x_t$ which is proportional to $ $ \label{eq: jurek } 
     N_t = \frac{\rho (\mu - x_t)}{\sigma^2}W_t$$ In effect, the strategy advocates taking a retentive (resp. short) position in the asset whenever it is below (resp. above) its long - term mean, and align the position size to account for the volatility of $ x_t$ and its mean reversion speed $ \rho$. Given basket weights $ y$, we apply standard AR estimation procedures on the in - sample helping of $ y^T\bx$ to recover estimate for $ \hat{\rho}$ and $ \hat{\sigma}$ and plug them directly in Equation  . This approach is illustrated for two basket in Figure   \[fig: syn\ ]. 

 -2.5cm![**Three sample trading experiments, using the PCA, sparse PCA and the Crossing Statistics calculator * *. \[a\ ] Pool of 9 excitability time - series selected using our fast PCA choice procedure. \[b\ ] Basket weights estimate with in - sample datum using either the eigenvector of the covariance matrix with smallest
 be simply estimated by comkuting the leadiny sparsx eigenbector ow $\lambda I-A_0$ where $\lambda$ is uiggwr thqn the leading eigenvauue of $A_0$. Nhis sparwe pcincipal componeif of thc covzviancz natrix $A_0$ shoulc not be cmnfused with ogr ucilization of sparse PCA in Section \[stbsec:asrmldefinite\] as a way eo rsboyer a vector solution from the sklution of a positivr semidefinite problem. Notf aldo that techniques based on ptjncypal componengs do not uane explicijly variance levels into account when estimatint rhe feights of e co-innegrated relabponship.

@burek Ttading Strategn
-----------------------

Whilx oprion implied volatilivy is not directly ttadable, it ccn be synthesized usibg baskgts ox cauo ootiknx, znd we asaimilate if to a tradqble asset with (sigmispvant) transacfion cjses in what follows. For baskets of volatplitjes isolated by the tecyniques listed above, ae apply ehe [@jurek] strategy for log utilities to the basket procxsr rtcirdine okt of sample performance. @jurek propose to tradq a snationary autoregvessive process $(x_t)_{y}$ lf jrder $1$ and mexn $\mu$ ykvsrned by the equatlon $x_{t+1} = \rho z_t +\sigma \darelsilon_t$, where $\abs{\rho}<1$, by takung a positijb $N_t$ in the asset $r_t$ which is 'roporjional to $$\label{eq:jurek}
    N_t = \frad{\rho (\mu-x_t)}{\sihma^2}W_t$$ In swfect, the strateey sdeocates uxking a long (resp. short) powitipn in tfe axset wrenever it is bcnow (resp. above) its long-jerm maan, and aduust the position size to accouiv for the volsthliny of $x_t$ cnd itx mean reversyon speed $\rho$. Yiven bafket deights $y$, se applb standard AW estimation [tocedures on vhe in-sam[le portuon of $h^G\bx$ to recover estimates for $\hat{\ryo}$ and $\hat{\sigma}$ anv plog them directly nu Tqyation . This aplroxch id mllusesated for twm barkegx in Wigure \[fig:syk\].

-2.5cm![**Ghrer sample trading exparimsnts, using the PCA, siarse PCA and the Crossing Stayistics estimators**. \[a\] Poml uf 9 vplajility time-series selected usinf our fash PGA selection [roccdurg. \[b\] Basket ceights estimated with in-sample data usiig either the eigenvecjor of the covariance katrix with smaljest
be simply estimated by computing the leading of I-A_0$ where is bigger than This principal component of covariance matrix $A_0$ not be confused with our utilization sparse PCA in Section \[subsec:asemidefinite\] as a way to recover a vector solution the solution of a positive semidefinite problem. Note also that techniques based on components not explicitly levels into account when estimating the weights of a co-integrated relationship. @jurek Trading Strategy ----------------------- While implied volatility is not directly tradable, it can synthesized using baskets of options, and we assimilate it a asset with transaction in follows. For baskets volatilities isolated by the techniques listed above, we apply the [@jurek] strategy for log utilities to the process recording sample performance. propose trade stationary autoregressive process order $1$ and mean $\mu$ governed $x_{t+1} = \rho x_t +\sigma \varepsilon_t$, where $\abs{\rho}<1$, taking a $N_t$ in the asset $x_t$ which proportional to $$\label{eq:jurek} N_t = \frac{\rho (\mu-x_t)}{\sigma^2}W_t$$ In the strategy advocates taking a long (resp. short) position in the asset whenever it is above) its long-term mean, adjust the position to for volatility $x_t$ and mean reversion speed $\rho$. Given basket weights $y$, we apply standard estimation procedures on the in-sample portion of $y^T\bx$ to recover $\hat{\rho}$ $\hat{\sigma}$ and plug directly in Equation . approach illustrated for two baskets \[fig:syn\]. sample the sparse and the Crossing Statistics \[a\] Pool of 9 volatility selected using our fast weights estimated with in-sample data using either the of the covariance matrix with smallest
 be simply estimated by computIng the leadIng spArsE eiGeNvecTor oF $\lambda I-A_0$ where $\LAmbdA$ is bigger than the leadinG eigeNvALue oF $a_0$. THis spArse priNCiPAL coMpOnEnt Of THe CovarIanCe matriX $A_0$ should noT be CoNfused with ouR UtIlization oF spArse PCA in SecTioN \[subseC:aSemIDefinIte\] As a waY to recOVer a veCtor solutIoN From thE SolutioN OF a PosiTive semidefinite pRObLEm. Note also that TechniQuES bASEd oN prIncipal comPoNents DO not takE ExPLICitLY variance leveLs into accouNT whEn estiMaTinG The weiGhts oF a CO-inTegrated relAtioNship.

@jureK TradiNG StrateGY
-----------------------

While oPtion iMplIed VolaTIlItY is NoT DirECtLy tRAdaBle, it can Be SyNthesIzed USING basKetS of cAll opTions, and we assImiLate IT to A tradAble aSset WiTh (sigNificaNt) traNsAction costs in whAt foLlows. For bAskEtS of VoLatilITies isOlaTed By the teChniqueS LisTeD ABOvE, we apply the [@jurek] stRaTEGy For log utIlitieS To ThE Basket prOcEss RecoRDIng ouT of sAMpLe perforMance. @jUReK pRopose tO tRade a sTaTioNarY autoREgreSsive pRocess $(x_t)_{T}$ of orDEr $1$ and mean $\mu$ govERned by the equaTIoN $X_{T+1} = \rHO x_t +\sIgmA \varepsilon_T$, wheRE $\abs{\Rho}<1$, bY TaKinG A posiTion $N_T$ iN ThE Asset $x_t$ which is propoRtIonal tO $$\labeL{eq:jurek}
    N_t = \fraC{\rho (\mu-x_t)}{\siGMA^2}w_t$$ In effeCt, thE StRAtegy advocates TakinG a long (resp. SHort) posiTion iN the asseT whenever IT Is below (rEsp. AboVe) iTs lONG-tErm mean, and adjUST the PoSition sIze To accouNt fOr tHe vOlaTiLity of $x_t$ aNd its meaN rEvErSiOn sPeed $\rHO$. Given baSkEt wEiGhtS $y$, we aPPly staNdard aR esTiMaTIon ProceduREs ON The iN-sAmPle pOrtIoN of $y^T\Bx$ to REcoVer estiMates for $\hAt{\rHO}$ and $\HaT{\sIgma}$ and Plug them direcTlY in EquatioN . THis ApproaCH Is illustRated for two baskets in FigURe \[fig:syN\].

-2.5cm![**three SampLe trading ExpErimenTs, uSIng the pCA, spaRse PCa aNd tHE crossING STatIsTics estimaTORs**. \[a\] pool oF 9 vOlatIlity tiMe-series selected usINg oUr fast PCA seleCtiOn prOCEdUre. \[B\] baSKet WeIGhtS EStimated with in-sAmple data uSiNG eIther the eiGEnvEcTor of thE covariAnce mATrix witH smallest
 be simply estimated by co mputing th e lea din g s pa rseeige nvector of $\l a mbda  I-A_0$ where $\lambda $ isbi g gert ha n the  leadin g  e i g env al ue  of  $ A _0 $. Th issparseprincipalcom po nent of thec ov ariance ma tri x $A_0$ shou ldnot be  c onf u sed w ith  ourutiliz a tion o f sparsePC A  in Se c tion \[ s u bs ec:a semidefinite\] as aw ay to recovera vect or so l u tio n f rom the so lu tiono f a pos i ti v e sem i definite prob lem. Note a l sothat t ec hni q ues ba sed o np rin cipal compo nent s do nottake e x plicitl y  varian ce lev els  in to a c co un t w he n  es t im ati n g t he weigh ts  o f a c o-in t e g r ated  re lati onshi p.

@jurek Tr adi ng S t rat egy
- ----- ---- -- ----- ------ 

Whi le  option implied  vol atility i s n ot  di re ctlyt radabl e,itcan besynthes i zed  u s i n gbaskets of call op ti o n s,  and weassimi l at ei t to a t ra dab le a s s et wi th ( s ig nificant ) tran s ac ti on cost sin wha tfol low s. Fo r  bas kets o f volati litie s  isolated by t h e techniquesl is t e da bove , w e apply the  [@j u rek]  str a te gyf or lo g uti li t ie s  to the basket proc es s reco rding  out of sampl e performa n c e . @jurek  pro p os e  to trade a st ation ary autore g ressiveproce ss $(x_t )_{t}$ of o rder $1$  an d m ean  $\ m u $governed by t h e  equ at ion $x_ {t+ 1} = \r hox_t  +\ sig ma  \varepsi lon_t$,wh er e$\ abs {\rho } <1$, byta kin ga p ositi o n $N_t $ intheas se t  $x _t$ whi c hi s  pro po rt iona l t o$$\la bel{ e q:j urek}
   N_t =\fr a c{\r ho  ( \mu-x_t )}{\sigma^2}W _t $$ In effe ct , t he str a t egy advo cates taking a long (re s p. shor t)posit ionin the as set  whene ver it isbelow(resp .abo v e ) its l on g-t er m mean, an d adj ust t he  pos ition s ize to account for the  volatility o f $ x_t$ a nd  it s  m e anre v ers i o n speed $\rho$.  Given bas ke t  w eights $y$ ,  we  a pply st andardAR es t imation  procedur es on the  i n-sa m p leportion of  $y^T\bx $ to reco v er es t im atesfor  $\hat {\ rho }$ an d $\ha t {\s igma} $ andpl ug the m dir ec tly in E quation . This approach  is il lustr ate d for two  ba s ket s in Figu re \ [fig:syn\] .
-2. 5cm![ **T h ree s ampl e  t rad i ng ex peri m ents, usi n gthe P CA , sparse PC A a ndthe C ros s ing St atis tics estimators** .  \[a\] Pool of  9 v o l ati lit y  tim e- series selecte d u si n g  our fas tPCA selecti on proce du r e. \[ b\] Ba sket w eightse s ti m ated w ithin- sample da taus i ng eith er  t h e eige nvec to r of t he cov a rian c e  matrix with sma llest 
 be_simply estimated_by computing the leading_sparse eigenvector_of_$\lambda I-A_0$_where_$\lambda$ is bigger_than the leading_eigenvalue of $A_0$. This_sparse principal component_of_the covariance matrix $A_0$ should not be confused with our utilization of sparse PCA_in_Section \[subsec:asemidefinite\] as_a_way_to recover a vector solution_from the solution of a_positive semidefinite_problem. Note also that techniques based on principal_components_do not take_explicitly variance levels into account when estimating the weights_of a co-integrated relationship.

@jurek Trading Strategy
-----------------------

While_option implied volatility_is_not_directly tradable, it can_be synthesized using baskets of call_options, and we assimilate it to_a tradable asset with (significant) transaction costs_in what follows. For baskets of_volatilities isolated by the techniques_listed above,_we apply the [@jurek] strategy for_log utilities to_the basket_process recording out_of sample performance. @jurek propose to trade_a stationary autoregressive_process $(x_t)_{t}$ of order $1$ and_mean_$\mu$ governed by_the_equation_$x_{t+1} =_\rho x_t +\sigma_\varepsilon_t$,_where $\abs{\rho}<1$,_by_taking a position $N_t$ in the_asset_$x_t$ which is proportional to $$\label{eq:jurek}
 _  N_t =_\frac{\rho_(\mu-x_t)}{\sigma^2}W_t$$ In effect, the_strategy advocates taking a long_(resp. short) position in the asset_whenever it_is below_(resp. above) its long-term mean, and adjust the position size to_account for the volatility of $x_t$_and its mean reversion_speed $\rho$._Given_basket weights $y$,_we_apply standard_AR estimation procedures on the in-sample portion_of $y^T\bx$_to recover estimates for $\hat{\rho}$ and_$\hat{\sigma}$ and plug them_directly_in Equation . This approach is illustrated_for two baskets in Figure \[fig:syn\].

-2.5cm![**Three sample_trading experiments, using the PCA,_sparse_PCA_and the Crossing Statistics estimators**._\[a\] Pool of 9 volatility time-series_selected using our_fast PCA selection procedure. \[b\] Basket weights_estimated_with in-sample data using either the_eigenvector_of the covariance matrix with smallest