and Sun (14) and Due and Singlen (19) apply he maximum likelihd mehd while Singh (15), and Lngsa and Schwarz (12) respecively emply he hreesage leas s

Size: px
Start display at page:

Download "and Sun (14) and Due and Singlen (19) apply he maximum likelihd mehd while Singh (15), and Lngsa and Schwarz (12) respecively emply he hreesage leas s"

Transcription

1 A MONTE CARLO FILTERING APPROACH FOR ESTIMATING THE TERM STRUCTURE OF INTEREST RATES Akihik Takahashi 1 and Seish Sa 2 1 The Universiy f Tky, Kmaba, Megur-ku, Tky Japan 2 The Insiue f Saisical Mahemaics, Minami-Azabu, Mina-ku, Tky Japan (March, 2000) Absrac. We develp new mehdlgy fr esimain f he erm srucure f ineres raes based n a Mne Carl lering apprach. As an example, we apply he mehd LIBORs (Lndn InerBank Oered Raes) and ineres raes swaps in he Japanese inerbank marke. Key wrds and phrases: Generalized Sae Sapce Mdel, Mne Carl Filer, Term Srucure f Ineres Rae, Mne Carl Simulain, Term Srucure Mdel 1. Inrducin We prpse a new framewrk f he esimain f he erm srucure f ineres raes based n he sae space mdel. In paricular, ur knwledge, i is he rs applicain f he Mne Carl ler he esimain f he erm srucure mdels described by sae variables fllwing muli-dimensinal Markv prcesses. Fr example, ur mehd can be applied he erm srucure mdels based n he dynamic general equilibrium hery f Cx, Ingersll and Rss (15a,b)(CIR) which includes muli-facr CIR mdels used by Chen and Sc (13), Pearsn and Sun (14), Singh (15) and Due and Singlen (19), and he schasic vlailiy mdel develped by Lngsa and Schwarz (12). I is well-knwn ha a leas w sae variables are necessary explain he dynamics f he erm srucure in he real wrld. Chen and Sc (13), Pearsn and Sun (14), Singh (15), and Due and Singlen (19) cncluded ha ne-facr mdels are n enugh describe he variain f he erm srucure by analyzing reasury r swap markes in he Unied Saes. Hwever, muli-facr mdels fen esed in he analysis are n necessarily amng he bes candidaes. Fr insance, in he muli-facr CIR mdel where he sp ineres rae is described by he sum f several sae variables independenly fllwing square-r prcesses, he inuiive inerpreain f he sae variables is n clear, and smeimes he esimaed parameers d n saisfy he nnnegaive cndiin f he sp rae. One f he reasns why he mdels which d n explain he daa very well are fen emplyed in empirical analyses is ha hey allw analyic sluins f zer cupn bnds' prices. I is mainly due he limiain f he mehds applied he esimain. In paricular, Chen and Sc (13), Pearsn 1

2 and Sun (14) and Due and Singlen (19) apply he maximum likelihd mehd while Singh (15), and Lngsa and Schwarz (12) respecively emply he hreesage leas square mehd wih he principal cmpnen analysis and he generalized mmen mehd(gmm) wih GARCH. Hwever, i is subsanially dicul apply hse mehds wihu analyic sluins f he zer-cupn bnds' prices. Mrever, exising researches end replace unbservable sae variables such as he sp rae and he vlailiy by sme bservable variables (Chan e al. (12), Lngsa and Schwarz (12)), bu hey may subsanially suer frm he measuremen errrs. While sme f hem ake he measuremen errrs in accun explicily, he ways f he cnsiderain are n naural and smewha ad hc. (Chen and Sc (13), Due and Singlen (19)) We prpse a Mne Carl lering apprach based n he generalized sae space mdel vercme he prblems f exising researches. The sae space mdel cnsiss f he sysem mdel describing he prcesses f sae variables and he bservain mdel represening he funcinal relain beween he sae variables and he bservainal daa in he real wrld, which implies ha he mehd can be naurally applied he esimain f he erm srucure mdels based n he Markv sae prcesses. I is als pssible inrduce measuremen errrs in he general way wihu any bias. Mrever, he Mne Carl ler can be applied much brader class f he erm srucure mdels, especially even he mdels in which he zer-cupn bnds' prices can n be analyically bained. The paper is rganized as fllws. In he nex chaper, we will rs summarize he sae space mdel and he erm srucure mdels based n Markv sae prcesses. Then, we will clarify he relain beween ineres rae mdels as well as bservainal daa and he sae space mdel. Nex, we will give a cncree algrihm f he Mne Carl ler applied he empirical analysis. In chaper hree, we will es he validiy f he Mne Carl ler using LIBORs (Lndn InerBank Oered Raes) and ineres rae swaps in he Japanese inerbank marke. In chaper fur, we will make cncluding remarks. 2. The Esimain f he Term Srucure Based n he Sae Space Mdeling We explain in his chaper he esimain mehd fr he sae variables and he parameers f erm srucure mdels based n he sae space mdeling. Firs, we give he general frm f sae space mdels. (See Kiagawa and Gersh(16) fr he deail.) A sae space mdel cnsiss f he fllwing sysem mdel and he bservain mdel. Tha is, 8 >< Y = F (Y, ;v ) sysem mdel (2.1) >: Z = H(Y ;u ) bservain mdel 2

3 where Y, Z and dene a N dimensinal sae vecr, a M dimensinal bservain vecr a ime and he ime inerval f bservainal daa respecively while v and u dene he sysem nise and he bservainal nise whse densiy funcins are given respecively by q(v) and (u). F and H are in general nn-linear funcins f R N R N 7! R N and R N R M 7! R M, and he iniial sae vecr Y 0 is assumed be a randm variable whse densiy funcin is given by p 0 (Y ). Befre describing hw sae space mdels are uilized in ur analysis, we briey explain erm srucure mdels based n Markv sae prcesses. We rs assume ha he ime hrizn f ecnmic aciviies is [0;T ], T < 1. Given he lered prbabiliy space (;F;fF g;p), we suppse ha a k dimensinal vecr f sae variables dened by Y fllws a k dimensinal Markv prcess, (2.2) dy = (Y; )d + S(Y; )db ; where B is he n dimensinal sandard Brwnian min under he lered prbabiliy space, and (Y; ) and S(Y; ) dene real-valued funcins f R k [0;T ] 7! R k, and R k [0;T ] 7! R kn, respecively. Suppse als ha he insananeus shr-erm ineres rae a ime dened by r() and a zer cupn bnd's price a wih he mauriy T dened by P (; T ) are sme funcins f Y where 2 [0;T ] and T 2 [; T ]. We ne ha he se f he zer cupn bnds' prices fp (; T )g T 2[0;T ] represens he erm srucure f ineres raes a ime. Then, based n he arbirage-free argumen f nancial ecnmics, P (; T ) saises a parial dierenial equain(pde) 1 2 race(ss0 P YY )+[, ]P Y + P, rp =0; wih he erminal bundary cndiin, P (T;T)=1. Here, (Y; ) denes s called he risk premium which is a funcin f Y and, and fr insance, i can be deermined by he general equilibrium asse pricing hery f ecnmics such as CIR(15a). Mrever, i is well knwn ha he sluin f his PDE is represened by he cndiinal expecain given infrmain a ime, (2.3) R P (; T )=E Q [e, T r(u)du ]; where E Q [] denes he cndiinal expecain perar a ime under s called he risk-neural prbabiliy measure Q. I is als knwn ha under he measure Q, he vecr f sae variables Y fllws a schasic dierenial equain, (2.4) dy = f(y; ), (Y; )gd + S(Y; )db where B denes he n dimensinal sandard Brwnian min under he measure Q. Nex, we clarify hw he sae space mdels can be applied he esimain f erm srucure mdels. When is sucienly small, he Euler apprximain he equain 3

4 (2:2) can be used fr he sysem mdel, Y = F (Y, ;v ). Tha is, (2.5) Y = Y, + (Y, ;, ) + S(Y, ;, )v p where he sysem nise v fllws he N dimensinal sandard nrmal disribuin. Of curse, he her apprximain schemes culd be applied he descreizain f (2:2). (Fr example, see secin D and nes f chaper 11 in Due(16).) Mrever, when Y is explicily slved given Y, as in he case f a linear schasic dierenial equain, i is beer use ha represenain. Tha is, in (2:2), suppse Y is represened by a linear schasic dierenial equain: (2.6) dy =(AY + )d + SdB where and A dene R k -valued funcins f he ime parameer and k k cnsan marices respecively, and S denes a k n cnsan marix s ha = SS 0 is psiive semidenie. Then, given Y,, Y is slved by (2.7) Y = e A Y, + Z, e (,s)a ds + v : Here, v fllws he nrmal disribuin wih he mean zer and he cvariance marix where Z 0 e sa e sa0 ds: In his case, he sysem mdel is given by (2:7) and he densiy funcin q(v) f he sysem nise v is ha f he nrmal disribuin. In he bservain mdel, he bservain vecr a ime, Z is expressed by n( 1) unis f zer-cupn bnds' prices and he nise vecr u. Tha is, Z = H(fP (; + T i )g n i=1 ;u ): In he equain abve, P (; + T i ) can be evaluaed by he cmpuain f he equain (2:3) under he prcess (2:4). In addiin, we assume hereafer ha he bservain mdel is f he frm, (2.8) Z = H(fP (; + T i )g n i=1)+u ; u N(0; u ); where he densiy funcin (u) f he nise vecr u is given by ha f he mulidimensinal nrmal disribuin wih he mean 0, and he variance-cvariance marix u. Here, u denes a M M diagnal marix wih psiive diagnal elemens where M is dened be he dimensin f he bservainal daa a each ime. LIBORs and ineres rae swap raes a ime in he inerbank marke are he ypical examples f he bservain vecr Z. In his case, H() is deermined by he hereical relain 4

5 beween LIBORs/ineres rae swap raes and he zer-cupn bnds' prices. Tha is, LIBOR wih he erm T n dened by L (T n ) and swap raes wih he erm T n dened by S (T n ) are expressed by he zer-cupn bnds' prices as fllws. (2.9) (2.10)! 1 1 L (T n ) = P (; + T n ), 1 T n S (T n ) = 1 P, P (; + T n) n i=1 P (; + i); where denes he inerval f cash ws and T n = n. Fr example, =0:5 is sandard in he Japanese yen swap marke. We will explain he acual daa f LIBORs and swap raes used fr he empirical analysis in he nex chaper. Frm he discussin abve, I can be seen ha he framewrk f he sae space mdel is naurally applied he esimain prblem f he erm srucure f ineres raes. Finally, we disucuss abu ur esimain mehd in deail. We ne ha he sandard Kalman ler can n be applied he esimain as bh he sysem mdel and he bservain mdel described abve are generally nn-linear, and hence we uilize he Mne Carl ler. While several appraches are prpsed fr he Mne Carl ler (see Duce, Bara, and Duvau(15), Durbin, and Kpman(19), Grdn, Salmnd, and Smih(13), Tanizaki(13), fr insance.), we adp he apprach develped by Kiagawa(16). In he fllwing, we describe he uline f he algrihm f he Mne Carl ler applied he empirical analysis in he nex chaper. Firs, we summarize he nain fllwing Kiagawa(16). p(y jz, ), called \ne sep ahead predicin" denes he cndiinal densiy funcin f Y given Z, where is he inerval f ime series. p(y jz ), called \ler" denes he cndiinal densiy funcin f Y given Z. fp (1) ; ;p (m) g and ff (1) ; ;f (m) g represen he vecrs f he realizain f m rials f Mne Carl frm p(y jz, ) and p(y jz ), respecively. Then, if we se ff (1) 0 ; ;f (m) 0 g as he realizain f Mne Carl frm p 0 (Y ), he densiy funcin f he iniial sae vecr Y 0, he algrihm f he Mne Carl ler is as fllws. [The summary f he algrihm f he Mne Carl ler] (i) Generae he iniial sae vecr ff (1) 0 ; ;f (m) 0 g. (ii) Apply he fllwing seps (a)(d) each ime = 0; ; 2; ; (T, );T where T denes he nal ime pin fhedaa. (a) Generae he sysem nise v (j), j =1; ;m accrding he densiy funcin q(v). (b) Cmpue fr each j =1; ;m p (j) = F (f (j), ;v(j) ): 5

6 (c) Evaluae n[x;0; u ], he densiy funcin f N(0; u )ax = Z,H(p (j) ),j = 1; ;m and dene hse as (j), j = 1; ;m. Here, H() represen fr insance, he equain (2:9) and (2:10) which are regarded as funcins f he sae vecr Y. The prices f zer-cupn bnds in hse equains are cmpued hrugh he equain (2:3) by using he prcess (2:4), and if i is n analyically slvable, sme numerical mehd such as Mne Carl simulain is implemened. g. Mre pre- g wih (d) Implemen resampling f ff (1) ; ;f (m) cisely, implemen resampling f each f (i) he prbabiliy Prb.(f (i) = p (j) jz )= (j) P m k=1 (k) g frm fp (1) ; ;p (m), i =1; ;mfrm fp (1) ; ;p (m) ; j =1; ;m; i =1; ;m: The esimain f unknwn parameers is based n he maximum likelihd mehd. If denes he vecr represening whle unknwn parameers, he lg-likelihd l() is given by l() =p(z ; ;Z T j) = T k=1 p(z kjz ; ;Z (k,1) ;) where p(z jz 0 ) = p 0 (Z ). The lg-likelihd l() is cmpued apprximaely wihin he framewrk f he Mne Carl ler by l() = T X k=1 lg m X j=1 (j) k 1 T A, lg m: Then, maximize l() wih respec bain he maximum likelihd esimar ^. Fr pimizain, grid search and a self-rganizing mehd are applied. (See Kiagawa(19) fr deails f a self-rganizing sae-space mdel.) Finally, we uilize AIC(Akaike's Infrmain Crierin) as a crierin selec he erm srucure mdels if here are several candidaes. mdel. Tha is, he mdel wih he smaller AIC can be regarded as he beer 3. An Applicain he Time Series f he Japanese Ineres Raes In his chaper, we examine he validiy f ur mehd using he ime series f ineres raes in he Japaneses inerbank marke. As an example f ineres rae mdels, we use Hull and Whie(14) in which he dynamics f a sae vecr Y can be represened by a linear schasic dierenial equain (2:6), and he resuling sysem mdel is described by he equain (2:7); given, he sysem mdel can be expressed as a linear mdel in he sae space seing: (3.1) Y = FY, + ()+v ; 6

7 where Y = (Y i ());i = 1; 2 denes a w dimensinal sae vecr, () is a R 2 -valued funcin f he ime parameer, F is a 2 2 marix wih cnsan elemens, and v fllws a w dimensinal nrmal disribuin wih he mean zer. They als assume ha he sp rae r is expressed as a funcin f Y 1, r = g(y 1 ) where g() is sme realvalued funcin. Fr a funcinal frm f r = g(y 1 ), if we ake g(y 1 ) = Y 1, we allw negaive ineres raes because f he nrmaliy f he sp rae while we can bain an analyic sluin f P (; T ), which subsanially reduce he cmpuainal burden. In paricular, he mdel implies relaively high prbabiliy f he negaive ineres raes in such lw ineres raes envirnmen f he recen Japan and his seems inapprpriae. Hence, we specify g(y 1 ) based n Hull(19) (chaper 21, pp.588) such ha g(y 1 ) = Y 1 fr Y 1 ", and g(y 1 ) = "e (Y 1,") " fr Y 1 < ", where " is sme predeermined psiive cnsan. Clearly, g() is psiive, mnnically increasing, and lim Y1!,1 g(y 1 ) = 0. We nex deermine he bservain mdel as he equain (2:8) wih he equains (2:9) and (2:10). Mrever, he zer-cupn bnds' prices in H() f (2:8) are cmpued by he equain, (3.2) P (; T )=E Q [e,r T g(y 1 (u))du jy ()]: We ne ha his is n analyically slvable and hence is evaluaed by sme numerical mehd such as Mne Carl simulains. Nex, we explain he resul f esimains fr he Japanese inerbank marke. The daa used fr he analysis is summarized as fllws. he perid and he frequency f he daa: daily daa f 19/1/1-19/7/22(662 series) Japanese yen LIBOR: six-mnh, welve-mnh Japanese yen swap raes: w-year, hree-year, fur-year, ve-year, seven-year, enyear Figures 1(1)(3) shw he bservainal daa fr he perid f he analysis; (1) and (2) respecively shw he series f LIBORs and hse f swap raes while (3) shws he spread beween en-year and w-year. We apply he algrihm f he Mne Carl ler described in he previus chaper he esimain. The dimensin f he bservainal daa is M =8,he equains (2:9) and (2:10) are applied H() in he equain (2:8), and zer-cupn bnds' prices in he equain (3:2) are cmpued by using Mne Carl simulains. The numbers f rials in he Mne Carl ler and in he Mne Carl simulains are 5000 and 300 respecively where he randm numbers are generaed by ran2 in Numerical Recipes (Secnd Ediin), and he aniheic variables mehd is uilized in he Mne Carl simulains. Mrever, we implemen parallel cmpuain 7

8 wih 40 prcessrs by IBM RS/6000-SP. In addiin, we suppse " be 0:0005(5 basis pin) fr r = g(y 1 ). Firs we invesigae he ime series f he esimaed facrs. Figures 2(1)(2) shw he relain beween he esimaed facrs Y i, i = 1; 2 bained by averaging m = 5000 samples f lers and he frward raes cmpued frm he daa. The crrelain beween he level f each facr and he crrespnding frward rae is als lised belw each gure. Frm he graphical bservain, here seems be srng relain beween Y 2 and enyear frward rae which is dened be he raes wih he erm 0.5 year saring frm 9.5 years frward, and he variain f Y 1 is similar ha f six-mnh LIBOR. We nex see he ing f he mdel he daa. The ime series f each bservain and crrespnding esimae are shwn in gures 3(1)(8), and he explanain pwer f each esimae is lised abve crrespnding gure. We ne ha he esimaes f LIBORs and swap raes are cmpued by Mne Carl simulains based n esimaed parameers and esimaed facrs, where he number f rials f he simulains are They shw ha he mdel s very well he swap raes while i des n LIBORs: The explanain pwers are mre han.5% fr all he swap raes while hey are 39.4% and 77.3% fr six-mnh and welve-mnh LIBORs respecively. Finally, gures 4(1)(4) shw he bservains and esimaes f erm srucures a fur daes (/6/24, /12/4, /9/30, /5/14), and hese indicae ha he shapes f he erm srucure f he mdel are very clse he real nes excep fr he LIBORs a ember 4h f 19 as well as a 14h f 19. Hence, we can cnclude ha he mdel explains he variains f he swap raes very well while i des n explain hse f LIBORs. This implies ha anher facr may be necessary imprve he ing LIBORs. Hence, we nex implemen he case in which he sae vecr Y = (Y i ()), i =1; 2; 3 is hree dimensinal, and has he same frm f (3:1). As in he previus case, we rs invesigae he ime series f he esimaed facrs. Figures 5(1)(3) shw he relains beween he esimaed facrs Y i, i =1; 2; 3 and he frward raes as well as heir crrelains. In his case, here is much clearer relain beween facrs and frward raes han in he previus case. Tha is, here seems exis srng relains beween Y 1 and six-mnh LIBOR, and Y 2 and he w--en-year frward raes' spread as well as beween Y 3 and en-year frward rae. The ime series f each bservain and crrespnding esimae as well as he explanain pwers fr his case are shwn in gures 6(1)(8). Clearly, he explanain pwers are mre han.5% fr all he raes and mre han % excep fr welve-mnh LIBOR and w-year swap rae. Especially, we ne ha he mdel remarkably imprves he ing LIBORs. Figures 7(1)(4) shw he bservains and esimaes f he erm srucures a fur daes(/6/24, /12/4, /9/30, /5/14), which implies ha he mdel can replicaes he real erm srucures including LIBORs very clsely. 8

9 Finally, we ne ha AIC in his case(aic= ) is mre imprved han in he previus case(aic= ), and hus, we can cnclude ha he mdel explains he variains f all he raes fairly well. Frm he abve analysis f he Japanese inerbank marke, we cnrm ha ur apprach uilizing he Mne Carl ler is valid fr he esimain f muli-facr erm srucure mdels which is usually cnsidered ugh ask. 4. Cncluding Remarks We develp a new framewrk fr he empirical analysis f he erm srucure f ineres raes based n he generalized sae space mdel. As an example, we apply he Mne Carl ler he ime series f LIBORs and f ineres rae swaps in he Japanese inerbank marke, and cnrm he validiy f ur mehd. References Chan, K.C., G.A. Karlyi, F.A. Lngsa, and A.B. Sanders(12), \An Empirical Cmparisn f Alernaive Mdels f he Shr-Term Ineres Rae," The Jurnal f Finance 47, pp Chen, R., and L. Sc(13), \Maximum Likelihd Esimain fr a Mulifacr Equilibrium Mdel f he Term Srucure f Ineres Raes," Jurnal f Fixed Incme, ember, pp Cx, J.C., J.E. Ingersll, and S.A. Rss(15a), \An Inerempral General Equilibrium Mdel Asse Prices," Ecnmerica 53, pp Cx, J.C., J.E. Ingersll, and S.A. Rss(15b), \A Thery f he Term Srucure f Ineres Raes," Ecnmerica 53, pp Duce, A., E. Bara, and P. Duvau(15), \A Mne Carl Apprach Recursive Bayesian Sae Esimain," Prceedings IEEE Signal Prcessing/Ahs Wrkshp n Higher Order Saisics, Girna, Spain. Due, D.(16), Dynamic Asse Pricing Thery(secnd ediin), Princen Universiy Press. Due, D., and K.J. Singlen(19), \An Ecnmeric Mdel f he Term Srucure f Ineres-Rae Swap Yields," The Jurnal f Finance, 52, ember pp Durbin, J., and S.J. Kpman(19), \Mne Carl Maximum Likelihd Esimain fr Nn-Gaussian Sae Space Mdels," Bimerika 84, pp Grdn, N., D.J. Salmnd, and A.F.M. Smih(13), \Nvel Apprach Nnlinear/Nn-Gaussian Bayesian Sae Esimain," IEE Prceedings-F, 140, N.2, pp Hull, J.(19), Opins, Fuures, and Oher Derivaives (furh ediin), Prenice-Hall. Hull, J., and A. Whie(14), \Numerical Prcedures fr Implemening Term Srucure Mdels II:Tw- Facr Mdels," Jurnal f Derivaives 2, Winer pp Kiagawa, G.(16), \Mne Carl Filer and Smher fr Nn-Gaussian Nnlinear Sae Space Mdels," Jurnal f Cmpuainal and Graphical Saisics 5, N.1, pp Kiagawa, G.(19), \A Self-Organizing Sae-Space Mdel," Jurnal f he American Saisical Assciain, 93, N.443, pp Kiagawa, G. and W. Gersch(16), Smhness Prir Analysis f Time Series, Lecure Nes in Saisics N. 116, Springer-Verlag. Lngsa,F., and E.S. Schwarz(12), \Ineres Rae Vlailiy and he Term Srucure:A Tw-Facr General Equilibrium Mdel," The Jurnal f Finance, 47, pp

10 Pearsn, N.D., and T.S. Sun(14), \Expliing he Cndiinal Densiy in Esimaing he Term Srucure: An Applicain he Cx, Ingersll, and Rss Mdel," The Jurnal f Finance 49, pp Press, W.H., B.P. Flannery, S.A. Teuklsky and W.T. Veerling (12), Numerical Recipes in C (2nd Ediin), Cambridge Universiy Press. Singh, M.K.(15)\Esimain f Mulifacr Cx, Ingersll, and Rss Term Srucure Mdel," Jurnal f Fixed Incme, ember,pp Tanizaki, H.(13), Nnlinear Filers, Lecure Nes in Ecnmics and Mahemaical Sysems, N. 400, Springer-Verlag. 10

11 (1) LIBORs (3) Spread (Swap10Y - Swap2Y) LIBOR-1Y LIBOR-6M (BP) Aug Nv Jul Mar Aug Nv Jul (2) SWAP raes SWAP-10Y SWAP-7Y SWAP-5Y SWAP-4Y SWAP-3Y SWAP-2Y Aug Nv Fig. 1. Jul Mar ime series f bserved LIBORs and swap raes (1) Libr6M(slid), Y1(ds) (Scale:lef) (Scale:righ) (2) Frward10Y(slid), Y2(ds) (Scale:lef) (Scale:righ) Crr= 0.90 Crr= 0. (3) Frward 2Y-10Y(slid), Y1(ds) (Scale:lef) (Scale:righ) Crr= 0.67 Fig. 2. esimaed facrs (w-facr case) 11

12 (1) Libr6M (2) Libr1Y R2= 39.4 % R2= 77.3 % Observains (slid line), Esimaes (ds line) Observains (slid line), Esimaes (ds line) (3) Swap2Y (4) Swap3Y R2=.7 % R2=.5 % Observains (slid line), Esimaes (ds line) Observains (slid line), Esimaes (ds line) (5) Swap4Y (6) Swap5Y R2=.6 % R2=.7 % Observains (slid line), Esimaes (ds line) Observains (slid line), Esimaes (ds line) (7) Swap7Y (8) Swap10Y R2=.8 % R2=.6 % h Observains (slid line), Esimaes (ds line) i Observains (slid line), Esimaes (ds line) Fig. 3. bservains and esimaes (w-facr case); he explanain pwer (R2) is dened by max 1, he variance f residuals ; he variance f bservains 0. 12

13 (1) 24,19 (2). 4, (Yr) Observain(), Esimaes() (Yr) Observain(), Esimaes() (3). 30, 19 (4) 14, (Yr) (Yr) Observain(), Esimaes() Observain(), Esimaes() Fig. 4. he erm srucures(w-facr case) 13

14 (1) Libr6M(slid), Y1(ds) (Scale:lef) (Scale:righ) (2) Frward10Y(slid), Y3(ds) (Scale:lef) (Scale:righ) Crr= 0.90 Crr= 0. (3) Frward 2Y-10Y(slid), Y2(ds) (Scale:lef) (Scale:righ) Crr= 0.90 Fig. 5. esimaed facrs (hree-facr case) 14

15 (1) Libr6M (2) Libr1Y R2=.6 % R2=.6 % Observains (slid line), Esimaes (ds line) Observains (slid line), Esimaes (ds line) (3) Swap2Y (4) Swap3Y R2=.6 % R2=.8 % Observains (slid line), Esimaes (ds line) Observains (slid line), Esimaes (ds line) (5) Swap4Y (6) Swap5Y R2= % R2=.9 % Observains (slid line), Esimaes (ds line) Observains (slid line), Esimaes (ds line) (7) Swap7Y (8) Swap10Y R2=.9 % R2=.6 % Observains (slid line), Esimaes (ds line) Observains (slid line), Esimaes (ds line) Fig. 6. bservains and esimaes (hree-facr case) 15

16 (1) 24,19 (2). 4, (Yr) Observain(), Esimaes() (Yr) Observain(), Esimaes() (3). 30, 19 (4) 14, (Yr) (Yr) Observain(), Esimaes() Observain(), Esimaes() Fig. 7. he erm srucures (hree-facr case) 16

independenly fllwing square-r prcesses, he inuiive inerpreain f he sae variables is n clear, and smeimes i seems dicul nd admissible parameers fr whic

independenly fllwing square-r prcesses, he inuiive inerpreain f he sae variables is n clear, and smeimes i seems dicul nd admissible parameers fr whic A MONTE CARLO FILTERING APPROACH FOR ESTIMATING THE TERM STRUCTURE OF INTEREST RATES Akihik Takahashi 1 and Seish Sa 2 1 The Universiy f Tky, 3-8-1 Kmaba, Megur-ku, Tky 153-8914 Japan 2 The Insiue f Saisical

More information

Brace-Gatarek-Musiela model

Brace-Gatarek-Musiela model Chaper 34 Brace-Gaarek-Musiela mdel 34. Review f HJM under risk-neural IP where f ( T Frward rae a ime fr brrwing a ime T df ( T ( T ( T d + ( T dw ( ( T The ineres rae is r( f (. The bnd prices saisfy

More information

AP Physics 1 MC Practice Kinematics 1D

AP Physics 1 MC Practice Kinematics 1D AP Physics 1 MC Pracice Kinemaics 1D Quesins 1 3 relae w bjecs ha sar a x = 0 a = 0 and mve in ne dimensin independenly f ne anher. Graphs, f he velciy f each bjec versus ime are shwn belw Objec A Objec

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 15 10/30/2013. Ito integral for simple processes

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 15 10/30/2013. Ito integral for simple processes MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.7J Fall 13 Lecure 15 1/3/13 I inegral fr simple prcesses Cnen. 1. Simple prcesses. I ismery. Firs 3 seps in cnsrucing I inegral fr general prcesses 1 I inegral

More information

An application of nonlinear optimization method to. sensitivity analysis of numerical model *

An application of nonlinear optimization method to. sensitivity analysis of numerical model * An applicain f nnlinear pimizain mehd sensiiviy analysis f numerical mdel XU Hui 1, MU Mu 1 and LUO Dehai 2 (1. LASG, Insiue f Amspheric Physics, Chinese Academy f Sciences, Beijing 129, China; 2. Deparmen

More information

The Components of Vector B. The Components of Vector B. Vector Components. Component Method of Vector Addition. Vector Components

The Components of Vector B. The Components of Vector B. Vector Components. Component Method of Vector Addition. Vector Components Upcming eens in PY05 Due ASAP: PY05 prees n WebCT. Submiing i ges yu pin ward yur 5-pin Lecure grade. Please ake i seriusly, bu wha cuns is wheher r n yu submi i, n wheher yu ge hings righ r wrng. Due

More information

10.7 Temperature-dependent Viscoelastic Materials

10.7 Temperature-dependent Viscoelastic Materials Secin.7.7 Temperaure-dependen Viscelasic Maerials Many maerials, fr example plymeric maerials, have a respnse which is srngly emperaure-dependen. Temperaure effecs can be incrpraed in he hery discussed

More information

GMM Estimation of the Number of Latent Factors

GMM Estimation of the Number of Latent Factors GMM Esimain f he Number f aen Facrs Seung C. Ahn a, Marcs F. Perez b March 18, 2007 Absrac We prpse a generalized mehd f mmen (GMM) esimar f he number f laen facrs in linear facr mdels. he mehd is apprpriae

More information

PRINCE SULTAN UNIVERSITY Department of Mathematical Sciences Final Examination Second Semester (072) STAT 271.

PRINCE SULTAN UNIVERSITY Department of Mathematical Sciences Final Examination Second Semester (072) STAT 271. PRINCE SULTAN UNIVERSITY Deparmen f Mahemaical Sciences Final Examinain Secnd Semeser 007 008 (07) STAT 7 Suden Name Suden Number Secin Number Teacher Name Aendance Number Time allwed is ½ hurs. Wrie dwn

More information

Productivity changes of units: A directional measure of cost Malmquist index

Productivity changes of units: A directional measure of cost Malmquist index Available nline a hp://jnrm.srbiau.ac.ir Vl.1, N.2, Summer 2015 Jurnal f New Researches in Mahemaics Science and Research Branch (IAU Prduciviy changes f unis: A direcinal measure f cs Malmquis index G.

More information

5.1 Angles and Their Measure

5.1 Angles and Their Measure 5. Angles and Their Measure Secin 5. Nes Page This secin will cver hw angles are drawn and als arc lengh and rains. We will use (hea) represen an angle s measuremen. In he figure belw i describes hw yu

More information

Visco-elastic Layers

Visco-elastic Layers Visc-elasic Layers Visc-elasic Sluins Sluins are bained by elasic viscelasic crrespndence principle by applying laplace ransfrm remve he ime variable Tw mehds f characerising viscelasic maerials: Mechanical

More information

GAMS Handout 2. Utah State University. Ethan Yang

GAMS Handout 2. Utah State University. Ethan Yang Uah ae Universiy DigialCmmns@UU All ECAIC Maerials ECAIC Repsiry 2017 GAM Handu 2 Ehan Yang yey217@lehigh.edu Fllw his and addiinal wrs a: hps://digialcmmns.usu.edu/ecsaic_all Par f he Civil Engineering

More information

Convex Stochastic Duality and the Biting Lemma

Convex Stochastic Duality and the Biting Lemma Jurnal f Cnvex Analysis Vlume 9 (2002), N. 1, 237 244 Cnvex Schasic Dualiy and he Biing Lemma Igr V. Evsigneev Schl f Ecnmic Sudies, Universiy f Mancheser, Oxfrd Rad, Mancheser, M13 9PL, UK igr.evsigneev@man.ac.uk

More information

An Introduction to Wavelet Analysis. with Applications to Vegetation Monitoring

An Introduction to Wavelet Analysis. with Applications to Vegetation Monitoring An Inrducin Wavele Analysis wih Applicains Vegeain Mniring Dn Percival Applied Physics Labrary, Universiy f Washingn Seale, Washingn, USA verheads fr alk available a hp://saff.washingn.edu/dbp/alks.hml

More information

Lecture 4 ( ) Some points of vertical motion: Here we assumed t 0 =0 and the y axis to be vertical.

Lecture 4 ( ) Some points of vertical motion: Here we assumed t 0 =0 and the y axis to be vertical. Sme pins f erical min: Here we assumed and he y axis be erical. ( ) y g g y y y y g dwnwards 9.8 m/s g Lecure 4 Accelerain The aerage accelerain is defined by he change f elciy wih ime: a ; In analgy,

More information

Kinematics Review Outline

Kinematics Review Outline Kinemaics Review Ouline 1.1.0 Vecrs and Scalars 1.1 One Dimensinal Kinemaics Vecrs have magniude and direcin lacemen; velciy; accelerain sign indicaes direcin + is nrh; eas; up; he righ - is suh; wes;

More information

Motion Along a Straight Line

Motion Along a Straight Line PH 1-3A Fall 010 Min Alng a Sraigh Line Lecure Chaper (Halliday/Resnick/Walker, Fundamenals f Physics 8 h ediin) Min alng a sraigh line Sudies he min f bdies Deals wih frce as he cause f changes in min

More information

Return and Volatility Spillovers Between Large and Small Stocks in the UK

Return and Volatility Spillovers Between Large and Small Stocks in the UK eurn and Vlailiy Spillvers Beween Large and Small Scks in he UK ichard D. F. Harris Xfi Cenre fr Finance and Invesmen Universiy f Exeer, UK Aniru Pisedasalasai Deparmen f Accuning, Finance and Infrmain

More information

Section 12 Time Series Regression with Non- Stationary Variables

Section 12 Time Series Regression with Non- Stationary Variables Secin Time Series Regressin wih Nn- Sainary Variables The TSMR assumpins include, criically, he assumpin ha he variables in a regressin are sainary. Bu many (ms?) ime-series variables are nnsainary. We

More information

THE DETERMINATION OF CRITICAL FLOW FACTORS FOR NATURAL GAS MIXTURES. Part 3: The Calculation of C* for Natural Gas Mixtures

THE DETERMINATION OF CRITICAL FLOW FACTORS FOR NATURAL GAS MIXTURES. Part 3: The Calculation of C* for Natural Gas Mixtures A REPORT ON THE DETERMINATION OF CRITICAL FLOW FACTORS FOR NATURAL GAS MIXTURES Par 3: The Calculain f C* fr Naural Gas Mixures FOR NMSPU Deparmen f Trade and Indusry 151 Buckingham Palace Rad Lndn SW1W

More information

Coherent PSK. The functional model of passband data transmission system is. Signal transmission encoder. x Signal. decoder.

Coherent PSK. The functional model of passband data transmission system is. Signal transmission encoder. x Signal. decoder. Cheren PSK he funcinal mdel f passand daa ransmissin sysem is m i Signal ransmissin encder si s i Signal Mdular Channel Deecr ransmissin decder mˆ Carrier signal m i is a sequence f syml emied frm a message

More information

Stability of the SDDRE based Estimator for Stochastic Nonlinear System

Stability of the SDDRE based Estimator for Stochastic Nonlinear System 26 ISCEE Inernainal Cnference n he Science f Elecrical Engineering Sabiliy f he SDDRE based Esimar fr Schasic Nnlinear Sysem Ilan Rusnak Senir Research Fellw, RAFAEL (63, P.O.Bx 225, 322, Haifa, Israel.;

More information

Numerical solution of some types of fractional optimal control problems

Numerical solution of some types of fractional optimal control problems Numerical Analysis and Scienific mpuing Preprin Seria Numerical sluin f sme ypes f fracinal pimal cnrl prblems N.H. Sweilam T.M. Al-Ajmi R.H.W. Hppe Preprin #23 Deparmen f Mahemaics Universiy f Husn Nvember

More information

Practical Considerations when Estimating in the Presence of Autocorrelation

Practical Considerations when Estimating in the Presence of Autocorrelation CS-BIGS (): -7 008 CS-BIGS hp://www.benley.edu/csbigs/vl-/jaggia.pdf Pracical Cnsiderains when Esimaing in he Presence f Aucrrelain Sanjiv Jaggia Orfalea Cllege f Business, Cal Ply, USA Alisn Kelly-Hawke

More information

Ramsey model. Rationale. Basic setup. A g A exogenous as in Solow. n L each period.

Ramsey model. Rationale. Basic setup. A g A exogenous as in Solow. n L each period. Ramsey mdel Rainale Prblem wih he Slw mdel: ad-hc assumpin f cnsan saving rae Will cnclusins f Slw mdel be alered if saving is endgenusly deermined by uiliy maximizain? Yes, bu we will learn a l abu cnsumpin/saving

More information

i-clicker Question lim Physics 123 Lecture 2 1 Dimensional Motion x 1 x 2 v is not constant in time v = v(t) acceleration lim Review:

i-clicker Question lim Physics 123 Lecture 2 1 Dimensional Motion x 1 x 2 v is not constant in time v = v(t) acceleration lim Review: Reiew: Physics 13 Lecure 1 Dimensinal Min Displacemen: Dx = x - x 1 (If Dx < 0, he displacemen ecr pins he lef.) Aerage elciy: (N he same as aerage speed) a slpe = a x x 1 1 Dx D x 1 x Crrecin: Calculus

More information

ON THE COMPONENT DISTRIBUTION COEFFICIENTS AND SOME REGULARITIES OF THE CRYSTALLIZATION OF SOLID SOLUTION ALLOYS IN MULTICOMPONENT SYSTEMS*

ON THE COMPONENT DISTRIBUTION COEFFICIENTS AND SOME REGULARITIES OF THE CRYSTALLIZATION OF SOLID SOLUTION ALLOYS IN MULTICOMPONENT SYSTEMS* METL 006.-5.5.006, Hradec nad Mravicí ON THE OMPONENT DISTRIUTION OEFFIIENTS ND SOME REGULRITIES OF THE RYSTLLIZTION OF SOLID SOLUTION LLOYS IN MULTIOMPONENT SYSTEMS* Eugenij V.Sidrv a, M.V.Pikunv b, Jarmír.Drápala

More information

American Society for Quality

American Society for Quality American Sciey fr Qualiy Nnparameric Esimain f a Lifeime Disribuin When Censring Times Are Missing Auhr(s): X. Jan Hu, Jerald F. Lawless, Kazuyuki Suzuki Surce: Technmerics, Vl. 4, N. 1 (Feb., 1998), pp.

More information

Machine Learning for Signal Processing Prediction and Estimation, Part II

Machine Learning for Signal Processing Prediction and Estimation, Part II Machine Learning fr Signal Prceing Predicin and Eimain, Par II Bhikha Raj Cla 24. 2 Nv 203 2 Nv 203-755/8797 Adminirivia HW cre u Sme uden wh g really pr mark given chance upgrade Make i all he way he

More information

RAPIDLY ADAPTIVE CFAR DETECTION BY MERGING INDIVIDUAL DECISIONS FROM TWO-STAGE ADAPTIVE DETECTORS

RAPIDLY ADAPTIVE CFAR DETECTION BY MERGING INDIVIDUAL DECISIONS FROM TWO-STAGE ADAPTIVE DETECTORS RAPIDLY ADAPIVE CFAR DEECION BY MERGING INDIVIDUAL DECISIONS FROM WO-SAGE ADAPIVE DEECORS Analii A. Knnv, Sung-yun Chi and Jin-a Kim Research Cener, SX Engine Yngin-si, 694 Krea kaa@ieee.rg; dkrein@nesx.cm;

More information

Microwave Engineering

Microwave Engineering Micrwave Engineering Cheng-Hsing Hsu Deparmen f Elecrical Engineering Nainal Unied Universiy Ouline. Transmissin ine Thery. Transmissin ines and Waveguides eneral Sluins fr TEM, TE, and TM waves ; Parallel

More information

Module 4. Analysis of Statically Indeterminate Structures by the Direct Stiffness Method. Version 2 CE IIT, Kharagpur

Module 4. Analysis of Statically Indeterminate Structures by the Direct Stiffness Method. Version 2 CE IIT, Kharagpur Mdle Analysis f Saically Indeerminae Srcres by he Direc Siffness Mehd Versin CE IIT, Kharagr Lessn The Direc Siffness Mehd: Temerare Changes and Fabricain Errrs in Trss Analysis Versin CE IIT, Kharagr

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility Saisics 441 (Fall 214) November 19, 21, 214 Prof Michael Kozdron Lecure #31, 32: The Ornsein-Uhlenbeck Process as a Model of Volailiy The Ornsein-Uhlenbeck process is a di usion process ha was inroduced

More information

Revelation of Soft-Switching Operation for Isolated DC to Single-phase AC Converter with Power Decoupling

Revelation of Soft-Switching Operation for Isolated DC to Single-phase AC Converter with Power Decoupling Revelain f Sf-Swiching Operain fr Islaed DC Single-phase AC Cnverer wih wer Decupling Nagisa Takaka, Jun-ichi Ih Dep. f Elecrical Engineering Nagaka Universiy f Technlgy Nagaka, Niigaa, Japan nakaka@sn.nagakau.ac.jp,

More information

Lecture 3: Resistive forces, and Energy

Lecture 3: Resistive forces, and Energy Lecure 3: Resisive frces, and Energy Las ie we fund he velciy f a prjecile ving wih air resisance: g g vx ( ) = vx, e vy ( ) = + v + e One re inegrain gives us he psiin as a funcin f ie: dx dy g g = vx,

More information

if N =2 J, obtain analysis (decomposition) of sample variance:

if N =2 J, obtain analysis (decomposition) of sample variance: Wavele Mehds fr Time Series Analysis Eamples: Time Series X Versus Time Inde Par VII: Wavele Variance and Cvariance X (a) (b) eamples f ime series mivae discussin decmpsiin f sample variance using waveles

More information

Fractional Order Disturbance Observer based Robust Control

Fractional Order Disturbance Observer based Robust Control 201 Inernainal Cnference n Indusrial Insrumenain and Cnrl (ICIC) Cllege f Engineering Pune, India. May 28-30, 201 Fracinal Order Disurbance Observer based Rbus Cnrl Bhagyashri Tamhane 1, Amrua Mujumdar

More information

Unit-I (Feedback amplifiers) Features of feedback amplifiers. Presentation by: S.Karthie, Lecturer/ECE SSN College of Engineering

Unit-I (Feedback amplifiers) Features of feedback amplifiers. Presentation by: S.Karthie, Lecturer/ECE SSN College of Engineering Uni-I Feedback ampliiers Feaures eedback ampliiers Presenain by: S.Karhie, Lecurer/ECE SSN Cllege Engineering OBJECTIVES T make he sudens undersand he eec negaive eedback n he llwing ampliier characerisics:

More information

CHAPTER 7 CHRONOPOTENTIOMETRY. In this technique the current flowing in the cell is instantaneously stepped from

CHAPTER 7 CHRONOPOTENTIOMETRY. In this technique the current flowing in the cell is instantaneously stepped from CHAPTE 7 CHONOPOTENTIOMETY In his echnique he curren flwing in he cell is insananeusly sepped frm zer sme finie value. The sluin is n sirred and a large ecess f suppring elecrlye is presen in he sluin;

More information

Forward guidance. Fed funds target during /15/2017

Forward guidance. Fed funds target during /15/2017 Forward guidance Fed funds arge during 2004 A. A wo-dimensional characerizaion of moneary shocks (Gürkynak, Sack, and Swanson, 2005) B. Odyssean versus Delphic foreign guidance (Campbell e al., 2012) C.

More information

Impact Switch Study Modeling & Implications

Impact Switch Study Modeling & Implications L-3 Fuzing & Ordnance Sysems Impac Swich Sudy Mdeling & Implicains Dr. Dave Frankman May 13, 010 NDIA 54 h Annual Fuze Cnference This presenain cnsiss f L-3 Crprain general capabiliies infrmain ha des

More information

SMKA NAIM LILBANAT KOTA BHARU KELANTAN. SEKOLAH BERPRESTASI TINGGI. Kertas soalan ini mengandungi 7 halaman bercetak.

SMKA NAIM LILBANAT KOTA BHARU KELANTAN. SEKOLAH BERPRESTASI TINGGI. Kertas soalan ini mengandungi 7 halaman bercetak. Name : Frm :. SMKA NAIM LILBANAT 55 KOTA BHARU KELANTAN. SEKOLAH BERPRESTASI TINGGI PEPERIKSAAN PERCUBAAN SPM / ADDITIONAL MATHEMATICS Keras ½ Jam ½ Jam Unuk Kegunaan Pemeriksa Arahan:. This quesin paper

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

CAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank

CAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank CAUSAL INFERENCE Technical Track Sessin I Phillippe Leite The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Phillippe Leite fr the purpse f this wrkshp Plicy questins are causal

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

An recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes

An recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,

More information

Bayesian Dynamic Factor Analysis of a Simple Monetary DSGE Model

Bayesian Dynamic Factor Analysis of a Simple Monetary DSGE Model WP//29 Bayesian Dynamic Facr Analysis f a Simple Mneary DSGE Mdel Maxym Kryshk 20 Inernainal Mneary Fund WP//29 IMF Wrking Paper IMF Insiue Bayesian Dynamic Facr Analysis f a Simple Mneary DSGE Mdel Prepared

More information

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) > Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);

More information

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA Mdelling f Clck Behaviur Dn Percival Applied Physics Labratry University f Washingtn Seattle, Washingtn, USA verheads and paper fr talk available at http://faculty.washingtn.edu/dbp/talks.html 1 Overview

More information

Large-scale Distance Metric Learning with Uncertainty

Large-scale Distance Metric Learning with Uncertainty i Large-scale Disance Meric Learning wih Uncerainy Qi Qian Jiasheng Tang Ha Li Shenghu Zhu Rng Jin Alibaba Grup, Bellevue, WA, 98004, USA {qi.qian, jiasheng.js, liha.lh, shenghu.zhu, jinrng.jr}@alibaba-inc.cm

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Finite time L 1 Approach for Missile Overload Requirement Analysis in Terminal Guidance

Finite time L 1 Approach for Missile Overload Requirement Analysis in Terminal Guidance Chinese Jurnal Aernauics 22(29) 413-418 Chinese Jurnal Aernauics www.elsevier.cm/lcae/cja Finie ime L 1 Apprach r Missile Overlad Requiremen Analysis in Terminal Guidance Ji Dengga*, He Fenghua, Ya Yu

More information

Successive ApproxiInations and Osgood's Theorenl

Successive ApproxiInations and Osgood's Theorenl Revisa de la Unin Maemaica Argenina Vlumen 40, Niimers 3 y 4,1997. 73 Successive ApprxiInains and Osgd's Therenl Calix P. Caldern Virginia N. Vera de Seri.July 29, 1996 Absrac The Picard's mehd fr slving

More information

Efficient and Fast Simulation of RF Circuits and Systems via Spectral Method

Efficient and Fast Simulation of RF Circuits and Systems via Spectral Method Efficien and Fas Simulain f RF Circuis and Sysems via Specral Mehd 1. Prjec Summary The prpsed research will resul in a new specral algrihm, preliminary simular based n he new algrihm will be subsanially

More information

Tom BLASINGAME Texas A&M U. Slide 1

Tom BLASINGAME Texas A&M U. Slide 1 Fundamenal Fl Lecure 7 The Diffusiviy Equain fr Mulihase Fl Slide 1 Fundamenal Fl Lecure 7 The Diffusiviy Equain fr Mulihase Fl Slide Fundamenal Fl Lecure 7 The Diffusiviy Equain fr Mulihase Fl Slide 3

More information

447. Assessment of damage risk function of structural components under vibrations

447. Assessment of damage risk function of structural components under vibrations 447. Assessmen f damage risk funcin f srucural cmnens under virains J. Dulevičius, A. Žiliukas Kaunas Universiy f Technlgy, Kesuci s. 27, LT-4432 Kaunas, Lihuania e-mail: jnas.dulevicius@ku.l, ananas.ziliukas@ku.l

More information

initially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur

initially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur Cdewrd Distributin fr Frequency Sensitive Cmpetitive Learning with One Dimensinal Input Data Aristides S. Galanpuls and Stanley C. Ahalt Department f Electrical Engineering The Ohi State University Abstract

More information

Online Influence Maximization under Independent Cascade Model with Semi-Bandit Feedback

Online Influence Maximization under Independent Cascade Model with Semi-Bandit Feedback Online Influence Maximizain under Independen Cascade Mdel wih Semi-Bandi Feedback Zheng Wen Adbe Research zwen@adbe.cm Michal Valk SequeL eam, INRIA Lille - Nrd Eurpe michal.valk@inria.fr Branislav Kven

More information

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The

More information

Physics Courseware Physics I Constant Acceleration

Physics Courseware Physics I Constant Acceleration Physics Curseware Physics I Cnsan Accelerain Equains fr cnsan accelerain in dimensin x + a + a + ax + x Prblem.- In he 00-m race an ahlee acceleraes unifrmly frm res his p speed f 0m/s in he firs x5m as

More information

ELEG 205 Fall Lecture #10. Mark Mirotznik, Ph.D. Professor The University of Delaware Tel: (302)

ELEG 205 Fall Lecture #10. Mark Mirotznik, Ph.D. Professor The University of Delaware Tel: (302) EEG 05 Fall 07 ecure #0 Mark Mirznik, Ph.D. Prfessr The Universiy f Delaware Tel: (3083-4 Email: mirzni@ece.udel.edu haper 7: apacirs and Inducrs The apacir Symbl Wha hey really lk like The apacir Wha

More information

INFLUENCE OF WIND VELOCITY TO SUPPLY WATER TEMPERATURE IN HOUSE HEATING INSTALLATION AND HOT-WATER DISTRICT HEATING SYSTEM

INFLUENCE OF WIND VELOCITY TO SUPPLY WATER TEMPERATURE IN HOUSE HEATING INSTALLATION AND HOT-WATER DISTRICT HEATING SYSTEM Dr. Branislav Zivkvic, B. Eng. Faculy f Mechanical Engineering, Belgrade Universiy Predrag Zeknja, B. Eng. Belgrade Municipal DH Cmpany Angelina Kacar, B. Eng. Faculy f Agriculure, Belgrade Universiy INFLUENCE

More information

The Impact of Nonresponse Bias on the Index of Consumer Sentiment. Richard Curtin, Stanley Presser, and Eleanor Singer 1

The Impact of Nonresponse Bias on the Index of Consumer Sentiment. Richard Curtin, Stanley Presser, and Eleanor Singer 1 The Impac f Nnrespnse Bias n he Index f Cnsumer Senimen Richard Curin, Sanley Presser, and Eleanr Singer 1 Inrducin A basic ene f survey research is he abslue preference fr high respnse raes. A lw respnse

More information

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

PHY305F Electronics Laboratory I. Section 2. AC Circuit Basics: Passive and Linear Components and Circuits. Basic Concepts

PHY305F Electronics Laboratory I. Section 2. AC Circuit Basics: Passive and Linear Components and Circuits. Basic Concepts PHY305F Elecrnics abrary I Secin ircui Basics: Passie and inear mpnens and ircuis Basic nceps lernaing curren () circui analysis deals wih (sinusidally) ime-arying curren and lage signals whse ime aerage

More information

A Matrix Representation of Panel Data

A Matrix Representation of Panel Data web Extensin 6 Appendix 6.A A Matrix Representatin f Panel Data Panel data mdels cme in tw brad varieties, distinct intercept DGPs and errr cmpnent DGPs. his appendix presents matrix algebra representatins

More information

21.9 Magnetic Materials

21.9 Magnetic Materials 21.9 Magneic Maerials The inrinsic spin and rbial min f elecrns gives rise he magneic prperies f maerials è elecrn spin and rbis ac as iny curren lps. In ferrmagneic maerials grups f 10 16-10 19 neighbring

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

Money in OLG Models. 1. Introduction. Econ604. Spring Lutz Hendricks

Money in OLG Models. 1. Introduction. Econ604. Spring Lutz Hendricks Mne in OLG Mdels Ecn604. Spring 2005. Luz Hendricks. Inrducin One applicain f he mdels sudied in his curse ha will be pursued hrughu is mne. The purpse is w-fld: I prvides an inrducin he ke mdels f mne

More information

Physics 111. Exam #1. September 28, 2018

Physics 111. Exam #1. September 28, 2018 Physics xam # Sepember 8, 08 ame Please read and fllw hese insrucins carefully: Read all prblems carefully befre aemping slve hem. Yur wrk mus be legible, and he rganizain clear. Yu mus shw all wrk, including

More information

The Buck Resonant Converter

The Buck Resonant Converter EE646 Pwer Elecrnics Chaper 6 ecure Dr. Sam Abdel-Rahman The Buck Resnan Cnverer Replacg he swich by he resnan-ype swich, ba a quasi-resnan PWM buck cnverer can be shwn ha here are fur mdes f pera under

More information

Review of HAARP Experiment and Assessment of Ionospheric Effects

Review of HAARP Experiment and Assessment of Ionospheric Effects Third AL PI ympsium Kna, Hawaii Nvember 9-3, 009 Review f HAARP Experimen and Assessmen f Inspheric Effecs T. L. Ainswrh, Y. Wang, J.-. Lee, and K.-. Chen Naval Research Labrary, Washingn DC, UA CRR, Nainal

More information

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate. Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since

More information

The 37th International Physics Olympiad Singapore. Experimental Competition. Wednesday, 12 July, Sample Solution

The 37th International Physics Olympiad Singapore. Experimental Competition. Wednesday, 12 July, Sample Solution The 37h Inernainal Physics Olypiad Singapre Experienal Cpeiin Wednesday, July, 006 Saple Sluin Par a A skech f he experienal seup (n required) Receiver Raing able Gnieer Fixed ar Bea splier Gnieer Mvable

More information

6 th International Conference on Trends in Agricultural Engineering 7-9 September 2016, Prague, Czech Republic

6 th International Conference on Trends in Agricultural Engineering 7-9 September 2016, Prague, Czech Republic THEORETICAL INVESTIGATIONS OF MINERAL FERTILISER DISTRIBTION BY MEANS OF AN INCLINED CENTRIFGAL TOOL V. Bulgakv 1, O. Adamchuk, S. Ivanvs 3 1 Nainal niversiy Lie and Envirnmenal Sciences kraine Nainal

More information

Index-based Most Similar Trajectory Search

Index-based Most Similar Trajectory Search Index-based Ms Similar rajecry Search Elias Frenzs, Ksas Grasias, Yannis hedridis Labrary f Infrmain Sysems eparmen f Infrmaics Universiy f Piraeus Hellas echnical Repr Series UNIPI-ISL-R-6- Nvember 6

More information

Poincaré s Equations for Cosserat Media: Application to Shells

Poincaré s Equations for Cosserat Media: Application to Shells Pincaré s Equains fr Cssera Media: Applicain Shells Frédéric Byer, Federic Renda T cie his versin: Frédéric Byer, Federic Renda. Pincaré s Equains fr Cssera Media: Applicain Shells. Jurnal f Nnlinear Science,

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Estimation of Poses with Particle Filters

Estimation of Poses with Particle Filters Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU

More information

OBJECTIVES OF TIME SERIES ANALYSIS

OBJECTIVES OF TIME SERIES ANALYSIS OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging

More information

Index-based Most Similar Trajectory Search

Index-based Most Similar Trajectory Search Index-based Ms Similar rajecry Search Elias Frenzs Ksas Grasias Yannis hedridis ep. f Infrmaics, Universiy f Piraeus, Greece ep. f Infrmaics, Universiy f Piraeus, Greece ep. f Infrmaics, Universiy f Piraeus,

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Subject: Turbojet engines (continued); Design parameters; Effect of mass flow on thrust.

Subject: Turbojet engines (continued); Design parameters; Effect of mass flow on thrust. 16.50 Leure 19 Subje: Turbje engines (ninued; Design parameers; Effe f mass flw n hrus. In his haper we examine he quesin f hw hse he key parameers f he engine bain sme speified perfrmane a he design ndiins,

More information

Soccer Player Tracking across Uncalibrated Camera Streams

Soccer Player Tracking across Uncalibrated Camera Streams EEE nernainal rkshp n Visual Surveillance and erfrmance Evaluain f Tracking and Surveillance ETS 3 n cnjuncin wih V Ocber 3 Nice France. Sccer layer Tracking acrss Uncalibraed amera Sreams Jinman Kang

More information

Acta Scientiarum. Technology ISSN: Universidade Estadual de Maringá Brasil

Acta Scientiarum. Technology ISSN: Universidade Estadual de Maringá Brasil Aca cieniarum. Technlgy IN: 86-2563 eduem@uem.br Universidade Esadual de Maringá Brasil hang, Hsu Yang A mehdlgy fr analysis f defecive pipeline by inrducing sress cncenrain facr in beam-pipe finie elemen

More information

Answers: ( HKMO Heat Events) Created by: Mr. Francis Hung Last updated: 21 September 2018

Answers: ( HKMO Heat Events) Created by: Mr. Francis Hung Last updated: 21 September 2018 nswers: (009-0 HKMO Hea Evens) reaed by: Mr. Francis Hung Las updaed: Sepember 08 09-0 Individual 6 7 7 0 Spare 8 9 0 08 09-0 8 0 0.8 Spare Grup 6 0000 7 09 8 00 9 0 0 Individual Evens I In hw many pssible

More information

Financial Econometrics Kalman Filter: some applications to Finance University of Evry - Master 2

Financial Econometrics Kalman Filter: some applications to Finance University of Evry - Master 2 Financial Economerics Kalman Filer: some applicaions o Finance Universiy of Evry - Maser 2 Eric Bouyé January 27, 2009 Conens 1 Sae-space models 2 2 The Scalar Kalman Filer 2 21 Presenaion 2 22 Summary

More information

Business Cycles. Approaches to business cycle modeling

Business Cycles. Approaches to business cycle modeling Business Cycles 73 Business Cycles Appraches business cycle mdeling Definiin: Recurren paern f dwnswings and upswings: Acrss many indusries Wih cmmn paern f c-mvemen amng majr variables Oupu Emplymen Invesmen

More information

Math 2214 Solution Test 1A Spring 2016

Math 2214 Solution Test 1A Spring 2016 Mah 14 Soluion Tes 1A Spring 016 sec Problem 1: Wha is he larges -inerval for which ( 4) = has a guaraneed + unique soluion for iniial value (-1) = 3 according o he Exisence Uniqueness Theorem? Soluion

More information

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1 Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies

More information

Analysis on the Stability of Reservoir Soil Slope Based on Fuzzy Artificial Neural Network

Analysis on the Stability of Reservoir Soil Slope Based on Fuzzy Artificial Neural Network Research Jurnal f Applied Sciences, Engineering and Technlgy 5(2): 465-469, 2013 ISSN: 2040-7459; E-ISSN: 2040-7467 Maxwell Scientific Organizatin, 2013 Submitted: May 08, 2012 Accepted: May 29, 2012 Published:

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Fuzzy Spaces for Neutrosophic Information Representation

Fuzzy Spaces for Neutrosophic Information Representation Fzzy Spaces r Nersphic Inrmain Represenain Vasile Parasc Tarm Inrmain Technlgy Bchares Rmania e-mail: parasc.v@gmail.cm Absrac. The paper presens sme seps r mli-valed represenain nersphic inrmain. These

More information

Time-dependent behaviour of inhomogeneous Restructures: application to long term analysis of R C arch and arch-frame bridges

Time-dependent behaviour of inhomogeneous Restructures: application to long term analysis of R C arch and arch-frame bridges Timedependen behaviur f inhmgeneus Resrucures: applicain lng erm analysis f R C arch and archframe bridges E Mla* Pliecnic di Milan aly F Pigni Pliecnic di Milan aly 26h Cnference n OUR WORLD N CONCRETE

More information

1 Solutions to selected problems

1 Solutions to selected problems 1 Soluions o seleced problems 1. Le A B R n. Show ha in A in B bu in general bd A bd B. Soluion. Le x in A. Then here is ɛ > 0 such ha B ɛ (x) A B. This shows x in B. If A = [0, 1] and B = [0, 2], hen

More information

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11. 1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be

More information

The lower limit of interval efficiency in Data Envelopment Analysis

The lower limit of interval efficiency in Data Envelopment Analysis Jurnal f aa nelpmen nalysis and ecisin Science 05 N. (05) 58-66 ailable nline a www.ispacs.cm/dea lume 05, Issue, ear 05 ricle I: dea-00095, 9 Pages di:0.5899/05/dea-00095 Research ricle aa nelpmen nalysis

More information