The periodogram of fractional processes

Size: px
Start display at page:

Download "The periodogram of fractional processes"

Transcription

1 The periodogram of fracioal processes Carlos Velasco y Deparameo de Ecoomía Uiversidad Carlos III de Madrid December 4, 2005 Absrac We aalyze asympoic properies of he discree Fourier rasform ad he periodogram of ime series obaied hrough (rucaed) liear lerig of saioary processes. The class of lers coais he fracioal di erecig operaor ad is coe cies decay a a algebraic rae, implyig log-rage depedece properies for he lered processes whe he degree of iegraio is posiive. These iclude fracioal ime series which are osaioary for ay value of he memory parameer ( 6= 0) ad possibly osaioary redig ( 0:5). We cosider boh fracioal di erecig or iegraio of wealy depede ad log memory saioary ime series. The resuls obaied for he momes of he Fourier rasform ad he periodogram a Fourier frequecies i a degeeraig bad aroud he origi are weaer compared o he saioary orucaed case for > 0; bu su cie for he aalysis of parameric ad semiparameric memory esimaes. They are applied o he sudy of he properies of he log-periodogram regressio esimae of he memory parameer for Gaussia processes, for which asympoic ormaliy could o be showed usig previous resuls. However oly cosisecy ca be showed for he redig cases, 0:5 < : Several deredig ad iiializaio mechaisms are sudied ad oly local codiios o specral desiies of saioary ipu series ad rasfer fucios of lers are assumed. MSC: Primary: 62M5; Secodary: 62M0, 60G8 Keywords: Discree Fourier rasform; Log-rage depedece; Log memory; Nosaioary series; Log-periodogram regressio; Asympoic ormaliy. This paper is a revised versio of he worig paper "Log-periodogram iferece for osaioary iegraed ad di ereced processes". y Correspodece address: Deparameo de Ecoomía, Uiversidad Carlos III de Madrid, Calle Madrid 26, Geafe (Madrid), Spai. Fax: Carlos.Velasco@uc3m.es. Research fuded by he Spaish Miiserio de Educació y Ciecia, Ref. o. SEJ /ECON.

2 Iroducio Mos models for osaioary redig ad log memory ime series are based o some uderlyig wealy depede saioary processes. Thus for saioary log memory or log-rage depede processes i is ofe assumed ha x = ( L) d " ; 0 < d < 0:5; = ; 2; : : : ; where he fracioal di erece ler ( formal expasio, 6= ; 2; : : : L) d is de ed i erms of he lag operaor L hrough he ( L) d := =0 ' (d) L ; ' (d) = ( + d) (d) ( + ) ; where (x) = R 0 z x e z dz is he Gamma fucio such ha (d) = for d = 0; ; 2; : : :. Here ", = 0; ; : : : ; is assumed o be covariace saioary ad wealy depede accordig o some sadard de iio, which usually eails ha he process " has a posiive bouded specral desiy a zero frequecy. For isace, Hosig (98) ad Grager ad Joyeux (980) origially proposed ha " is a ARMA process. The x is ermed iegraed of order d, is covariace saioary wih ie variace uder d < 0:5, ad is specral desiy has a pole (zero) a zero frequecy if d > 0 (d < 0). Tredig observaios ca be modelled by parial sums (Sowell (990); Hurvich ad Ray (995); Velasco (999)), y = y 0 + x ; = ; 2; : : : ; (2) = leadig o di erece-saioary models where he icreme series ( L)y = x is saioary iegraed of order d, 0:5 d < 0:5, saisfyig. The we say ha y is iegraed of order d, 0:5 d < :5. Usig successive parial sums i is sraighforward o de e higher order iegraed series for ay d 0:5. Aleraively, i is possible o cosider (e.g. Robiso ad Mariucci (200), Phillips (999)) ha a process is geeraed by a rucaed lerig as = ' ; > 0; = ; 2; : : : ; (3) =0 so all he pas wealy depede saioary iovaios ; 0, are igored. This coveio maes esseial he dae of sar of observaio. However his framewor ca easily be geeralized seig a warmig up period where he i ow of iformaio ca begi before we acually observe he process. The lered process ; hough wih ie variace for xed, is o-saioary for ay value of 6= 0: However if < 0:5, i coverges i mea square as! o he covariace saioary x obaied i for he same sequece of iovaios = ", = ; : : : ;. Usig Sirlig s formula, we d ha ' as!, where meas ha he raio of lefad righ-had sides eds o oe. Therefore, whe 0:5 he variace of ca grow wihou limi wih ad is osaioary log-rage depede i he sese of Heyde ad Yag (997). Of course, boh ideas ca be combied o de e fracioal iegraed processes, achievig complee geeraliy i he geeraio of ime series wih log rage depede behaviour. 2

3 The log-rage properies of he processes -(2) ad (3) are described by he memory parameers d ad, ad uder regulariy codiios ad appropriaely ormalized, such processes coverge o di ere versios of fracioal Browia moio of parameers d > 0:5 ad > 0:5 respecively (see Mariucci ad Robiso (999) for a discussio). The memory parameers ca be esimaed by a variey of mehods uder geeral assumpios. Mai focus i he lieraure has bee paid o saioary log memory series (e.g. Fox ad Taqqu (986), Robiso (995a, b)), hough di erecesaioary models have bee also cosidered (e.g. Velasco (999), Velasco ad Robiso (2000)). Usig he oher aleraive, Phillips (999) cosidered asympoic properies of he discree Fourier rasform (DFT) of give by (3), where is a liear process wih coe cies i he Wold decomposiio associaed o a uiformly smooh specral desiy, exedig ideas of Phillips ad Solo (992) i erms of a Barle s approximaio o he DFT of he liear iid iovaios. Robiso (2005) cosidered bouds bewee he DFT of ad is saioary versio x give by whe = d 2 [ 0:5; 0:5) ad " =, ad sudied he e ec of aperig whe hese series are possibly furher fracioally iegraed o achieve osaioary memory levels. These resuls are used o derive he properies of di ere parameer esimaes which are liear fucios of he periodogram, ad whose aalysis had bee coduced previously for (di erece) saioary processes. I his paper we geeralize ad improve he previous resuls o he momes of he DFT ad he periodogram for frequecy domai iferece o fracioal series give by (3) followig a di ere roue. I a similar spiri o Robiso (995a), we sudy direcly he DFT of he observed series ad se our local codiios i erms of he specral desiy of he wealy depede iovaios ad of a geeral osaioary log rage depede ler as i Robiso ad Mariucci (200). We also cosider i his paper several exesios of model (3). We cosider simulaeously series wih egaive memory ( < 0), which are releva for saisical iferece o fracioally di ereced daa; processes wih lers iiialized a a poi di ere from zero; ad fracioal di erecig ad iegraio of saioary log memory ime series wih = x give, e.g. by wih 0 < d < 0:5. Our resuls are weaer ha he oes ha ca be obaied i he orucaed case for boh < 0:5 ad 0:5 (cf. Robiso (995a) ad Velasco (999), respecively). However, hey are su cie o usify valid large sample log-periodogram iferece o ; < 0:5 (Gewee ad Porer- Huda (983), Robiso (995a)) for Gaussia processes, hough we ca oly obai cosisecy resuls whe 0:5 <. The resuls allow for saioary log memory iovaios, log-periodogram regressio esimaig he overall memory of he resulig ime series i his case. Previous resuls i Phillips (999) ad Robiso (2005) could o usify such asympoic properies. Our local codiios o he specral desiy of are also weaer ad more geeral ha hose of Kim ad Phillips (999) who sudied his semiparameric esimae whe i a relaed se-up followig Phillips (999) aalysis. The resuls obaied i his paper are also useful for oher problems ad, for isace, he arrow-bad While esimae ivesigaed i Robiso (995b) for saioary log memory liear processes has bee aalyzed i Marmol ad Velasco (2004) uder (3) for liear ad < 0:5. The res of he paper is orgaized as follows. Nex secio aalyzes he covariace properies of he DFT of fracioally iegraed series, while Secio 3 cosiders he log-periodogram esimae of : Secio 4 discusses he iiializaio of he process ad di ere deredig mechaisms. The case whe is a log-memory (saioary) process is sudied i Secio 5. Some echical lemmas ad proofs are coaied i a appedix. 3

4 2 Discree Fourier rasform of fracioally iegraed ime series We follow he same oaio as i Robiso ad Mariucci (200) bu i a uivariae framewor o simplify he preseaio. We cosider liear lers more geeral ha he fracioal ( L), wih similar rs order asympoic properies. We cosider simulaeously he cases of posiive ad egaive memory parameer. The ler wih posiive memory parameer iduces posiive auocorrelaio for log lags ad is he releva case for mos daa i pracice, hough egaives values are liely o occur if some deredig mechaism, such as di erecig, has bee applied previously. We suppose ha he observed sequece is give by = = where he ler coe cies ; is rasfer fucio = ; = ; : : : ; ; =0 ad he process saisfy he followig Assumpio. exp (i) Assumpio.. f (0) g 2, where ; < ; is he class of sequeces f ; = 0; ; : : :g, such ha = ( = 0); ad + ; as! 0; (4) = O( ); as! 0:.2. (0) = P =0 = 0; ad 2 as! 0, < < 0:.3. is a zero mea covariace saioary process wih specral desiy f ; posiive ad coiuous i a ierval aroud = 0. I is easy o show ha he sequece ' give by he fracioal ler ( L) belogs o he class : I fac, he class is more geeral ha ha de ed i erms of he fracioal iegraio ler sice we allow for coribuios of a smaller order of iegraio ; > 0; which ca be of ieres i some applicaios regardig coiegraed ime series ad volailiy sequeces (see e.g. he survey i Velasco (2005)). This possibiliy will be e ecively accoued for i Assumpio 3.2. For < 0; he class could be de ed i erms of di ereces of lers belogig o ( + ): Thus = (+) (+) ; = ; 2; : : : ; 0 = ; as whe we ae fracioal di ereces of a ie legh vecor observed i = ; : : : ;. For > 0; he asympoic behavior of whe! 0, 4

5 is deduced i Lemma 2 i he Appedix, usig Assumpio., ad described by Assumpio.2 for < < 0: The codiios o imposed by Assumpio.3 are o resricive as we leave all he specral desiy f uparamerized excep a zero frequecy, where i is smooh, ad oly require iegrabiliy for covariace saioariy. We furher relax his codiio i Secio 5 ad cosider saioary log-rage depede whose specral desiy has sigulariies a = 0. We are maily cocered wih he asympoic properies of he DFT of ; = ; : : : ; ; w ( ) = (2) =2 = exp(i ); a he se of Fourier frequecies = 2=; i a degeeraig bad aroud he origi. These are he releva frequecies o describe log-rage properies of. To reproduce resuls parallel o Theorem 2 of Robiso (995a) or Theorem of Velasco (999), for saioary ad di erece saioary processes, respecively, we cosider a (possibly o-iegrable) bu saioary specrum for ; f = f 2 ; (see e.g. Velasco (999) for di erece-saioary processes). The specrum f is o-iegrable for 0:5 because we show i Lemma 2 (for > 0) ha f f (0) 2 ; as! 0; (5) wih 0 < f (0) < : This specrum describes he secod order dyamics properies of ad would be he referece arge for log-periodogram regressios because i equals he limi of he expecaio of he periodogram, as happes for (di erece) saioary process, see also Solo (992). The specrum f is he limi, as!, of he ime-varyig specral desiy f () = f 2 ; = exp(i); = which is direcly relaed o he expecaio of he periodogram I = w 2 : To see his we ca use ha E [w ( )w ( )] = 2 ; f d; (6) =0 see (8.) i Robiso ad Mariucci (200), where = [ ; ] ad ; = ( ; ) ( ; ); (; ) = exp (i( + )) ( = ): The we ca chec ha 2 ; d = ( ) 2 = ad, if f is smooh eough a = ; we ca approximae he periodogram expecaio E [I ( )] = E [w ( )w ( )] by f () ( ): Here, he fucio ; plays he oi role of 5

6 Feer s erel, K = (2) P = exp(i)2, ad of he fracioal rasfer fucio ' 2 = P =0 '(d) exp(i) 2 = exp(i) 2d i he (di erece) saioary case -(2) for d < :5, sice i his case E [I yy ( )] = K ( )f "" ' 2 d: The erel K eds o a Dirac s dela fucio as! ; so he above quaiy eds o f yy ( ) = f "" ( )'( ) 2 : The sudy of he properies of (6) based o he erel (; ) geeralizes usual saioary Fourier aalysis echiques for ime depede liear lers, possibly wih algebraic decay o zero. Now we impose some exra smoohess codiios o f a = 0 hrough ad f : Assumpio f is boudedly di ereiable for 2 ( ; ) ; some > 0: 2.2. is di ereiable i ( ; ) f0g, wih d d = O( ) as! 0: This assumpio is similar o Assumpio 2 of Theorem 2 of Robiso (995a), where he covariace marix of he DFT of saioary log memory series is aalyzed. I does o impose exra smoohess compared o his case ad i implies ha (d=d) log f = O( ) as! 0, i view of Lemma 2 ( > 0) ad Assumpio.2 ( < 0). We give ow a equivale resul for he DFT of a Fourier frequecies, whose proof ca be foud i he Appedix. Theorem Uder Assumpios ad 2, <, for < < m; m=! 0 as! ;.. E [w ( )w ( )] = f ( ) + O f ( ) + (log ) =0:5o ; > 0; = f ( ) + O f ( ) + log (log ) = 0:5o ; < 0:.2. E [w ( )w ( )] ; E [w ( )w ( )] ; E [w ( )w ( )] q = O f ( )f ( ) + (=) log + (log ) =0:5o ; > 0; q = O f ( )f ( ) + (=) log (log ) = 0:5o ; < 0: pf For = 0 all he bouds are replaced by O ( )f ( ) log ; : The bouds obaied i his heorem for he ormalized periodogram expecaio f ( )E [I ( )] he saioary oes, O( log ); of osaioary daa a Fourier frequecies are di ere for ay from 0:5 < < 0:5; give by Robiso (995a) ad repored here 6

7 for he = 0 case. They are also di ere from he di erece saioary resul, O( 2( ) log ); 0:5 < (see Velasco, 999), bu depede o he degree of iegraio as well. For isace, he ew coribuio i whe > 0 arises from he emporal ihomogeeiy of ; = ; : : : ; ; whe he ime ivaria specral desiy f is compared wih he ime-varyig specral desiy f () : However, whe 0:5 < < 0 he resuls i Theorem are he same as i he saioary case, up o a logarihmic facor, ad, as expeced, he e ecs of rucaio seem o be alleviaed i case of egaive memory. Robiso (2005) obaied for fracioal series as i (3) ha h f ( )E w ( ) w x ( ) 2i = O (log ) = 0:5 (7) for 2 [ 0:5; 0:5) whe x is give by for d ad he same sequece of iovaios. Usig his resul ad a riagle iequaliy i his equaio (C.4), ogeher wih Robiso (995a, Theorem 2(i)), we ca obai ha for 2 (0; 0:5) ; f ( )E [I ( )] = + O =2 : This resul is improved by he correspodig boud O( ) i our Theorem for ay < 0:5. This improveme is impora o allow for he asympoic ormaliy of he log-periodogram regressio esimae of ; for < 0:5; bu oly cosisecy ca be obaied whe 0:5; see ex secio. Usig he same argumes i he proof of Theorem ad cosiderig oly he coribuio from 2 i (2), we ca obai uder Assumpios ad 2 ha he lef had side of (7) is O + log for 2 (0; 0:5) ; see (24), which geeralizes Robiso s (2005) resul o log memory lers i he class by meas of local assumpios o he smoohess of f ad a = 0: Phillips (999) provides some probabiliy bouds o he di erece w ( ) e i w ( ); > 0; for fracioally iegraed a each frequecy ; eiher xed or i a degeeraig bad. Apar from some correcio erms, similar o he oes discussed i Secio 4 below, which are O p =2 whe < 0:5; he remaider is O p, which would correspod o a similar boud o he oe i our Theorem., > 0. Taig hose bouds as uiform, some properies of di ere esimaes are deduced i Phillips (999), ad some ad-hoc modi caios of he usual semiparameric procedures of memory esimaio ca be proposed, see also Phillips ad Kim (999). However, oe ha for he aalysis of he properies of saisics which are oliear fucios of he periodogram of a observed fracioally iegraed series, such as he log periodogram memory esimae sudied i ex secio, he aalysis of di ere approximaios of he DFT of daa is o su cie. Uiform bouds o he periodogram expecaios of uder Gaussiaiy are required i his case because i is o eough o corol he (liear) disace wih he DFT of he saioary versio of he origial process, cf. (7). For he aalysis of he periodogram properies a all frequecies i [0; ], icludig hose xed for which is of he same order as, we eed o mae Assumpio 2 uiform for ay 2 [0; ]: The he resuls of Theorem will hold ad full parameric esimaes for osaioary redig daa could be sudied as i Velasco ad Robiso (2000). I ex secio we cocerae o a paricular semiparameric esimae uder our local codiios. 7

8 3 The log-periodogram regressio esimae of The semiparameric esimae of he memory parameer based o he log-periodogram regressio esimae was proposed by Gewee ad Porer-Huda (983). Robiso (995a) showed ha he periodogram of saioary log memory series is asympoically ubiased ad ucorrelaed, as wih shor memory, whe evaluaed a he harmoic frequecies for growig wih sample size. This is he basis o wrie he logarihm of (5) as a liear regressio model wih approximaely homoscedasic ad ucorrelaed errors, log I ( ) log f (0) 2 log + log I ( ) ; = ; : : : ; m; (8) f ( ) where m is small compared wih : The log-periodogram regressio esimae b is he leas squares esimae of ; 0 m b (r r) 2 A = m = (r r) log I ( ); where he log-periodogram is he depede variable ad he regressor is r = 2 log, r = m P m = r. The asympoic properies of b were aalyzed rigourously for muliple saioary Gaussia series ( 0:5 < < 0:5) by Robiso (995a). He cosidered a poolig of coribuios from adace frequecies o achieve e ciecy gais ad showed ha he periodograms a he very rs frequecies do o have ice asympoic saisical properies, so followig Küsch (987) he proposed o exclude he rs ` frequecies from he regressio. Hurvich, Deo ad Brodsy (998), HDB heceforh, have show uder some addiioal codiios ha rimmig of very low frequecies may o be ecessary for he aalysis of he asympoic properies of b. HDB oly cosider fracioal processes wih wealy depede iovaios which possess a specral desiy wih hree bouded derivaives a = 0, ad he regressor i (8) which arises aurally for fracioal processes, z = log(4 si 2 ( =2)); as origially proposed by Gewee ad Porer-Huda. However, we geeralize heir resuls uder he followig assumpio. We say ha a fucio g (x) is Hölder, 0 < ; i a ierval of he origi if g (x) g (0) Cx uiformly for x 2 ( ; ) ; for some > 0 ad C < : Assumpio 3 Furher o Assumpio 2, we se 3.. For some 2 (; 3], for i a ierval of he origi, f is di ereiable wih f 0 Hölder( ), < 2; f is wice di ereiable wih f 00 Hölder( 2), 2 < For some 2 (0; 3], 2 = 2 ( + O( )) as! 0, 0 < 2; 2 = 2 ( + A 2 + O( )) as! 0, 2 < 3: Assumpio 3 imposes a rae o he approximaio (5), beig Assumpio 3.2 similar o Robiso s (995a) Assumpio. The reaso for cosiderig values of ad larger ha 2 is o esimae explicily he bias of ^, hough cosisecy oly requires mi f; g > 0 ad he bes possible rae for he asympoic bias is already achieved by mif; g = 2. Thus for fracioal ime series, we 8

9 have ha 2 = (2 si =2) 2 = O( 3 ) ; as! 0; so = 3 ad followig Robiso (995a), for hese models we have ha f = f (0) 2 f 00 (0) + f (0) O( ) ; as! 0; 2 for > 2, usig ha f(0) 0 = 0. For example, he presece of addiive oise i of smaller order of iegraio ha would lead o values of smaller ha 2. Therefore, we ca adap he bias esimaio resuls of HDB, hough hey de ed he log-periodogram regressio esimae usig he fracioal regressors z. Theorem 2 Le 0:5 < < 0:5. Uder Assumpios, 2 ad 3, = mif; g 2 (2; 3], ad m + (m ) 2 log m! 0 as! ; (9) for Gaussia E(^) = f 00 (0) m 2 f (0) + A 2 ( + o) + m Var(^) = 2 24m + o(m ); where m := O(m log m) >0 + O(m log 2 m) <0 + O(m log 3 m) =0 : Proof of Theorem 2. Follows from HDB, us usig our Theorem isead of Theorem 2 or Robiso (995a). The all Lemmas of HDB hold uder he codiios of he heorem, us seig i heir Lemma 3 ha = +(=) log ; 0 < < 0:5; log m < m; m +m! 0 ad he correspodig modi caios for 0 from Theorem. The HDB s Lemma 4 follows ow from he appropriae modi caio of Lemma 4. of Robiso ad Mariucci (200). I Lemmas 6 ad 7 of HDB we obai ha E[" ] = where = O ; > 0; = O ; < 0; ad Var[" ] = 2 =6 + O( ); log m m; whereas i Lemma 8 of HDB, he boud is ow O(m log m); > 0; O(m log 2 m); < 0; uder (9). The codiio (m ) 2 log m! 0 ad our Hölder assumpio are used isead of hree bouded derivaives i HDB, while he codiio m log m! 0 i HDB is implied by (9), 3: If we use z isead of r i he log-periodogram regressio ad 2 = e i 2 ; he he erm A does o show up i Theorem 2. Whe 0:5; i is o possible o show he p m-cosisecy of ^; bu us adapig he mehods of he previous heorem we ca prove he followig resul which gives codiios o m for he cosisecy of ^ for ay value > 0 ad 0:5 < <, 6= 0:5. Theorem 3 Le 0:5 < <, 6= 0:5. Uder Assumpios, 2 ad 3, = mif; g 2 (0; 2], ad m + (m ) log m! 0 as! ; (0) 9

10 for Gaussia E(^) = O (m ) log m + m ; Var(^) = O m + O m 2( ) log 2 m : The badwidh codiios (9) ad 0 hold if m K a for 0 < a < : However we are o able o show he cosisecy whe = 0:5 ad 0:5 ad all Fourier frequecies ; : : : ; m are used because of he erms growig wih log ad i Theorem. However, if we wish o rim ou he rs ` frequecies i he log-periodogram regressio, for some ` growig slowly wih m bu faser ha log i he asympoics, he i is easy o adap HDB s argumes ad a similar resul o Theorem 3 is valid also for = 0:5 ad 0:5: Fially if 0:5 < < 0:5 we ca adap Robiso s (995a) ceral limi heorem as i HBD o avoid ay rimmig of low frequecies, ad usig a badwidh m K a such ha 0 < a < 2=(2 + ): Theorem 4 Le 0:5 < < 0:5. Uder Assumpios, 2 ad 3, = mif; g 2 (0; 2], ad m log 2 + m 2+ 2! 0 as! ; for Gaussia p m(^ )! d N(0; 2 24 ): Proof of Theorem 4. Follows as Theorem 2 of HDB from Robiso (995a), adapig for a geeral smoohess codiio o f i erms of. The aalysis of o-gaussia ad apered series ca also be pursued usig he ideas pu forward here (see Velasco (2000), Hurvich, Moulies ad Soulier (2002)). 4 Periodogram modi caios ad iiial codiios I his secio we discuss some modi caios of he periodogram proposed i he lieraure, which are aleraives o he usual deredig procedure of osaioary series cosisig o aig di ereces. We also brie y explore he implicaios of he de iio of fracioally iegraed processes i erms of a rucaed ler i he previous asympoic heory. Phillips (999) proposed he followig correcio for he DFT a each usig a represeaio of he DFT ha ivolves he las observaio (we assume o iiial codiios, 0 0 i his oaio), w ( ) := w ( ) + p e i 2 e ; 6= 0mod: i Phillips (999) ad Kim ad Phillips (999) moivaed such correcio i erms of improvig he properies of he DFT for redig processes wih memory aroud = ; ad i is suiable o deal wih he cases where, sice i his case he usual log periodogram regressio esimae is o cosise. I fac, his correcio ca be see as a liear deredig i he ime domai seig = ; = ; : : : ; ; 0

11 so = 0 ad w ( ) is acually he DFT of. If we se addiioally = 0; so ( is easy o chec ha e i w ( ) = w ( ) e ; i L) = 0; i for he icremes = ( L). Therefore, he aalysis of Secios 2 ad 3 is valid for osaioary series wih < 2 replacig by = because, accordig o (5), has memory ad we ca se := e i ( e i ) where = ( e i ). I his sese he rasformaio is equivale o aig rs di ereces, ad higher order correcios ca be developed for furher di erecig. However whe 0 < < ad he modi ed DFT is used isead of w ( ), he implici process has egaive memory (because ypically has) ad he properies of is DFT ad log-periodogram esimaes are immediae cosequeces of our previous resuls wih < 0. We ow explore he iiial codiios problem. Le assume for some T for = (T ) ; = = T : This icludes he se-ups of Robiso ad Mariucci (200), T = he followig model ; Phillips (999), T = 0; ad oher possibiliies wih T xed wih or wih T icreasig wih, e.g. T = b ; b > 0; which implies ha he srech of iiial codiios icreases wih he widow of observed daa. As far as T is xed, i is easy o chec i he Appedix ha all he resuls of Secios 2 ad 3 hold as if T = us modifyig he boud for 2 +T + 2 used i he proof of Theorem. We oly cosider he case where 0; sice egaive values of usually origiae from some di erecig process give a sample of size : Corollary The coclusios of Theorems -4 are valid for ay T 0, 0: I fac, if T is xed or O(), he asympoic bouds give i Theorem ca o be easily improved furher, bu if T is icreasig fas eough wih ; some improvemes seem possible. I paricular, he boud O( ) i Theorem, > 0; ca be muliplied by a facor (T ) which is o as! if = o(t ); 0 < < : Noe ha if T! as! (iiial codiios growig faser ha observed daa) he we are closer o he (di erece) saioary framewor of Hurvich ad Ray (995) ad Velasco (999), ad he coribuio of he pas iovaios o each observaio, P (T ) := P 0 = T = P T + =, ca be showed o be O p ( 2 ); 0 < < 0:5; or O p ((T + ) 2 ) for 0:5; uder some regulariy codiios o : Thus, assumig ha Cov[ ; ] = O( ); > ; we have

12 ha Var [P (T )] = T + T + = = Cov[ ; ] T + C T + = =+ T + C = T + 2( ) = T + C 2( ) ; = ad he claim follows. For redig observaios wih > 0:5; he ime-depede iiial codiios ca have a domia coribuio compared o he iformaio accumulaed i he curre periods ; 2; : : : ; ; as log as =T! 0: 5 Fracioal iegraio ad di erecig of saioary lograge depede ime series We cosider i his secio he siuaio where a ie observed srech, = ; : : : ;, of a saioary log memory ime series x, wih memory 0:5 < d < 0:5 is fracioally iegraed (di ereced) of order > 0 ( < 0). These are he basic operaios performed whe usig fracioal values of o boh asympoically saioary ad redig processes, where we use a esimae close o d such ha he lered series is ear saioary ad shor-rage depede. We oly cosider 0:5 < 0:5 sice he fracioal lers are cumulaive i he sese ha if 2, ; > 0 he =0 2 ( + ); 2, ad (cf. Robiso ad Mariucci s (200) Lemma 3.), as moivaed i Secio 2 whe preseig he lers wih egaive. Oherwise we ca rs iegrae or ae di ereces, ad he properies of he DFT of hese processes ca be obaied as a direc cosequece of he resuls of his secio. Assumpio 4 For f g 2 as i Assumpio, <, = x ; (2) = wih d d = O( ) as! 0; ad where x is covariace saioary wih specral desiy saisfyig f xx = O( 2d ); d < 0:5; ad d d f xx = O(f xx ) as! 0: 0:5 2

13 The class of processes x de ed Assumpio 4 is more geeral ha, sice we do o require x o be a fracioal series, oly o have a specral desiy wih a similar asympoic behaviour as! 0. Therefore, he se of fracioally iegraed processes described by Assumpio 4 is more geeral ha ha of Phillips (999) where oly uiformly smooh specral desiies f xx are cosidered, so x is ecessarily shor memory. Usig our previous de iio, i holds ha f = 2 f xx G 2(d+) as! 0; ad i order o esimae he expoe h := d + we aalyze he asympoic properies of he momes of he DFT of log memory -fracioally iegraed or di ereced ime series. I order o dered a observed series, he ideally d; so h 0; ad upper bouds o h should o be very resricive i pracical applicaios. The followig heorems ca be cosidered ypically for series wih d < 0 ( > 0) ad d > 0 ( < 0) respecively, for which h is small. The proof is i he Appedix. Theorem 5 Uder Assumpio 4, 0:5 < 0:5, d 2 [ 0:5; 0:5); h < 0:5; h = d + ; for < < m; m=! 0 as! ; 5.. E [w ( )w ( )] = f ( ) + O f ( ) ; > 0; = f ( ) + O f ( ) (log ) = 0:5 ; < 0: 5.2. E [w ( )w ( )] ; E [w ( )w ( )] ; E [w ( )w ( )] q = O f ( )f ( ) + (=) +d log o ; > 0; q h = O f ( )f ( ) (log ) = 0:5 + (=) +d log i ; < 0: pf For = 0; all he bouds are O ( )f ( ) log ; : Noe ha whe d = 0 we basically recover he bouds of Theorem, us chagig by i he secod par of his heorem, because Assumpio 4 replaces he bouded di ereiabiliy of f a he origi, so ow (d=d) log f xx = O : Furhermore, we also recover he resuls from Robiso s (995a) Theorem 2 for saioary log memory ime series ( = 0), excep for he erm whe > 0; so i is sill possible o reproduce all he resuls of HBD o he asympoic properies of he log-periodogram esimae of he memory h = + d for saioary Gaussia ime series wih memory d possibly di ere from zero, iegraed or di ereced -imes (by he ( L) operaor i = ; : : : ; ). To esimae he MSE ad obai he asympoic disribuio of ^h we iroduce his ew assumpio, which geeralizes Assumpio 3. We omi he proof because is similar o ha of previous resuls. Assumpio 5 Furher o Assumpio 4, we se f xx = 2d f xx ad 5.. For some 2 (; 3], 2 (0; ); > 0; f xx is di ereiable wih f 0 xx Hölder( ), < 2; f xx is wice di ereiable wih f 00 xx Hölder( 2), 2 < 3: 3

14 5.2. For some 2 (0; 3], 2 = 2 ( + O( )) as! 0, 0 < 2; 2 = 2 ( + A 2 + O( )) as! 0, 2 < 3: Theorem 6 Le < 0:5; h = d + 2 ( 0:5; 0:5). The he coclusios of Theorems 2 ad 4 hold for he log periodogram esimae ^h of Gaussia series (2) uder Assumpio 5 (replacig f by fxx), for 2 (2; 3] ad 2 (0; 2] respecively; = mi f; g ; ad he same choices of m: The heory developed i his secio does o iclude he siuaios where a iiial esimae bd is compued ad he sequece ^d ; = ; : : : ;, is obaied by fracioal di erecig. This operaio should approach he wea depede iovaios d if ^d coverges fas eough o he rue d; bu i geeral he limi migh deped o he disribuio of ^d. This problem also relaes o memory esimaio based o residuals, obaied eiher from deermiisic ime deredig or by sochasic deredig i coiegraed sysems, as is aalyzed i Marmol ad Velasco (200) ad Hassler, Marmol ad Velasco (2006) respecively. Similar echiques o hose of his paper ca also be used for he aalysis of log-rage depede ime series geeraed by di ere simulaio algorihms. 6 Appedix: Proofs of resuls I he followig lemma we collec some echical resuls abou ad o be used i he proof of he resuls coaied i his appedix ad which ca be of idepede ieres for Fourier aalysis wih osaioary lers. Noe ha he boud for improves upo (8.3) of Robiso ad Mariucci (200), wha is ey for may of our resuls. Lemma 3 below will exed hese resuls for < 0. I he sequel C deoes a geeric cosa ha may chage each ime is used. Lemma 2 Uder Assumpio, for f g 2 ; 0 < < ; as! 0; 2 2 ; C mi f ; g ; C mi ; ; where = ; ad (; ) C mi ; + mi f ; [ + ]g : Proof of Lemma 2. The bouds for 2, ad are a direc cosequece of Lemma 3.2 i Robiso ad Mariucci (200). To obai he bouds for (; ); rs we have ha, wih C ; (; ) = e i(+) ( ) C ( = while, usig C ; = ) C C + ; (3) = (; ) C = C : (4) Now, de ig D (x) = P = eix ad usig summaio by pars, (; ) = ( )D ( + ) e i D ( + ) (5) = 4

15 (compare o (8.2) i Robiso ad Mariucci (200), which we are o able o reproduce here, ad wih (8.3), which is a sraighforward cosequece of our resuls). The, usig ha D (x) mifx ; g; we have ha ( ca obai ha e i D ( + ) C + C + = )D ( + ) C mif ; g mif + ; g: Similarly we = = C + : (6) Fially, usig ha D (x) = e ix e ix = e ix we have ha he secod erm i (5) is e i(+) e i(+) f ( )g (7) ad is modulus is bouded by C + f + g ad he lemma follows from his boud ad (3), (4) ad (6). Lemma 3 Uder Assumpio, for f g 2, < < 0; C f + g; C mi ; ; ad (; ) C f + + g mif; + g. Proof of Lemma 3. Firs we have for he bouds of ; = = 0 e i = = < < 0; ha A = O( ); while also = = + 0 e i = a= = 2 A = O : The, because (0) = P =0 = 0 we have ha (0) = P = = O( ); ad also = (0) + O 0 ( ) for some 2 [0; ]: The, we ca chec ha 0 ( P ) = O =0 = + ; so = O + + = O ( + ) : Therefore, if ; he = O ( ) ; while if > ; + C + C ; so C f + g : For he aalysis of (; ) we ca use (5) ad (7), wih D (x) Cx ; obaiig ha (; ) C + + C + f + g C C C ; ad usig D (x) C ad (x) ; (; ) C + + C C + : 5

16 Proof of Theorem. We oly prove he bouds for he rs wo expecaios, sice he proof for he res is similar bu simpler. We rs deal wih he case > 0: Proceedig as i equaio (8.2) of Robiso ad Mariucci (200), ; d = = e i( ) p pe i( p) D mi(;p) ( ) = 2 = 2 = e i( ) a=0 a e ia b=0 b e ib e i(+ )( ) ( ) ( ): (8) = The E [w ( )w ( )] f ( ) = 2 ; ff f ( )g d (9) ( ) +f ( ) ( ) 2 ( ) 2 ; (20) = because R ;d = 2 P = ( ) 2 : Wriig 2 2 = ( ) + ( ) + 2 ; (2) we obai from Lemma 2 for > 0; 2 2 C + 2 mif; () 2( o ) g ; (22) (compare o discussio o p. 968 of Robiso ad Mariucci (200)) ad herefore ( ) 2 ( ) 2o (23) C o 2 ( ) + mif; ( ) 2( ) g 8 9 C < 2 : = 2( ) + + ( ) ; = = = O f ( ) + ( + log f = 0:5g) + 2( ) 2( )o (24) = O f ( ) + log f = 0:5g : (25) I view of (25), (20) is O f ( ) + log f = 0:5g ; > 0 Fix oe > 0 such ha f is boudedly di ereiable i he ierval ( wrie (9) as ; ). The we ca 2 ; ( + ) ff ( + ) f ( )g d; (26) 6

17 where ; ( + ) = ( ; ) ( ; + ); ad from Lemma 2, ( ; ) CL mi f ; [ + + ]g ; > 0; where L = ; < ; = ; ; so L mif ; g. Now we cosider di ere iervals iside [ ; ] for he iegral i (26). Firs we have ha < =2 C C < =2 < =2 usig ha R < =2 L d = O (log ), ad ( ; ) 2 d 2 L d = O f ( ) log ; =2 3 =2 =2 C ( ; ) 2 d 3 =2 =2 C + 2 d + C 3 =2 2 =2 3 =2 d =2 C =2 2 d + C 2 C 2 = O f ( ) ; whe < =2; while for =2 < < we obai ha =2 3 =2 =2 C ( ; ) 2 d 3 =2 =2 C mi 2 ; + 3 = L 2 d =2 C mi 2 ; =2 2 d + 2 C 2 (log ) =0:5 = O f ( ) 2 (log ) =0:5 = O f ( ) 2( ) (log ) =0:5 : Nex =2 C ( ; ) 2 d C L d = O 2 log = O f ( ) log ; while he coribuio from he ierval [ ; 3 =2] ca be deal wih similarly. 7

18 Fially, usig he iegrabiliy of f ; C C = ( ; ) 2 ff ( + ) + f ( )g d ! 2 O + f ( + )d 2 ff ( + ) + f ( )g d = O f ( ) : Now o boud E [w ( )w ( )] = (2) R ;f d we disiguish wo cases, values of ad close ad far away. If =2 <, usig ha P ei(+ )( ) = P ) = 0; ei( 6= (mod ), (2) R ;d is by (8) equal o e i(+ )( ) ( ) ( ) ( ) 2o + e i(+ )( ) ( ) 2 ( ) 2o = e i(+ )( ) ( ) f( ( ) ( )) + (( ) ( ))g + (( ) ( )) e i(+ )( ) ( ) + e i(+ )( ) ( ) 2 ( ) 2o : (27) Usig d d C (MVT), as! 0; as = O( ); =2 < ; ad he Mea Value Theorem ( ) ( ) + ( ) ( ) C ; ( ) ( ) C ; while usig summaio by pars e i(+ )( ) ( ) = ) ei(+)( a= a ( ) a+ ( ) a a= = C a = O( ): e i( ) (28) The, as = O( ); =2 < ; ad usig (22), we obai e i(+ )( ) ( ) ( ) C 2 ( ) h i + mif; () 2( ) g + 2 C + ( + log 0:5 ) = O qf ( )f ( ) + log 0:5 : (29) 8

19 Whe < =2 he 2 ; ad usig summaio by pars, 2 ; d = ( ) ( ) + ( ) + ( ) = a= e i(+ a)( ) (30) + ( ) ( ) e i(+ )( ) : (3) Now (3) is exacly zero whe 6= (mod ); ad wih P g; he modulus of he righ had side of (30) is bouded by a= ei(+ a)( ) C mif; C ( ) ( ) + ( ) + + ( ) ( ) + ( ) = f ( ) + ( )g C = = O ( ) =! = O + = O qf ( )f ( ) : (32) Therefore, from (29) ad (32) we obai ha, < < =2; 2 ; d = O qf ( )f ( ) + log 0:5 : Nex, wih ; = ( ; ) ( ; ), ( ; ) CL ( ) mif ; + g; followig Robiso (995a, p. 063), we obai ha E [w ( )w ( )] = 2 2 ; ff f ( )g d (33) ( + )=2 + 2 (+ )=2 ; ff f ( )g d (34) (+ )=2 2 ff ( ) f ( )g ; d (35) ( + ) + 2 ; ff f ( )g d: (36) 2 Now usig he di ereiabiliy of f i 2 ( ; ); (33) is bouded by sup f ( ; ) ( ; ) d ( + )=22 ( + )=2 C ( ) 2 ( + )=2 L ( ) d = O ( ) log : 9

20 Nex (34) is bouded by C 2 2 sup f 0 ( + )=2 (+ )=2 L ( ) d = O 2 log : Now (35) is bouded by, 2 ( ) sup f 0 = O 2 (+ )=2 L ( (+ )=2 ) d! ( ; ) ( ; ) d = O 2 log : Fially, cosider (36). If < 0:5; he 2 ( ; ) ( ; ) ff f ( )g d C sup f 0 L ( ) + + d C 2 + C 2 ; ad if 0:5 2 ( ; ) ( ; ) ff f ( )g d C sup f 0 L ( ) ^ C 2 (log ) =0:5 + C ( ) (log ) =0:5 C ( ) : + ^ o C 2 2 2( ) (log ) =0:5 + d The bouds for he remaiig iervals follow ow i a simpler way. The usig O 2 log = pf O ( )f ( )(=) log he heorem follows for > 0. We ow deal wih he case < 0: Followig he proof for > 0, usig (2) ad Lemma 2 ow, we have ha (23) is bouded by, 0:5 < < 0; 8 C < 2 + : = = 2 + = 2 + = 2 9 = ; = O 2 = O f ( ) : (37) ad if < 0:5; we obai ha (23) is O 2 + (log ) f=0:5g = O f ( ) (log ) f=0:5go : (38) 20

21 We have o boud also he followig iegrals, " # ; ( + ) ff ( + ) < <2 2 < f ( )g d: Now ; (+ ) = ( ; ) ( ; +); ( ; ) C f + + g L, so < C < ( ; ) 2 d C 2 d = O f ( ) ; < <2 C <2 ( ; ) 2 d C 2 L d = O f ( ) log ; <2 2 < C C 2 < 2 < ( ; ) 2 d L d = O = O f ( ) 2 2 ; ad ally, usig he iegrabiliy of f ; C C ( ; ) 2 ff ( + ) + f ( )g d 2 ff ( + ) + f ( )g d! = O + f ( + )d = O f ( ) 2 2 : For he proof of he secod saeme of he heorem, we have ha E [w ( )w ( )] = (2) R ;f d. Cosider rs =2 < : Usig d d C = O( ); =2 < ; ad he MVT, ( ) ( ) + ( ) ( ) C ( ); as! 0; as ( ) ( ) C ; 2

22 while by summaio by pars e i(+ )( ) ( ) = ) ei(+)( C a= a a ( ) a+ ( ) a= + C a= + a a = e i( ) = O( ): The, from (27), because = O( ); =2 < ; ad followig he mehods used o obai (37) ad (38), e i(+ )( ) ( ) ( ) C ( ) 2 ( ) + + ( ) 2 ( ) 2 o = O f ( ) log ( = 0:5) : Whe < =2 he 2 ; ad usig summaio by pars as i (30)-(3), we have ha (3) is zero whe 6= (mod ); ad usig Lemma 2 he modulus of (30) is bouded by C ( ) ( ) + ( ) + + ( ) ( ) + ( ) = ^ C = f ( ) + ( )g + = + f ( ) + ( )g C = + C = + f ( ) + ( )g C C log ( = 0:5) + C o = O f ( ) log ( = 0:5) : The, ; = ( ; ) ( ; ), ( ; ) C f + g mif + ; g; < 0; so followig Robiso (995a, p. 063), we are led o esimae he coribuio of (33) o (36) whe < < 0. Usig he di ereiabiliy of f i 2 ( ; ); (33) is bouded by sup f ( ; ) ( ; ) d ( + )=22 ( + )=2 2 C 2 L ( ( + )=2 ) d = O 2 log : 22

23 Nex (34) is bouded by C 2 ad (35) by sup f 0 ( + )=2 (+ )=2 2 ( ) sup f 0 = O (+ )=2 (! ( ; ) d = O 2 L d (+ )=2 ; ) d! = O 2 log ; ( ; ) ( ; ) d = O 2 log : Fially (36) is, usig he iegrabiliy of f ; ( 2 ( C ) ) ( ; ) ( ; ) ff f ( )g d f d = O : The usig O pf = O ( )f ( )() 2 O (=) 2 () 2 pf = O ( )f ( )() 2 ad O 2 = p f ( )f ( ) he heorem follows. Proof of Corollary. I is eough o show ha Lemma 2 holds. The proof for Lemma 3 is similar. To show ha he approximaios of +T + by have he same bouds as whe T = ; we have ow ha ; d = 2 e i(+ )( ) +T + ( ) +T + ( ): The +T + C mif ; (T + ) g ad T ++ C (T + ) C 0 ; T ++ C ; > 0; so we ca obai 2 +T + 2 = +T + T ++ ( ) + +T + ( ) T T + 2 ad ha, 0 < < ; 2 +T + 2 o C 2 [(T + )] mif; ((T + )) 2( ) g C 2 [] + [] 2( )o : The ;T ( ; ) = P ei(+) +T + ( ) = P +T +e i( +T +) D ( + ) sais es he same bouds as ( ; ) i Lemma ad he corollary follows. 23

24 Proof of Theorem 5. We oly prove agai he boud for he rs wo expecaios. We cosider agai rs he case wih > 0: Now E [w ( )w ( )] f ( ) is E [I ( )] f ( ) = 2 ; ff xx f xx ( )g d (39) +f xx ( ) 2 ; d ( ) 2 ; where he secod lie is O f ( ) ; > 0; as i Theorem. For he aalysis of (39) we cosider he same iervals of iegraio as i Theorem ad oly replace f by f xx ; oig ha f xx is di ereiable i (0; ); for some > 0; wih f 0 xx = O( 2d ) as! 0+. The all bouds i Theorem for (9), 0 < < ; should be muliplied by ; obaiig a boud for (39) of order O f ( ) log = O f ( ) log : The oly ierval ha requires special cosideraio is =2 3 =2 =2 C (f xx ( ) + f xx ( + )) ( ; ) 2 d C 2 3 =2 = d d + C 2 =2 3 =2 f xx ( ) 2 d C 2 2 2d + C 2 2d = O f ( ) : For he secod expecaio cosidered i he heorem we ca cosider agai he same decomposiio (33)-(36) replacig f by f xx. Now proceedig as i he proof of Theorem 2(c) i Robiso (995a), cosiderig large ad small ; he bouds should be adaped by muliplicaio pfxx by 2d + 2d ; so we obai he boud O ( )f xx ( ) 2 (=) d log = pf O ( )f ( )(=) +d log : The oly ierval ha requires furher sudy is 2 ( ; ) ( ; ) ff xx f xx ( )g d C =2 2d + + d +C =2 f xx ( ) C 2d 2 = O + + q f ( )f ( )(=) +d log d : We ow cosider he case < 0: For aalyzig he coribuio of (20) we muliply he boud for (23) i Theorem by f xx ( ); obaiig O f ( ) : For he erm (9), we have o cosider 24

25 di ere iervals of iegraio of (26), < =2 C fxx( 0 ) ( ; ) 2 d < =2 C 2 2d L d = O f ( ) ; < =2 =2 2 =2 (f xx ( ) + f xx ( + )) L 2 d C 2 2 C 2 2 2d+ + =2 2 2d d! C 2 2d = O f ( ) ; ad ally, he ypical erm for he remaiig iervals is give by 2 C f xx ( ) L 2 2 d 2 +C max 2 d 0:5 f xx ( + ) C f xx ( ) 2 + C d 0:5 = O f ( ) : 2 L 2 0:5 d 2 d 0:5 d 2 For he secod expecaio, rs we ae (38) muliplied by p f xx ( )f xx ( ) = O The i is possible o esimae (33) o (36) for d d : 0:5 < < 0 usig he same mehod of Robiso (995a, pp ), because wih he boud ( ; ) C f + g L ( +); < 0; we obai he desired resul muliplyig a improved boud ha ca be obaied collecig Robiso s resuls for = 0; O ( )f xx ( ) (=) d log = O ( )f xx ( ) log ; pfxx pfxx by he facor (=) ; < 0: We omi he deails. Refereces Gewee, J. ad S. Porer-Huda (983). The esimaio ad applicaio of log memory ime series models. Joural of Time Series Aalysis 4, Fox, R. ad M.S. Taqqu (986). Large-sample properies of parameer esimaes for srogly depede saioary Gaussia imes series. Aals of Saisics 4, Grager, C.W.J. ad R. Joyeux (980). A iroducio o log-memory ime series ad fracioal di erecig. Joural of Time Series Aalysis, Hassler, U., F. Marmol ad C. Velasco (2006). Residual periodogram iferece for log-ru relaioships. Joural of Ecoomerics 30,

26 Heyde, C.C. ad Y. Yag (997). O de ig log-rage depedece. Joural Applied Probabiliy 34, Hosig, J.R.M. (98). Fracioal di erecig. Biomeria 68, Hurvich, C.M. ad B.K. Ray (995). Esimaio of he memory parameer for osaioary or oiverible fracioally iegraed processes. Joural of Time Series Aalysis 6, Hurvich, C.M., R. Deo ad J. Brodsy (998). The mea square error of Gewee ad Porer-Huda s Esimaor of he memory parameer of a log memory ime series. Joural of Time Series Aalysis 9, Hurvich, C.M., E. Moulies ad P. Soulier (2002). The FEP esimaor for poeially osaioary liear ime series. Sochasic Processes ad heir Applicaios 97, Kim, C.S. ad P.C.B. Phillips (999). Log-periodogram regressio: The osaioary case. Prepri, Cowles Foudaio for Research i Ecoomics, Yale Uiversiy. Küsch, H.R. (987). Saisical aspecs of self-similar processes. Proceedigs s World Cogress of he Beroulli Sociey, VNU Sciece Press. Mariucci, D. ad P.M. Robiso (2000). Wea covergece of mulivariae fracioal processes. Sochasic Processes ad heir Applicaios 86, Marmol, F. ad C. Velasco (200). Tred saioariy versus log rage depedece i ime series aalysis. Joural of Ecoomerics 08, Marmol, F. ad C. Velasco (2004). Cosise esig of fracioal coiegraio. Ecoomerica 72, Phillips, P.C.B. (999). Discree Fourier rasforms of fracioal processes. Cowles Foudaio for Research i Ecoomics Discussio Paper 243, Yale Uiversiy. Phillips, P.C.B. ad V. Solo (992). Asympoics for liear processes. Aals of Saisics 20, Robiso, P.M. (995a). Log-periodogram regressio of ime series wih log rage depedece. Aals of Saisics 23, Robiso, P.M. (995b). Gaussia semiparameric esimaio of log rage depedece. Aals of Saisics 23, Robiso, P.M. ad D. Mariucci (200). Narrow-bad aalysis of osaioary processes. Aals of Saisics 29, Robiso, P.M. (2005). The disace bewee rival osaioary processes. Joural of Ecoomerics 28, Solo, V. (992). Irisic radom fucios ad he paradox of /f oise. Siam Joural of Applied Mahemaics 52, Sowell, F.B. (990). The fracioal ui roo disribuio. Ecoomerica 58, Velasco, C. (999). No-saioary log-periodogram regressio. Joural of Ecoomerics 9,

27 Velasco, C. (2000). No-Gaussia log-periodogram regressio. Ecoomeric Theory 6, Velasco, C. (2005). Semiparameric Esimaio of Log-Memory Models. Forhcomig i Palgrave Hadboo of Ecoomerics, Vol.. Ecoomeric Theory, K. Paerso ad T.C. Mills eds, Palgrave, MacMilla. Velasco, C. ad P.M. Robiso (2000). While Pseudo-Maximum Lielihood Esimaes of No- Saioary Time Series. Joural of he America Saisical Associaio 95,

STK4080/9080 Survival and event history analysis

STK4080/9080 Survival and event history analysis STK48/98 Survival ad eve hisory aalysis Marigales i discree ime Cosider a sochasic process The process M is a marigale if Lecure 3: Marigales ad oher sochasic processes i discree ime (recap) where (formally

More information

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS Opimal ear Forecasig Alhough we have o meioed hem explicily so far i he course, here are geeral saisical priciples for derivig he bes liear forecas, ad

More information

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory Liear Time-Ivaria Sysems (LTI Sysems) Oulie Basic Sysem Properies Memoryless ad sysems wih memory (saic or dyamic) Causal ad o-causal sysems (Causaliy) Liear ad o-liear sysems (Lieariy) Sable ad o-sable

More information

A Note on Prediction with Misspecified Models

A Note on Prediction with Misspecified Models ITB J. Sci., Vol. 44 A, No. 3,, 7-9 7 A Noe o Predicio wih Misspecified Models Khresha Syuhada Saisics Research Divisio, Faculy of Mahemaics ad Naural Scieces, Isiu Tekologi Badug, Jala Gaesa Badug, Jawa

More information

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP)

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP) ENGG450 Probabiliy ad Saisics for Egieers Iroducio 3 Probabiliy 4 Probabiliy disribuios 5 Probabiliy Desiies Orgaizaio ad descripio of daa 6 Samplig disribuios 7 Ifereces cocerig a mea 8 Comparig wo reames

More information

1 Notes on Little s Law (l = λw)

1 Notes on Little s Law (l = λw) Copyrigh c 26 by Karl Sigma Noes o Lile s Law (l λw) We cosider here a famous ad very useful law i queueig heory called Lile s Law, also kow as l λw, which assers ha he ime average umber of cusomers i

More information

Section 8 Convolution and Deconvolution

Section 8 Convolution and Deconvolution APPLICATIONS IN SIGNAL PROCESSING Secio 8 Covoluio ad Decovoluio This docume illusraes several echiques for carryig ou covoluio ad decovoluio i Mahcad. There are several operaors available for hese fucios:

More information

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition. ! Revised April 21, 2010 1:27 P! 1 Fourier Series David Radall Assume ha u( x,) is real ad iegrable If he domai is periodic, wih period L, we ca express u( x,) exacly by a Fourier series expasio: ( ) =

More information

Lecture 9: Polynomial Approximations

Lecture 9: Polynomial Approximations CS 70: Complexiy Theory /6/009 Lecure 9: Polyomial Approximaios Isrucor: Dieer va Melkebeek Scribe: Phil Rydzewski & Piramaayagam Arumuga Naiar Las ime, we proved ha o cosa deph circui ca evaluae he pariy

More information

Math 6710, Fall 2016 Final Exam Solutions

Math 6710, Fall 2016 Final Exam Solutions Mah 67, Fall 6 Fial Exam Soluios. Firs, a sude poied ou a suble hig: if P (X i p >, he X + + X (X + + X / ( evaluaes o / wih probabiliy p >. This is roublesome because a radom variable is supposed o be

More information

Lecture 15 First Properties of the Brownian Motion

Lecture 15 First Properties of the Brownian Motion Lecure 15: Firs Properies 1 of 8 Course: Theory of Probabiliy II Term: Sprig 2015 Isrucor: Gorda Zikovic Lecure 15 Firs Properies of he Browia Moio This lecure deals wih some of he more immediae properies

More information

Extremal graph theory II: K t and K t,t

Extremal graph theory II: K t and K t,t Exremal graph heory II: K ad K, Lecure Graph Theory 06 EPFL Frak de Zeeuw I his lecure, we geeralize he wo mai heorems from he las lecure, from riagles K 3 o complee graphs K, ad from squares K, o complee

More information

OLS bias for econometric models with errors-in-variables. The Lucas-critique Supplementary note to Lecture 17

OLS bias for econometric models with errors-in-variables. The Lucas-critique Supplementary note to Lecture 17 OLS bias for ecoomeric models wih errors-i-variables. The Lucas-criique Supplemeary oe o Lecure 7 RNy May 6, 03 Properies of OLS i RE models I Lecure 7 we discussed he followig example of a raioal expecaios

More information

Comparisons Between RV, ARV and WRV

Comparisons Between RV, ARV and WRV Comparisos Bewee RV, ARV ad WRV Cao Gag,Guo Migyua School of Maageme ad Ecoomics, Tiaji Uiversiy, Tiaji,30007 Absrac: Realized Volailiy (RV) have bee widely used sice i was pu forward by Aderso ad Bollerslev

More information

Extended Laguerre Polynomials

Extended Laguerre Polynomials I J Coemp Mah Scieces, Vol 7, 1, o, 189 194 Exeded Laguerre Polyomials Ada Kha Naioal College of Busiess Admiisraio ad Ecoomics Gulberg-III, Lahore, Pakisa adakhaariq@gmailcom G M Habibullah Naioal College

More information

Big O Notation for Time Complexity of Algorithms

Big O Notation for Time Complexity of Algorithms BRONX COMMUNITY COLLEGE of he Ciy Uiversiy of New York DEPARTMENT OF MATHEMATICS AND COMPUTER SCIENCE CSI 33 Secio E01 Hadou 1 Fall 2014 Sepember 3, 2014 Big O Noaio for Time Complexiy of Algorihms Time

More information

Inference of the Second Order Autoregressive. Model with Unit Roots

Inference of the Second Order Autoregressive. Model with Unit Roots Ieraioal Mahemaical Forum Vol. 6 0 o. 5 595-604 Iferece of he Secod Order Auoregressive Model wih Ui Roos Ahmed H. Youssef Professor of Applied Saisics ad Ecoomerics Isiue of Saisical Sudies ad Research

More information

Notes 03 largely plagiarized by %khc

Notes 03 largely plagiarized by %khc 1 1 Discree-Time Covoluio Noes 03 largely plagiarized by %khc Le s begi our discussio of covoluio i discree-ime, sice life is somewha easier i ha domai. We sar wih a sigal x[] ha will be he ipu io our

More information

Using Linnik's Identity to Approximate the Prime Counting Function with the Logarithmic Integral

Using Linnik's Identity to Approximate the Prime Counting Function with the Logarithmic Integral Usig Lii's Ideiy o Approimae he Prime Couig Fucio wih he Logarihmic Iegral Naha McKezie /26/2 aha@icecreambreafas.com Summary:This paper will show ha summig Lii's ideiy from 2 o ad arragig erms i a cerai

More information

Additional Tables of Simulation Results

Additional Tables of Simulation Results Saisica Siica: Suppleme REGULARIZING LASSO: A CONSISTENT VARIABLE SELECTION METHOD Quefeg Li ad Ju Shao Uiversiy of Wiscosi, Madiso, Eas Chia Normal Uiversiy ad Uiversiy of Wiscosi, Madiso Supplemeary

More information

On the Validity of the Pairs Bootstrap for Lasso Estimators

On the Validity of the Pairs Bootstrap for Lasso Estimators O he Validiy of he Pairs Boosrap for Lasso Esimaors Lorezo Campoovo Uiversiy of S.Galle Ocober 2014 Absrac We sudy he validiy of he pairs boosrap for Lasso esimaors i liear regressio models wih radom covariaes

More information

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods"

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods Suppleme for SADAGRAD: Srogly Adapive Sochasic Gradie Mehods" Zaiyi Che * 1 Yi Xu * Ehog Che 1 iabao Yag 1. Proof of Proposiio 1 Proposiio 1. Le ɛ > 0 be fixed, H 0 γi, γ g, EF (w 1 ) F (w ) ɛ 0 ad ieraio

More information

Comparison between Fourier and Corrected Fourier Series Methods

Comparison between Fourier and Corrected Fourier Series Methods Malaysia Joural of Mahemaical Scieces 7(): 73-8 (13) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Joural homepage: hp://eispem.upm.edu.my/oural Compariso bewee Fourier ad Correced Fourier Series Mehods 1

More information

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS PB Sci Bull, Series A, Vol 78, Iss 4, 2016 ISSN 1223-7027 CLOSED FORM EVALATION OF RESTRICTED SMS CONTAINING SQARES OF FIBONOMIAL COEFFICIENTS Emrah Kılıc 1, Helmu Prodiger 2 We give a sysemaic approach

More information

COS 522: Complexity Theory : Boaz Barak Handout 10: Parallel Repetition Lemma

COS 522: Complexity Theory : Boaz Barak Handout 10: Parallel Repetition Lemma COS 522: Complexiy Theory : Boaz Barak Hadou 0: Parallel Repeiio Lemma Readig: () A Parallel Repeiio Theorem / Ra Raz (available o his websie) (2) Parallel Repeiio: Simplificaios ad he No-Sigallig Case

More information

The Central Limit Theorem

The Central Limit Theorem The Ceral Limi Theorem The ceral i heorem is oe of he mos impora heorems i probabiliy heory. While here a variey of forms of he ceral i heorem, he mos geeral form saes ha give a sufficiely large umber,

More information

K3 p K2 p Kp 0 p 2 p 3 p

K3 p K2 p Kp 0 p 2 p 3 p Mah 80-00 Mo Ar 0 Chaer 9 Fourier Series ad alicaios o differeial equaios (ad arial differeial equaios) 9.-9. Fourier series defiiio ad covergece. The idea of Fourier series is relaed o he liear algebra

More information

CSE 241 Algorithms and Data Structures 10/14/2015. Skip Lists

CSE 241 Algorithms and Data Structures 10/14/2015. Skip Lists CSE 41 Algorihms ad Daa Srucures 10/14/015 Skip Liss This hadou gives he skip lis mehods ha we discussed i class. A skip lis is a ordered, doublyliked lis wih some exra poiers ha allow us o jump over muliple

More information

12 Getting Started With Fourier Analysis

12 Getting Started With Fourier Analysis Commuicaios Egieerig MSc - Prelimiary Readig Geig Sared Wih Fourier Aalysis Fourier aalysis is cocered wih he represeaio of sigals i erms of he sums of sie, cosie or complex oscillaio waveforms. We ll

More information

Mathematical Statistics. 1 Introduction to the materials to be covered in this course

Mathematical Statistics. 1 Introduction to the materials to be covered in this course Mahemaical Saisics Iroducio o he maerials o be covered i his course. Uivariae & Mulivariae r.v s 2. Borl-Caelli Lemma Large Deviaios. e.g. X,, X are iid r.v s, P ( X + + X where I(A) is a umber depedig

More information

A Note on Random k-sat for Moderately Growing k

A Note on Random k-sat for Moderately Growing k A Noe o Radom k-sat for Moderaely Growig k Ju Liu LMIB ad School of Mahemaics ad Sysems Sciece, Beihag Uiversiy, Beijig, 100191, P.R. Chia juliu@smss.buaa.edu.c Zogsheg Gao LMIB ad School of Mahemaics

More information

Stochastic Processes Adopted From p Chapter 9 Probability, Random Variables and Stochastic Processes, 4th Edition A. Papoulis and S.

Stochastic Processes Adopted From p Chapter 9 Probability, Random Variables and Stochastic Processes, 4th Edition A. Papoulis and S. Sochasic Processes Adoped From p Chaper 9 Probabiliy, adom Variables ad Sochasic Processes, 4h Ediio A. Papoulis ad S. Pillai 9. Sochasic Processes Iroducio Le deoe he radom oucome of a experime. To every

More information

Dynamic h-index: the Hirsch index in function of time

Dynamic h-index: the Hirsch index in function of time Dyamic h-idex: he Hirsch idex i fucio of ime by L. Egghe Uiversiei Hassel (UHassel), Campus Diepebeek, Agoralaa, B-3590 Diepebeek, Belgium ad Uiversiei Awerpe (UA), Campus Drie Eike, Uiversieisplei, B-260

More information

Stationarity and Unit Root tests

Stationarity and Unit Root tests Saioari ad Ui Roo ess Saioari ad Ui Roo ess. Saioar ad Nosaioar Series. Sprios Regressio 3. Ui Roo ad Nosaioari 4. Ui Roo ess Dicke-Fller es Agmeed Dicke-Fller es KPSS es Phillips-Perro Tes 5. Resolvig

More information

Moment Generating Function

Moment Generating Function 1 Mome Geeraig Fucio m h mome m m m E[ ] x f ( x) dx m h ceral mome m m m E[( ) ] ( ) ( x ) f ( x) dx Mome Geeraig Fucio For a real, M () E[ e ] e k x k e p ( x ) discree x k e f ( x) dx coiuous Example

More information

BRIDGE ESTIMATOR AS AN ALTERNATIVE TO DICKEY- PANTULA UNIT ROOT TEST

BRIDGE ESTIMATOR AS AN ALTERNATIVE TO DICKEY- PANTULA UNIT ROOT TEST The 0 h Ieraioal Days of Saisics ad Ecoomics Prague Sepember 8-0 06 BRIDGE ESTIMATOR AS AN ALTERNATIVE TO DICKEY- PANTULA UNIT ROOT TEST Hüseyi Güler Yeliz Yalҫi Çiğdem Koşar Absrac Ecoomic series may

More information

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER Maerials Physics ad Mechaics 3 (5) 36-4 Received: March 7 5 METHOD OF THE EQUIVAENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBEM FOR EASTIC DIFFUSION AYER A.V. Zemsov * D.V. Tarlaovsiy Moscow Aviaio Isiue

More information

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY U.P.B. Sci. Bull., Series A, Vol. 78, Iss. 2, 206 ISSN 223-7027 A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY İbrahim Çaak I his paper we obai a Tauberia codiio i erms of he weighed classical

More information

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi Liear lgebra Lecure #9 Noes This week s lecure focuses o wha migh be called he srucural aalysis of liear rasformaios Wha are he irisic properies of a liear rasformaio? re here ay fixed direcios? The discussio

More information

L-functions and Class Numbers

L-functions and Class Numbers L-fucios ad Class Numbers Sude Number Theory Semiar S. M.-C. 4 Sepember 05 We follow Romyar Sharifi s Noes o Iwasawa Theory, wih some help from Neukirch s Algebraic Number Theory. L-fucios of Dirichle

More information

SUMMATION OF INFINITE SERIES REVISITED

SUMMATION OF INFINITE SERIES REVISITED SUMMATION OF INFINITE SERIES REVISITED I several aricles over he las decade o his web page we have show how o sum cerai iiie series icludig he geomeric series. We wa here o eed his discussio o he geeral

More information

N! AND THE GAMMA FUNCTION

N! AND THE GAMMA FUNCTION N! AND THE GAMMA FUNCTION Cosider he produc of he firs posiive iegers- 3 4 5 6 (-) =! Oe calls his produc he facorial ad has ha produc of he firs five iegers equals 5!=0. Direcly relaed o he discree! fucio

More information

The Moment Approximation of the First Passage Time for the Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier

The Moment Approximation of the First Passage Time for the Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier America Joural of Applied Mahemaics ad Saisics, 015, Vol. 3, No. 5, 184-189 Available olie a hp://pubs.sciepub.com/ajams/3/5/ Sciece ad Educaio Publishig DOI:10.1691/ajams-3-5- The Mome Approximaio of

More information

Solutions to selected problems from the midterm exam Math 222 Winter 2015

Solutions to selected problems from the midterm exam Math 222 Winter 2015 Soluios o seleced problems from he miderm eam Mah Wier 5. Derive he Maclauri series for he followig fucios. (cf. Pracice Problem 4 log( + (a L( d. Soluio: We have he Maclauri series log( + + 3 3 4 4 +...,

More information

F D D D D F. smoothed value of the data including Y t the most recent data.

F D D D D F. smoothed value of the data including Y t the most recent data. Module 2 Forecasig 1. Wha is forecasig? Forecasig is defied as esimaig he fuure value ha a parameer will ake. Mos scieific forecasig mehods forecas he fuure value usig pas daa. I Operaios Maageme forecasig

More information

Skewness of Gaussian Mixture Absolute Value GARCH(1, 1) Model

Skewness of Gaussian Mixture Absolute Value GARCH(1, 1) Model Commuicaios for Saisical Applicaios ad Mehods 203, Vol. 20, No. 5, 395 404 DOI: hp://dx.doi.org/0.535/csam.203.20.5.395 Skewess of Gaussia Mixure Absolue Value GARCH(, Model Taewook Lee,a a Deparme of

More information

Time Series, Part 1 Content Literature

Time Series, Part 1 Content Literature Time Series, Par Coe - Saioariy, auocorrelaio, parial auocorrelaio, removal of osaioary compoes, idepedece es for ime series - Liear Sochasic Processes: auoregressive (AR), movig average (MA), auoregressive

More information

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x) 1 Trasform Techiques h m m m m mome : E[ ] x f ( x) dx h m m m m ceral mome : E[( ) ] ( ) ( x) f ( x) dx A coveie wa of fidig he momes of a radom variable is he mome geeraig fucio (MGF). Oher rasform echiques

More information

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws. Limi of a fucio.. Oe-Sided..... Ifiie limis Verical Asympoes... Calculaig Usig he Limi Laws.5 The Squeeze Theorem.6 The Precise Defiiio of a Limi......7 Coiuiy.8 Iermediae Value Theorem..9 Refereces..

More information

Solution. 1 Solutions of Homework 6. Sangchul Lee. April 28, Problem 1.1 [Dur10, Exercise ]

Solution. 1 Solutions of Homework 6. Sangchul Lee. April 28, Problem 1.1 [Dur10, Exercise ] Soluio Sagchul Lee April 28, 28 Soluios of Homework 6 Problem. [Dur, Exercise 2.3.2] Le A be a sequece of idepede eves wih PA < for all. Show ha P A = implies PA i.o. =. Proof. Noice ha = P A c = P A c

More information

B. Maddah INDE 504 Simulation 09/02/17

B. Maddah INDE 504 Simulation 09/02/17 B. Maddah INDE 54 Simulaio 9/2/7 Queueig Primer Wha is a queueig sysem? A queueig sysem cosiss of servers (resources) ha provide service o cusomers (eiies). A Cusomer requesig service will sar service

More information

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract A ieresig resul abou subse sums Niu Kichloo Lior Pacher November 27, 1993 Absrac We cosider he problem of deermiig he umber of subses B f1; 2; : : :; g such ha P b2b b k mod, where k is a residue class

More information

Actuarial Society of India

Actuarial Society of India Acuarial Sociey of Idia EXAMINAIONS Jue 5 C4 (3) Models oal Marks - 5 Idicaive Soluio Q. (i) a) Le U deoe he process described by 3 ad V deoe he process described by 4. he 5 e 5 PU [ ] PV [ ] ( e ).538!

More information

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form Joural of Applied Mahemaics Volume 03, Aricle ID 47585, 7 pages hp://dx.doi.org/0.55/03/47585 Research Aricle A Geeralized Noliear Sum-Differece Iequaliy of Produc Form YogZhou Qi ad Wu-Sheg Wag School

More information

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions.

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions. Iera. J. Mah. & Mah. Si. Vol. 6 No. 3 (1983) 559-566 559 ASYMPTOTIC RELATIOHIPS BETWEEN TWO HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS TAKA KUSANO laculy of Sciece Hrosh llersy 1982) ABSTRACT. Some asympoic

More information

Affine term structure models

Affine term structure models /5/07 Affie erm srucure models A. Iro o Gaussia affie erm srucure models B. Esimaio by miimum chi square (Hamilo ad Wu) C. Esimaio by OLS (Adria, Moech, ad Crump) D. Dyamic Nelso-Siegel model (Chrisese,

More information

A Generalized Cost Malmquist Index to the Productivities of Units with Negative Data in DEA

A Generalized Cost Malmquist Index to the Productivities of Units with Negative Data in DEA Proceedigs of he 202 Ieraioal Coferece o Idusrial Egieerig ad Operaios Maageme Isabul, urey, July 3 6, 202 A eeralized Cos Malmquis Ide o he Produciviies of Uis wih Negaive Daa i DEA Shabam Razavya Deparme

More information

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs America Joural of Compuaioal Mahemaics, 04, 4, 80-88 Published Olie Sepember 04 i SciRes. hp://www.scirp.org/joural/ajcm hp://dx.doi.org/0.436/ajcm.04.4404 Mea Square Coverge Fiie Differece Scheme for

More information

Exercise 3 Stochastic Models of Manufacturing Systems 4T400, 6 May

Exercise 3 Stochastic Models of Manufacturing Systems 4T400, 6 May Exercise 3 Sochasic Models of Maufacurig Sysems 4T4, 6 May. Each week a very popular loery i Adorra pris 4 ickes. Each ickes has wo 4-digi umbers o i, oe visible ad he oher covered. The umbers are radomly

More information

The analysis of the method on the one variable function s limit Ke Wu

The analysis of the method on the one variable function s limit Ke Wu Ieraioal Coferece o Advaces i Mechaical Egieerig ad Idusrial Iformaics (AMEII 5) The aalysis of he mehod o he oe variable fucio s i Ke Wu Deparme of Mahemaics ad Saisics Zaozhuag Uiversiy Zaozhuag 776

More information

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3 Ieraioal Joural of Saisics ad Aalysis. ISSN 48-9959 Volume 6, Number (6, pp. -8 Research Idia Publicaios hp://www.ripublicaio.com The Populaio Mea ad is Variace i he Presece of Geocide for a Simple Birh-Deah-

More information

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM Available a hp://pvamu.edu/aam Appl. Appl. Mah. ISSN: 193-9466 Vol. 5, No. Issue (December 1), pp. 575 584 (Previously, Vol. 5, Issue 1, pp. 167 1681) Applicaios ad Applied Mahemaics: A Ieraioal Joural

More information

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE Mohia & Samaa, Vol. 1, No. II, December, 016, pp 34-49. ORIGINAL RESEARCH ARTICLE OPEN ACCESS FIED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE 1 Mohia S. *, Samaa T. K. 1 Deparme of Mahemaics, Sudhir Memorial

More information

Order Determination for Multivariate Autoregressive Processes Using Resampling Methods

Order Determination for Multivariate Autoregressive Processes Using Resampling Methods joural of mulivariae aalysis 57, 175190 (1996) aricle o. 0028 Order Deermiaio for Mulivariae Auoregressive Processes Usig Resamplig Mehods Chaghua Che ad Richard A. Davis* Colorado Sae Uiversiy ad Peer

More information

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4) 7 Differeial equaios Review Solve by he mehod of udeermied coefficies ad by he mehod of variaio of parameers (4) y y = si Soluio; we firs solve he homogeeous equaio (4) y y = 4 The correspodig characerisic

More information

Approximately Quasi Inner Generalized Dynamics on Modules. { } t t R

Approximately Quasi Inner Generalized Dynamics on Modules. { } t t R Joural of Scieces, Islamic epublic of Ira 23(3): 245-25 (22) Uiversiy of Tehra, ISSN 6-4 hp://jscieces.u.ac.ir Approximaely Quasi Ier Geeralized Dyamics o Modules M. Mosadeq, M. Hassai, ad A. Nikam Deparme

More information

Solutions to Problems 3, Level 4

Solutions to Problems 3, Level 4 Soluios o Problems 3, Level 4 23 Improve he resul of Quesio 3 whe l. i Use log log o prove ha for real >, log ( {}log + 2 d log+ P ( + P ( d 2. Here P ( is defied i Quesio, ad parial iegraio has bee used.

More information

State and Parameter Estimation of The Lorenz System In Existence of Colored Noise

State and Parameter Estimation of The Lorenz System In Existence of Colored Noise Sae ad Parameer Esimaio of he Lorez Sysem I Eisece of Colored Noise Mozhga Mombeii a Hamid Khaloozadeh b a Elecrical Corol ad Sysem Egieerig Researcher of Isiue for Research i Fudameal Scieces (IPM ehra

More information

King Fahd University of Petroleum & Minerals Computer Engineering g Dept

King Fahd University of Petroleum & Minerals Computer Engineering g Dept Kig Fahd Uiversiy of Peroleum & Mierals Compuer Egieerig g Dep COE 4 Daa ad Compuer Commuicaios erm Dr. shraf S. Hasa Mahmoud Rm -4 Ex. 74 Email: ashraf@kfupm.edu.sa 9/8/ Dr. shraf S. Hasa Mahmoud Lecure

More information

A note on deviation inequalities on {0, 1} n. by Julio Bernués*

A note on deviation inequalities on {0, 1} n. by Julio Bernués* A oe o deviaio iequaliies o {0, 1}. by Julio Berués* Deparameo de Maemáicas. Faculad de Ciecias Uiversidad de Zaragoza 50009-Zaragoza (Spai) I. Iroducio. Le f: (Ω, Σ, ) IR be a radom variable. Roughly

More information

Online Supplement to Reactive Tabu Search in a Team-Learning Problem

Online Supplement to Reactive Tabu Search in a Team-Learning Problem Olie Suppleme o Reacive abu Search i a eam-learig Problem Yueli She School of Ieraioal Busiess Admiisraio, Shaghai Uiversiy of Fiace ad Ecoomics, Shaghai 00433, People s Republic of Chia, she.yueli@mail.shufe.edu.c

More information

Review Exercises for Chapter 9

Review Exercises for Chapter 9 0_090R.qd //0 : PM Page 88 88 CHAPTER 9 Ifiie Series I Eercises ad, wrie a epressio for he h erm of he sequece..,., 5, 0,,,, 0,... 7,... I Eercises, mach he sequece wih is graph. [The graphs are labeled

More information

Four equations describe the dynamic solution to RBC model. Consumption-leisure efficiency condition. Consumption-investment efficiency condition

Four equations describe the dynamic solution to RBC model. Consumption-leisure efficiency condition. Consumption-investment efficiency condition LINEARIZING AND APPROXIMATING THE RBC MODEL SEPTEMBER 7, 200 For f( x, y, z ), mulivariable Taylor liear expasio aroud ( x, yz, ) f ( x, y, z) f( x, y, z) + f ( x, y, z)( x x) + f ( x, y, z)( y y) + f

More information

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003 ODEs II, Suppleme o Lecures 6 & 7: The Jorda Normal Form: Solvig Auoomous, Homogeeous Liear Sysems April 2, 23 I his oe, we describe he Jorda ormal form of a marix ad use i o solve a geeral homogeeous

More information

On The Eneström-Kakeya Theorem

On The Eneström-Kakeya Theorem Applied Mahemaics,, 3, 555-56 doi:436/am673 Published Olie December (hp://wwwscirporg/oural/am) O The Eesröm-Kakeya Theorem Absrac Gulsha Sigh, Wali Mohammad Shah Bharahiar Uiversiy, Coimbaore, Idia Deparme

More information

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL JANUARY 29, 2013

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL JANUARY 29, 2013 LINEAR APPROXIMATION OF THE BASELINE RBC MODEL JANUARY 29, 203 Iroducio LINEARIZATION OF THE RBC MODEL For f( x, y, z ) = 0, mulivariable Taylor liear expasio aroud f( x, y, z) f( x, y, z) + f ( x, y,

More information

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation.

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation. ecure Phys 375 Eergy Desiy / Eergy Flu / oal Eergy i D Overview ad Moivaio: Fro your sudy of waves i iroducory physics you should be aware ha waves ca raspor eergy fro oe place o aoher cosider he geeraio

More information

Semiparametric and Nonparametric Methods in Political Science Lecture 1: Semiparametric Estimation

Semiparametric and Nonparametric Methods in Political Science Lecture 1: Semiparametric Estimation Semiparameric ad Noparameric Mehods i Poliical Sciece Lecure : Semiparameric Esimaio Michael Peress, Uiversiy of Rocheser ad Yale Uiversiy Lecure : Semiparameric Mehods Page 2 Overview of Semi ad Noparameric

More information

Additive nonparametric models with time variable and both stationary and nonstationary regressions

Additive nonparametric models with time variable and both stationary and nonstationary regressions Addiive oparameric models wih ime variable ad boh saioary ad osaioary regressios Chaohua Dog Oliver Lio The Isiue for Fiscal Sudies Deparme of Ecoomics, UCL cemmap workig paper CWP59/7 Addiive oparameric

More information

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are: 3. Iiial value problems: umerical soluio Fiie differeces - Trucaio errors, cosisecy, sabiliy ad covergece Crieria for compuaioal sabiliy Explici ad implici ime schemes Table of ime schemes Hyperbolic ad

More information

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION Malaysia Joural of Mahemaical Scieces 2(2): 55-6 (28) The Soluio of he Oe Species Loka-Volerra Equaio Usig Variaioal Ieraio Mehod B. Baiha, M.S.M. Noorai, I. Hashim School of Mahemaical Scieces, Uiversii

More information

HYBRID STOCHASTIC LOCAL UNIT ROOTS. Offer Lieberman and Peter C. B. Phillips. November 2017 COWLES FOUNDATION DISCUSSION PAPER NO.

HYBRID STOCHASTIC LOCAL UNIT ROOTS. Offer Lieberman and Peter C. B. Phillips. November 2017 COWLES FOUNDATION DISCUSSION PAPER NO. HYBRID STOCHASTIC LOCAL UNIT ROOTS By Offer Lieberma ad Peer C. B. Phillips November 7 COWLES FOUNDATION DISCUSSION PAPER NO. 33 COWLES FOUNDATION FOR RESEARCH IN ECONOMICS YALE UNIVERSITY Box 88 New Have,

More information

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL SEPTEMBER 17, 2013

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL SEPTEMBER 17, 2013 LINEAR APPROXIMATION OF THE BASELINE RBC MODEL SEPTEMBER 7, 203 Iroducio LINEARIZATION OF THE RBC MODEL For f( xyz,, ) = 0, mulivariable Taylor liear expasio aroud f( xyz,, ) f( xyz,, ) + f( xyz,, )( x

More information

Economics 8723 Macroeconomic Theory Problem Set 2 Professor Sanjay Chugh Spring 2017

Economics 8723 Macroeconomic Theory Problem Set 2 Professor Sanjay Chugh Spring 2017 Deparme of Ecoomics The Ohio Sae Uiversiy Ecoomics 8723 Macroecoomic Theory Problem Se 2 Professor Sajay Chugh Sprig 207 Labor Icome Taxes, Nash-Bargaied Wages, ad Proporioally-Bargaied Wages. I a ecoomy

More information

Applying the Moment Generating Functions to the Study of Probability Distributions

Applying the Moment Generating Functions to the Study of Probability Distributions 3 Iformaica Ecoomică, r (4)/007 Applyi he Mome Geerai Fucios o he Sudy of Probabiliy Disribuios Silvia SPĂTARU Academy of Ecoomic Sudies, Buchares I his paper, we describe a ool o aid i provi heorems abou

More information

C(p, ) 13 N. Nuclear reactions generate energy create new isotopes and elements. Notation for stellar rates: p 12

C(p, ) 13 N. Nuclear reactions generate energy create new isotopes and elements. Notation for stellar rates: p 12 Iroducio o sellar reacio raes Nuclear reacios geerae eergy creae ew isoopes ad elemes Noaio for sellar raes: p C 3 N C(p,) 3 N The heavier arge ucleus (Lab: arge) he ligher icomig projecile (Lab: beam)

More information

The Eigen Function of Linear Systems

The Eigen Function of Linear Systems 1/25/211 The Eige Fucio of Liear Sysems.doc 1/7 The Eige Fucio of Liear Sysems Recall ha ha we ca express (expad) a ime-limied sigal wih a weighed summaio of basis fucios: v ( ) a ψ ( ) = where v ( ) =

More information

Clock Skew and Signal Representation

Clock Skew and Signal Representation Clock Skew ad Sigal Represeaio Ch. 7 IBM Power 4 Chip 0/7/004 08 frequecy domai Program Iroducio ad moivaio Sequeial circuis, clock imig, Basic ools for frequecy domai aalysis Fourier series sigal represeaio

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 4 9/16/2013. Applications of the large deviation technique

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 4 9/16/2013. Applications of the large deviation technique MASSACHUSETTS ISTITUTE OF TECHOLOGY 6.265/5.070J Fall 203 Lecure 4 9/6/203 Applicaios of he large deviaio echique Coe.. Isurace problem 2. Queueig problem 3. Buffer overflow probabiliy Safey capial for

More information

Asymptotic statistics for multilayer perceptron with ReLu hidden units

Asymptotic statistics for multilayer perceptron with ReLu hidden units ESANN 8 proceedigs, Europea Symposium o Arificial Neural Neworks, Compuaioal Ielligece ad Machie Learig. Bruges (Belgium), 5-7 April 8, i6doc.com publ., ISBN 978-8758747-6. Available from hp://www.i6doc.com/e/.

More information

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme America Joural of Compuaioal ad Applied Maemaics, (6): 77-8 DOI:.59/.acam.6. Numerical Soluio of Parabolic Volerra Iegro-Differeial Equaios via Bacward-Euler Sceme Ali Filiz Deparme of Maemaics, Ada Mederes

More information

Available online at J. Math. Comput. Sci. 4 (2014), No. 4, ISSN:

Available online at   J. Math. Comput. Sci. 4 (2014), No. 4, ISSN: Available olie a hp://sci.org J. Mah. Compu. Sci. 4 (2014), No. 4, 716-727 ISSN: 1927-5307 ON ITERATIVE TECHNIQUES FOR NUMERICAL SOLUTIONS OF LINEAR AND NONLINEAR DIFFERENTIAL EQUATIONS S.O. EDEKI *, A.A.

More information

Research Article Testing for Change in Mean of Independent Multivariate Observations with Time Varying Covariance

Research Article Testing for Change in Mean of Independent Multivariate Observations with Time Varying Covariance Joural of Probabiliy ad Saisics Volume, Aricle ID 969753, 7 pages doi:.55//969753 Research Aricle Tesig for Chage i Mea of Idepede Mulivariae Observaios wih Time Varyig Covariace Mohamed Bouahar Isiue

More information

Supplementary Information for Thermal Noises in an Aqueous Quadrupole Micro- and Nano-Trap

Supplementary Information for Thermal Noises in an Aqueous Quadrupole Micro- and Nano-Trap Supplemeary Iformaio for Thermal Noises i a Aqueous Quadrupole Micro- ad Nao-Trap Jae Hyu Park ad Predrag S. Krsić * Physics Divisio, Oak Ridge Naioal Laboraory, Oak Ridge, TN 3783 E-mail: krsicp@orl.gov

More information

Math 2414 Homework Set 7 Solutions 10 Points

Math 2414 Homework Set 7 Solutions 10 Points Mah Homework Se 7 Soluios 0 Pois #. ( ps) Firs verify ha we ca use he iegral es. The erms are clearly posiive (he epoeial is always posiive ad + is posiive if >, which i is i his case). For decreasig we

More information

Fermat Numbers in Multinomial Coefficients

Fermat Numbers in Multinomial Coefficients 1 3 47 6 3 11 Joural of Ieger Sequeces, Vol. 17 (014, Aricle 14.3. Ferma Numbers i Muliomial Coefficies Shae Cher Deparme of Mahemaics Zhejiag Uiversiy Hagzhou, 31007 Chia chexiaohag9@gmail.com Absrac

More information

Calculus BC 2015 Scoring Guidelines

Calculus BC 2015 Scoring Guidelines AP Calculus BC 5 Scorig Guidelies 5 The College Board. College Board, Advaced Placeme Program, AP, AP Ceral, ad he acor logo are regisered rademarks of he College Board. AP Ceral is he official olie home

More information

Procedia - Social and Behavioral Sciences 230 ( 2016 ) Joint Probability Distribution and the Minimum of a Set of Normalized Random Variables

Procedia - Social and Behavioral Sciences 230 ( 2016 ) Joint Probability Distribution and the Minimum of a Set of Normalized Random Variables Available olie a wwwsciecedireccom ScieceDirec Procedia - Social ad Behavioral Scieces 30 ( 016 ) 35 39 3 rd Ieraioal Coferece o New Challeges i Maageme ad Orgaizaio: Orgaizaio ad Leadership, May 016,

More information

Introduction to the Mathematics of Lévy Processes

Introduction to the Mathematics of Lévy Processes Iroducio o he Mahemaics of Lévy Processes Kazuhisa Masuda Deparme of Ecoomics The Graduae Ceer, The Ciy Uiversiy of New York, 365 Fifh Aveue, New York, NY 10016-4309 Email: maxmasuda@maxmasudacom hp://wwwmaxmasudacom/

More information

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming* The Eighh Ieraioal Symposium o Operaios esearch ad Is Applicaios (ISOA 9) Zhagjiajie Chia Sepember 2 22 29 Copyrigh 29 OSC & APOC pp 33 39 Some Properies of Semi-E-Covex Fucio ad Semi-E-Covex Programmig*

More information