THE SMALL SAMPLE PROPERTIES OF TESTS OF THE EXPECTATIONS HYPOTHESIS: A MONTE CARLO INVESTIGATION E. GARGANAS, S.G.HALL

Size: px
Start display at page:

Download "THE SMALL SAMPLE PROPERTIES OF TESTS OF THE EXPECTATIONS HYPOTHESIS: A MONTE CARLO INVESTIGATION E. GARGANAS, S.G.HALL"

Transcription

1 ISSN THE SMALL SAMPLE PROPERTIES OF TESTS OF THE EXPECTATIONS HYPOTHESIS: A MONTE CARLO INVESTIGATION E. GARGANAS, S.G.HALL Taaka Busiess School Discussio Papers: TBS/DP04/6 Lodo: Taaka Busiess School,

2 The small sample properies of ess of he expecaios hypohesis: A Moe Carlo ivesigaio By E. Gargaas (Naioal Bak of Greece), S.G. Hall (Imperial College) ABSTRACT I his paper, we exed resuls from he fiace lieraure ha explores small-sample bias, due o persise variables, i ess of prese value asse pricig models. Usig Moe Carlo experimes we aalyze he fiie sample behaviour of a variey of ess of he expecaios hypohesis of he erm srucure, whe ieres raes are icreasigly persise. Resuls overwhelmigly idicae ha asympoic iferece ca be misleadig. I he Campbel ad Shiller (1991) vecor auoregressio (VAR) framework, he Wald es overrejecs he correc ull ad he slope of he correspodig ordiary leas squares (OLS) es is srogly biased. The disorios ca be exreme for recursive esimaors, such as, he Kalma filer ime varyig parameer model. Ocober 2004 JEL classificaio: C15; E43; G12; Keywords: Ieres raes; Expecaios hypohesis; Small sample bias; Kalma filer 2

3 1. Iroducio The icoclusive empirical evidece agais he expecaios hypohesis of he erm srucure of ieres raes has ofe led researchers o sugges ha he heory remais popular primarily o a priori grouds, as i sill seems o provide he mos iuiive explaaio for he erm srucure. I is pure form he expecaios hypohesis (EH) saes ha he log rae is a forecas of expeced fuure shor raes, assumig zero or a ime-ivaria risk premium ad raioal expecaios. A he same ime, hese empirical puzzles have brough abou he developme of a small variey of ieresig explaaios for he heory s rejecio. A recurrig issue i he rece empirical lieraure, has bee he adopio of a formal mehodology for excludig he uforuae case i which he EH is merely rejeced as a resul of saisical or ecoomeric disorios i small samples (ad correcig for hem). The paricular case of srog small sample bias due o persise ieres raes is a ypical problem ha applies o may markes, icludig ha of he U.S. Geerally, here are hree broad explaaios for empirical failures of he EH. Firs, he heory simply is rue. The key behavioural assumpios do o hold e.g. ages are irraioal or risk averse. This has ured a relaively small par of research away from US ieres raes for which he EH is usually rejeced (e.g. Campbel ad Shiller, 1991; Bekaer e al., 1997) owards oher, maily Europea, marke raes for several of which he EH is o rejeced (e.g., Hardouvelis, 1994; Gerlach ad Smes, 1997; Bekaer ad Hodrick, 2000). Secod, he EH ca be saved by icorporaig ime-varyig risk premia o he model specificaio (GARCH models) (Egle, e al., 1987). Third, he saisical ess hemselves lead o false rejecios due o heir poor small sample properies. This leads o a family of explaaios icludig peso problems wih regime shifs (e.g. Bekaer e al., 2001), learig behaviour ad, fially, persise variables (Makiw ad Shapiro, 1986; Bekaer e al., 1997). The laer is he explaaio ha we se ou o explore i his paper, usig arificial raher ha real world ieres raes, i a Moe Carlo simulaio evirome. 3

4 I a iflueial paper, Bekaer e al. (1997) show ha he small sample properies of sadard ess of he EH ca be surprisigly wised whe ieres raes are highly persise. I heir resuls, he wo well kow Campbel ad Shiller (1991) spread regressio specificaios ad he wo VAR mehodology regressios, all provide srogly biased coefficies ha deped o he persisece parameer of a AR(1) oe-period rae, eve for samples of 2000 observaios. The high dispersio of he small sample disribuio of he slope coefficie, makes i more difficul o rejec he EH i small samples. They coclude ha small sample iferece usig carefully desiged Moe Carlo simulaio experimes has esseially become a sie qua o for esig he EH i moder fiace research (e.g. Logsaf, 2000). I his paper, we coduc a series of Moe Carlo experimes o examie he small sample properies of he perfec foresigh spread regressio es, he variace bouds iequaliy ad four VAR mehodology saisics (Campbel ad Shiller, 1991), all i respose o a icreasig persisece parameer. Give ha he lieraure is o always coclusive o he properies of hese small sample disribuios, depedig o he daa geeraig assumpios (or he heoreical approximaios i closed form esimaes), we cofirm ha for mos of hese ess high persisece leads o a bias owards overrejecio of he EH whe asympoic iferece is used. The surprisig excepio is he variace bouds raio for which he EH is uderrejeced. We also fid a paricularly srog bias i he case of he Wald es of he VAR coefficie resricios ad i he slope coefficie of he regressio of he (VAR-cosruced) heoreical agais he acual spread. Cosiderig he furher issue of ieres rae mauriy horizo we fid ha i maers sigificaly, for all ess. Overall, our Moe Carlo resuls idicae ha he saisical size disorios are larger for some ess ha for ohers. This would help explai why i is ofe he case ha i acual markes differe ess may be more favourable or less favourable o he EH. The paper is orgaised as follows. Secio 2 briefly oulies ess of he EH model. Secio 3 explais why we may expec small sample bias o be paricularly srog for mos of hese ess whe ieres raes are highly persise. We discuss he lieraure fidigs ad focus o he properies of he slope coefficie i he Campbel ad Shiller perfec foresigh spread regressio. Secio 4 explores Moe Carlo evidece regardig he small sample behaviour of 6 sadard EH ess ad a Kalma filer model. A sesiiviy aalysis is performed wih respec o a chagig persisece coefficie i a AR(1) daa geeraig process for he shor rae. The las secio he summarises he key fidigs wih some closig remarks. 4

5 2. Empirical ess of he EH 2.1 The perfec foresigh spread regressio ad volailiy ess The mai erm srucure relaioship implied by he EH usig coiuously compouded ieres raes ad a logarihmic approximaio for pure discou bods is: R = 1 1 Εr + c i = 0 (1) where he shor rae, wih oe period o mauriy, is deoed by r, he log rae wih (>1) periods o mauriy is represeed by R ad c is a cosa erm premium. Equaio (1) he saes ha he log rae, a cosa, hrough ime, erm premium, c. R, is a weighed average of expeced fuure shor ieres raes, r, plus I ca be show ha he fudameal EH equaio (1) may be rasformed hrough subracio of he shor rae, r, from boh sides, io he followig expressio (seig he erm premium c o zero, for simpliciy): i 1 (2) R r = 1 Ε r + i i= 1 where S = R r is he acual spread bewee he log ad he shor rae. Campbell ad Shiller * S = 1 i r + i = 1 (1991) defie he variable ( ) i 1 i expressio (2) as he spread prediced by he model if ages have perfec foresigh. The perfec foresigh spread (PFS) is hus cosruced as a weighed average of chages i fuure shor raes, usig ex-pos realisaios of ieres raes from he sample. Expressio (2) he saes ha, accordig o he EH, he acual spread, S, is he expeced PFS, Ε S *. This formulaio (2) is ormally adoped whe ieres raes are believed o be o-saioary I(1) processes. 5

6 Addig o eiher of hese expressios ((1),(2)) he raioal marke expecaios (RE) assumpio ha he shor rae is equal o is forecas plus a whie - oise error erm: r = Ε r + η (3) has led o a umber of ess of he EH proposed i he lieraure. Oe of he OLS regressio ess suggesed by Campbell ad Shiller (1991) ha we cosider here is: S * = α + βs + ε (4) We expec β = 1 ad α = 0 uder he ull of he joi pure expecaios hypohesis (EH wih zero erm premium c) wih he RE assumpio (hereby simply referred o as EH). The error erm is serially correlaed ad follows a MA(-2) process ad a GMM correcio o he variace covariace marix of he coefficies ca be used. We also cosider he correspodig variace boud iequaliy, reformulaed by Campbel ad Shiller (1991) i erms of S ad S * as: Var ( S ) Var( S ) (5) The variace (or sadard deviaio) of he ieres rae spread may o exceed ha of he PFS whe he EH is rue. Equivalely, he raio of he correspodig sadard deviaios (or variaces) mus be less ha 1. I is geerally expeced ha he wo ess (regressio (4) ad variace boud (5)) mus give similar ifereces. 2.2 The Campbell - Shiller VAR mehodology Campbell ad Shiller (1991) show ha he fudameal EH equaio (1) may sill be esed by replacig he RE assumpio wih a auoregressive scheme for predicig he perfec foresigh spread (i expressio (2)), as log as he acual spread S is icluded i he forecasig variables. 6

7 Assumig, firs, ha ad r S are saioary sochasic processes, vecor [ ] S r Z is also = a saioary sochasic process. I is he assumed ha a p-h order vecor auoregressio (VAR) represeaio of Z may be used o forecas he PFS. This bivariae VAR(p) may be wrie i compaio form (as a 1 s order VAR): z ω (6) + 1 = Az + +1 where, A, is a square (2p x 2p) marix of coefficies ad z is a (2p x 1) vecor of regressors, [ ] r r r S S... S z = -1 p+ 1 1 p The compaio form VAR may he be used o forecas (chages i) fuure shor raes as required by he EH equaio (2), by defiig wo (2p x 1) selecio vecors e1 ad e2 such ha S = e1 z ad r = e2 z. These vecors choose he wo key variables of ieres, S ad r, hrough oe eleme equal o 1, which is he 1 s eleme, for e1, ad he p+1 s eleme (row) for e2, wih all remaiig (2p-1) elemes beig zero. Usig he he VAR from equaio (6) o forecas oe period ahead: Ε Ε r 1 = e2 z 1 = e2 + + Az (7) More geerally, he chai rule of forecasig for muli period expecaios usig he VAR implies: E r e2 A z i + i = (8) The geeralised VAR forecas for equaio (2): Ε r + i from equaio (8) ca he be replaced i he EH 1 * i i S = Ε S = 1 e2 A z i = 1 (9) 7

8 The forecas of he PFS usig he VAR, i he righ had side of equaio (9), is he heoreical spread, S. Therefore, if he EH (2) is rue, he acual spread should be equal o he heoreical spread, e 1 z = S. From (9), Campbell ad Shiller (1991) derive he followig resricios o he VAR parameers whe he EH is rue: 1 f ( a) = e1 e2 A I 1 1 ( I A )( I A) ( I A) = 0 (10) The (2p x 1) o-liear resricios, f ( a), expressed i (10) deped oly o he esimaed parameers, Â, of he uresriced VAR model (6). A Wald es may he be used (as i Campbel ad Shiller, 1987) o es heir validiy: W = f 1 ( a) { [ f ( a) ]} f ( a) var (11) The variace of he resricios, var [ f ( a) ] VAR coefficies, Σ:, depeds o he variace-covariace marix of he var = (12) [ f ( a) ] f ( a) Σ f ( a) a a where f a (a) is he firs derivaive of he resricios evaluaed a he coefficie esimaes (from A). I addiio o he Wald es, he VAR mehodology employs a small lis of oher measures of proximiy bewee S ad S, ha ca be used o evaluae he EH, as i suggess hrough (9) ha he wo series mus be equal. Besides he iformal use of a graph of he wo series, saisical measures iclude a raio of sadard deviaios, he correlaio saisic, ad a OLS regressio es: ( S ) sd( S ) = 1 sd (13) ( S, S ) = 1 ρ = Corr (14) 8

9 S = α + βs + ν (15) where we expec α=0 ad β=1. 3. Small sample bias ad ieres rae persisece Cosiderig agai he previous es saisics, i ca be show ha early all of hem suffer from small sample bias ha becomes quie severe i some cases, depedig o he ime series properies of ieres raes ad of course he size of he sample. All OLS based ess icludig he VAR regressio ess poeially provide biased coefficie esimaes (Bekaer e al., 1997). Similarly, he variace bouds family of ess (Flavi, 1983) also eds o disor iferece i samples of ypical size. This secio focuses o he bias i he slope coefficie of he PFS regressio es i (4). Briefly, i his case, he OLS coefficie esimaes are biased due o he presece of a sochasic ad o ecessarily exogeous regressor (S ). Firs, we model he oe period ieres rae as a firs order auoregressive (AR(1)) process. Besides he desirable propery of greaer aalyical racabiliy (i he EH framework), he AR(1) model usually describes well ieres rae daa for acual ecoomies. Thus, he oe period ieres rae, r, is: r = µ + ρ 1 + ε (16) r where ρ 1 ad he error erm, ε, is a ormally ad idepedely disribued process wih mea zero ad a cosa variace, σ ε 2. Whe he EH (+RE) is rue, he log rae ad he spread are proporioal o r, hece, hey oo are persise. I greaer deail, whe µ = 0 (i 16) ad he AR(1) is subsiued ou o forecas fuure oe-period shor raes i he fudameal EH equaio (1) he log rae is: R 1 1 j = r 1 ρ = j = 0 1 ( 1 ρ ) r ( ρ ) (17) 9

10 Accordigly, he spread, S, is also aalogous o r : ( 1 ρ ) 1 S = R r = 1 r ( ) ( ) = η ρ, r 1 ρ (18) 1 where ( ) ( 1 ρ ) η ρ, = 1. ( 1 ρ ) Usig he AR(1) represeaio for r, expressio (18) implies ha he spread is persise ad he persisece parameer is, agai, ρ: S ( ρ )( ρr 1 + ε ) = ρη( ρ, ) r 1 + η( ρ, ) ε = ρs ν = 1 η, + (19) The PFS regressio (4) ad expressio (19) esseially form he followig sysem of equaios: S S * = α + βs = ξ + ρs 1 + u + v 2 where u ~ i.i.d(0, σ ) 2 where v ~ i.i.d(0, σ ) u v (20) (21) The secod equaio implies ha S is o fixed i repeaed samples, violaig a OLS requireme for he esimaio of he PFS equaio (20). Moreover, he error erms, u ad v, are coemporaeously correlaed wih covariace σ uv. Sadard ecoomeric heory ells us ha if S, he idepede variable i he PFS regressio, is sochasic he his is o ecessarily problemaic for he ypical properies of he OLS esimaor. The Gauss-Markov heorem may sill hold provided he error process of he regressor is idepede of he error process of he regressio iself. However, eve i ha case, he small sample properies of he OLS esimaor may be poor. 1 Neverheless, i his sysem, here is also a correlaio (σ uv ) bewee he disurbaces of he idepede variable (S ) ad he disurbaces of he PFS regressio (20) iself ad herefore he Gauss-Markov heorem o loger applies. I he paricular example, his is because he acual spread iheris o oly he persisece bu also he sochasic compoe of he shor rae, ε, 10

11 from he AR(1) model (16). The PFS o he lef had side of (20) also coais ha compoe as i icludes r +1 (from EH (2)). Therefore, he OLS coefficie esimaes, ˆ α ad ˆ β, are cosise bu biased, i samples of ypical size. Moreover he samplig disribuio of his small-sample esimaor is o sadard. This problem is well kow i he fiace lieraure, as i applies o forecasig equaios ofe used i ess of asse pricig heories, icludig he PFS regressio (e.g., Sambaugh 1986, 1999; Makiw ad Shapiro, 1986). Sambaugh (1986) has show ha if he sochasic regressor is also persise (as S is here) he he bias i esimaig he OLS slope coefficie depeds o he bias i esimaig he persisece coefficie, ρ, as follows: Ε ( ˆ σ uv β β ) = ( ˆ ρ ρ ) σ Ε 2 v (22) The bias i esimaig he auocorrelaio coefficie (ρ) of a auoregressive process is egaive (dowward) (Marrio ad Pope, 1954). The closer he rue coefficie, ρ, is o 1, he more severe his ype of bias becomes. Moreover, i he PFS regressio, he covariace, σ uv, is egaive so he slope bias (i (22)) is posiive (upward). Bekaer e al. (1997), illusrae ha i he case of he PFS regressio he slope bias expressio (22) will implicily ivolve esimaio of higher order auocorrelaio coefficies (due o he fuure shor raes ivolved i cosrucig he PFS) as well. They use Kedal s (1954) aalyical approximaios for he bias i higher order auocorrelaio coefficies o reach a more specific expressio (ha 22) of he bias i βˆ, i he PFS regressio. The bias is expressed as: 1 1 ˆ 1 (1 ) Ε( β ρ 1) = θ j = θ 0 j η (1 ρ) (1 ρ ) (23) ( ρ, ) j= 1 1 where ( ) ( 1 ρ ) j= 1 η ρ, = 1 ad θ ( 1 ρ) j is Kedall s approximaio o he esimaed auocorrelaio coefficie bias: 1 E.g. see pp i Cuhberso e al. (1992). 11

12 ( ρ ρ ) ( 1+ ρ ) 1 θ j = E ˆ j j = ( ) 2 (24) T 1 ( ) ( j ) j 1 ρ + jρ 1 ρ where j is he order of auocorrelaio. Hece, i his case: Ε ( ˆ β 1) = f (, T, ρ). The esimaed βˆ is posiively biased, above 1. The bias icreases wih higher ρ, smaller T ad smaller. The bias i expressio (23) is paricularly srog for very high auocorrelaio coefficies, ρ, akig values bewee 0.90 ad Bekaer e al. (1997) show ha he empirical disribuio is highly dispersed ad skewed o he righ. I akes almos as may as 2,000 observaios for symmery. Cosiderig, ow, he VAR ess ad saisics, Bekaer e al. (1997) show ha he VAR OLS coefficies are all biased, depedig o persisece. The bias passes hrough o early all he VAR-based ess ad merics. I summary, hey show ha: a) o a firs order approximaio, Corr spread ( S, S ) = 1, he coefficie of correlaio bewee he spread S ad he heoreical S is ubiased for small samples ad b) he sadard deviaio raio ( S ) sd( ) sd is S biased (upwards) as ρ rises, exceedig is heoreical value of 1 2. Bekaer ad Hodrick (2001) show ha he small sample disribuio of he Wald es of he EH resricios for he VAR coefficies is severely disored. The Wald saisic eds o overrejec he EH (biased upwards). 4. Moe Carlo experimes 4.1 Moe Carlo Srucure The resuls discussed i he previous secio show ha small sample iferece is esseial for esig he EH, especially whe ieres rae daa display a high degree of persisece. Give ha asympoic iferece is misleadig i fiie samples, well-desiged Moe Carlo simulaio (MCS) experimes ca be used o correc for small sample bias (e.g., Logsaf, 2000). Aalyical values from closed form esimaes ca eiher be edious o reach (i he case of he VAR) or of doubful validiy for high ρ. Moreover, esimaes of he bias i he coefficies are o sufficie for saisical iferece give also he skewess of he small sample disribuios. 2 Campbel ad Shiller (1991) also repor Moe Carlo evidece for he small sample disribuios of he VAR saisics. They fid a smaller bias, o large eough o overur he rejecio of he EH for heir daa. 12

13 We use Moe Carlo simulaio experimes o explore he sesiiviy of he small sample disribuios of he differe EH ess, described i secio 2, o chagig persisece levels. The AR(1) process for r i (16) is ceral o he Moe Carlo aalysis, as we are heoreically ieresed i he resposiveess of he shape of hese disribuios o persisece, ρ. However, researchers wishig o es acual daa migh boosrap heir esimaes ad correc he daa geeraig process for he bias i ρ, while also adopig a richer represeaio of realiy wih a VAR or VAR-GARCH model for he shor rae (e.g., Logsaf, 2000; Bekaer e al., 1997). We geerae daa uder he ull of EH (=EH+RE). This requires us o use he fudameal EH equaio (1) o geerae he log rae R (for =12), which i ur requires us o subsiue expeced fuure oe period raes. We, firs, geerae (T+200) observaios of he oe period rae, r, usig he AR(1) model i equaio (16) wih ormally disribued radom errors ad he recursively subsiue for fuure shor raes (RE) i he log rae model (EH). For he geeraed series, we he discard he firs 200 observaios o miimise he ifluece of iiial codiios. The acual spread we cosruc he PFS ( S 12 S is he he differece bewee he wo ieres rae series, R 12 ad r, ad *12 ) (from he EH equaio (2)) as he weighed average of chages i he forwarded shor raes, subsiuig recursively he AR(1) model. Havig geeraed series for he PFS ad he acual spread, complyig wih he EH, we esimae i each experime: 1) he PFS OLS regressio es (4) ad 2) i a ew roud of Moe Carlo experimes, he sadard deviaio raio of he wo series (5). I a hird roud of experimes, wih he same process for daa geeraio, ad afer calculaig he acual spread 12 S ad r we esimae a VAR(1) model ad use he coefficie esimaes o cosruc he heoreical spread, S. We he esimae i each replicaio he correspodig OLS regressio es (15) ad he VAR merics (correlaio ad sadard deviaio raio). For all hese experimes, we derive he small sample disribuios of he coefficies ad saisics, replicaig he exercise 5,000 imes. We repea he 5,000 Moe Carlo replicaios for arificial samples of (T=) 124, 524, 2,000 ad 20,000 observaios. Aiheic variaes are used i all experimes wih he excepio of he variace boud es for which we repea 10,000 imes (isead of 5,000) he Moe Carlo rials. For he AR(1) parameers i (16) we assume fixed values for µ = 1.24 ad σ 2 ε = 2.28, wih ρ chagig. We also validae he code by ryig he same se of parameer values for he AR(1) process as hose used i Bekaer e al. (1997) i heir similar experime. For he slope of he PFS regressio es, we fid bea=

14 agais heir 1.51 wih a sadard deviaio s.d.=0.62 agais 0.65 (furher similariies are documeed i figure 1 ad figure 2). 14

15 4.2 Simulaio resuls The perfec foresigh spread OLS es Table 1 repors Moe Carlo resuls (for 5,000 rials) ha oulie he small sample disribuio of he slope coefficie β i he PFS regressio es from experimes i which (AR(1)) ieres raes (from equaio (16)) were geeraed, as described i he previous secio, for differe ρ. We also explore larger sample properies of he es for samples of 524, 2,000 ad 20,000 observaios. Table 1. OLS slope coefficie β for PFS 12: ρ ad sample size sesiiviy Pael A, repors Moe Carlo resuls for differe samples ad for a broad rage of ρ bewee 0 ad 1. For each ρ, he mea slope coefficie β ad he sadard deviaio aroud his esimae are repored. The firs value repored is from OLS esimaio ha icludes a cosa erm, while he value udereah deoes OLS esimaio wihou a cosa erm. Pael B, repors resuls for very high values of ρ ad a sample of 124 observaios, for hree quarile levels of he small sample disribuio (besides he mea ad he sadard deviaio). Pael A Larger samples: Number of Observaios (T) = ,000 20,000 Mea βˆ Mea βˆ Mea βˆ Mea βˆ ρ (s.e.) (s.e.) (s.e.) (s.e.) ρ Mea βˆ (s.e.) Pael B High Persisece T=124 Quariles (0.044) (0.044) (0.015) (0.016) (0.007) (0.007) (0.002) (0.002) (0.226) (0.226) (0.053) (0.056) (0.023) (0.023) (0.011) (0.011) (0.004) (0.004) (0.416) (0.395) (0.078) (0.081) (0.038) (0.038) (0.019) (0.019) (0.006) (0.006) (0.689) (0.628) (0.146) (0.150) (0.077) (0.078) (0.041) (0.040) (0.013) (0.013) (2.305) (1.673) (0.159) (0.158) (0.079) (0.080) (0.041) (0.041) (0.013) (0.013) (0.159) (0.158) *Noe: 1s row for each value of he auocorrelaio coefficie, ρ, refers o sadard OLS value of he slope ad he secod row refers o he resriced-cosa OLS value of he slope. 15

16 Overall, he resuls are i accordace wih he aalyic approximaio 3 ad Moe Carlo evidece i Bekaer e al. (1997). For very high ρ, above 0.85, (i Pael B) we fid sigifica posiive bias ad disorios of he small sample disribuio of he esimaor. However, for lower degrees of persisece (Pael A), a egaive bias is also prese. For example, he aalyical approximaio esimae for ρ = 0.75 ad a sample of 124 observaios is β = 1.08, which implies a very small posiive bias (agais he heoreical value of 1). The mea esimae i our Moe Carlo resuls is 0.98, suggesig a very small egaive bias i he esimae of he slope coefficie, β. The bias is prese for every ρ i pael A of he able, for samples of 124 ad 524 observaios. I is due o small sample iflueces, oher ha ha from he persisece of ieres raes, perhaps also o be expeced. The resuls i Pael A of able 1 sugges wo paers as we icrease ρ ad he sample size, T. As ρ varies from 0 o 0.75 he mea slope, βˆ, (ad he sadard deviaio aroud i) rises from 0.92 o 0.98 for 124 observaios ad similarly for samples up o 2000 observaios 4. Thus as auocorrelaio rises, he implied posiive bias i esimaig β (described i Bekaer e al., 1997) rises, parly offseig he egaive bias foud i small samples. Cosiderig ow larger samples, he value of he slope β approximaes he heoreical value of 1, cofirmig he cosisecy propery of he OLS esimaor. I is srikig ha eve for a subsaially large sample of 2000 observaios he mea slope, which is very close o 1, sill slighly rises wih ρ, a he hird decimal poi. The slope disribuio collapses exacly o 1, wih a very small sadard deviaio, oly for a sample of 20,000 observaios. Pael B of able 1, repors resuls for high persisece, wih ρ bewee 0.90 ad 0.99 ad a sample of 124 observaios. The posiive bias ow becomes large (e.g. mea βˆ =1.63 for ρ = 0.96) ad he small sample disribuio is widely dispersed (sadard deviaio = 0.69). We fid ha 97.5% of he slope coefficie esimaes for ρ = 0.96 are below 3.11, implyig ha if a slope coefficie is esimaed o be as high as 3 i acual daa for a marke of he same degree of ieres rae persisece, he PEH should sill o be rejeced, usig his able. The resuls sugges ha, usig small sample, raher ha sadard asympoic, iferece, i is more difficul o rejec he ull 3 This is proposiio 3 i Bekaer e al. (1997), saed i equaio (23) here. 4 We also repor he ui roo case, for ρ = 1, which behaves differely. The aalyical approximaio of he bias refers oly o values smaller ha 1. 16

17 hypohesis whe ieres rae persisece is high due o he large dispersio of he small sample disribuio 5. Figure 1. PFS es slope coefficie: small sample agais asympoic disribuio Desiy curves (scaled) for small sample vs Asympoic Normal disribuio 0.04 Relaive frequecy Value of he slope coefficie (β ) small sample db N(1,0.0924) Noe: MCS resuls for sample size T=124 observaios ad ρ= NW-s.e. = Figure 2. Moe Carlo desiy curves for larger samples 0.25 ρ = Relaive Frequecy Value of slope coefficie (β) We repeaed all he above experimes excludig his ime he cosa erm from he OLS regressios. We fid ha for high ρ, he mea esimae of β, i he small sample disribuio (i able 1) ad he media (which is o repored i he able) are closer o 1 (see figure 3) ha hey 5 See Bekaer e al,

18 were i sadard OLS. The bias is subsaially reduced hough o elimiaed. However, he sadard deviaio of he small sample disribuio is agai high. Therefore, for high values of ρ, he slope coefficie bias is almos elimiaed whe he cosa erm is excluded from OLS regressios bu i is sill more difficul o rejec he ull hypohesis (EH), ha i is i large samples. The disicio bewee he wo esimaios is o impora for lower ρ bewee 0 ad 0.75 (pael A). Figure 3. Shif i he mea of he small sample disribuio: zero OLS cosa erm 0.03 Icludig vs excludig OLS cosa (124 obs) Relaive frequecy Value of he slope coefficie (β ) exclude c iclude c Table 2 repors Moe Carlo evidece for he sesiiviy of he small sample disribuio across differe spread horizos 6. The spread is calculaed bewee he (EH-cosruced) log rae for = 12, 9, 6, 3, 2 ad he AR(1) oe period rae. The sample cosiss of 124 observaios ad ρ is fixed o 0.93 (µ = 1.24). 6 For a sudy of marke efficiecy ess across differe horizos see Campbell (2001). 18

19 Table 2. OLS es for PFS across differe spreads: ρ=0.93, T=124. Spread mauriy, Mea β (s.e.) (aalyic al) - raio agais 1 Quariles: media (0.416) (1.415) (0.464) (1.429) (0.517) (1.443) (0.589) (1.459) (0.628) (1.465) Noe: Aalyical refers o he heoreical approximaio, usig proposiio 3 i Bekaer e al. (1997). Resuls i able 2 show ha as he mauriy () of he log rae i he spread S becomes shorer, he mea βˆ (ad he media) from he MCS experimes ad he sadard deviaio ed o rise. I he case of S 2, he mea Moe Carlo esimae of βˆ is 1.45 (compared o he aalyical esimae of 1.47), agais 1.27 foud for S 12. These resuls are broadly i agreeme wih he aalyical approximaios: as rises ( ˆ β 1 ) E (i equaio (23)) becomes smaller. Therefore, i all he previous Moe Carlo resuls i able 1, for which S 12 was used as a bechmark case, he bias i βˆ is udersaed whe oe refers o oher spreads. 19

20 4.2.2 Variace bouds I he ex roud of Moe Carlo experimes we explore he small sample disribuio of he * variace boud raio (i (5)), d. ( S ) s. d. ( S ). s, for he same se of values for chagig ρ ad T. We use agai he AR(1) oe period rae o geerae he log rae uder he ull hypohesis (EH), for 10,000 replicaios 7. Table 3: Volailiy raio for chagig ρ ad T * The raio VR= s d. ( S ) s. d. ( S ). is calculaed over 10,000 rials, for =12. Pael A, covers a wide rage of values for ρ, bewee 0 ad 1 for samples of 124, 524, 2,000 ad 20,000 observaios. Pael B cosiders oly high ρ for a small sample of 124 observaios. ρ Pael A Larger samples: Number of Observaios (T) = ,000 20,000 Mea VR (s.e.) % of raios >1 (ou of all rials) Mea VR (s.e.) Mea VR (s.e.) Mea VR (s.e.) Pael B High Persisece Number of Observaios (T)= 124 ρ Mea VR (s.e.) Quarile s (0.025) 45.8% (0.008) (0.004) (0.001) (0.077) (0.029) 15.7% (0.012) (0.006) (0.002) (0.085) (0.040) 3.4% (0.018) (0.009) (0.003) (0.078) (0.063) 0.6% (0.030) (0.016) (0.005) (0.077) (0.018) 0.0% (0.005) (0.001) (0.000) (0.018) We double he umber of replicaios o 10,000, i his case, as he aiheics echique does o make a differece for calculaig sadard deviaios. 20

21 We fid ha he mea volailiy raio esimaed over all rials, for a sample of 124 observaios ad for ρ = 0, is exacly equal o 1 ad he 97.5% quarile (which is o repored i he able) is a 1.06, which violaes he boud. This seems surprisig a firs, give ha he heoreical value of 1 is he exac upper boud ad ha daa were geeraed uder he ull. The violaio of he boud is observed for higher ρ as well. This is he ype of violaio of he boud explaied i J. Flavi s kow paper (1983), hough for ieres rae, raher ha spread, volailiy raios. I small samples he raio is biased upwards edig o violae he boud ad rejec he EH oo ofe, whe he ull hypohesis (EH) is rue. However, i our Moe Carlo resuls, as ρ rises, he mea raio eds o fall. For ρ = 0.75 ad a sample of 124 observaios, he mea raio is equal o 0.79 agais 0.91 for ρ = 0.5 ad correspodigly he 97.5% quarile level is foud a 0.94 agais 1.01 (o i he able). The sadard error aroud he mea Moe Carlo esimae also rises wih ρ. Therefore, i opposie direcio o he upward small sample bias here is also a dowward bias o he esimaed raio, which is sroger whe persisece rises. The persisece effec is more evide i pael B (of able 3) for ρ > The secod colum of he able (3) shows he opposig biases more clearly. Whe ρ is low here are freque violaios of he boud (e.g., a ρ = 0, 46% of he raios are above 1). As persisece rises, however, i is foud ha violaios fall owards 0, as a resul of he domiaig, dowward, persisece bias foud i hese experimes. For larger samples, biases ed o be elimiaed. However, he bias due o higher persisece sill exiss give ha he raio eds o fall whe persisece is higher. Eve for a very large sample of 20,000 observaios, a ρ = 0.50, he mea raio is 0.93, fallig o 0.77 a ρ = Overall, large sample resuls do o collapse o a specific value, ulike he case of OLS, where he esimaor was foud cosise. Ulike he OLS case, he variace bouds raio is o expeced heoreically o coverge o a specific value ad i mus oly saisfy he iequaliy. 8 The ui roo sceario, ρ = 1, is agai a special case ad he mea raio dramaically falls o

22 Figure 4. Moe Carlo hisogram: Volailiy raio ρ = 0.93, 124 observaios Table 4. Volailiy raios across differe spreads: ρ =0.93, T=124 Spread mauriy, 12 Mea VR (s.e.) (0.085) media (0.076) (0.063) (0.046) (0.035) Table 4, repors resuls for differe mauriy horizos,, of he spreads S. The level of persisece is fixed a ρ = 0.93 i all experimes, for samples of 124 observaios. Similarly o our previous resuls for he OLS regressio es, whe he log rae mauriy,, becomes shorer, he bias is larger ad he mea raio for S 2 falls o 0.17 agais 0.47 for S 12. Agai, he resuls repored i able 3 ad referrig o S 12, udersae he persisece bias for shorer mauriy spreads. 22

23 4.2.3 VAR ess ad merics I his roud of Moe Carlo experimes we explore he properies of he VAR model. Usig he same parameer values, we geerae r as a AR(1) process ad he log rae accordig o he EH, as before. I each rial, we geerae daa uder he ull, we esimae by OLS a bivariae VAR(1) model: S a S a r u (25) + 1 = , + 1 r a S a r u (26) + 1 = , + 1 We he calculae he heoreical spread ad he associaed saisics of proximiy o he acual spread: i) he sadard deviaio raio, sd(s )/sd(s ), ii) he correlaio saisic, Corr(S,S ), iii) he slope coefficie from he OLS regressio of S agais S ad iv) he Wald es of resricios for he VAR coefficies, A, implied by he EH (i equaio (10)) ad heir (Newey-Wes) robus variace covariace marix. Table (5) repors resuls describig he small-sample disribuio of sd(s )/sd(s ) ad is sesiiviy o chagig values of ρ ad sample size, T. I accordace wih he aalyic approximaio i Bekaer e al. (1997), he raio sd(s )/sd(s ) is biased upwards i a small sample, whe ρ is risig. For ρ up o 0.75, i pael A, he mea raio, sd(s )/sd(s ), is oly slighly above he heoreical value of 1 ad risig wih ρ. As he sample size, T, rises, for each ρ, he mea raio falls closer o 1. The bias is more proouced for ρ bewee 0.85 ad 0.96, i pael B, as he mea sadard deviaio raio rises from 1.01 o The small sample disribuio becomes more dispersed, wih sadard deviaios aroud he mea raios correspodigly risig from 0.18 o For ρ = 0.96, heoreical sadards of closeess o 1 would say ha sd(s )/sd(s ) = 1.6, is o close eough, however, uder small-sample iferece usig able 5, 1.6 is he mea raio, (wih a high sadard deviaio of 0.66) ad he EH is correcly suppored. 23

24 Table 5: VAR merics: sd(s )/sd(s ) for PFS 12: ρ ad T sesiiviy Pael A Larger samples: Number of Observaios (T) = ,000 Pael B High Persisece Number of Observaios (T)= 124 ρ Mea raio (s.e.) Mea raio (s.e.) Mea raio (s.e.) ρ Mea raio (s.e.).975 (1-righ ail) (0.011) (0.005) (0.028) (0.177) (0.018) (0.010) (0.005) (0.353) (0.035) (0.019) (0.098) (0.663) (0.094) (0.053) (0.017) (1.785) Table 6. VAR merics: sd(s )/sd(s) across spreads, ρ =0.93, T=124 Spread mauriy, 12 9 Mea raio (s.e.) (0.353) (0.447) media (0.528) (0.618) (0.6513)

25 Table 6 preses resuls for he sadard deviaio raio across differe mauriy horizos. The mea raio (ad sadard deviaio aroud i) rises as he log rae mauriy,, of he spread, S, shores. This resul is o i accordace wih he correspodig aalyic approximaio i Bekaer e al. (1997), accordig o which his saisic does o (o a firs approximaio) deped o he mauriy horizo,. However, a similar paer ca be observed i heir Moe Carlo resuls. I he same Moe Carlo exercise we esimae he correlaio coefficie, Corr(S,S ), ad fid ha i is ubiased wih respec o risig persisece. As expeced from heory, i collapses almos perfecly o is heoreical value (=1), for all ρ. Moreover, i coras wih heory, we fid ha he mea Corr(S,S ) varies o a small degree across mauriy horizos,. We do o repor hese resuls i a able, as ay variaios are oo small. Lookig ow a he mai saisical es of he VAR mehodology, he Wald es of he VAR cross-equaio resricios, able 7 shows resuls from he same experimes (for varyig ρ ad sample sizes, T). Give ha here are o racable aalyic approximaios for he heoreical bias i his saisic, wih respec o persisece, we oly rely o resuls from Moe Carlo evidece. Similarly o Bekaer ad Hodrick (2001), we fid ha he small sample properies of he Wald es of he EH resricios are very poor. 25

26 Table 7. VAR ess: Wald es for PFS 12: ρ ad T sesiiviy Pael A Larger samples: Number of Observaios (T) = ,000 ρ 0.95 quarile (1-righ ail) ρ Mea (s.e.) Pael B High Persisece T= (8.510) (13.434) (9.895) (10.935)

27 Table 8. VAR ess: Wald es across spreads, ρ =0.93, T=124 The firs row of he able quoes criical values from he χ 2 (2) disribuio, from sadard saisical ables. Spread S Mea (s.e.) quariles media χ 2 (2 d.f.) Na S 2 S (3.470) (3.992) S (5.091) S (6.924) S (11.521) Overall, we fid a proouced upward bias i he mea Wald es saisic, esimaed over 5,000 Moe Carlo rials, ad a more dispersed disribuio, as ρ rises, eve i larger samples. For fixed ρ, he bias falls, as T rises. More paricularly, whe persisece, ρ, is low ad samples are large he Wald es saisic is sill close o he criical values suggesed by he χ 2 disribuio, i pael A (able 7). For example, for large T=2,000, whe ρ =0.5, he 95% quarile level is a 6.21 agais he (χ 2 ) criical value of However, whe ρ is raised o 0.75, for he same T, he 95% quarile level rises o Pael B (i able 7) shows ha for very high values of ρ ad for a sample of 124 observaios, he small sample disribuio of he Wald saisic is cosiderably disored. Uder small sample iferece ad for ρ = 0.96, he EH should oly be rejeced for Wald es values higher ha (agais he sadard iferece value of 5.99), a he 95% quarile. Hece usig he Wald es uder asympoic iferece he EH eds o be over-rejeced i small samples. Nex, we explore i Moe Carlo resuls i able 8, differeces across chagig mauriy horizos,, for he spread, fixig ρ = 0.93, for a sample of 124 observaios. The bias i he mea Wald es saisic is higher for loger horizo spreads ad he 95% quarile level rises o for S 12 agais oly 9.67, for S 2. Therefore, ulike all previous ess, S 12 which was used i all he 27

28 experimes repored i able 7, oversaes he degree of persisece bias ha would apply o shorer mauriy spreads. We examie las he OLS regressio es of he heoreical spread S agais S i Moe Carlo resuls repored i able 9. I should be oed ha durig hese experimes, OLS esimaio failed o a grea exe i large samples. Therefore, our resuls focus o a small sample of 124 observaios, for differe ρ. Table 9. OLS βˆ from regressio of S agais S, T= 124, high ρ ρ Mea β (s.e.) - raio agais 1 quariles: media (1-righ ail) 0.05 (lef ail) (0.011) (0.036) (0.095) (0.178) (0.354) (0.665) (2.521) The resuls show ha a posiive bias i he slope coefficie is prese ad he small sample disribuio is highly dispersed. The mea (ad he media) βˆ rises above 1 as ρ, rises owards The sadard deviaio aroud βˆ also rises. Focusig o high ρ, able 9 shows a more proouced slope coefficie bias ha is eve larger ha ha recorded for βˆ i he PFS es, i able 1. Table 9 shows ha, for a real marke i which ρ is ear 0.92, if we fid βˆ = 1.90 he he EH mus sill o be rejeced, usig small sample 28

29 iferece, a he 97.5% quarile level. Moreover, i is ieresig o oe ha i larger samples he mea βˆ foud here coverges o 1 more quickly ha i did i he PFS regressios. Figure 5. Moe Carlo resuls: OLS βˆ (S agais S ) hisograms for high ad low ρ Table 10. OLS βˆ from regressio of S agais S, across spreads, T=124 Spread S Mea β (s.e.) -raio agais 1 quariles: media (1-righ ail) 0.05 (lef ail) S (0.642) S 3 S (0.614) (0.528) S 9 S (0.448) (0.354)

30 Table 10 preses resuls across mauriy horizos for differe S. Similarly o our resuls for he PFS regressio es i able 2, he small sample bias i βˆ eds o fall whe he log rae mauriy horizo,, is larger. For S 2, he mea βˆ is 1.48 ad i falls for larger (e.g., for S 12 βˆ = 1.26). Thus, he EH may be rejeced, a he 97.5% level, for βˆ = 2.9 for some spreads (e.g., for S 9 ) while i may o be rejeced, for he same value, for oher spreads (e.g. S 2 ) The Kalma filer-ime varyig parameer model Oe of he explaaios for possible empirical failures of he EH is ha especially whe here are regime chages, parameers may be usable raher ha cosa as i is assumed by OLS. I learig models, isead of raioally kowig a cosa parameer model, ages are updaig he model s parameers as hey gai more iformaio hrough ime. I his las roud of Moe Carlo experimes we cosider he properies of a ime varyig parameer model (TVP). We fid ha TVP models are a ieresig case of exreme small sample bias. Alhough hese models help elimiae he bias from regime shifs, hey are cosiderably more exposed o persisece bias ha he OLS esimaor. We geerae daa uder he ull (EH), exacly as i he previous experimes, for a AR(1) shor rae ad a log rae cosruced accordig o he EH. We he esimae, i each Moe Carlo rial, a ime varyig parameer model 9 of he same PFS equaio (4) ha was esimaed by OLS i secio (4.2.1). For a more horough descripio of he model, see Kim ad Nelso (1999). The wo mai equaios of he sae-space form are: S * = α + β S + ε (Measureme equaio) (3.23) α 1 = β 0 0 α 1 β 1 1 η1 + η 2 (Trasiio equaio) (3.24) where S * is agai he perfec foresigh spread (PFS). Give ha he EH is rue, we geerally expec ha ˆ 2 2 = 1ad ˆ α = 0, for mos. We esimae he wo hyper-parameers σ ad σ β ad he variace of he error erm i he measureme equaio, σ 2 ε. Table 11 repors resuls from 5,000 Moe Carlo replicaios for a sample of 124 observaios ad for ρ risig from 0 o η 1 η 2 9 We use Gauss code wrie by Kim ad Nelso (1999) for Kalma filer esimaio. 30

31 Table 11. Time varyig parameer PFS regressio, S 12 : Sesiiviy o ρ, T=124 The ime varyig cosa ad slope esimae are averaged over all replicaios ad all observaios (= T). ρ Mea σ 2 ε ( PFS equaio ) ( -raio) Mea σ 2 η1 ( α equaio ) (- raio) Mea σ 2 η2 ( β equaio ) (- raio) Mea α (s.d.) (- raio) Mea β (s.d.) (- raio) 0 1*10-8 (-0.004) *10-7 (0.001) *10-10 (-0.001) (9*10-5 ) (4*10-3 ) *10-9 (0.001) *10-8 (0.008) (14.581) (25.234) (19.160) (14.459) (13.000) (12.551) (11.220) 6*10-10 (0.009) 8*10-6 (0.824) 4*10-4 (3.698) (4.204) (1.341) (2.194) (4.071) (0.641) (0.002) (1.187) (-0.011) (2.103) (-0.003) (3.045) (0.007) (4.548) (-0.223) (6.500) (-0.154) (10.392) (-0.593) (0.020) (2.277) (0.032) (4.294) (0.065) (5.936) (0.131) (5.834) (0.232) (7.371) (0.549) (6.632) (3.094) (4.350) Cosiderig overall resuls from Kalma filer esimaio, for all ρ, he average σ 2 ε (over 5000 replicaios) is ear 0. The average hyper-parameer ha correspods o α, σ 2 η1, is relaively 2 large, implyig a very high sigal o oise raio ( σ σ 2 η1 ε ), compared o σ 2 η2, which correspods o β. Neverheless, he laer is also large eough o imply a high sigal o oise raio, for all ρ (excep a ρ = 0). Turig ow o he mea levels ad pahs for he wo TVPs, αˆ ad βˆ (i able 11), implied by he Kalma filer algorihm ad by he hyperparameer esimaes, we fid he same paer of sesiiviy o ρ as i he case of he OLS es of he EH (i secio 4.2.1). The ime varyig cosa erm αˆ is o average ear 0 bu he mea ime varyig slope coefficie βˆ rises from 1 owards much higher values, as ρ rises. The sadard deviaio aroud hese mea esimaes for αˆ ad βˆ rises as ρ also rises. More imporaly, he resuls sugges ha he upward bias i i he ime varyig parameer model, is impressively srog compared o ha foud for OLS esimaio i secio (e.g. for ρ = 0.96, βˆ = 4.64 agais 1.63 i OLS resuls). βˆ, 31

32 This is o surprisig give ha, i Kalma filerig he sample used for esimaio is esseially smaller ha 124 observaios. I he basic filer 10 algorihm, iformaio uil ime is used, which is before he ed of he sample, T(=124), ad i is gradually updaed a every ime icreme. As a resul, he sample size is smaller ha 124 a each. Moreover, he Kalma filer algorihm ca be ierpreed as a weighed leas squares esimaor, aachig smaller weighs o observaios ha are farher back i ime. Therefore, eve whe he las observaio T is reached, he ime varyig coefficie is esimaed usig a smaller (ha 124) widow of observaios. Kalma filer esimaes are more comparable o hose from a OLS esimaor ha uses a sample smaller ha 124 observaios. I a ime varyig parameer framework, he average values ad he able are oly iformaive whe complemeed by he correspodig plos of average coefficie esimaes for every. The properies of he esimaed ime varyig parameer model are bes summarised i he fourwidow plo i figure 6. The upper wo widows display he plos of he forecas error ad is associaed codiioal variace. The forecas error 11 is defied as: ν = ˆ ˆ ˆ β S ad he codiioal variace of he forecas error is 1 * y y 1 = S α ( ) f ν. 1 = Ε 1 10 Basic filer is used here o refer o simple filerig as opposed o smoohig which would use all iformaio up o T, he ed of he sample. 11 The firs few observaios display a higher forecas error, uil observaio 13. This is because a high uceraiy coefficie (=500,000), aached o he ages iiial guesses of he TVPs, is assumed i he code. 32

33 Figure 6. αˆ ad βˆ i deviaios from heoreical value (2 s.e. bads), for ρ =0.94 The lower wo widows i figure 6, display he plos of coefficies as deviaios from heir heoreical values ( αˆ =0 ad βˆ =1) (i.e. esimaes of he bias) 12. Figure 6 shows ha for ρ = 0. 94, βˆ - 1= 2.4, wihi a small bad of 2 sadard errors, while hypohesis would have bee rejeced based o sadard iferece. 13 αˆ = 0. The ull 12 The sadard error bads display he sadard deviaio of he disribuio obaied from he Moe Carlo rials. 13 We also fid ha ρ=0.99 is a limiig case. The resulig pah of he ime varyig coefficies is uusual ad he properies of he ime varyig model are o loger robus. The ime varyig slope coefficie, βˆ, rapidly rises over ime o 15. The TVP model is o applicable whe r is a ear- radom walk ad he EH is rue. 33

34 Figure 7. βˆ deviaio from 1, for 2000 observaios, ρ =0.93 TVP b for 2000 observaios bea: deviaio from observaio umber We also ivesigae larger samples of 524 ad of 2,000 observaios for his model ad fid ha he average Moe Carlo ˆ β 1 for samples of 524 ad 2000 observaios (figure 7) ad for a fixed ρ = is always equal o 2 as he sample size icreases. Ulike OLS, he TVP model is o a cosise esimaor, give ha, a each, he Kalma filer effecively uses a widow of observaios for he esimaio of βˆ, so ha raisig T will o elimiae he bias. 5. Coclusios We have used Moe Carlo Simulaio experimes o explore he small sample disribuio for a umber of ess of he EH. Our resuls cofirm ha he small sample disribuio for all hese saisics is more ad more disored i respose o risig persisece i ieres raes ad he spread. For all hese ess, we also fid ha mauriy horizo maers for small sample bias ad i mos cases, he laer is smaller for ess wih loger horizo spreads (S 12 ). I accordace wih resuls i Bekaer e al. (1997) we fid srog posiive bias i he OLS slope coefficie of he PFS es ad a small sample disribuio ha has greaer dispersio ad is skewed o he righ. The slope bias is almos elimiaed whe we ru he same regressios wihou a cosa erm. * Oe of he mos srikig resuls is ha he variace boud raio, s d. ( S ) s. d. ( S )., is biased dowwards, for high ρ, hece i eds o favour he EH oo ofe. This is i coras wih he 34

35 OLS es ha eds o rejec he EH oo ofe ad wih previous Moe Carlo evidece i he variace bouds lieraure, ha deals hough wih differe formulaios of he boud (Flavi, 1983). Furher Moe Carlo evidece shows ha all he VAR mehodology saisics are biased wih he excepio of Corr(S,S ). Small sample iferece would be preferable i he case of he sadard deviaio raio ha exhibis wide dispersio i is small sample disribuio. However, he case of he Wald es saisic sads ou, give ha asympoic iferece is pracically redered ieffecive whe persisece is high. The Wald es srogly eds o over-rejec he ull hypohesis i small samples, as ρ rises. We also look a he OLS slope coefficie i he regressio of S agais S, which, o our kowledge, is o cosidered elsewhere, ad fid a relaively high posiive bias. The Kalma filer ime varyig parameer model is o a ypical es of he EH, however, i belogs o a broader class of learig models ha may be more appropriae for esig he EH i samples which display impora regime shifs. A he same ime, i provides a appealig example of more exreme small sample bias (agais ha foud for OLS) as demosraed by our Moe Carlo evidece. Moreover, hese resuls are releva o a wider class of ecoomeric echiques ha may be applied o esig he EH or oher heoreical models. These would iclude weighed leas squares esimaio, recursive esimaio ad rollig widow esimaes, all of which effecively use a smaller widow of observaios ha OLS (for he same sample). The resuls preseed i his paper aim o provide a guide o praciioers wishig o es he EH i markes ha display high ieres rae auocorrelaio. Their coefficie esimaes are more likely o be comparable o he quarile levels documeed i our ables ha o he sadard criical values of asympoic iferece. Tha beig said, Moe Carlo simulaio would ehace esimaio ad iferece if ew criical values were used, derived from daa geeraed wih parameer calibraio o he specific marke. Fially, our resuls also serve as a heoreical ivesigaio of he properies of hese ess. I is hoped ha Moe Carlo evidece of differeces i he exe of small sample bias amog aleraive ess ca help explai he discrepacies foud i he empirical lieraure whe differe ess are employed. Towards he same direcio, a ivesigaio of he differeces i he power of hese ess, wih daa geeraio uder a specified aleraive o he ull hypohesis, would provide a possible exesio o our resuls. 35

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

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

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

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

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

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

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

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

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

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

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

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

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

Fresnel Dragging Explained

Fresnel Dragging Explained Fresel Draggig Explaied 07/05/008 Decla Traill Decla@espace.e.au The Fresel Draggig Coefficie required o explai he resul of he Fizeau experime ca be easily explaied by usig he priciples of Eergy Field

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

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

FOR 496 / 796 Introduction to Dendrochronology. Lab exercise #4: Tree-ring Reconstruction of Precipitation

FOR 496 / 796 Introduction to Dendrochronology. Lab exercise #4: Tree-ring Reconstruction of Precipitation FOR 496 Iroducio o Dedrochroology Fall 004 FOR 496 / 796 Iroducio o Dedrochroology Lab exercise #4: Tree-rig Recosrucio of Precipiaio Adaped from a exercise developed by M.K. Cleavelad ad David W. Sahle,

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

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

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

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

λ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

Stationarity and Error Correction

Stationarity and Error Correction Saioariy ad Error Correcio. Saioariy a. If a ie series of a rado variable Y has a fiie σ Y ad σ Y,Y-s or deeds oly o he lag legh s (s > ), bu o o, he series is saioary, or iegraed of order - I(). The rocess

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

A Two-Level Quantum Analysis of ERP Data for Mock-Interrogation Trials. Michael Schillaci Jennifer Vendemia Robert Buzan Eric Green

A Two-Level Quantum Analysis of ERP Data for Mock-Interrogation Trials. Michael Schillaci Jennifer Vendemia Robert Buzan Eric Green A Two-Level Quaum Aalysis of ERP Daa for Mock-Ierrogaio Trials Michael Schillaci Jeifer Vedemia Rober Buza Eric Gree Oulie Experimeal Paradigm 4 Low Workload; Sigle Sessio; 39 8 High Workload; Muliple

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

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

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

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

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

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

Development of Kalman Filter and Analogs Schemes to Improve Numerical Weather Predictions

Development of Kalman Filter and Analogs Schemes to Improve Numerical Weather Predictions Developme of Kalma Filer ad Aalogs Schemes o Improve Numerical Weaher Predicios Luca Delle Moache *, Aimé Fourier, Yubao Liu, Gregory Roux, ad Thomas Warer (NCAR) Thomas Nipe, ad Rolad Sull (UBC) Wid Eergy

More information

Chapter Chapter 10 Two-Sample Tests X 1 X 2. Difference Between Two Means: Different data sources Unrelated. Learning Objectives

Chapter Chapter 10 Two-Sample Tests X 1 X 2. Difference Between Two Means: Different data sources Unrelated. Learning Objectives Chaper 0 0- Learig Objecives I his chaper, you lear how o use hypohesis esig for comparig he differece bewee: Chaper 0 Two-ample Tess The meas of wo idepede populaios The meas of wo relaed populaios The

More information

Electrical Engineering Department Network Lab.

Electrical Engineering Department Network Lab. Par:- Elecrical Egieerig Deparme Nework Lab. Deermiaio of differe parameers of -por eworks ad verificaio of heir ierrelaio ships. Objecive: - To deermie Y, ad ABD parameers of sigle ad cascaded wo Por

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

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

6/10/2014. Definition. Time series Data. Time series Graph. Components of time series. Time series Seasonal. Time series Trend

6/10/2014. Definition. Time series Data. Time series Graph. Components of time series. Time series Seasonal. Time series Trend 6//4 Defiiio Time series Daa A ime series Measures he same pheomeo a equal iervals of ime Time series Graph Compoes of ime series 5 5 5-5 7 Q 7 Q 7 Q 3 7 Q 4 8 Q 8 Q 8 Q 3 8 Q 4 9 Q 9 Q 9 Q 3 9 Q 4 Q Q

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

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

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

Local Influence Diagnostics of Replicated Data with Measurement Errors

Local Influence Diagnostics of Replicated Data with Measurement Errors ISSN 76-7659 Eglad UK Joural of Iformaio ad Compuig Sciece Vol. No. 8 pp.7-8 Local Ifluece Diagosics of Replicaed Daa wih Measureme Errors Jigig Lu Hairog Li Chuzheg Cao School of Mahemaics ad Saisics

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

A Bayesian Approach for Detecting Outliers in ARMA Time Series

A Bayesian Approach for Detecting Outliers in ARMA Time Series WSEAS RASACS o MAEMAICS Guochao Zhag Qigmig Gui A Bayesia Approach for Deecig Ouliers i ARMA ime Series GUOC ZAG Isiue of Sciece Iformaio Egieerig Uiversiy 45 Zhegzhou CIA 94587@qqcom QIGMIG GUI Isiue

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

Pure Math 30: Explained!

Pure Math 30: Explained! ure Mah : Explaied! www.puremah.com 6 Logarihms Lesso ar Basic Expoeial Applicaios Expoeial Growh & Decay: Siuaios followig his ype of chage ca be modeled usig he formula: (b) A = Fuure Amou A o = iial

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

Page 1. Before-After Control-Impact (BACI) Power Analysis For Several Related Populations. Richard A. Hinrichsen. March 3, 2010

Page 1. Before-After Control-Impact (BACI) Power Analysis For Several Related Populations. Richard A. Hinrichsen. March 3, 2010 Page Before-Afer Corol-Impac BACI Power Aalysis For Several Relaed Populaios Richard A. Hirichse March 3, Cavea: This eperimeal desig ool is for a idealized power aalysis buil upo several simplifyig assumpios

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

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

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

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

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

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

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

Lecture 8 April 18, 2018

Lecture 8 April 18, 2018 Sas 300C: Theory of Saisics Sprig 2018 Lecure 8 April 18, 2018 Prof Emmauel Cades Scribe: Emmauel Cades Oulie Ageda: Muliple Tesig Problems 1 Empirical Process Viewpoi of BHq 2 Empirical Process Viewpoi

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

Detection of Level Change (LC) Outlier in GARCH (1, 1) Processes

Detection of Level Change (LC) Outlier in GARCH (1, 1) Processes Proceedigs of he 8h WSEAS I. Cof. o NON-LINEAR ANALYSIS, NON-LINEAR SYSTEMS AND CHAOS Deecio of Level Chage () Oulier i GARCH (, ) Processes AZAMI ZAHARIM, SITI MERIAM ZAHID, MOHAMMAD SAID ZAINOL AND K.

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

Convergence theorems. Chapter Sampling

Convergence theorems. Chapter Sampling Chaper Covergece heorems We ve already discussed he difficuly i defiig he probabiliy measure i erms of a experimeal frequecy measureme. The hear of he problem lies i he defiiio of he limi, ad his was se

More information

14.02 Principles of Macroeconomics Fall 2005

14.02 Principles of Macroeconomics Fall 2005 14.02 Priciples of Macroecoomics Fall 2005 Quiz 2 Tuesday, November 8, 2005 7:30 PM 9 PM Please, aswer he followig quesios. Wrie your aswers direcly o he quiz. You ca achieve a oal of 100 pois. There are

More information

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models Oulie Parameer esimaio for discree idde Markov models Juko Murakami () ad Tomas Taylor (2). Vicoria Uiversiy of Welligo 2. Arizoa Sae Uiversiy Descripio of simple idde Markov models Maximum likeliood esimae

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

t = s D Overview of Tests Two-Sample t-test: Independent Samples Independent Samples t-test Difference between Means in a Two-sample Experiment

t = s D Overview of Tests Two-Sample t-test: Independent Samples Independent Samples t-test Difference between Means in a Two-sample Experiment Overview of Te Two-Sample -Te: Idepede Sample Chaper 4 z-te Oe Sample -Te Relaed Sample -Te Idepede Sample -Te Compare oe ample o a populaio Compare wo ample Differece bewee Mea i a Two-ample Experime

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

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1)

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1) Fourier Series Iroducio I his secio we will sudy periodic sigals i ers o heir requecy is said o be periodic i coe Reid ha a sigal ( ) ( ) ( ) () or every, where is a uber Fro his deiiio i ollows ha ( )

More information

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i) Mah PracTes Be sure o review Lab (ad all labs) There are los of good quesios o i a) Sae he Mea Value Theorem ad draw a graph ha illusraes b) Name a impora heorem where he Mea Value Theorem was used i he

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

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

A Probabilistic Nearest Neighbor Filter for m Validated Measurements.

A Probabilistic Nearest Neighbor Filter for m Validated Measurements. A Probabilisic Neares Neighbor iler for m Validaed Measuremes. ae Lyul Sog ad Sag Ji Shi ep. of Corol ad Isrumeaio Egieerig, Hayag Uiversiy, Sa-og 7, Asa, Kyuggi-do, 45-79, Korea Absrac - he simples approach

More information

ASSESSING GOODNESS OF FIT

ASSESSING GOODNESS OF FIT ASSESSING GOODNESS OF FIT 1. Iroducio Ofe imes we have some daa ad wa o es if a paricular model (or model class) is a good fi. For isace, i is commo o make ormaliy assumpios for simpliciy, bu ofe i is

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

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

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

Institute of Actuaries of India

Institute of Actuaries of India Isiue of cuaries of Idia Subjec CT3-robabiliy ad Mahemaical Saisics May 008 Eamiaio INDICTIVE SOLUTION Iroducio The idicaive soluio has bee wrie by he Eamiers wih he aim of helig cadidaes. The soluios

More information

Problems and Solutions for Section 3.2 (3.15 through 3.25)

Problems and Solutions for Section 3.2 (3.15 through 3.25) 3-7 Problems ad Soluios for Secio 3 35 hrough 35 35 Calculae he respose of a overdamped sigle-degree-of-freedom sysem o a arbirary o-periodic exciaio Soluio: From Equaio 3: x = # F! h "! d! For a overdamped

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 LINEAR APPROXIMATION OF THE BASELINE RBC MODEL FEBRUARY, 202 Iroducio For f(, y, z ), mulivariable Taylor liear epasio aroud (, yz, ) f (, y, z) f(, y, z) + f (, y, z)( ) + f (, y, z)( y y) + f (, y, z)(

More information

Relationship between education and GDP growth: a mutivariate causality analysis for Bangladesh. Abstract

Relationship between education and GDP growth: a mutivariate causality analysis for Bangladesh. Abstract Relaioship bewee educaio ad GDP growh: a muivariae causaliy aalysis for Bagladesh Tariq Saiful Islam Deparme of Ecoomics, Rajshahi Uiversiy Md Abdul Wadud Deparme of Ecoomics, Rajshahi Uiversiy Qamarullah

More information

in insurance : IFRS / Solvency II

in insurance : IFRS / Solvency II Impac es of ormes he IFRS asse jumps e assurace i isurace : IFRS / Solvecy II 15 h Ieraioal FIR Colloquium Zürich Sepember 9, 005 Frédéric PNCHET Pierre THEROND ISF Uiversié yo 1 Wier & ssociés Sepember

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

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

Using GLS to generate forecasts in regression models with auto-correlated disturbances with simulation and Palestinian market index data

Using GLS to generate forecasts in regression models with auto-correlated disturbances with simulation and Palestinian market index data America Joural of Theoreical ad Applied Saisics 04; 3(: 6-7 Published olie December 30, 03 (hp://www.sciecepublishiggroup.com//aas doi: 0.648/.aas.04030. Usig o geerae forecass i regressio models wih auo-correlaed

More information

Statistical Estimation

Statistical Estimation Learig Objecives Cofidece Levels, Iervals ad T-es Kow he differece bewee poi ad ierval esimaio. Esimae a populaio mea from a sample mea f large sample sizes. Esimae a populaio mea from a sample mea f small

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

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

O & M Cost O & M Cost

O & M Cost O & M Cost 5/5/008 Turbie Reliabiliy, Maieace ad Faul Deecio Zhe Sog, Adrew Kusiak 39 Seamas Ceer Iowa Ciy, Iowa 54-57 adrew-kusiak@uiowa.edu Tel: 39-335-5934 Fax: 39-335-5669 hp://www.icae.uiowa.edu/~akusiak Oulie

More information

INVESTMENT PROJECT EFFICIENCY EVALUATION

INVESTMENT PROJECT EFFICIENCY EVALUATION 368 Miljeko Crjac Domiika Crjac INVESTMENT PROJECT EFFICIENCY EVALUATION Miljeko Crjac Professor Faculy of Ecoomics Drsc Domiika Crjac Faculy of Elecrical Egieerig Osijek Summary Fiacial efficiecy of ivesme

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

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

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

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

Distribution of Estimates

Distribution of Estimates Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion

More information

Optimization of Rotating Machines Vibrations Limits by the Spring - Mass System Analysis

Optimization of Rotating Machines Vibrations Limits by the Spring - Mass System Analysis Joural of aerials Sciece ad Egieerig B 5 (7-8 (5 - doi: 765/6-6/57-8 D DAVID PUBLISHING Opimizaio of Roaig achies Vibraios Limis by he Sprig - ass Sysem Aalysis BENDJAIA Belacem sila, Algéria Absrac: The

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

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

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

BAYESIAN ESTIMATION METHOD FOR PARAMETER OF EPIDEMIC SIR REED-FROST MODEL. Puji Kurniawan M

BAYESIAN ESTIMATION METHOD FOR PARAMETER OF EPIDEMIC SIR REED-FROST MODEL. Puji Kurniawan M BAYESAN ESTMATON METHOD FOR PARAMETER OF EPDEMC SR REED-FROST MODEL Puji Kuriawa M447 ABSTRACT. fecious diseases is a impora healh problem i he mos of couries, belogig o doesia. Some of ifecious diseases

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

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