Computer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures

Size: px
Start display at page:

Download "Computer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures"

Transcription

1 MPRA Munich Personal RePEc Archive Compuer Simulaes he Effec of Inernal Resricion on Residuals in Linear Regression Model wih Firs-order Auoregressive Procedures Mei-Yu Lee Deparmen of Applied Finance, Yuanpei Universiy 2014 Online a hps://mpra.ub.uni-muenchen.de/60362/ MPRA Paper No , posed 4 December :08 UC

2 Journal of Saisical and Economeric Mehods, vol.3, no.3, 2014, 1-22 ISSN: (prin), (online) Scienpress Ld, 2014 Compuer Simulaes he Effec of Inernal Resricion on Residuals in Linear Regression Model wih Firs-order Auoregressive Procedures Mei-Yu Lee 1 Absrac his paper demonsraes he impac of paricular facors such as a non-normal error disribuion, consrains of he residuals, sample size, he muli-collinear values of independen variables and he auocorrelaion coefficien on he disribuions of errors and residuals. his explains how residuals increasingly end o a normal disribuion wih increased linear consrains on residuals from he linear regression analysis mehod. Furhermore, reduced linear requiremens cause he shape of he error disribuion o be more clearly shown on he residuals. We find ha if he errors follow a normal disribuion, hen he residuals do as well. However, if he errors follow a U-quadraic disribuion, hen he residuals have a mixure of he error disribuion and a normal disribuion due o he ineracion of linear requiremens and sample size. hus, increasing he consrains on he residual from more independen variables causes he residuals o follow a normal 1 Deparmen of Applied Finance, Yuanpei Universiy, Hsinchu, aiwan. mylee@gm.ypu.edu.w Aricle Info: Received : June 14, Revised : July 27, Published online : Augus 15, 2014.

3 2 Inernal Resricion on Residuals disribuion, leading o a poor esimaor in he case where errors have a non-normal disribuion. Only when he sample size is large enough o eliminae he effecs of hese linear requiremens and muli-collineariy can he residuals be viewed as an esimaor of he errors. Mahemaics Subjec Classificaion: 37M10; 62M10; 37M05; 68U20 Keywords: ime series; Auoregressive model; Compuer simulaion; Non-normal disribuion 1 Inroducion I is reasonable o quesion why business sudies always use linear regression models, bu he engineering and qualiy managemen fields do no rely on such models. Researchers face a residual disribuion differen o he error disribuion, and his diverges from he resul of Box and Pierce (1970). Box and Pierce suppose ha residuals ha have a good fi should be he rue errors and can be regarded as esimaors of he errors in an auoregressive process. hus, some residuals ha are viewed as good esimaors of errors include he Durbin-Wason es saisic (Durbin and Wason, 1950, 1951) and he LaGrange Muliplier es saisic (Berusch and Pagan, 1980) in he linear regression model wih an auoregressive error process. In general, residuals are viewed as he highes represenaion of errors and are combined as an esimaor ha is used o esimae he properies of errors such as he serial correlaion of errors or he error disribuion. In he lieraure, Yule (1921) firs discusses he problem of serial correlaion, and hen Roos (1936) provides basic soluions regarding how independen variables are independen by he use of choosing lagged ime and how he rend and flucuaion can be grabbed. Box and Pierce (1970) also noe ha residual auocorrelaion can be approximaed

4 Mei-Yu Lee 3 o he linear ransformaion of error auocorrelaion and possesses a normal disribuion. herefore, he residual disribuion plays an imporan role residuals can be used o es he iniially assumed error disribuion. For example, he errors have a normal disribuion if he residuals are idenified as having a normal disribuion. However, ime series daa may have a non-normal error disribuion which conradics an assumpion in he linear regression model. hen, he auoregressive error procedure conradics he assumpion of he linear regression model regarding he independence of errors. his paper aims o explain ha he disribuions of he residuals are separae from he disribuions of he errors by using a probabiliy simulaion approach. ha is, we run a compuer simulaion wih a random number able wih a uniform disribuion beween 0 and 1 o invesigae he following hree poins: (1) if he errors are disribued normally, hen he residuals also have a normal disribuion, (2) if he errors have a non-normal disribuion, hen he residuals also have a non-normal disribuion, and (3) as he sample size becomes larger, he residual disribuion approaches o he error disribuion, ha is, he law of large numbers gradually has influence over he disribuion of he residuals. he second poin is based on an example of he U-quadraic disribued errors. If he argumen of Box and Pierce were righ, hen he residuals should have a U-quadraic disribuion, and no a normal disribuion. However, he calculaed residuals appear o have a normal disribuion, conradicing Box and Pierce. his paper presens compuer simulaion resuls showing ha he normaliy of he residuals resuls from he number of independen variables ha deermine he consrains of he residuals, X εˆ = 0, wih 1 plus he number of independen variables. he consrains of he residuals, X εˆ = 0, are referred o as he linear requiremens of he linear regression model and is number of consrains is he same as he degree of freedom (Lee, 2014a, 2014b, 2014c). he hird poin refers o he change in he residual disribuion wih a fixed number of independen variables and a U-quadraic error disribuion when he sample sizes become larger.

5 4 Inernal Resricion on Residuals he paper is srucured as follows. Secion 2 describes he model seing and simulaion procedure. Secion 3 gives he resuls of he hree cases where he error has a normal disribuion, he error has a U-quadraic disribuion, and he sample size is changed. Secion 4 concludes his paper. 2 Model and simulaion procedures Consider a linear regression model wih k independen variables and sample sizes, as Y = X β + ε ( 1) ( k ) ( k 1) ( 1) on he condiions of E( ε) = 0, E( X ε) = 0 procedure is and he firs-order auoregressive ε + +1 = ρε µ where = 1,2,, 1 and < 1 Yˆ = Xβˆ, hen he residuals are εˆ = ε X X 1 ˆ = and ρ. he esimaors are β ( X X) X Y 1 ( X ) ε, 1 ( X) X ε = I X( X X) and linear requiremen are X εˆ = 0 where Y = Yˆ + εˆ (Balagi 2011). When Yˆ is normal disribuion, he regression coefficien poin esimaors also are normal disribuion and εˆ is normal disribuion, Ŷ approximaes o normal disribuion, if k is enough large. We calculae he mean-square-error (MSE) as ˆ = from ( ε) 0 ˆˆ MSE = ( Yˆ Xβˆ ) ( Yˆ Xβˆ ) E and E( ε ε ) = ( X X) MSE ˆ ˆ 1. hus, he wo main facors ha

6 Mei-Yu Lee 5 affec he probabiliy disribuion of he residuals are he assumpion of he error disribuion and he linear requiremen of X εˆ = 0, which has k+1 consrains. Moreover, he sample size, muli-collinear values of independen variables and he auocorrelaion coefficien affec he disribuion of he residuals. However, as he disribuion of he residuals in he above model is difficul o formulae, his paper only uses compuer simulaions o show how facors affec he relaion beween he errors and he residuals. 2.1 he simulaor mehod he sampling disribuion of a es saisic may or may no be known. In paricular, some sample disribuion of he es saisic canno be ransferred using radiional mahemaical echniques such as calculus or Mone Carlo mehods. he concep of he Mone Carlo mehod is a good simulaion mehod, bu he use of coninuous daa isn possible in our compuer program. o manage he coninuous normal disribuion (or U-quadraic disribuion), we run a compuer simulaion using a sofware program ha can work wih any probabiliy disribuion ransformaion. he probabiliy heory of his paper has been creaed using he basic concep of a probabiliy disribuion simulaor, and he funcions of coninuous random variables can be ransferred from a uniform disribuion wih parameers of 0 and 1, X ~ U(0,1). hus, we can generae he daa and compue he coefficiens and images. 2 he compuer simulaion is based on he following seps. 2 he sofware program is named as Whie model I ha can be download from hp://goo.gl/oudpsp. he disribuions of he h error and residual, he disribuions of sum of he errors and residuals can be simulaed by he sofware. he simulaion echnology is from C.C.C. Ld. (hp://psccc.com.w/en/produc). And he U-quadraic disribuion formula can reference a hp://psccc.com.w/uploads/files/probabiliy/1/chaper_one_02.pdf.

7 6 Inernal Resricion on Residuals (1) Generae daa from a random number able of U( 0,1). Each value of he U-quadraic disribuion can be obained when he value is from he inverse funcion of a Normal disribuion (or U-quadraic disribuion). (2) Collec he values whose number maches he number of errors. Clearly, hese values are i.i.d. Normal disribuion (or U-quadraic disribuion). ha is he se values of µ from he serial correlaion model where E 2 ( µ ) = 0, Var( µ ) = σ, 1,...,, = and or U-quadraic disribuion, ε = ρ ε + µ, = 1,2,..., 1, ρ < When ρ is known, ε + 1 can be found. (3) he residuals follow he poin esimae requiremens of he linear model. Y = β 0 + β1 X 1, β k X k, + ε, 1,...,, = If he number and values of independen variables are known, hen ˆ β ˆ β,..., ˆ 0, 1 β k can be esimaed, whereas εˆ is consrained by X εˆ = 0. Meanwhile, he simulaor obeys he linear model mehod and creaes he residual values ( εˆ ). hus, he simulaion process includes: Sep 1: Giving he inercep and slope value, and he daa se of independen variables. Sep 2: Using he simulaion mehod o ge he error daa se from he probabiliy disribuion wih sample size,. Sep 3: According o he linear model, compuing he daa se of dependen variables: Y = Xβ + ε.

8 Mei-Yu Lee 7 Sep 4: Calculaing he poin esimaor values of he regression coefficien and geing he esimaed values of dependen variables: Yˆ = Xβˆ. Sep 5: Calculaing he daa se of residuals. εˆ = Y Xβˆ. ε is simulaed by he Sep 1 and 2, and εˆ is simulaed by he Sep 1 o 5. ε and εˆ are simulaed oally imes, and hen generaed values, respecively, o form he frequency disribuion ha can reach o he real ε and εˆ disribuions. 3 Main Resuls 3.1 he errors follow normal disribuion Secion 3.1 discusses he condiion ha he disribuion of he errors is a normal disribuion wih 6 independen variables, 15 samples and he auocorrelaion coefficien of he errors is zero. hus, he firs column of Figure 1 is he disribuion of he errors, which is a sandard normal disribuion, and he second column of Figure 1 illusraes he shape and coefficien of he 7 h residual disribuion where he coefficiens of he mean, skewness, and kurosis represen a normal disribuion and are he same as he error disribuion. he residuals can be viewed as an esimaor of he errors because Figure 1 guaranees ha he disribuion of he residuals is he same as he disribuion of he errors. hus, Figure 1 suppors he conclusion of Box and Pierce (1970) when he errors are normally disribued. he reason is as follows: he residuals are a combinaion of he errors in he linear regression model, ha is ( β) ε ε ˆ = Y Yˆ = X β ˆ +, and he auoregressive procedure akes he errors oward he auocorrelaion wih each oher wihou changing he propery beween he residuals and he errors. A he same ime, he addiive propery of a normal

9 8 Inernal Resricion on Residuals disribuion impacs he linear combinaion of he errors. hus, he residuals show he normally disribued propery of he errors. f(w1),f(w1), W1=error(1) f(w17),f(w17),w17=residual(7) Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Figure 1: he errors and residuals disribuions as he errors are normal disribuion Anoher reason is ha he mahemaical formula of he normal disribuion includes sin and cos funcions, which are cyclical funcions. herefore, he residuals follow he normal disribuion when he error disribuion is a normal disribuion. On he oher hand, if he error disribuion has no properies of sin or cos funcions, such as logisic, uniform, U-quadraic, or an exponenial disribuion, hen he errors have no cyclical propery and may no be shaped as per he normal disribuion, bu like oher disribuions insead. his paper confirms he proposiion below. Proposiion 1 he residuals are normally disribued when he errors are normally disribued because a normal disribuion wih a sin and cos funcion form is subjeced o he addiive propery.

10 Mei-Yu Lee 9 However, i is hard o promise ha he errors will always be normally disribued when researchers do no have he populaion daa. Samples should always be esed o classify wha disribuions hey follow. 3.2 he errors follow U-quadraic disribuion he paper gives an example of non-normal error disribuion ha he errors, were U-quadraic disribuion, hen he compuer simulaion offers evidence abou he disribuions of he errors and residuals on he condiion of 6 independen variables, 15 samples, he 1 lagged period, variance of error is 1 and he auocorrelaion coefficien of he errors is zero. Lee (2014c) supposed ha he values of independen variables have serious impac on he residuals, hus he paper simulaneously discusses he disribuions of he residuals a wo cases where he values of independen variables are separaely wih low and high muli-collinearliy in he linear regression model wih he firs-order auoregressive error procedure. he hird column of Figure 2 presens he disribuion of he firs residual ha is generaed from he independen variables wih he populaion correlaion coefficien is Figure 2 shows ha he firs column is he error disribuion which is U-quadraic disribuion, he second column is he disribuion of he firs residual wih low muli-collineariy and he hird column is he disribuion of he firs residual wih high muli-collineariy. he disribuions of he firs residual in he second and hird columns are as similar as normal disribuion while he disribuion of he errors is U-quadraic disribuion. he residual disribuions in Figure 2 are differen from he error disribuion, hus, he residuals canno be regarded as an esimaor of he errors when he errors are non-normal disribuion, moreover, he serial correlaion es for auocorrelaion of he errors is no suiable o use he mahemaic combinaion of he residuals because he difference beween he µ,

11 10 Inernal Resricion on Residuals disribuions of he errors and residuals. he researchers should firs invesigae he disribuion of he daa o classify wha disribuion he daa is or hey always obain he resul ha he errors follow normal disribuion from he residuals. f(w1),f(w1), W1=error(1) f(w11),f(w11),w11=residual(1) wih low muli-collineariy f(w11),f(w11),w11=residual(1) wih high muli-collineariy Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Figure 2: he error and residual disribuions when error follows U-quadraic disribuion (=15) he main difference beween Figure 1 and Figure 2 is he assumpion of he error disribuion, bu he residuals show a normal disribuion in Figure 2. here mus be some special facors no ye discovered, independen of he error disribuion which has he propery of sin and cos funcions or he propery of addiion from he normal disribuion. Lee (2014b) discovered ha he number of independen variables from 1 o 6 causes he residual disribuion o end owards a normal disribuion in he linear regression model wih an auoregressive procedure. herefore, he paper supposes ha errors are no resriced, bu ha he residuals are resriced by X εˆ = 0 which has k+1 consrains. his linear requiremen of he residuals in he linear regression model disors he error disribuion away

12 Mei-Yu Lee 11 from he original disribuion o a normal disribuion. he linear requiremen of he residuals can produce wo sources in ANOVA able when he number of independen variables is k and he sample size is. he firs source is he sum of squares in regression (SSR) whose degree of freedom is k. he second source is he sum of squares of error (SSE) whose degree of freedom is -k-1. hese wo sources decide he shape of he residual disribuion. As k becomes larger, he residuals are resriced by more consrains, similar o a linear combinaion of random variables. If more random variables are added o he linear combinaion, he new random variable ends oward he normal disribuion. he residuals ha are resriced by more consrains also represen similar saes such ha he disribuion of he residuals ends oward a normal disribuion. If we consider he collinear values of independen variables, hen he collinear effec affecs he convergence o a normal disribuion when he error disribuion is a U-quadraic disribuion. he hree columns of Figure 2 show ha he residuals wih convergence o normaliy are weakened by high muli-collinear values of independen variables when he number of independen variables is fixed. he coefficiens of Figure 2 express he variance of he residuals, S ( εˆ ). MSE increases as he values of independen variables change from low o high muli-collineariy. he shape of he disribuion in he hird column of Figure 2 shows ha he disribuion of he firs residual has a fla region around is mean and has more lef-skewness and less cenralizaion. hus, he high muli-collineariy causes a big problem wih he residuals, and so he disribuion of he firs residual is no he same as he shape of he disribuion in he second column of Figure 2. he muli-collineariy of he independen variables resuls in he calculaion and combinaion of he residuals becoming more complex. Moreover, he variance and disribuion of he residuals are boh disurbed by he collineariy when he errors follow a non-normal disribuion. When he errors are asymmeric, he residuals

13 12 Inernal Resricion on Residuals are delayed in heir convergence o normaliy. o address his, he linear regression model adds more independen variables o accelerae he residuals convergence o normaliy because more and more independen variables can resric he residuals and weaken he collinear effec. Meanwhile, he degree of freedom also decreases. 3 his paper proposes a second proposiion as follows. Proposiion 2. (1) A larger number of independen variables, k, brings faser convergence o normaliy on he residuals when he errors have a non-normal disribuion. (2) Higher muli-collinear values of independen variables cause slower convergence o normaliy for he residuals. (3) Higher muli-collinear values and a smaller number of independen variables cause he disribuion of he residuals o follow neiher a non-normal disribuion nor he error disribuion, bu a mixure of he error disribuion and a normal disribuion according o he consrains of he residuals. he muli-collinear propery and consrains on he residuals have opposie effecs, simulaneously disurbing he disribuion of he residuals which becomes a mixed disribuion beween a normal disribuion and he error disribuion. 3.3 Sample size is changed he sample size effec plays a very imporan role in ime series models ha can represen he law of large numbers. hus, he residuals may be regarded as a good esimaor of he errors. he simulaion case only changes he sample size from 9 3 In fac, he mulicollinear case of he compuer simulaion implies ha only he number of independen variables is more han 20 and he degree of freedom is more han 2, hen he residuals and coefficiens of regression model wih firs-order auoregressive error process will have normal disribued poin esimaors when he low mulicollieariy of independen variables exiss and he errors follow independenly idenically disribuion wih symmeric a zero.

14 Mei-Yu Lee 13 o 107 on he condiion of 6 independen variables. he populaion variance is 1, he auocorrelaion coefficien is zero, he values of independen variables are from he fron of he daa se, and he simulaed seing follows secion 3.2. he residuals gradually revealed he propery of he errors when he sample sizes increased from 9 o 107. Wihou loss of generalizaion, his paper only shows he shapes and coefficiens of disribuions for he able -1, which are in Appendix A. h 0.5 residual in Figure 3 and =14, he 7 h residual =60, he 30 h residual Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Figure 3: he h 0.5 residual s disribuion h Figure 3 explores he disribuions of he 0.5 residual a = 14 and 60. he second column of Figure 3 has a less regular disribuion shape han he firs column a = 14. Meanwhile, he shape of he 30 h residual a = 60 ends oward a U-quadraic disribuion more han he disribuion of he 7 h residual a = 14. h Figure 3 clearly shows ha he larger sample sizes cause he 0.5 residual o end owards he error disribuion. hus, he sample size effec can limi he effec

15 14 Inernal Resricion on Residuals of he linear requiremens on he residuals. A he same ime, i allows he disribuion of he residuals o represen more properies of he errors. he coefficiens in able A-1 show ha he h 0.5 residual has a mean of around zero, skewed coefficiens, larger variances, and decreasing kurosis when he sample sizes increased from 9 o 107. In comparison wih he firs column in Figure 2, he coefficiens in able A-1 approach he coefficiens in he firs column in Figure 2 when he sample sizes increased, especially for he variance and kurosis coefficiens. he coefficiens ha change wih differen sample sizes also show ha he disribuion of he residuals can be an esimaor of he errors if he sample size effec is sufficienly larger han he effec of consrains on he residuals. he ineracion beween he sample size and he linear requiremen h causes he differen shape and coefficiens of he disribuion of he 0.5 residual. he residuals become more cenralized when he linear requiremen has more consrains. Neverheless, he residuals are affeced by he error disribuion and are more likely o end oward he error disribuion when he number of consrains is no large enough. In oher words, he residual disribuion is vasly differen from he error disribuion because he residuals are affeced by he linear requiremen of linear regression models, he error disribuion, sample sizes, and he values of independen variables. I is generally assumed ha errors are normally disribued, which is a symmeric disribuion. However, he residuals can be an esimaor of he errors when k+1 is greaer han 20, ranges from 23 o 23+(k-19), and here is lile or no muli-collineariy. 3.4 Auocorrelaion coefficien of he errors is 0.7 he above saemens are on he condiion of zero auocorrelaion of he errors. However, nonzero auocorrelaion cases are usually seen in he daa. We discuss he case where he auocorrelaion coefficien of he errors is 0.7 and oher

16 Mei-Yu Lee 15 condiions are as he same as secion 3.2. here are also wo cases of low and high muli-collinear values of independen variables. We show he shapes and coefficiens of disribuion of firs and 15 h errors and residuals in Figure 4 and Figure 5, respecively. he 1s error he 1s residual Low muli-collineariy he 1s residual High muli-collineariy Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Figure 4: he firs error and residual disribuions when error follows U-quadraic disribuion and auocorrelaion coefficien is 0.7 he disribuion of he firs error follows a U-quadraic disribuion, bu he firs residuals ha are affeced by he degree of muli-collineariy show differen shapes for he disribuion of he residuals in Figure 4. he disribuion of he firs residuals canno represen he propery of he firs error and greaer muli-collinear values of independen variables cause more cenralized disribuions of he firs residuals. Anoher imporan view from coefficiens in Figure 4 is ha he coefficiens are similar o each oher among he hree columns, bu he populaion disribuion of he firs column and sampling disribuions of he second and hird columns have compleely differen shapes from each oher. On he oher hand, he coefficiens

17 16 Inernal Resricion on Residuals may highligh quesionable resuls ha are used o pass hypohesis esing when researchers only invesigae he means and variances of residuals. Comparing he hree columns, he linear requiremen and muli-collineariy push he firs residual owards cenralizaion. Higher muli-collinear values of independen variables induce more cenralized residuals. his sudy also ran he disribuions of he 15 h error and residuals which are divided ino wo pars including low and high muli-collineariy in Figure 5. 4 he disribuion of he 15 h error is no he same as he U-quadraic disribuion, bu is more cenralized as per he normal disribuion because of he nonzero auocorrelaion coefficien. Meanwhile, he disribuions of he 15 h residual are more similar o he disribuion of he 15 h error in Figure 5 han he firs residual in Figure 4. he means and variances a he firs and hird columns in Figure 5 highligh ha he 15 h residual is very similar o he 15 h error. However, he diagrams and kurosis coefficiens show ha here is a vas difference beween he 15 h residual and he 15 h error. Comparing Figure 4 and Figure 5, he coefficiens of he errors show similar coefficiens, excep for he kurosis coefficiens. he diagrams of he firs column from Figure 4 and Figure 5 are from he U-quadraic shape o he opposie of he U-quadraic shape because of he non-zero auocorrelaion coefficien. he diagrams of low muli-collinear residuals show bulging from he firs residual o he 15 h residual, and he diagrams of high muli-collinear residuals are similar. his highlighs ha he nonzero auocorrelaion coefficien and muli-collineariy inerac on he residuals. Higher muli-collineariy decreases he effec of he auocorrelaion coefficien on he residuals 4 We also simulae he case of -0.7 auocorrelaion coefficien in Appendix B.

18 Mei-Yu Lee 17 he 15 h error he 15 h residual Low muli-collineariy he 15 h residual High muli-collineariy Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : Figure 5: he 15h error and residual disribuions when error follows U-quadraic disribuion and auocorrelaion coefficien is Conclusion he purpose of his paper is o explain why residuals canno perfecly represen errors. his paper confirms ha residuals ha require specific condiions can be viewed as an esimaor of he errors, bu i may no be appropriae o assume a normal disribuion because of he properies of he daa. However, if he errors are normally disribued, he residuals can be a good esimaor of he errors because of he properies of a normal disribuion. his paper also shows, via compuer simulaion resuls, how a non-normal disribuion, linear requiremens, muli-collineariy, sample size, and he auocorrelaion coefficien affec he disribuions of he errors and he residuals. Firs, residuals are resriced, bu errors are no. Hence, he values of he residuals are consrained by a linear requiremen from he regression analysis, and his causes he residuals o no perfecly represen errors. Second, his paper supposes ha he linear requiremen, sample sizes, and he values of independen variables

19 18 Inernal Resricion on Residuals inerac in he linear regression model. hus: (1) he consrains of he linear requiremens are high enough ha he residuals follow a normal disribuion when he sample sizes are fixed, regardless of he error erm assumed. (2) If he error erm is assumed o follow a normal disribuion, hen he residuals follow a normal disribuion. (3) Larger sample sizes resul in he residuals revealing properies of he error erm when he linear requiremen is fixed. References [1] B.H. Balagi, Economerics, Fifh ediion, Springer: New York, [2] C.F. Roos, Annual Survey of Saisical echniques: he Correlaion and Analysis of ime Series--Par II, Economerica, 4(4), (1936), [3] G.E. Box, D.A. Pierce, Disribuion of Residual Auocorrelaions in Auoregressive-inegraed Moving Average ime Series Models, Journal of he American Saisical Associaion, 65(332), (1970), [4] G.U. Yule, On he ime-correlaion Problem, wih Especial Reference o he Variae-Difference Correlaion Mehod, Journal of he Royal Saisical Sociey, 84(4), (1921), [5] J. Durbin, G.S. Wason, esing for Serial Correlaion in Leas Squares Regression: I, Biomerika, 37(3/4), (1950), [6] J. Durbin, G.S. Wason, esing for Serial Correlaion in Leas Squares Regression. II, Biomerika, 38(1/2), (1951), [7] M.Y. Lee, he Paern of R-Square in Linear Regression Model wih Firs-Order Auoregressive Error Process and Bayesian propery: Compuer Simulaion, Journal of Accouning & Finance Managemen Sraegy, 9(1), (2014a). [8] M.Y. Lee, Limiing Propery of Durbin-Wason es Saisic, manuscrip, (2014b).

20 Mei-Yu Lee 19 [9] M.Y. Lee, he Conflic of Residual and Error Simulaed in Linear Regression Model wih AR(1) Error Process, manuscrip, (2014c). [10].S. Breusch, L.G. Godfrey, A review of recen work on esing for auocorrelaion in dynamic economic models," Universiy of Souhampon, (1980).

21 20 Inernal Resricion on Residuals Appendix A. Secion 3.3 shows how he change of sample size affecs he disribuion of each residual. However, each differen sample size case has residuals, he paper shows he coefficiens of he 0.5 h residual whose 0.5 h is half. he coefficiens of able A-1 almos have he same means of he 0.5 h residual and canno le readers know wha he difference among hose disribuions of he 0.5 h residual from differen sample sizes. hus, he paper pu he coefficiens of he 0.5 h residual in Appendix and he graphs of disribuion of he 0.5 h residual in he main conen. able A-1. he coefficiens of he 0.5 h residual in differen sample sizes h Coefficiens Mahemaical Mean: Variance : S.D. : Skewed Coef.: Kurosis Coef.: MAD : Range : Median : IQR : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : MAD : Range : Median : IQR : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : MAD : Range : Median : IQR : h Coefficiens Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef.: MAD : Range : Median : IQR : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef.: MAD : Range : Median : IQR : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : MAD : Range : Median : IQR : h Coefficiens Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : MAD : Range : Median : IQR : Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef. : MAD : Range : Median : IQR : h Coefficiens Mahemaical Mean: Variance : S.D. : Skewed Coef. : Kurosis Coef.: MAD : Range : Median : IQR : Mahemaical Mean: Variance : S.D. : Skewed Coef.: Kurosis Coef.: MAD : Range : Median : IQR : Mahemaical Mean: Variance : S.D. : Skewed Coef.: Kurosis Coef.: MAD : Range : Median : IQR :

22 Mei-Yu Lee 21 Appendix B. he paper also simulaes he siuaion ha he auocorrelaion of he errors is he residuals have he smaller mean and variance and more negaive skewness on he condiion of high mulicollineariy. he 1 s error he 1 s residual Low muli-collineariy he 1 s residual High muli-collineariy Figure B-1. Mahemaical Mean: Mahemaical Mean: Mahemaical Mean: Variance : Variance : Variance : S.D. : S.D. : S.D. : Skewed Coef. : Skewed Coef. : Skewed Coef. : Kurosis Coef. : Kurosis Coef. : Kurosis Coef. : he firs error and residual disribuions when error follows U-quadraic disribuion and auocorrelaion coefficien is -0.7 Wih comparison of he auocorrelaion coefficiens of he errors ha are 0.7 and -0.7, he means of he firs error and he 15 h error have less means when auocorrelaion coefficien is Second, he probabiliy disribuion of he firs error in Figure B-1 is as he same as in he lef side of Figure 4, so does he 15 h error in Figure B-2 and Figure 5. hird, he low muli-collineariy siuaion shows ha he firs residual has larger variance and less cenralizaion in Figure B-1 han in Figure 4. However, he 15 h residual is more cenralized in Figure B-2 han in Figure 5. Finally, he high muli-collineariy siuaion explores ha he firs residual in Figure B-1 has larger variance and more cenralizaion han in Figure 4,. However, he 15 h residual has smaller variance and less cenralizaion in Figure B-2 han in Figure 5.

23 22 Inernal Resricion on Residuals he 15 h error he 15 h residual Low muli-collineariy he 15 h residual High muli-collineariy Figure B-2. Mahemaical Mean: Mahemaical Mean: Mahemaical Mean: Variance : Variance : Variance : S.D. : S.D. : S.D. : Skewed Coef. : Skewed Coef. : Skewed Coef. : Kurosis Coef. : Kurosis Coef. : Kurosis Coef. : he 15 h error and residual disribuions when error follows U-quadraic disribuion and auocorrelaion coefficien is -0.7

The Effect of Nonzero Autocorrelation Coefficients on the Distributions of Durbin-Watson Test Estimator: Three Autoregressive Models

The Effect of Nonzero Autocorrelation Coefficients on the Distributions of Durbin-Watson Test Estimator: Three Autoregressive Models EJ Exper Journal of Economi c s ( 4 ), 85-9 9 4 Th e Au h or. Publi sh ed by Sp rin In v esify. ISS N 3 5 9-7 7 4 Econ omics.e xp erjou rn a ls.com The Effec of Nonzero Auocorrelaion Coefficiens on 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

DEPARTMENT OF STATISTICS

DEPARTMENT OF STATISTICS A Tes for Mulivariae ARCH Effecs R. Sco Hacker and Abdulnasser Haemi-J 004: DEPARTMENT OF STATISTICS S-0 07 LUND SWEDEN A Tes for Mulivariae ARCH Effecs R. Sco Hacker Jönköping Inernaional Business School

More information

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

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

More information

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

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

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In

More information

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

Vehicle Arrival Models : Headway

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

More information

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model:

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model: Dynamic Economeric Models: A. Auoregressive Model: Y = + 0 X 1 Y -1 + 2 Y -2 + k Y -k + e (Wih lagged dependen variable(s) on he RHS) B. Disribued-lag Model: Y = + 0 X + 1 X -1 + 2 X -2 + + k X -k + e

More information

GMM - Generalized Method of Moments

GMM - Generalized Method of Moments GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

4.1 Other Interpretations of Ridge Regression

4.1 Other Interpretations of Ridge Regression CHAPTER 4 FURTHER RIDGE THEORY 4. Oher Inerpreaions of Ridge Regression In his secion we will presen hree inerpreaions for he use of ridge regression. The firs one is analogous o Hoerl and Kennard reasoning

More information

Wednesday, November 7 Handout: Heteroskedasticity

Wednesday, November 7 Handout: Heteroskedasticity Amhers College Deparmen of Economics Economics 360 Fall 202 Wednesday, November 7 Handou: Heeroskedasiciy Preview Review o Regression Model o Sandard Ordinary Leas Squares (OLS) Premises o Esimaion Procedures

More information

How to Deal with Structural Breaks in Practical Cointegration Analysis

How to Deal with Structural Breaks in Practical Cointegration Analysis How o Deal wih Srucural Breaks in Pracical Coinegraion Analysis Roselyne Joyeux * School of Economic and Financial Sudies Macquarie Universiy December 00 ABSTRACT In his noe we consider he reamen of srucural

More information

The General Linear Test in the Ridge Regression

The General Linear Test in the Ridge Regression ommunicaions for Saisical Applicaions Mehods 2014, Vol. 21, No. 4, 297 307 DOI: hp://dx.doi.org/10.5351/sam.2014.21.4.297 Prin ISSN 2287-7843 / Online ISSN 2383-4757 The General Linear Tes in he Ridge

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

Econ Autocorrelation. Sanjaya DeSilva

Econ Autocorrelation. Sanjaya DeSilva Econ 39 - Auocorrelaion Sanjaya DeSilva Ocober 3, 008 1 Definiion Auocorrelaion (or serial correlaion) occurs when he error erm of one observaion is correlaed wih he error erm of any oher observaion. This

More information

Stationary Time Series

Stationary Time Series 3-Jul-3 Time Series Analysis Assoc. Prof. Dr. Sevap Kesel July 03 Saionary Time Series Sricly saionary process: If he oin dis. of is he same as he oin dis. of ( X,... X n) ( X h,... X nh) Weakly Saionary

More information

Generalized Least Squares

Generalized Least Squares Generalized Leas Squares Augus 006 1 Modified Model Original assumpions: 1 Specificaion: y = Xβ + ε (1) Eε =0 3 EX 0 ε =0 4 Eεε 0 = σ I In his secion, we consider relaxing assumpion (4) Insead, assume

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

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

More information

Testing the Random Walk Model. i.i.d. ( ) r

Testing the Random Walk Model. i.i.d. ( ) r he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp

More information

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

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

More information

A unit root test based on smooth transitions and nonlinear adjustment

A unit root test based on smooth transitions and nonlinear adjustment MPRA Munich Personal RePEc Archive A uni roo es based on smooh ransiions and nonlinear adjusmen Aycan Hepsag Isanbul Universiy 5 Ocober 2017 Online a hps://mpra.ub.uni-muenchen.de/81788/ MPRA Paper No.

More information

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Asymmery and Leverage in Condiional Volailiy Models Michael McAleer WORKING PAPER

More information

Regression with Time Series Data

Regression with Time Series Data Regression wih Time Series Daa y = β 0 + β 1 x 1 +...+ β k x k + u Serial Correlaion and Heeroskedasiciy Time Series - Serial Correlaion and Heeroskedasiciy 1 Serially Correlaed Errors: Consequences Wih

More information

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H. ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models

More information

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Y 0.4Y 0.45Y Y to a proper ARMA specification. HG Jan 04 ECON 50 Exercises II - 0 Feb 04 (wih answers Exercise. Read secion 8 in lecure noes 3 (LN3 on he common facor problem in ARMA-processes. Consider he following process Y 0.4Y 0.45Y 0.5 ( where

More information

1. Diagnostic (Misspeci cation) Tests: Testing the Assumptions

1. Diagnostic (Misspeci cation) Tests: Testing the Assumptions Business School, Brunel Universiy MSc. EC5501/5509 Modelling Financial Decisions and Markes/Inroducion o Quaniaive Mehods Prof. Menelaos Karanasos (Room SS269, el. 01895265284) Lecure Noes 6 1. Diagnosic

More information

The Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form

The Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form Chaper 6 The Simple Linear Regression Model: Reporing he Resuls and Choosing he Funcional Form To complee he analysis of he simple linear regression model, in his chaper we will consider how o measure

More information

Time series Decomposition method

Time series Decomposition method Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,

More information

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H. ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen

More information

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model

More information

Forecasting optimally

Forecasting optimally I) ile: Forecas Evaluaion II) Conens: Evaluaing forecass, properies of opimal forecass, esing properies of opimal forecass, saisical comparison of forecas accuracy III) Documenaion: - Diebold, Francis

More information

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange

More information

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.

More information

Dynamic Models, Autocorrelation and Forecasting

Dynamic Models, Autocorrelation and Forecasting ECON 4551 Economerics II Memorial Universiy of Newfoundland Dynamic Models, Auocorrelaion and Forecasing Adaped from Vera Tabakova s noes 9.1 Inroducion 9.2 Lags in he Error Term: Auocorrelaion 9.3 Esimaing

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics

More information

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013 Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Distribution of Least Squares

Distribution of Least Squares Disribuion of Leas Squares In classic regression, if he errors are iid normal, and independen of he regressors, hen he leas squares esimaes have an exac normal disribuion, no jus asympoic his is no rue

More information

GINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad,

GINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad, GINI MEAN DIFFEENCE AND EWMA CHATS Muhammad iaz, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. E-Mail: riaz76qau@yahoo.com Saddam Akbar Abbasi, Deparmen of Saisics, Quaid-e-Azam Universiy

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

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

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

More information

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance Lecure 5 Time series: ECM Bernardina Algieri Deparmen Economics, Saisics and Finance Conens Time Series Modelling Coinegraion Error Correcion Model Two Seps, Engle-Granger procedure Error Correcion Model

More information

Solutions to Odd Number Exercises in Chapter 6

Solutions to Odd Number Exercises in Chapter 6 1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Solutions: Wednesday, November 14

Solutions: Wednesday, November 14 Amhers College Deparmen of Economics Economics 360 Fall 2012 Soluions: Wednesday, November 14 Judicial Daa: Cross secion daa of judicial and economic saisics for he fify saes in 2000. JudExp CrimesAll

More information

Cointegration and Implications for Forecasting

Cointegration and Implications for Forecasting Coinegraion and Implicaions for Forecasing Two examples (A) Y Y 1 1 1 2 (B) Y 0.3 0.9 1 1 2 Example B: Coinegraion Y and coinegraed wih coinegraing vecor [1, 0.9] because Y 0.9 0.3 is a saionary process

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims Problem Se 5 Graduae Macro II, Spring 2017 The Universiy of Nore Dame Professor Sims Insrucions: You may consul wih oher members of he class, bu please make sure o urn in your own work. Where applicable,

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

1 Differential Equation Investigations using Customizable

1 Differential Equation Investigations using Customizable Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for

More information

A New Unit Root Test against Asymmetric ESTAR Nonlinearity with Smooth Breaks

A New Unit Root Test against Asymmetric ESTAR Nonlinearity with Smooth Breaks Iran. Econ. Rev. Vol., No., 08. pp. 5-6 A New Uni Roo es agains Asymmeric ESAR Nonlineariy wih Smooh Breaks Omid Ranjbar*, sangyao Chang, Zahra (Mila) Elmi 3, Chien-Chiang Lee 4 Received: December 7, 06

More information

CHAPTER 17: DYNAMIC ECONOMETRIC MODELS: AUTOREGRESSIVE AND DISTRIBUTED-LAG MODELS

CHAPTER 17: DYNAMIC ECONOMETRIC MODELS: AUTOREGRESSIVE AND DISTRIBUTED-LAG MODELS Basic Economerics, Gujarai and Porer CHAPTER 7: DYNAMIC ECONOMETRIC MODELS: AUTOREGRESSIVE AND DISTRIBUTED-LAG MODELS 7. (a) False. Economeric models are dynamic if hey porray he ime pah of he dependen

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

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

More information

Measurement Error 1: Consequences Page 1. Definitions. For two variables, X and Y, the following hold: Expectation, or Mean, of X.

Measurement Error 1: Consequences Page 1. Definitions. For two variables, X and Y, the following hold: Expectation, or Mean, of X. Measuremen Error 1: Consequences of Measuremen Error Richard Williams, Universiy of Nore Dame, hps://www3.nd.edu/~rwilliam/ Las revised January 1, 015 Definiions. For wo variables, X and Y, he following

More information

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,

More information

Properties of Autocorrelated Processes Economics 30331

Properties of Autocorrelated Processes Economics 30331 Properies of Auocorrelaed Processes Economics 3033 Bill Evans Fall 05 Suppose we have ime series daa series labeled as where =,,3, T (he final period) Some examples are he dail closing price of he S&500,

More information

Appendix to Creating Work Breaks From Available Idleness

Appendix to Creating Work Breaks From Available Idleness Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember

More information

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1 Chaper 5 Heerocedasic Models Inroducion o ime series (2008) 1 Chaper 5. Conens. 5.1. The ARCH model. 5.2. The GARCH model. 5.3. The exponenial GARCH model. 5.4. The CHARMA model. 5.5. Random coefficien

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

OBJECTIVES OF TIME SERIES ANALYSIS

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

More information

The equation to any straight line can be expressed in the form:

The equation to any straight line can be expressed in the form: Sring Graphs Par 1 Answers 1 TI-Nspire Invesigaion Suden min Aims Deermine a series of equaions of sraigh lines o form a paern similar o ha formed by he cables on he Jerusalem Chords Bridge. Deermine he

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If

More information

Empirical Process Theory

Empirical Process Theory Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued

More information

FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA

FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA N. Okendro Singh Associae Professor (Ag. Sa.), College of Agriculure, Cenral Agriculural Universiy, Iroisemba 795 004, Imphal, Manipur

More information

Modeling and Forecasting Volatility Autoregressive Conditional Heteroskedasticity Models. Economic Forecasting Anthony Tay Slide 1

Modeling and Forecasting Volatility Autoregressive Conditional Heteroskedasticity Models. Economic Forecasting Anthony Tay Slide 1 Modeling and Forecasing Volailiy Auoregressive Condiional Heeroskedasiciy Models Anhony Tay Slide 1 smpl @all line(m) sii dl_sii S TII D L _ S TII 4,000. 3,000.1.0,000 -.1 1,000 -. 0 86 88 90 9 94 96 98

More information

Department of Economics East Carolina University Greenville, NC Phone: Fax:

Department of Economics East Carolina University Greenville, NC Phone: Fax: March 3, 999 Time Series Evidence on Wheher Adjusmen o Long-Run Equilibrium is Asymmeric Philip Rohman Eas Carolina Universiy Absrac The Enders and Granger (998) uni-roo es agains saionary alernaives wih

More information

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

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

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j = 1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

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

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

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

Failure of the work-hamiltonian connection for free energy calculations. Abstract

Failure of the work-hamiltonian connection for free energy calculations. Abstract Failure of he work-hamilonian connecion for free energy calculaions Jose M. G. Vilar 1 and J. Miguel Rubi 1 Compuaional Biology Program, Memorial Sloan-Keering Cancer Cener, 175 York Avenue, New York,

More information

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Stochastic Model for Cancer Cell Growth through Single Forward Mutation Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com

More information

Hypothesis Testing in the Classical Normal Linear Regression Model. 1. Components of Hypothesis Tests

Hypothesis Testing in the Classical Normal Linear Regression Model. 1. Components of Hypothesis Tests ECONOMICS 35* -- NOTE 8 M.G. Abbo ECON 35* -- NOTE 8 Hypohesis Tesing in he Classical Normal Linear Regression Model. Componens of Hypohesis Tess. A esable hypohesis, which consiss of wo pars: Par : a

More information

Analyze patterns and relationships. 3. Generate two numerical patterns using AC

Analyze patterns and relationships. 3. Generate two numerical patterns using AC envision ah 2.0 5h Grade ah Curriculum Quarer 1 Quarer 2 Quarer 3 Quarer 4 andards: =ajor =upporing =Addiional Firs 30 Day 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 andards: Operaions and Algebraic Thinking

More information

Modeling Economic Time Series with Stochastic Linear Difference Equations

Modeling Economic Time Series with Stochastic Linear Difference Equations A. Thiemer, SLDG.mcd, 6..6 FH-Kiel Universiy of Applied Sciences Prof. Dr. Andreas Thiemer e-mail: andreas.hiemer@fh-kiel.de Modeling Economic Time Series wih Sochasic Linear Difference Equaions Summary:

More information

(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum.

(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum. January 01 Final Exam Quesions: Mark W. Wason (Poins/Minues are given in Parenheses) (15) 1. Suppose ha y follows he saionary AR(1) process y = y 1 +, where = 0.5 and ~ iid(0,1). Le x = (y + y 1 )/. (11)

More information

Stat 601 The Design of Experiments

Stat 601 The Design of Experiments Sa 601 The Design of Experimens Yuqing Xu Deparmen of Saisics Universiy of Wisconsin Madison, WI 53706, USA December 1, 2016 Yuqing Xu (UW-Madison) Sa 601 Week 12 December 1, 2016 1 / 17 Lain Squares Definiion

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Chapter 4. Truncation Errors

Chapter 4. Truncation Errors Chaper 4. Truncaion Errors and he Taylor Series Truncaion Errors and he Taylor Series Non-elemenary funcions such as rigonomeric, eponenial, and ohers are epressed in an approimae fashion using Taylor

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

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

More information