A ROBUST NON-LINEAR MULTIVARIATE KALMAN FILTER FOR ARBITRAGE IDENTIFICATION IN HIGH FREQUENCY DATA
|
|
- Bartholomew Armstrong
- 6 years ago
- Views:
Transcription
1 A ROBUST NON-LINEAR MULTIVARIATE KALMAN FILTER FOR ARBITRAGE IDENTIFICATION IN HIGH FREQUENCY DATA P. J. BOLLAND AND J. T. CONNOR London Busness School Deparmen of Decson Scence Sussex Place, Regens Park London NW 4SA Phone -(+44) FAX - (+44) E-Mal - pbolland@lbs.lon.ac.uk E-Mal - jconnor@lbs.lon.ac.uk ABSTRACT We presen a mehodology for modellng real world hgh frequency fnancal daa. The mehodology copes wh he errac arrval of daa and s robus o addve oulers n he daa se. Arbrage prcng relaonshps are formulaed no a lnear sae space represenaon. Arbrage opporunes volae hese prcng relaonshps and are analogous o mulvarae addve oulers. Robus denfcaon/flerng of arbrage opporunes n he daa s accomplshed by Kalman flerng. The sae space model used o descrbe he prcng relaonshps s general enough o handle boh lnear and non-lnear models. The recursve Kalman equaons are adaped o fler ck daa, cope wh he errac arrval of observaons and produce esmaes of all he arbrage prces on every me sep. We demonsrae he mehodology wh a robus neural nework fler appled o foregn exchange rangular arbrage. Tck daa from hree markes s used: $/DM, /$, /DM The fler produces esmaes of he arbrage prce for all exchange raes on every second, ncreasng boh he speed and effcency of arbrage denfcaon. KEYWORDS: Arbrage, Foregn Exchange, Mulvarae Kalman Fler, Neural Nework, Oulers, Robus, Tck Daa.. Inroducon Arbrage s a fundamenal mechansm for achevng effcency n he fnancal markes (Ross 976). An arbrage opporuny occurs when a prce dscrepancy exss beween wo or more hghly relaed asses. The opporuny can be exploed by buyng he under prced asse and sellng he over prced asse, producng a prof whou ncurrng any rsk. Msprcng s rapdly correced n hghly compeve markes (Frenkel and Levch 975,977), herefore arbrage raders need rapd denfcaon, fas ransacons and low ransacon coss. Many arbrage
2 relaonshps have been denfed n he fnancal markes. Our mehodology can be appled o any sysem of lnear arbrage prcng relaonshps. Secon. descrbes he rangular foregn exchange arbrage we use o demonsrae he mehodology. Prevous sudes of arbrage denfcaon have manly been lmed o examnng daly daa and so mgh mss many of avalable nraday opporunes. Sudes ha have examned nraday daa (Rhee and Chang 992) have been lmed o examnng only a mnue fracon of he daa because of he need for smulaneous observaons. The mehodology we presen allows arbrage opporunes o be denfed wh rregular (non-smulaneous) observaons. Irregular mes seres presens a serous problem o convenonal modellng mehodologes. Several mehodologes for dealng wh errac daa have been suggesed n he leraure. Muller e al 990, sugges mehods of lnear nerpolaon beween errac observaons o oban a regular homogenous mes seres. Oher auhors (Ghysels and Jasak 995) have favoured non-lnear me deformaon ( busness-me or ck-me ), however hs mehodology has no smple equvalen for mulvarae seres. The mehodology we presen descrbes he dynamcs of fundamenal underlyng arbrage saes whch are observed as errac nosy exchange raes. We rea he errac arrval as a mssng daa problem. The Kalman fler descrbed s dscree, as he daa s only provde n quansed me seps (.e. seconds), however he mehodology could be exended o connuous me problems wh he Kalman-Bucy fler (Medch 969). The sae space represenaon descrbed n secon 2 allows us model he sysem a he maxmum resoluon of he avalable daa (Reuers daa quoed by he second). Convenonal modellng mehodologes may also be napproprae for modellng ck daa as he dsrbuon s ofen heavy aled (Dacorogna 995). Fnancal daa, especally quoaons, are prone o daa corrupon and oulers. Chung 99, dscovered 0.25% of he MMI fuures quoes were ousde of he daly hgh and low and are herefore serous daa corrupon s. Secon 3 deals our robus mehodology whch s smlar o ha descrbed by Masrelez and Marn 975 and Marn and Vandaele 983. The sae space represenaon s capable of ncorporang boh lnear and nonlnear models. The esmaon of he models s performed usng an E.M. algorhm descrbed n secon 4, whch was nroduced by Dempser, Lard and Rubn 977 o esmae parameers of models when some of he daa s mssng.
3 The mehodology we presen s suable for real world fnancal daa and ncreases boh he speed and effcency of arbrage denfcaon. We demonsrae he fler on $/, $/DM, /DM daa from , he resuls of he esmaon and flerng are shown n secon 5.. Foregn Exchange Arbrage We examne foregn exchange rangulaon s for arbrage opporunes (alhough he same mehodology can be appled o many varees of arbrage). In he absence of ransacons coss and bd-ask spread he followng equlbrum relaonshps mus hold for currency raes, EX(0,) EX(,2) EX(2,0) = EX(0,) EX(,2) EX(2,3) EX(3,0) = EX(0,) EX(,2) EX(2,3)......EX(m,0) =, () where EX(,j) represens he spo rae for currency j when expressed n uns of currency. If he equlbrum relaonshps n Eq.() hold, hen a sngle counres m exchange raes can be used o produce esmaes of all he cross raes, EX(, j) = EX(0, j) / EX(0,), n hs paper he US Dollar s used as he base currency. Takng logarhms of he rangular relaonshps, allows he cross raes o be expressed as log(ex(,2)) = log(ex(0,2)) - log(ex(0,)) log(ex(, j)) = log(ex(0, j)) - log(ex(0,)) log(ex(m -, m)) = log(ex(0,m)) - log(ex(0,m -)). (2) If he addve relaonshps of Eq.(2) are volaed, an arbrage opporuny exss where rskless, profable ransacons can occur. Volaons of he rangular relaonshps are analogous o an ouler n he daa se, he larger he msprcng he larger he ouler. When marke frcon' s are ncluded (ransacon coss and bd-ask spread) slgh msprcng can occur whn small bands around he arbrage prce. In he followng secon he rangular currency relaonshps are encoded whn a
4 sae space form and a mulvarae Kalman fler s used o denfy any sgnfcan volaons of Eq.(2). 2. Space Represenaon of FX Relaonshps The mehodology we presen below descrbes how oulers can be robusly denfed/flered n mulvarae non-lnear daa. In hs applcaon he oulers ha he Kalman fler denfes are suaons n whch an arbrage opporunes exs. The Kalman fler has been adaped o fler ck daa and o updae he esmaes of he exchange raes every me sep. The Kalman fler used s general enough o handle boh lnear and non-lnear models. For non-lnear models a pon-wse lnearzaon s performed o predc he Kalman fler s sae changes, and o updae he recursve esmaes of he error predcon covarance (Connor, Marn, Alas 994). The parameers of he models used n he Kalman fler are robusly esmaed from cleaned daa, descrbed n secon 4. The observaon vecor z n he sae space model represens he logarhm of each exchange rae observed. If all possble exchange raes (cks) are observed n a gven second hen, z s gven by z z z z z z, (3) ( m m m m T z = ( 0, ) (, 0, 2 ) (,, 0, ) (,, 2 ) (, ) (, ),,,, ) (, j where z ) = log(ex(, j)). Usually only a subse of Eq.(3) are observed. The elemens of z come n wo prncple groups : The log of he m exchange raes for he base currency (0,j), ( 0, ) ( 0, 2) ( 0, m) z, z,, z = Log ( Base Raes). (4) The log of he correspondng cross raes (,j), (, 2) (, m) ( 2, 3) ( m, m) z,, z, z,, z = Log ( Cross Raes). (5) The exchange rae msprcng problem can be formulaed no a famlar sae space model, x = f ( x ) + e, (6) z = Hx + v. (7)
5 The sae vecor x n Eq.(6) represens he log of he arbrage value of he m exchange raes for he base currency ($) as well as he auo-regressve srucure of he sysem. The sae ranson vecor f ( x ) n Eq.(6) represens he sysem dynamcs ha may be lnear or non-lnear (n he case of a lnear sysem f ( x ) s smply he sae ranson marx Φ ). The observaon marx, H n Eq.(7), exracs he base raes and uses he logarhmc arbrage equaons o esmaes he cross raes The sysem descrbed n Eq.(6) and Eq.(7) have wo ypes of drvng nose, e he sae nose and v he observaon nose. The sae nose e represens he varaon due o he exchange raes underlyng arbrage dynamcs. The observaon nose v has wo componens v = u + w, he frs componen u, represens he varaon caused by he ransacon coss and bd ask spread, whch allow he prce o move freely whn unprofable bounds. The second componen w, represens he addve oulers whn he daa (wheher hey are daa corrupon's or marke msprce anomales). The sae ranson vecor f ( x ) n Eq.(6), can be descrbed as a non-lnear mulvarae auo-regressve (NMAR) process for each of he base currency s exchange raes. The mulvarae auoregressve process, NMAR( p, p 2,, p m ) s defned ( ) ( ) ( ) by, 2 2 m m x = f ( x,, x, x,, x,, x,, x ) + ε ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) p ( 2 ) p ( m) p m m x = f ( x,, x, x,, x,, x,, x ) + ε ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( 2) ( ) p ( 2 ) p ( m) p ( m) ( m) ( ) ( ) ( 2) ( 2) ( m) ( m) ( m) x = f ( x,, x, x,, x,, x,, x ) + ε, ( ) p ( 2) p ( m) p ( ) ( 2) ( m) where x, x,, x are he log exchange raes for he base currency, and p () s he number of auoregressve erms for he h exchange rae. In Eq.(8) f ( ) s denoe non-lnear funcons governng each ndvdual exchange rae. The sae vecor x, he sae ranson vecor f ( x ) and he dsurbance vecor e n Eq.(6) are defned as follows: (8) x x x x x x m, (9) ( ) ( ) ( 2) ( 2) ( m) ( m) T x = (,, ( ), p +,, ( 2),, p +,, ( ) ) p + ( ) ( ) ( ) f ( x ) = ( f ( x ), x,, x, ( ) p + ( 2) ( 2) ( 2) f ( x ), x,, x,, ( 2) p + ( m) ( m) ( m) T, f ( x ), x,, x m ), ( ) p + (0)
6 ( ) ( 2) ( m) T e = ( ε,,, ε,,, ε,, ), () ( j where ε ), he random sae nose assocaed wh he j h exchange rae, appears n he + p j = ( ) poson n he dsurbance vecor e (where p ( 0) = 0 ). The observaon marx H n Eq.(7) exracs he base raes and he cross raes from he sae vecor x. Each of he rows of H relae o a specfc exchange rae, he rows are defned as follows, ( T m T T m T m m T H = [ h 0, ) (, ) (, ) (, ) (, ),, h 0, h 2,, h,, h ]. (2) For base currency exchange rae (0,j) : ( 0, j ) h [ ], (3) = ( ) ( j ) ( j ) ( m) x p x p x p x p so each h ( 0, j) exracs he base currency rae x j from x. For cross currency exchange raes (,j) : (, j) h = [ 0 0 ( ) 0 0 ( ) 0 ( j ) 0 0 ( m) ], (4) (, j) so each h h (, j ) j x = x x ). x p x p x p x p esmaes he cross rae (,j) usng he addve log relaonshps (.e. For regular (evenly spaced) me seres all he raes would be observed on every me sep. Tck daa, however, requres a mehodology capable of modellng rregular me seres. On any gven second only currences for whch a rade (ck) has occurred ener he observaon vecor and only he rows of he observaon marx whch correspond o an acual observaon are used o updae he flerng equaons. The observaon nose vecor v, n Eq.(7) also reconfgures s dmensonaly o correspond o he acual rades ha occur. Ths gves rse o several possble suaons : No observaons on any marke, = [ NULL], H = [ NULL], v = [ NULL]. z
7 The observaon vecor z, he observaon marx H and he observaon nose v are se o null. One or more markes produce observaons, ( 0, ) ( j, k ) e.g. base rae z and cross rae z are raded. z z = ( 0, ) h ( z j, k ), H = ( 0, ) j k v, v (, ) = h v ( 0, ) ( j, k ). All he markes produce observaons n one second, ( 0, ) z ( 0, ) h ( 0, m) z ( 0, m) h (, 2) (, 2) z z h =, H =, v (, m) (, m) z h ( z m, m ) ( m, m) h = v v v ( 0, ) ( 0, m) (, 2) v v (, m) ( m, m ). Expandng and conracng he observaon equaon n hs way, allows he sae space model o cope wh he errac arrval of ck daa and mmedaely ncorporae all new nformaon o updae he sae esmaes for all exchange raes. The mehodology produces an esmae of he saes (exchange raes), an esmae of he assocaed predcon error covarance, as well as he predcons of fuure saes a every second, regardless of any ck beng observed. 3. Robus Kalman Fler The underlyng saes x n Eq.(6), are unknown, hey are esmaed by a robus Kalman fler (Kalman 96). Usng robus mehodologes proecs he modellng procedure from serous performance degradaon caused by ll condoned daa. The recursve non-gaussan Kalman fler equaons as descrbed by Masrelez 975 and Marn and Vandaele 983 are dealed below. The robus one sep ahead predced sae vecor x and he predced observaon vecor z are gven by, x ( ~ = f x ), (5) z = H x. (6)
8 where ~ x s he flered sae vecor. The dsrbuon erms e and v n Eq.(6) and Eq.(7) are assumed o be zero mean, serally ndependen and muually ndependen, however no assumpons abou her dsrbuons are made. The covarance marces of e and v are denoed by Q = E( e e ) and R = E( v v ) respecvely. The modellng mehodology we employ assumes ha he nose covarance marces reman consan over me. For fnancal daa hs assumpon may be nvald (Ruz 994). The mehodology we presen can be exend o ncorporae sochasc volaly, see Harvey Ruz and Shephard 992. In an effor o lm he mpac of sochasc volaly he esmaon of Q and R was made usng a rollng wndow, see secon 4. The robusly flered esmae of he sae vecor ~ x and he predcon error covarance marx M are defned by he followng recursve relaonshps, ~ x = x + M H g ( z ), (7) M + = Φ P Φ + Q, (8) P = M M H G z H M, (9) ( ) where z = z z s he nnovaons vecor (he observed resdual), g ( z ) s he score funcon of he nnovaons wh componens, { g ( z )} { z Z } p = [ p{ z Z }], (20) ( z ) and G ( z ) s defned as he dfferenal of he score funcon, wh elemens, { ( z )} G j { g ( z )} = ( z ) j. (2) For a non-lnear sae space model he sae ranson marx Φ n Eq.(8) s esmaed by a pon-wse lnearzaon of he non-lnear model. The elemens of Φ are he paral dervaves of f evaluaed abou he robusly flered esmaes of he sae vecor ~ x, f ( x ) Φ j =, x ~ j = x x. (22)
9 In he sandard Kalman fler he densy funcon for he nnovaons s assumed o be Gaussan. In order o oban robusness we assume z has a symmerc heavyaled densy funcon. The score funcon g ( z ) for a Gaussan nnovaon process s lnear. For a heavy aled densy g ( z ) s gven by a non-lnear gan funcon ha lms he nfluence of large nnovaons. There s some laude gven n choce of g ( z ) n he above equaons. The score funcons for he Gaussan dsrbuon, Huber s leas nformave dsrbuon (Huber 98) and he Hampel re-descendng funcon are shown n fgure (for he one-dmensonal case). g(innovaons) Score Funcons Innovaons Gaussan Huber Hampel Fgure : Score Funcons: Gaussan, Huber, Hampel. For he case of Gaussan nnovaons he score funcon g ( z ) s gven by, g ( z ) = ( H M H + R ) z. (23) The dervaons of g ( z ) for he n-dmensonal Huber and Hampel denses s gven n Bolland and Connor 995. The sze of he nnovaons s he crcal value whch deermnes wheher he observaon s an ouler (arbrage opporuny). The magnude of he ouler s 2 defned by r = ( z z ) ( z z ) where s he covarance marx of he nnovaons. The measuremen r allows us o se a defnon for an ouler, so ha only msprcng of suffcen magnude o allow for profable rades are denfed. 4. Model Specfcaon and Esmaon
10 In order o produce he robus Kalman fler, esmaes of f ( x ), Q and R need o be obaned. The sae ranson funcon f ( x ) can be approxmaed wh many dfferen modellng mehodologes. There s a large body of emprcal evdence o sugges ha he domnan srucure n f ( x ) wll be a mean reverng process (Fama 965). The mean reverson could be an arfcal arefac of he daa. Roll 984, showed ha bd-ask bounce nduced srong negave auo-correlaon n fnancal daa. Tme seres of marke prces conan boh bd and ask prces so f no new nformaon arrves he rue value remans consan, any observed varaon s caused by he dfference n bd and ask prce. Bouncng beween bd and ask prces gves rse o a srong negave auo-correlaon shown n fgure 2. Ask Prce Value Bd Prce - + Fgure 2: Bd-Ask Bounce. To lm he mpac of he bd-ask bounce, md-prces were modelled. The mean reverson could also be an arefac of he prce quansaon (prces quoed n dscree uns). The saes of he sysem descrbed are he arbrage values of he exchange raes. The frs dfferences of he sae were aken o produce a saonary seres. Predcng sae changes raher han here levels requres a slgh re-formulaon of he f ( x ) s n Eq.(0). The sae ransons are formed by wo componens, a ( ) random walk componen, x (he prevous sae), plus he sae changes ( ) ( x x 2 ), so Eq.(0) becomes, d ( ) ( ) ( ) ( ) f ( x ) = ( x + d ( x x ), x,, x, 2 ( ) p + ( 2) ( 2) ( 2) ( 2) x + d ( x x ), x,, x ( 2 ),, 2 p + ( m) ( m) ( m) ( m) T, x + d ( x x ), x,, x ( m) ), 2 p + (24) ( where d ) ( x x 2 ) represens he NMAR srucure of he sae changes.
11 Tck daa for hree currences $/DM, /$, /DM ( ) was used o demonsrae he mehodology. The fundamenal saes are herefore he arbrage values of $/DM, /$, so x = (log($ / DM ),log($ / )). The saes were esmaed by several dfferen models: Nave random walk model (no fler), ~ x ~ x = Random walk Kalman fler, ~ x = ~ x + M H g ( z ) Lnear MAR(,), ~ x = f ( ~ x ) + M H g ( z ) Neural Nework NMAR(n,n), ~ x ( ~,, ~ = f x x n ) + M H g ( z ) To produce an esmaon se for an MAR(n,n) we requre daa examples wh n consecuve sae changes across all saes. Due o he errac naure of he ck daa consrucng an esmaon daa se whou mssng observaons for wh more han one lag proved dffcul. An nal lnear MAR(,) model was esmaed from he daa where cks occurred on consecuve seconds for boh he /$ and /DM exchange raes. Usng he MAR(,) as an esmae of f ( x ), he robus Kalman fler was appled o he whole daa se produced a flered daa se. Ths daa se conaned he flered values of he acual observed sae changes and he esmaed sae changes produced by he smple lnear model. An MAR(n,n) esmaon se could now be consruced from he flered daa, wh he esmaed sae changes fllng he problemac mssng observaons. To avod he problems of nonsaonary n he dynamcs over he whole daa se (wo years), he parameers were esmaed on a rollng wndow of one radng week, and re-esmaed daly. 4. Esmaon of Model Parameers An esmaon maxmsaon (EM) algorhm s employed o esmae he neural nework parameers, denoed λ, of Eq.(0) or equvalenly Eq.(24). The EM algorhm, see Dempser, Lard, and Rubn (977), s he sandard approach when esmang model parameers wh mssng daa. The EM algorhm has been used n he neural nework communy before, see for example Jordan and Jacobs, and Connor, Marn, and Alas (994). The E-Sep, gven n secon 4., esmaes he mssng daa. Wh he esmaed mssng daa assumed o be rue, he parameers are hen chosen by way of
12 maxmsng he lkelhood. Ths procedure s erave wh new parameer esmaes gvng rse o new esmaes of mssng daa whch n urn gve rse newer parameer esmaes. 4. E-Sep Durng he esmaon sep, he mssng daa, namely he x, v, and e, of Eq.(6) and Eq.(7) mus be esmaed. Ths s accomplshed usng he robus Kalman fler of secon 3. The esmaed mssng daa s denoed ~ x, ~ v, and ~ e. 4.2 M-Sep The robus lkelhood for he sysem defned by Eq.(6) and Eq.(7) s gven by where p{ z Z } N { z Z } l( λ ) = p = (25), s a funcon of r 2, he magnude of he nnovaons 2 r = z j z j Z j z j z j Z ( ~ (, )) ( ~ T λ (, λ )) and = H M H + R s he covarance marx. Ths has been derved by De Jong 988 for he case where nal sae esmaes and nose varances are consdered. The log lkelhood s defned by N { z Z } L( λ ) = log( p ) j= The parameers o maxmse he lkelhood, arg max N j λ = log( p{ z Z }). (27) λ = Smple graden descen s used o maxmse (27), hs s done by gnorng he dependence of on λ N L λ p{ } λ ( ) log( z Z ) ( ~,, ~ p, ) = H x x x (28) λ z λ = Nong ha { g( z )} { z Z } (26) log( p ) = from Eq.(20) and usng Eq.(7), ( z ) g( z ) H M ~ x x ( ~ x,, ~ = x, λ ) (29) ( ) p whch can be used wh Eq.(28) o ge he graden n erms of clean daa L ( λ ) N λ ( ~ p p λ ) λ ( ~,, ~, ) ( ~,, ~, ) = x x x x M x x x λ = whch s equvalen o dong back-propagaon on clean daa. (30)
13 Ouler Observaons 2 When an ouler, z, s observed r wll be very large causng he erm 2 g. ( r ) lm he conrbuon of z o he lkelhood n (26). As menoned n secon 3, he predcon varance reflecng greaer uncerany n he fundamenal raes. Fuure conrbuons o he lkelhood wll be effeced by hs ncreased uncerany. Unvarae Case If only a sngle quoe s observed a a gven me, he assocaed graden s gven by L ( λ) z ( λ) = 2( z z ( λ)) szz (3)) λ λ whch f z corresponds o a fundamenal dollar denomnaed rae s he same as found n he unvarae case frs explored n Connor, Marn and Alas (994). Cross Currency If he quoe corresponds o a cross currency, relaed o he dollar denomnaed raes (, j) ( ) ( j) by z = x x, he graden wll be composed of conrbuons from all he predcors relaed hrough H M ( k ) L ( λ) (, j) (, j) x ( λ) = 2( z z ( λ)) ( M k,, M k, j, ). (32) λ λ k Several Quoes of he Same Currency If several quoes are observed for one of he underlyng dollar denomnaed raes, he learnng algorhm smplfes grealy. The predcon covarance marx and correspondng nverse are gven by: a b b b c d d d b a b b d c d d = =, (33) b b a b d d c d b b b a d d d c where c and d are derved from ac + ( N ) bd = and ad + bc + ( N 2) bd = 0 whch allows he reducon of Eq.(33) o he smpler form N L ( λ) ( ) x ( λ) = ( x x ( λ)) k N λ 2 (34) N λ = o
14 where k = c + ( N ) d. Snce several quoes are avalable, he addve nose s N smoohed ou and one has more confdence n he average of he quoes, N ( ) x, han any of he quoes would be gven alone. Ths added confdence s N = expressed n erms of he sronger graden n () where k N > k for N >. No Observaons Due o he errac naure of ck daa, ofen here wll be no observaons durng a gven perod. For hs mssng daa, here wll be no conrbuon o he lkelhood gven n (26). Bu as n he case wh exreme oulers, he uncerany n fuure predcons, wll grow and effec he lkelhood of fuure observaons. 4.2 Esmaon of Nose Varances An erave procedure was appled o produce esmaes of Q and R. The observaon nose covarance R, has wo componens, u he marke frcon's (bd-ask spread, ransacon coss ec.) and w he addve oulers (prcng anomales). As flerng procedure nerpolaes hrough addve oulers, he esmae of R s only dependen on he frs componen u, R = E( u u ). Inal esmae of observaon nose covarance R, and he sae nose covarance Q, were produced by maxmsng he lkelhood of he MAR(,) sae space model. These esmaes of R and Q were refned durng he neural nework esmaon process by repeaed applcaon of maxmum lkelhood. 5. Resuls The ck daa was obaned from Reuers, Aprl 993-Aprl 995, he md-prce was modelled. The nal lnear MAR(,) model produced had a sae change ransons marx gven by, ( ) ( ) ( ) x.. x x f ( x ) = ( ) ( ) ( ) x x x 2 (35) The nal daa se from whch he MAR(,) parameer s were esmaed, was consruced from regons of consecuve cks (500 regons of consecuve pons coverng approxmaely wo weeks radng were used). The srengh of he dagonal erms demonsrae he meanng reverng. Ths model was used n he Kalman fler
15 o produce a flered daa se whou mssng observaons, from whch subsequen models were produced. Table show he flerng resuls for he neural nework Kalman fler and he nave random walk hypohess. The Kalman fler produces superor resuls for boh he mean squared error (MSE) and he robus medan absolue devaons (MAD). Table : Model Comparson. Model MSE MAD RW Model 2.69 (0-4 ) (No Fler) Neural Nework Kalman Fler.5 (0-4 ) Fgures 3 and 4 demonsrae he Kalman fler denfyng oulers. In Fgure 3 he flered esmaed saes for he $/DM exchange rae are represened by he sold lne, he acual observed rades by crcles, and rades occurrng on he oher exchange raes by vercal lnes. The flered saes represen he esmae of he rue arbrage value (no he bd or ask values). The effec of marke frcon s has been ncorporaed no he esmae R, so flered saes are always whn he acual bd and ask values observed. The Kalman fler denfes an ouler and uses pure predcon for he esmae of he $/DM rae a me and does no follow he spurous prce movemen. The msprcng s flered and classfed as an ouler ( r = 4. 65) by he robus algorhm presened here. Fgure 3 also demonsraes how new nformaon occurrng on any of he exchange raes s mmedaely ncorporaed n he robus esmaes of all he saes. A me , he vercal doed lne ndcaes he observaon of a DM/ rade, he esmae of he $/ sae s nsanly updae o ncorporae he affec of he rse n he DM/.
16 Robus Idenfcaon of Arbrage Opporunes $/DM r = Tme Fgure 3: Idenfcaon of Marke Msprcng. $/DM Tcks $/DM Saes BP/$ Tcks DM/BP Tcks Table 2: Tck Daa (bold) and Esmaed Raes (normal). $/DM Tcks /$ Tcks /Dm Tcks $/DM Esmae /$ Esmae /DM Esmae SD (0-5 ) Ln($/D) SD (0-5 ) Ln( /$) Table 2 demonsraes he ably of he Kalman fler o deal wh he errac arrval of ck daa. Bold fon exchange raes represen seconds when cks are observed, and normal fon represens he Kalman fler esmaes. When ncomplee observaon vecors are observed he Kalman fler uses he mulvarae auoregressve srucure o esmae he unseen raes. In seconds where no observaons are observed he fler uses pure predcon o esmae he mssng currency raes. The las wo columns n able 2 show he esmaed predcon sandard devaons
17 for he log of he saes. Agan, bold fon ndcaes he seconds where cks where observed, and normal fon ndcaes he recursve esmaes of he Kalman fler. The predcon error sandard devaons grows seadly n perods where no cks are observed, and collapses o he one second sae predcon error sandard devaon for he predcon based on new nformaon (cks). $/DM Robus Idenfcaon of Arbrage Opporunes 5:2:58 5:3:06 5:3:5 5:3:24 5:3:32 5:3:4 5:3:49 5:3:58 5:4:07 $/DM $/DM :23:3 0:23:48 0:24:06 0:24:23 0:24:40 0:24:58 5:4:24 5:4:4 5:4:59 5:5:6 5:5:33 5:5:50 $/DM :07:2 5:07:29 5:07:47 5:08:04 5:08:2 5:08:38 Fgure 4: Idenfcaon of Marke Msprcng.
18 6. Concluson We have shown he effecveness of one flerng approach for denfyng arbrage opporunes on currency ck daa. The mehodology s deally sued o he poor daa qualy of he fnancal markes (errac arrval and addve oulers). The Kalman fler produces rae esmaes every second wheher or no any cks are acually observed, hs ncreases boh he speed and effcency of denfyng arbrage opporunes. I s sraghforward o exend he above analyss o many more exchange raes and cross raes, ncreasng he possbly of fndng msprcng. In addon hs mehodology could be readly appled o all forms of arbrage whch are descrbed by smlar ses of prce relaonshps. 7. References Bolland, P. J. and Connor, J. T., "Mulvarae Non-Lnear Kalman Flers", Techncal Repor, London Busness School, June 995. Connor, J. T., Marn, R. D., and Alas, L. E. Recurren Neural Neworks and Robus Tme Seres Predcon, IEEE Transacons on Neural Neworks, pp , March 994. Chung, P. Y., A Transacons Daa Tes of Sock Index Fuures Marke Effcency and Index Arbrage Profably, Journal of Fnance, Vol. 46, pp , Dec. 99. Dacorogna, M. M., Prce Behavour and Models for Hgh Frequency Daa n Fnance, Proceedngs of he NNCM conference, London, England, Oc. -3, 995. Dempser, A. P., Lard, N. M., and Rubn, D. B. Maxmum lkelhood from ncomplee daa va he EM algorhm. Journal of he Royal Sascal Socey, B, 39, -38. (977). De Jong, P., The Maxmum Lkelhood of a Sae Space Model, Bomerca, Vol. 75,, pp 65-69, 988. Fama, E. F., The Behavour of Sock Marke Prces, Journal of Busness, Vol. 38, pp 34-05, Jan 965.
19 Frenkel, J. A., Levch, R. M. Covered Ineres Arbrage: Unexploed Profs?, Journal of Polcal Economy, Vol. 83, pp , 975. Frenkel, J. A., Levch, R. M. Transacon Coss and Ineres Arbrage: Tranqul Versus Turbulen Perods, Journal of Polcal Economy, Vol. 85, pp , 977. Ghysels, E., Jasak, J., Sochasc Volaly and Tme Deformaon: An Applcaon of Tradng Volume and Leverage Effecs, Proceedngs of he HFDF-I conference, Zurch, Swzerland, March 29-3, Vol., pp -4, 995, Harvey, A. C., Ruz, E., Shephard N. G., Mulvarae Sochasc Varance Models, Fnancal Markes Group dscusson paper, London School of Economcs, 992. Huber, P. J. Robus Sascs, Wley, New York, 98. Jordan, M.I. and Jacobs, R.A. Herarchcal mxures of expers and he EM algorhm, submed o Neural Compuaon. Kalman, R. E. A new approach o lnear flerng and predcon problems., Trans. ASME J. Basc Eng. Seres D, Vol. 82, pp , Mar 96. Marn, R. D., Samarov, A., Vandaele W. Robus Mehods for Tme Seres, In Zellner, A. Appled Tme Seres Analyss of Economc Daa, U.S. Bureau of he Census, Washngon, pp 53-69, 983. Masrelez, C. J. Approxmae Non-Gaussan Flerng wh Lnear Sae and Observaon Relaons, IEEE Transacons on Auomac Conrol, pp. 07-0, Feb Medch, J. S. Sochasc Opmal Lnear Esmaon and Conrol, McGraw-Hll, New York, 969. Muller, U. A., Dacorogna, M. M., Oslen, R. B., Pce, O. V., Schwarz, M., Morgenegg, C., Sascal Sudy of Foregn Exchange Raes, Emprcal Evdence of a Prce Change Scalng Law, and Inraday Analyss, Journal of Bankng and Fnance, Vol. 4, , 990.
20 Rhee, S. G., Chang, R. P. Inra Day Arbrage Opporunes n Foregn Exchange and Eurocurrency Markes", Journal of Fnance, Vol. 47, pp , Mar Roll, R. "A Smple Implc Measure of he Effecve Bd-Ask Spread n an Effcen Marke", Journal of Fnance, Vol. 39, pp 27-40, Sep Ross, S. A. The Arbrage Theory Of Capal Asse Prcng, Journal of Economc Theory, Vol. 3, pp , Dec Ruz, E., Quas-maxmum Lkelhood Esmaon of Sochasc Volaly Models, Journal of Economercs, Vol. 63, pp , 993.
John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany
Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy
More informationRELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA
RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA Mchaela Chocholaá Unversy of Economcs Braslava, Slovaka Inroducon (1) one of he characersc feaures of sock reurns
More informationVariants of Pegasos. December 11, 2009
Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4
CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped
More information( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model
BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng
More informationDynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005
Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc
More informationV.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS
R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon
More informationRobustness Experiments with Two Variance Components
Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference
More informationJanuary Examinations 2012
Page of 5 EC79 January Examnaons No. of Pages: 5 No. of Quesons: 8 Subjec ECONOMICS (POSTGRADUATE) Tle of Paper EC79 QUANTITATIVE METHODS FOR BUSINESS AND FINANCE Tme Allowed Two Hours ( hours) Insrucons
More informationNew M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)
Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor
More informationTHE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS
THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he
More informationVolatility Interpolation
Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local
More informationTime-interval analysis of β decay. V. Horvat and J. C. Hardy
Tme-nerval analyss of β decay V. Horva and J. C. Hardy Work on he even analyss of β decay [1] connued and resuled n he developmen of a novel mehod of bea-decay me-nerval analyss ha produces hghly accurae
More informationSolution in semi infinite diffusion couples (error function analysis)
Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of
More informationData Collection Definitions of Variables - Conceptualize vs Operationalize Sample Selection Criteria Source of Data Consistency of Data
Apply Sascs and Economercs n Fnancal Research Obj. of Sudy & Hypoheses Tesng From framework objecves of sudy are needed o clarfy, hen, n research mehodology he hypoheses esng are saed, ncludng esng mehods.
More informationAnalysis And Evaluation of Econometric Time Series Models: Dynamic Transfer Function Approach
1 Appeared n Proceedng of he 62 h Annual Sesson of he SLAAS (2006) pp 96. Analyss And Evaluaon of Economerc Tme Seres Models: Dynamc Transfer Funcon Approach T.M.J.A.COORAY Deparmen of Mahemacs Unversy
More informationComparison of Supervised & Unsupervised Learning in βs Estimation between Stocks and the S&P500
Comparson of Supervsed & Unsupervsed Learnng n βs Esmaon beween Socks and he S&P500 J. We, Y. Hassd, J. Edery, A. Becker, Sanford Unversy T I. INTRODUCTION HE goal of our proec s o analyze he relaonshps
More informationAdvanced time-series analysis (University of Lund, Economic History Department)
Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng
More informationLecture 6: Learning for Control (Generalised Linear Regression)
Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson
More informationEcon107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)
Econ7 Appled Economercs Topc 5: Specfcaon: Choosng Independen Varables (Sudenmund, Chaper 6 Specfcaon errors ha we wll deal wh: wrong ndependen varable; wrong funconal form. Ths lecure deals wh wrong ndependen
More informationTime Scale Evaluation of Economic Forecasts
CENTRAL BANK OF CYPRUS EUROSYSTEM WORKING PAPER SERIES Tme Scale Evaluaon of Economc Forecass Anons Mchs February 2014 Worng Paper 2014-01 Cenral Ban of Cyprus Worng Papers presen wor n progress by cenral
More informationRobust and Accurate Cancer Classification with Gene Expression Profiling
Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem
More informationCHAPTER 10: LINEAR DISCRIMINATION
CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g
More informationStandard Error of Technical Cost Incorporating Parameter Uncertainty
Sandard rror of echncal Cos Incorporang Parameer Uncerany Chrsopher Moron Insurance Ausrala Group Presened o he Acuares Insue General Insurance Semnar 3 ovember 0 Sydney hs paper has been prepared for
More informationIn the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!
ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal
More informationMachine Learning Linear Regression
Machne Learnng Lnear Regresson Lesson 3 Lnear Regresson Bascs of Regresson Leas Squares esmaon Polynomal Regresson Bass funcons Regresson model Regularzed Regresson Sascal Regresson Mamum Lkelhood (ML)
More informationCHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING
CHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING 4. Inroducon The repeaed measures sudy s a very commonly used expermenal desgn n oxcy esng because no only allows one o nvesgae he effecs of he oxcans,
More informationDepartment of Economics University of Toronto
Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of
More informationF-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction
ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or
More informationForecasting customer behaviour in a multi-service financial organisation: a profitability perspective
Forecasng cusomer behavour n a mul-servce fnancal organsaon: a profably perspecve A. Audzeyeva, Unversy of Leeds & Naonal Ausrala Group Europe, UK B. Summers, Unversy of Leeds, UK K.R. Schenk-Hoppé, Unversy
More informationSingle-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method
10 h US Naonal Congress on Compuaonal Mechancs Columbus, Oho 16-19, 2009 Sngle-loop Sysem Relably-Based Desgn & Topology Opmzaon (SRBDO/SRBTO): A Marx-based Sysem Relably (MSR) Mehod Tam Nguyen, Junho
More informationAdditive Outliers (AO) and Innovative Outliers (IO) in GARCH (1, 1) Processes
Addve Oulers (AO) and Innovave Oulers (IO) n GARCH (, ) Processes MOHAMMAD SAID ZAINOL, SITI MERIAM ZAHARI, KAMARULZAMMAN IBRAHIM AZAMI ZAHARIM, K. SOPIAN Cener of Sudes for Decson Scences, FSKM, Unvers
More informationCH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal
More informationFall 2010 Graduate Course on Dynamic Learning
Fall 200 Graduae Course on Dynamc Learnng Chaper 4: Parcle Flers Sepember 27, 200 Byoung-Tak Zhang School of Compuer Scence and Engneerng & Cognve Scence and Bran Scence Programs Seoul aonal Unversy hp://b.snu.ac.kr/~bzhang/
More informationBayesian Inference of the GARCH model with Rational Errors
0 Inernaonal Conference on Economcs, Busness and Markeng Managemen IPEDR vol.9 (0) (0) IACSIT Press, Sngapore Bayesan Inference of he GARCH model wh Raonal Errors Tesuya Takash + and Tng Tng Chen Hroshma
More informationLecture VI Regression
Lecure VI Regresson (Lnear Mehods for Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure VI: MLSC - Dr. Sehu Vjayakumar Lnear Regresson Model M
More informationOn One Analytic Method of. Constructing Program Controls
Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna
More informationNotes on the stability of dynamic systems and the use of Eigen Values.
Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon
More informationComb Filters. Comb Filters
The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of
More informationBernoulli process with 282 ky periodicity is detected in the R-N reversals of the earth s magnetic field
Submed o: Suden Essay Awards n Magnecs Bernoull process wh 8 ky perodcy s deeced n he R-N reversals of he earh s magnec feld Jozsef Gara Deparmen of Earh Scences Florda Inernaonal Unversy Unversy Park,
More information2. SPATIALLY LAGGED DEPENDENT VARIABLES
2. SPATIALLY LAGGED DEPENDENT VARIABLES In hs chaper, we descrbe a sascal model ha ncorporaes spaal dependence explcly by addng a spaally lagged dependen varable y on he rgh-hand sde of he regresson equaon.
More informationUNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION
INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he
More informationClustering (Bishop ch 9)
Cluserng (Bshop ch 9) Reference: Daa Mnng by Margare Dunham (a slde source) 1 Cluserng Cluserng s unsupervsed learnng, here are no class labels Wan o fnd groups of smlar nsances Ofen use a dsance measure
More informationAn introduction to Support Vector Machine
An nroducon o Suppor Vecor Machne 報告者 : 黃立德 References: Smon Haykn, "Neural Neworks: a comprehensve foundaon, second edon, 999, Chaper 2,6 Nello Chrsann, John Shawe-Tayer, An Inroducon o Suppor Vecor Machnes,
More informationFI 3103 Quantum Physics
/9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure
More informationChapter 8 Dynamic Models
Chaper 8 Dnamc odels 8. Inroducon 8. Seral correlaon models 8.3 Cross-seconal correlaons and me-seres crosssecon models 8.4 me-varng coeffcens 8.5 Kalman fler approach 8. Inroducon When s mporan o consder
More informationTHEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that
THEORETICAL AUTOCORRELATIONS Cov( y, y ) E( y E( y))( y E( y)) ρ = = Var( y) E( y E( y)) =,, L ρ = and Cov( y, y ) s ofen denoed by whle Var( y ) f ofen denoed by γ. Noe ha γ = γ and ρ = ρ and because
More informationOnline Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading
Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng
More information( ) () we define the interaction representation by the unitary transformation () = ()
Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger
More informationWiH Wei He
Sysem Idenfcaon of onlnear Sae-Space Space Baery odels WH We He wehe@calce.umd.edu Advsor: Dr. Chaochao Chen Deparmen of echancal Engneerng Unversy of aryland, College Par 1 Unversy of aryland Bacground
More informationShould Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth
Should Exac Index umbers have Sandard Errors? Theory and Applcaon o Asan Growh Rober C. Feensra Marshall B. Rensdorf ovember 003 Proof of Proposon APPEDIX () Frs, we wll derve he convenonal Sao-Vara prce
More informationKernel-Based Bayesian Filtering for Object Tracking
Kernel-Based Bayesan Flerng for Objec Trackng Bohyung Han Yng Zhu Dorn Comancu Larry Davs Dep. of Compuer Scence Real-Tme Vson and Modelng Inegraed Daa and Sysems Unversy of Maryland Semens Corporae Research
More informationNPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management
P age NPTEL Proec Economerc Modellng Vnod Gua School of Managemen Module23: Granger Causaly Tes Lecure35: Granger Causaly Tes Rudra P. Pradhan Vnod Gua School of Managemen Indan Insue of Technology Kharagur,
More informationOrdinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s
Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class
More informationPolymerization Technology Laboratory Course
Prakkum Polymer Scence/Polymersaonsechnk Versuch Resdence Tme Dsrbuon Polymerzaon Technology Laboraory Course Resdence Tme Dsrbuon of Chemcal Reacors If molecules or elemens of a flud are akng dfferen
More informationMethods for the estimation of missing values in time series
Edh Cowan Unversy Research Onlne Theses: Docoraes and Masers Theses 6 Mehods for he esmaon of mssng values n me seres Davd S. Fung Edh Cowan Unversy Recommended Caon Fung, D. S. (6). Mehods for he esmaon
More informationMODELING TIME-VARYING TRADING-DAY EFFECTS IN MONTHLY TIME SERIES
MODELING TIME-VARYING TRADING-DAY EFFECTS IN MONTHLY TIME SERIES Wllam R. Bell, Census Bureau and Donald E. K. Marn, Howard Unversy and Census Bureau Donald E. K. Marn, Howard Unversy, Washngon DC 0059
More informationOutline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model
Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon
More information. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.
Lnear Algebra Lecure # Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
Ths documen s downloaded from DR-NTU, Nanyang Technologcal Unversy Lbrary, Sngapore. Tle A smplfed verb machng algorhm for word paron n vsual speech processng( Acceped verson ) Auhor(s) Foo, Say We; Yong,
More informationBayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance
INF 43 3.. Repeon Anne Solberg (anne@f.uo.no Bayes rule for a classfcaon problem Suppose we have J, =,...J classes. s he class label for a pxel, and x s he observed feaure vecor. We can use Bayes rule
More informationRelative controllability of nonlinear systems with delays in control
Relave conrollably o nonlnear sysems wh delays n conrol Jerzy Klamka Insue o Conrol Engneerng, Slesan Techncal Unversy, 44- Glwce, Poland. phone/ax : 48 32 37227, {jklamka}@a.polsl.glwce.pl Keywor: Conrollably.
More information( ) [ ] MAP Decision Rule
Announcemens Bayes Decson Theory wh Normal Dsrbuons HW0 due oday HW o be assgned soon Proec descrpon posed Bomercs CSE 90 Lecure 4 CSE90, Sprng 04 CSE90, Sprng 04 Key Probables 4 ω class label X feaure
More informationLi An-Ping. Beijing , P.R.China
A New Type of Cpher: DICING_csb L An-Png Bejng 100085, P.R.Chna apl0001@sna.com Absrac: In hs paper, we wll propose a new ype of cpher named DICING_csb, whch s derved from our prevous sream cpher DICING.
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu
More informationHigh frequency analysis of lead-lag relationships between financial markets de Jong, Frank; Nijman, Theo
Tlburg Unversy Hgh frequency analyss of lead-lag relaonshps beween fnancal markes de Jong, Frank; Nman, Theo Publcaon dae: 1995 Lnk o publcaon Caon for publshed verson (APA): de Jong, F. C. J. M., & Nman,
More informationA HIERARCHICAL KALMAN FILTER
A HIERARCHICAL KALMAN FILER Greg aylor aylor Fry Consulng Acuares Level 8, 3 Clarence Sree Sydney NSW Ausrala Professoral Assocae, Cenre for Acuaral Sudes Faculy of Economcs and Commerce Unversy of Melbourne
More information5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)
5h Inernaonal onference on Advanced Desgn and Manufacurng Engneerng (IADME 5 The Falure Rae Expermenal Sudy of Specal N Machne Tool hunshan He, a, *, La Pan,b and Bng Hu 3,c,,3 ollege of Mechancal and
More information[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5
TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as
More informationGraduate Macroeconomics 2 Problem set 5. - Solutions
Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K
More informationInfluence Diagnostics in a Bivariate GARCH Process
Influence Dagnoscs n a Bvarae GARCH Process an Qu Jonaan Dark Xbn Zang Deparmen of Economercs and Busness Sascs Monas Unversy Caulfeld Eas VIC 345 Ausrala Marc 6 Absrac: In s paper we examne nfluence dagnoscs
More informationIntroduction to Boosting
Inroducon o Boosng Cynha Rudn PACM, Prnceon Unversy Advsors Ingrd Daubeches and Rober Schapre Say you have a daabase of news arcles, +, +, -, -, +, +, -, -, +, +, -, -, +, +, -, + where arcles are labeled
More informationReactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times
Reacve Mehods o Solve he Berh AllocaonProblem wh Sochasc Arrval and Handlng Tmes Nsh Umang* Mchel Berlare* * TRANSP-OR, Ecole Polyechnque Fédérale de Lausanne Frs Workshop on Large Scale Opmzaon November
More informationSampling Procedure of the Sum of two Binary Markov Process Realizations
Samplng Procedure of he Sum of wo Bnary Markov Process Realzaons YURY GORITSKIY Dep. of Mahemacal Modelng of Moscow Power Insue (Techncal Unversy), Moscow, RUSSIA, E-mal: gorsky@yandex.ru VLADIMIR KAZAKOV
More informationAdaptive Multivariate Statistical Process Control for Monitoring Time-varying Processes 1
Adapve Mulvarae Sascal Process Conrol for Monorng me-varyng Processes Sang Wook Cho *, Elane. Marn *, A. Julan Morrs *, and In-eum Lee + * Cenre of Process Analycs and Conrol echnology, School of Chemcal
More informationTSS = SST + SSE An orthogonal partition of the total SS
ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally
More information2.1 Constitutive Theory
Secon.. Consuve Theory.. Consuve Equaons Governng Equaons The equaons governng he behavour of maerals are (n he spaal form) dρ v & ρ + ρdv v = + ρ = Conservaon of Mass (..a) d x σ j dv dvσ + b = ρ v& +
More information. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.
Mah E-b Lecure #0 Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons are
More informationDYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń 2008
DYNAMIC ECONOMETRIC MODELS Vol. 8 Ncolaus Coperncus Unversy Toruń 008 Monka Kośko The Unversy of Compuer Scence and Economcs n Olszyn Mchał Perzak Ncolaus Coperncus Unversy Modelng Fnancal Tme Seres Volaly
More informationMath 128b Project. Jude Yuen
Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally
More informatione-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov
June 7 e-ournal Relably: Theory& Applcaons No (Vol. CONFIDENCE INTERVALS ASSOCIATED WITH PERFORMANCE ANALYSIS OF SYMMETRIC LARGE CLOSED CLIENT/SERVER COMPUTER NETWORKS Absrac Vyacheslav Abramov School
More informationLet s treat the problem of the response of a system to an applied external force. Again,
Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem
More informationNormal Random Variable and its discriminant functions
Noral Rando Varable and s dscrnan funcons Oulne Noral Rando Varable Properes Dscrnan funcons Why Noral Rando Varables? Analycally racable Works well when observaon coes for a corruped snle prooype 3 The
More informationFall 2009 Social Sciences 7418 University of Wisconsin-Madison. Problem Set 2 Answers (4) (6) di = D (10)
Publc Affars 974 Menze D. Chnn Fall 2009 Socal Scences 7418 Unversy of Wsconsn-Madson Problem Se 2 Answers Due n lecure on Thursday, November 12. " Box n" your answers o he algebrac quesons. 1. Consder
More informationThe Impact of SGX MSCI Taiwan Index Futures on the Volatility. of the Taiwan Stock Market: An EGARCH Approach
The Impac of SGX MSCI Tawan Index Fuures on he Volaly of he Tawan Sock Marke: An EGARCH Approach Phlp Hsu, Asssan Professor, Deparmen of Fnance, Naonal Formosa Unversy, Tawan Yu-Mn Chang, Asssan Professor,
More informationOMXS30 Balance 20% Index Rules
OMX30 Balance 0% ndex Rules Verson as of 30 March 009 Copyrgh 008, The NADAQ OMX Group, nc. All rghs reserved. NADAQ OMX, The NADAQ ock Marke and NADAQ are regsered servce/rademarks of The NADAQ OMX Group,
More informationSignal Diffusion Mapping: Optimal Forecasting with Time Varying Lags.
Absrac Sgnal Dffuson Mappng: Opmal Forecasng wh Tme Varyng Lags. Paul Gaskell Unversy of Souhampon, WAIS Frank McGroary Unversy of Souhampon, Fnance and Bankng Group Thanasss Tropans Unversy of Souhampon,
More informationEndogeneity. Is the term given to the situation when one or more of the regressors in the model are correlated with the error term such that
s row Endogeney Is he erm gven o he suaon when one or more of he regressors n he model are correlaed wh he error erm such ha E( u 0 The 3 man causes of endogeney are: Measuremen error n he rgh hand sde
More informationTools for Analysis of Accelerated Life and Degradation Test Data
Acceleraed Sress Tesng and Relably Tools for Analyss of Acceleraed Lfe and Degradaon Tes Daa Presened by: Reuel Smh Unversy of Maryland College Park smhrc@umd.edu Sepember-5-6 Sepember 28-30 206, Pensacola
More informationComputing Relevance, Similarity: The Vector Space Model
Compung Relevance, Smlary: The Vecor Space Model Based on Larson and Hears s sldes a UC-Bereley hp://.sms.bereley.edu/courses/s0/f00/ aabase Managemen Sysems, R. Ramarshnan ocumen Vecors v ocumens are
More informationFactor models with many assets: strong factors, weak factors, and the two-pass procedure
Facor models wh many asses: srong facors, weak facors, and he wo-pass procedure Sanslav Anaolyev CERGE-EI and ES Anna Mkusheva MI Augus 07 PRELIMIARY AD ICOMPLEE. PLEASE, DO O DISRIBUE! Absrac hs paper
More informationExistence and Uniqueness Results for Random Impulsive Integro-Differential Equation
Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal
More informationDEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL
DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL Sco Wsdom, John Hershey 2, Jonahan Le Roux 2, and Shnj Waanabe 2 Deparmen o Elecrcal Engneerng, Unversy o Washngon, Seale, WA, USA
More informationAdvanced Machine Learning & Perception
Advanced Machne Learnng & Percepon Insrucor: Tony Jebara SVM Feaure & Kernel Selecon SVM Eensons Feaure Selecon (Flerng and Wrappng) SVM Feaure Selecon SVM Kernel Selecon SVM Eensons Classfcaon Feaure/Kernel
More informationHEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD
Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,
More informationJ i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.
umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal
More informationMachine Learning 2nd Edition
INTRODUCTION TO Lecure Sldes for Machne Learnng nd Edon ETHEM ALPAYDIN, modfed by Leonardo Bobadlla and some pars from hp://www.cs.au.ac.l/~aparzn/machnelearnng/ The MIT Press, 00 alpaydn@boun.edu.r hp://www.cmpe.boun.edu.r/~ehem/mle
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ
More information