See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/303842450

Time Series using Stata (Oscar Torres-Reyna version)

Code · December 2007

CITATIONS READS

0 2,942

1 author:

Aymen Ammari
INSEEC
21 PUBLICATIONS   86 CITATIONS   

SEE PROFILE

All content following this page was uploaded by Aymen Ammari on 08 June 2016.

The user has requested enhancement of the downloaded file.

If you have a format like ‘date1’ type Date variable
-----STATA 10.x/11.x:
gen datevar = date(date1,"DMY", 2012)
format datevar %td /*For daily data*/
-----STATA 9.x:
gen datevar = date(date1,"dmy", 2012)
format datevar %td /*For daily data*/

If you have a format like ‘date2’ type

-----STATA 10.x/11.x:
gen datevar = date(date2,"MDY", 2012)
format datevar %td /*For daily data*/
------STATA 9.x:
gen datevar = date(date2,"mdy", 2012)
format datevar %td /*For daily data*/

If you have a format like ‘date3’ type

------STATA 10.x/11.x:
tostring date3, gen(date3a)
gen datevar=date(date3a,"YMD")
format datevar %td /*For daily data*/
------STATA 9.x:
tostring date3, gen(date3a)
gen year=substr(date3a,1,4)
gen month=substr(date3a,5,2)
gen day=substr(date3a,7,2)
destring year month day, replace
gen datevar1 = mdy(month,day,year)
format datevar1 %td /*For daily data*/

If you have a format like ‘date4’ type

See http://www.princeton.edu/~otorres/Stata/
2
PU/DSS/OTR
 Date variable (cont.)
If the original date variable is string (i.e. color red):

gen week= weekly(stringvar,"wy")

NOTE: Remember to format the date variable accordingly. After creating it type:
gen month= monthly(stringvar,"my")
gen quarter= quarterly(stringvar,"qy") format datevar %t? /*Change ‘datevar’ with your date variable*/
gen half = halfyearly(stringvar,"hy") Change “?” with the correct format: w (week), m (monthly), q (quarterly), h (half), y (yearly).
gen year= yearly(stringvar,"y")

If the components of the original date are in different numeric variables (i.e. color black):

gen daily = mdy(month,day,year)

gen week = yw(year, week) NOTE: Remember to format the date variable accordingly. After creating it type:
gen month = ym(year,month) format datevar %t? /*Change ‘datevar’ with your date variable*/
gen quarter = yq(year,quarter)
Change “?” with the correct format: w (week), m (monthly), q (quarterly), h (half), y (yearly).
gen half = yh(year,half-year)

To extract days of the week (Monday, Tuesday, etc.) use the function dow()

gen dayofweek= dow(date)

Replace “date” with the date variable in your dataset. This will create the variable ‘dayofweek’ where 0 is ‘Sunday’, 1 is
‘Monday’, etc. (type help dow for more details)

To specify a range of dates (or integers in general) you can use the tin() and twithin() functions. tin() includes the
first and last date, twithin() does not. Use the format of the date variable in your dataset.

/* Make sure to set your data as time series before using tin/twithin */

tsset date

regress y x1 x2 if tin(01jan1995,01jun1995)
regress y x1 x2 if twithin(01jan2000,01jan2001)

3
PU/DSS/OTR
 Date variable (example)
Time series data is data collected over time for a single or a group of variables. For this kind of data the first thing
to do is to check the variable that contains the time or date range and make sure is the one you need: yearly,
monthly, quarterly, daily, etc.
The next step is to verify it is in the correct format. In the example below the time variable is stored in “date”
but it is a string variable not a date variable. In Stata you need to convert this string variable to a date variable.*
A closer inspection of the variable, for the years 2000 the format changes, we need to create a new variable with
a uniform format. Type the following:

use http://dss.princeton.edu/training/tsdata.dta
gen date1=substr(date,1,7)
gen datevar=quarterly(date1,"yq")
format datevar %tq
browse date date1 datevar

For more details type

help date

*Data source: Stock & Watson’s companion materials

 From daily/monthly date variable to quarterly

use "http://dss.princeton.edu/training/date.dta", clear

*Quarterly date from daily date

gen datevar=date(date2,"MDY", 2012) /*Date2 is a string date variable*/

format datevar %td
gen quarterly = qofd(datevar)
format quarterly %tq

*Quarterly date from monthly date

gen month = month(datevar)

gen day=day(datevar)
gen year=year(datevar)
gen monthly = ym(year,month)
format monthly %tm

gen quarterly1 = qofd(dofm(monthly))

format quarterly1 %tq

browse date2 datevar quarterly monthly quarterly1

 From daily to weekly and getting yearly
use "http://dss.princeton.edu/training/date.dta", clear
gen datevar = date(date2, "MDY", 2012)
format datevar %td
gen year= year(datevar)
gen w = week(datevar)
gen weekly = yw(year,w)
format weekly %tw

browse

********************************************************

* From daily to yearly

gen year1 = year(datevar)
* From quarterly to yearly
gen year2 = yofd(dofq(quarterly))
* From weekly to yearly
gen year3 = yofd(dofw(weekly))
 Setting as time series: tsset

Once you have the date variable in a ‘date format’ you need to declare your data as time series in order to
use the time series operators. In Stata type:
tsset datevar

. tsset datevar
time variable: datevar, 1957q1 to 2005q1
delta: 1 quarter

If you have gaps in your time series, for example there may not be data available for weekends. This
complicates the analysis using lags for those missing dates. In this case you may want to create a continuous
time trend as follows:

gen time = _n

Then use it to set the time series:

tsset time

In the case of cross-sectional time series type:

sort panel date

by panel: gen time = _n
xtset panel time
6
PU/DSS/OTR
 Filling gaps in time variables

Use the command tsfill to fill in the gap in the time series. You need to tset, tsset or xtset the data
before using tsfill. In the example below:

tset quarters

tsfill

Type help tsfill for more details.

7
PU/DSS/OTR
 Subsetting tin/twithin

With tsset (time series set) you can use two time series commands: tin (‘times in’, from a to b) and
twithin (‘times within’, between a and b, it excludes a and b). If you have yearly data just include the years.

. list datevar unemp if tin(2000q1,2000q4)

datevar unemp

173. 2000q1 4.033333

174. 2000q2 3.933333
175. 2000q3 4
176. 2000q4 3.9

. list datevar unemp if twithin(2000q1,2000q4)

datevar unemp

174. 2000q2 3.933333

175. 2000q3 4

/* Make sure to set your data as time series before using tin/twithin */

tsset date

regress y x1 x2 if tin(01jan1995,01jun1995)

regress y x1 x2 if twithin(01jan2000,01jan2001)
8
PU/DSS/OTR
 Merge/Append
See

http://dss.princeton.edu/training/Merge101.pdf

PU/DSS/OTR
 Lag operators (lag)
Another set of time series commands are the lags, leads, differences and seasonal operators.
It is common to analyzed the impact of previous values on current ones.
To generate values with past values use the “L” operator

generate unempL1=L1.unemp
generate unempL2=L2.unemp
list datevar unemp unempL1 unempL2 in 1/5
. generate unempL1=L1.unemp
(1 missing value generated)

. generate unempL2=L2.unemp
(2 missing values generated)

. list datevar unemp unempL1 unempL2 in 1/5

datevar unemp unempL1 unempL2

1. 1957q1 3.933333 . .
2. 1957q2 4.1 3.933333 .
3. 1957q3 4.233333 4.1 3.933333
4. 1957q4 4.933333 4.233333 4.1
5. 1958q1 6.3 4.933333 4.233333

In a regression you could type:

regress y x L1.x L2.x
or regress y x L(1/5).x
10
PU/DSS/OTR
 Lag operators (forward)
To generate forward or lead values use the “F” operator
generate unempF1=F1.unemp
generate unempF2=F2.unemp
list datevar unemp unempF1 unempF2 in 1/5

. generate unempF1=F1.unemp
(1 missing value generated)

. generate unempF2=F2.unemp
(2 missing values generated)

. list datevar unemp unempF1 unempF2 in 1/5

datevar unemp unempF1 unempF2

1. 1957q1 3.933333 4.1 4.233333

2. 1957q2 4.1 4.233333 4.933333
3. 1957q3 4.233333 4.933333 6.3
4. 1957q4 4.933333 6.3 7.366667
5. 1958q1 6.3 7.366667 7.333333

In a regression you could type:

regress y x F1.x F2.x
or regress y x F(1/5).x 11
PU/DSS/OTR
 Lag operators (difference)
To generate the difference between current a previous values use the “D” operator
generate unempD1=D1.unemp /* D1 = y t – yt-1 */
generate unempD2=D2.unemp /* D2 = (y t – yt-1) – (yt-1 – yt-2) */
list datevar unemp unempD1 unempD2 in 1/5

. generate unempD1=D1.unemp
(1 missing value generated)

. generate unempD2=D2.unemp
(2 missing values generated)

. list datevar unemp unempD1 unempD2 in 1/5

datevar unemp unempD1 unempD2

1. 1957q1 3.933333 . .
2. 1957q2 4.1 .1666665 .
3. 1957q3 4.233333 .1333332 -.0333333
4. 1957q4 4.933333 .7000003 .5666671
5. 1958q1 6.3 1.366667 .6666665

In a regression you could type:

regress y x D1.x
12
PU/DSS/OTR
 Lag operators (seasonal)
To generate seasonal differences use the “S” operator
generate unempS1=S1.unemp /* S1 = y t – yt-1 */
generate unempS2=S2.unemp /* S2 = (y t – yt-2) */
list datevar unemp unempS1 unempS2 in 1/5

. generate unempS1=S1.unemp
(1 missing value generated)

. generate unempS2=S2.unemp
(2 missing values generated)

. list datevar unemp unempS1 unempS2 in 1/5

datevar unemp unempS1 unempS2

1. 1957q1 3.933333 . .
2. 1957q2 4.1 .1666665 .
3. 1957q3 4.233333 .1333332 .2999997
4. 1957q4 4.933333 .7000003 .8333335
5. 1958q1 6.3 1.366667 2.066667

In a regression you could type:

regress y x S1.x
13
PU/DSS/OTR
 Correlograms: autocorrelation
To explore autocorrelation, which is the correlation between a variable and its previous values,
use the command corrgram. The number of lags depend on theory, AIC/BIC process or
experience. The output includes autocorrelation coefficient and partial correlations coefficients
used to specify an ARIMA model.
corrgram unemp, lags(12)

. corrgram unemp, lags(12)

-1 0 1 -1 0 1
LAG AC PAC Q Prob>Q [Autocorrelation] [Partial Autocor]

1 0.9641 0.9650 182.2 0.0000

2 0.8921 -0.6305 339.02 0.0000
3 0.8045 0.1091 467.21 0.0000
4 0.7184 0.0424 569.99 0.0000
5 0.6473 0.0836 653.86 0.0000
6 0.5892 -0.0989 723.72 0.0000
7 0.5356 -0.0384 781.77 0.0000
8 0.4827 0.0744 829.17 0.0000
9 0.4385 0.1879 868.5 0.0000
10 0.3984 -0.1832 901.14 0.0000
11 0.3594 -0.1396 927.85 0.0000
12 0.3219 0.0745 949.4 0.0000

AC shows that the PAC shows that the Box-Pierce’ Q statistic tests Graphic view of AC Graphic view of PAC
correlation between the correlation between the the null hypothesis that all which shows a slow which does not show
current value of unemp and current value of unemp and correlation up to lag k are decay in the trend, spikes after the second
its value three quarters ago its value three quarters ago equal to 0. This series show suggesting non- lag which suggests that
is 0.8045. AC can be use to is 0.1091 without the effect significant autocorrelation stationarity. all other lags are mirrors
define the q in MA(q) only of the two previous lags. as shown in the Prob>Q See also the ac of the second lag. See the
in stationary series PAC can be used to define value which at any k are less command. pac command.
the p in AR(p) only in than 0.05, therefore
stationary series rejecting the null that all lags
are not autocorrelated. 14
PU/DSS/OTR
 Correlograms: cross correlation
The explore the relationship between two time series use the command xcorr. The graph below shows the correlation
between GDP quarterly growth rate and unemployment. When using xcorr list the independent variable first and the
dependent variable second. type

xcorr gdp unemp, lags(10) xlabel(-10(1)10,grid) . xcorr gdp unemp, lags(10) table

-1 0 1
LAG CORR [Cross-correlation]
Cross-correlogram
1.00

1.00
-10 -0.1080
-9 -0.1052
Cross-correlations of gdp and unemp

-8 -0.1075
-7 -0.1144
0.50

0.50
-6 -0.1283
-5 -0.1412
-4 -0.1501
-3 -0.1578
0.00

0.00
-2 -0.1425
-1 -0.1437
0 -0.1853
-0.50

-0.50
1 -0.1828
2 -0.1685
3 -0.1177
4 -0.0716
-1.00

-1.00 5 -0.0325
-10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 6 -0.0111
Lag 7 -0.0038
8 0.0168
9 0.0393
10 0.0419

At lag 0 there is a negative immediate correlation between GDP growth rate and unemployment. This means that a drop
in GDP causes an immediate increase in unemployment.

15
PU/DSS/OTR
 Correlograms: cross correlation
xcorr interest unemp, lags(10) xlabel(-10(1)10,grid)

Cross-correlogram
1.00

1.00
. xcorr interest unemp, lags(10) table
Cross-correlations of interest and unemp

-1 0 1
0.50

0.50
LAG CORR [Cross-correlation]

-10 0.3297
-9 0.3150
0.00

0.00
-8 0.2997
-7 0.2846
-6 0.2685
-0.50

-0.50
-5 0.2585
-4 0.2496
-3 0.2349
-1.00

-1.00
-2 0.2323
-10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10
Lag -1 0.2373
0 0.2575
1 0.3095
2 0.3845
Interest rates have a positive effect on future level of 3 0.4576
unemployment, reaching the highest point at lag 8 (four 4 0.5273
quarters or two years). In this case, interest rates are positive 5 0.5850
correlated with unemployment rates eight quarters later. 6 0.6278
7 0.6548
8 0.6663
9 0.6522
10 0.6237

16
PU/DSS/OTR
 Lag selection
Too many lags could increase the error in the forecasts, too few could leave out relevant information*.
Experience, knowledge and theory are usually the best way to determine the number of lags needed. There
are, however, information criterion procedures to help come up with a proper number. Three commonly used
are: Schwarz's Bayesian information criterion (SBIC), the Akaike's information criterion (AIC), and the
Hannan and Quinn information criterion (HQIC). All these are reported by the command ‘varsoc’ in Stata.
. varsoc gdp cpi, maxlag(10)

Selection-order criteria
Sample: 1959q4 - 2005q1 Number of obs = 182

lag LL LR df p FPE AIC HQIC SBIC

0 -1294.75 5293.32 14.25 14.2642 14.2852

1 -467.289 1654.9 4 0.000 .622031 5.20098 5.2438 5.30661
2 -401.381 131.82 4 0.000 .315041 4.52067 4.59204 4.69672*
3 -396.232 10.299 4 0.036 .311102 4.50804 4.60796 4.75451
4 -385.514 21.435* 4 0.000 .288988* 4.43422* 4.56268* 4.7511
5 -383.92 3.1886 4 0.527 .296769 4.46066 4.61766 4.84796
6 -381.135 5.5701 4 0.234 .300816 4.47401 4.65956 4.93173
7 -379.062 4.1456 4 0.387 .307335 4.49519 4.70929 5.02332
8 -375.483 7.1585 4 0.128 .308865 4.49981 4.74246 5.09836
9 -370.817 9.3311 4 0.053 .306748 4.4925 4.76369 5.16147
10 -370.585 .46392 4 0.977 .319888 4.53391 4.83364 5.27329

Endogenous: gdp cpi

Exogenous: _cons

When all three agree, the selection is clear, but what happens when getting conflicting results? A paper from
the CEPR suggests, in the context of VAR models, that AIC tends to be more accurate with monthly data,
HQIC works better for quarterly data on samples over 120 and SBIC works fine with any sample size for
quarterly data (on VEC models)**. In our example above we have quarterly data with 182 observations,
HQIC suggest a lag of 4 (which is also suggested by AIC).
* See Stock & Watson for more details and on how to estimate BIC and SIC
** Ivanov, V. and Kilian, L. 2001. 'A Practitioner's Guide to Lag-Order Selection for Vector Autoregressions'. CEPR Discussion Paper no. 2685. London, Centre for Economic Policy
Research. http://www.cepr.org/pubs/dps/DP2685.asp. 17
PU/DSS/OTR
 Unit roots
Having a unit root in a series mean that there is more than one trend in the series.

. regress unemp gdp if tin(1965q1,1981q4)

Source SS df MS Number of obs = 68

F( 1, 66) = 19.14
Model 36.1635247 1 36.1635247 Prob > F = 0.0000
Residual 124.728158 66 1.88982058 R-squared = 0.2248
Adj R-squared = 0.2130
Total 160.891683 67 2.4013684 Root MSE = 1.3747

unemp Coef. Std. Err. t P>|t| [95% Conf. Interval]

gdp -.4435909 .1014046 -4.37 0.000 -.6460517 -.2411302

_cons 7.087789 .3672397 19.30 0.000 6.354572 7.821007

. regress unemp gdp if tin(1982q1,2000q4)

Source SS df MS Number of obs = 76

F( 1, 74) = 3.62
Model 8.83437339 1 8.83437339 Prob > F = 0.0608
Residual 180.395848 74 2.43778172 R-squared = 0.0467
Adj R-squared = 0.0338
Total 189.230221 75 2.52306961 Root MSE = 1.5613

unemp Coef. Std. Err. t P>|t| [95% Conf. Interval]

gdp .3306551 .173694 1.90 0.061 -.0154377 .6767479

_cons 5.997169 .2363599 25.37 0.000 5.526211 6.468126

18
PU/DSS/OTR
 Unit roots
Unemployment rate.
line unemp datevar

12 10
Unemployment Rate
6 84

1957q1 1969q1 1981q1 1993q1 2005q1

datevar

19
PU/DSS/OTR
 Unit root test

The Dickey-Fuller test is one of the most commonly use tests for stationarity. The null
hypothesis is that the series has a unit root. The test statistic shows that the unemployment
series have a unit root, it lies within the acceptance region.

One way to deal with stochastic trends (unit root) is by taking the first difference of the variable
(second test below).
. dfuller unemp, lag(5)

Augmented Dickey-Fuller test for unit root Number of obs = 187

Interpolated Dickey-Fuller
Unit root Test 1% Critical 5% Critical 10% Critical
Statistic Value Value Value

Z(t) -2.597 -3.481 -2.884 -2.574

MacKinnon approximate p-value for Z(t) = 0.0936

. dfuller unempD1, lag(5)

Augmented Dickey-Fuller test for unit root Number of obs = 186

Interpolated Dickey-Fuller
Test 1% Critical 5% Critical 10% Critical
No unit root Statistic Value Value Value

Z(t) -5.303 -3.481 -2.884 -2.574

MacKinnon approximate p-value for Z(t) = 0.0000 20

PU/DSS/OTR
 Testing for cointegration
Cointegration refers to the fact that two or more series share an stochastic trend (Stock &
Watson). Engle and Granger (1987) suggested a two step process to test for cointegration (an
OLS regression and a unit root test), the EG-ADF test.

regress unemp gdp Run an OLS regression

predict e, resid Get the residuals

dfuller e, lags(10) Run a unit root test on the residuals.

Augmented Dickey-Fuller test for unit root Number of obs = 181

Interpolated Dickey-Fuller
Unit root* Test 1% Critical 5% Critical 10% Critical
Statistic Value Value Value

Z(t) -2.535 -3.483 -2.885 -2.575

MacKinnon approximate p-value for Z(t) = 0.1071

Both variables are not cointegrated

See Stock & Watson for a table of critical values for the unit root test and the theory behind.

21
*Critical value for one independent variable in the OLS regression, at 5% is -3.41 (Stock & Watson) PU/DSS/OTR
 Granger causality: using OLS
If you regress ‘y’ on lagged values of ‘y’ and ‘x’ and the coefficients of the lag of ‘x’ are
statistically significantly different from 0, then you can argue that ‘x’ Granger-cause ‘y’, this is,
‘x’ can be used to predict ‘y’ (see Stock & Watson -2007-, Green -2008).

1 . regress unemp L(1/4).unemp L(1/4).gdp

Source SS df MS Number of obs = 188

F( 8, 179) = 668.37
Model 373.501653 8 46.6877066 Prob > F = 0.0000
Residual 12.5037411 179 .069853302 R-squared = 0.9676
Adj R-squared = 0.9662
Total 386.005394 187 2.06419997 Root MSE = .2643

unemp Coef. Std. Err. t P>|t| [95% Conf. Interval]

unemp
L1. 1.625708 .0763035 21.31 0.000 1.475138 1.776279
L2. -.7695503 .1445769 -5.32 0.000 -1.054845 -.484256
L3. .0868131 .1417562 0.61 0.541 -.1929152 .3665415
L4. .0217041 .0726137 0.30 0.765 -.1215849 .1649931

gdp
L1. .0060996 .0136043 0.45 0.654 -.0207458 .0329451
L2. -.0189398 .0128618 -1.47 0.143 -.0443201 .0064405
L3. .0247494 .0130617 1.89 0.060 -.0010253 .0505241
L4. .003637 .0129079 0.28 0.778 -.0218343 .0291083

_cons .1702419 .096857 1.76 0.081 -.0208865 .3613704

2 . test L1.gdp L2.gdp L3.gdp L4.gdp

( 1) L.gdp = 0
(
(
2)
3)
L2.gdp = 0
L3.gdp = 0
You cannot reject the null hypothesis that all
( 4) L4.gdp = 0 coefficients of lag of ‘x’ are equal to 0.
F( 4, 179) = 1.67 Therefore ‘gdp’ does not Granger-cause
Prob > F = 0.1601
‘unemp’.
. 22
PU/DSS/OTR
 Granger causality: using VAR
The following procedure uses VAR models to estimate Granger causality using the command
‘vargranger’

1 . quietly var unemp gdp, lags(1/4)

2 . vargranger

Granger causality Wald tests

Equation Excluded chi2 df Prob > chi2

unemp gdp 6.9953 4 0.136

unemp ALL 6.9953 4 0.136

gdp unemp 6.8658 4 0.143

gdp ALL 6.8658 4 0.143

The null hypothesis is ‘var1 does not Granger-cause var2’. In both cases, we cannot reject
the null that each variable does not Granger-cause the other

23
PU/DSS/OTR
 Chow test (testing for known breaks)
The Chow test allows to test whether a particular date causes a break in the regression coefficients. It is named after Gregory
Chow (1960)*.

Step 1. Create a dummy variable where 1 if date > break date and 0 <= break date. Below we’ll test whether the first quarter of
1982 causes a break in the regression coefficients.

tset datevar
Change “tq” with the correct date format: tw (week), tm (monthly), tq (quarterly), th (half), ty (yearly) and the
gen break = (datevar>tq(1981q4)) corresponding date format in the parenthesis

Step 2. Create interaction terms between the lags of the independent variables and the lag of the dependent variables. We will
assume lag 1 for this example (the number of lags depends on your theory/data)

generate break_unemp = break*l1.unemp

generate break_gdp = break*l1.gdp

Step 3. Run a regression between the outcome variables (in this case ‘unemp’) and the independent along with the interactions
and the dummy for the break.

reg unemp l1.unemp l1.gdp break break_unemp break_gdp

Step 4. Run an F-test on the coefficients for the interactions and the dummy for the break

test break break_unemp break_gdp

The null hypothesis is no break. If the p-value
is < 0.05 reject the null in favor of the
. test break break_unemp break_gdp alternative that there is a break. In this
example, we fail to reject the null and
( 1) break = 0 conclude that the first quarter of 1982 does
( 2) break_unemp = 0 not cause a break in the regression
( 3) break_gdp = 0 coefficients.
F( 3, 185) = 1.14
Prob > F = 0.3351 * See Stock & Watson for more details 24
PU/DSS/OTR
 Testing for unknown breaks
The Quandt likelihood ratio (QLR test –Quandt,1960) or sup-Wald statistic is a modified version of the Chow test used to identify
break dates. The following is a modified procedure taken from Stock & Watson’s companion materials to their book Introduction
to Econometrics, I strongly advise to read the corresponding chapter to better understand the procedure and to check the critical
values for the QLR statistic. Below we will check for breaks in a GDP per-capita series (quarterly).

/* Replace the words in bold with your own variables, do not change anything else*/
/* The log file ‘qlrtest.log’ will have the list for QLR statistics (use Word to read it)*/
/* See next page for a graph*/
/* STEP 1. Copy-and-paste-run the code below to a do-file, double-check the quotes (re-type them if necessary)*/

log using qlrtest.log

tset datevar
sum datevar
local time=r(max)-r(min)+1
local i = round(`time'*.15)
local f = round(`time'*.85)
local var = "gdp"
gen diff`var' = d.`var'
gen chow`var' = .
gen qlr`var' = .
set more off
while `i'<=(`f') {
gen di = (_n > `i')
cap gen d_`var'1 = di*l1.`var'
cap gen d_`var'2 = di*l2.`var'
cap gen d_`var'3 = di*l3.`var'
cap gen d_`var'4 = di*l4.`var'
qui reg diff`var' L(1/4).diff`var' di, r
qui test di
sca chow = r(F)
cap replace chow`var' = r(F) in `i'
qui reg diff`var' L(1/4).diff`var' di d_`var'1 d_`var'2 d_`var'3 d_`var'4, r
qui test di d_`var'1 d_`var'2 d_`var'3 d_`var'4
sca qlr = r(F)
cap replace qlr`var' = r(F) in `i'
dis "`i' " %tq datevar[`i'] " " %8.3f chow " " %8.3f qlr
drop di d_`var'1 d_`var'2 d_`var'3 d_`var'4
local i = `i' + 1
}
25
NOTE: change %tq if your date variable is other than quarterly. Use %td for daily, %tm for monthly, %ty for yearly PU/DSS/OTR
 Testing for unknown breaks: graph
/* Replace the words in bold with your own variables, do not change anything else*/
/* The code will produce the graph shown in the next page*/
/* The critical value 3.66 is for q=5 (constant and four lags) and 5% significance*/
/* STEP 2. Copy-and-paste-run the code below to a do-file, double-check the quotes (re-type them if necessary)*/

sum qlr`var'
local maxvalue=r(max)
gen maxdate=datevar if qlr`var'==`maxvalue'
local maxvalue1=round(`maxvalue',0.01)
local critical=3.66 /*Replace with the appropriate critical value (see Stock & Watson)*/
sum datevar
local mindate=r(min)
sum maxdate
local maxdate=r(max)
gen break=datevar if qlr`var'>=`critical' & qlr`var'!=.
dis "Below are the break dates..."
list datevar qlr`var' if break!=.
levelsof break, local(break1)
twoway tsline qlr`var', title(Testing for breaks in GDP per-capita (1957-2005)) ///
xlabel(`break1', angle(90) labsize(0.9) alternate) ///
yline(`critical') ytitle(QLR statistic) xtitle(Time) ///
ttext(`critical' `mindate' "Critical value 5% (`critical')", placement(ne)) ///
ttext(`maxvalue' `maxdate' "Max QLR = `maxvalue1'", placement(e))

datevar qlrgdp

129. 1989q1 3.823702

131. 1989q3 6.285852
132. 1989q4 6.902882
133. 1990q1 5.416068
134. 1990q2 7.769114

135. 1990q3 8.354294

136. 1990q4 5.399252
137. 1991q1 8.524492
138. 1991q2 6.007093
139. 1991q3 4.44151

140. 1991q4 4.946689

141. 1992q1 3.699911
160. 1996q4 3.899656
161. 1997q1 3.906271 26
PU/DSS/OTR
 Testing for unknown breaks: graph (cont.)

Testing for breaks in GDP per-capita (1957-2005)

Max QLR = 8.52
8
6
QLR statistic
4

Critical value 5% (3.66)

2
0

1989q1

1989q4
1990q2
1990q4
1991q2
1991q4

1996q4
1989q3
1990q1
1990q3
1991q1
1991q3
1992q1

1997q1
Time

27
PU/DSS/OTR
 Time Series: white noise
White noise refers to the fact that a variable does not have autocorrelation. In Stata use the
wntestq (white noise Q test) to check for autocorrelation. The null is that there is no serial
correlation (type help wntestq for more details):

. wntestq unemp

Portmanteau test for white noise

Portmanteau (Q) statistic = 1044.6341

Prob > chi2(40) = 0.0000

. wntestq unemp, lags(10)

Serial correlation
Portmanteau test for white noise

Portmanteau (Q) statistic = 901.1399

Prob > chi2(10) = 0.0000

If your variable is not white noise then see the page on correlograms to see the order of the
autocorrelation.

28
PU/DSS/OTR
 Time Series: Testing for serial correlation
Breush-Godfrey and Durbin-Watson are used to test for serial correlation. The null in both tests
is that there is no serial correlation (type help estat dwatson, help estat dubinalt
and help estat bgodfrey for more details).
. regress unempd gdp

Source SS df MS Number of obs = 192

F( 1, 190) = 1.62
Model .205043471 1 .205043471 Prob > F = 0.2048
Residual 24.0656991 190 .126661574 R-squared = 0.0084
Adj R-squared = 0.0032
Total 24.2707425 191 .12707195 Root MSE = .3559

unempd1 Coef. Std. Err. t P>|t| [95% Conf. Interval]

gdp -.015947 .0125337 -1.27 0.205 -.0406701 .0087761

_cons .0401123 .036596 1.10 0.274 -.0320743 .1122989

. estat dwatson

Durbin-Watson d-statistic( 2, 192) = .7562744

. estat durbinalt

Durbin's alternative test for autocorrelation

lags(p) chi2 df Prob > chi2

1 118.790 1 0.0000

H0: no serial correlation

. estat bgodfrey
Serial correlation
Breusch-Godfrey LM test for autocorrelation

lags(p) chi2 df Prob > chi2

1 74.102 1 0.0000

H0: no serial correlation 29

PU/DSS/OTR
 Time Series: Correcting for serial correlation
Run a Cochrane-Orcutt regression using the prais command (type help prais for more
details)

. prais unemp gdp, corc

Iteration 0: rho = 0.0000

Iteration 1: rho = 0.9556
Iteration 2: rho = 0.9660
Iteration 3: rho = 0.9661
Iteration 4: rho = 0.9661

Cochrane-Orcutt AR(1) regression -- iterated estimates

Source SS df MS Number of obs = 191

F( 1, 189) = 2.48
Model .308369041 1 .308369041 Prob > F = 0.1167
Residual 23.4694088 189 .124176766 R-squared = 0.0130
Adj R-squared = 0.0077
Total 23.7777778 190 .125146199 Root MSE = .35239

unemp Coef. Std. Err. t P>|t| [95% Conf. Interval]

gdp -.020264 .0128591 -1.58 0.117 -.0456298 .0051018

_cons 6.105931 .7526023 8.11 0.000 4.621351 7.590511

rho .966115

Durbin-Watson statistic (original) 0.087210

Durbin-Watson statistic (transformed) 0.758116
30
PU/DSS/OTR
Useful links / Recommended books

• DSS Online Training Section http://dss.princeton.edu/training/

• UCLA Resources to learn and use STATA http://www.ats.ucla.edu/stat/stata/
• DSS help-sheets for STATA http://dss/online_help/stats_packages/stata/stata.htm
• Introduction to Stata (PDF), Christopher F. Baum, Boston College, USA. “A 67-page description of Stata, its key
features and benefits, and other useful information.” http://fmwww.bc.edu/GStat/docs/StataIntro.pdf
• STATA FAQ website http://stata.com/support/faqs/
• Princeton DSS Libguides http://libguides.princeton.edu/dss
Books
• Introduction to econometrics / James H. Stock, Mark W. Watson. 2nd ed., Boston: Pearson Addison
Wesley, 2007.
• Data analysis using regression and multilevel/hierarchical models / Andrew Gelman, Jennifer Hill.
Cambridge ; New York : Cambridge University Press, 2007.
• Econometric analysis / William H. Greene. 6th ed., Upper Saddle River, N.J. : Prentice Hall, 2008.
• Designing Social Inquiry: Scientific Inference in Qualitative Research / Gary King, Robert O.
Keohane, Sidney Verba, Princeton University Press, 1994.
• Unifying Political Methodology: The Likelihood Theory of Statistical Inference / Gary King, Cambridge
University Press, 1989
• Statistical Analysis: an interdisciplinary introduction to univariate & multivariate methods / Sam
Kachigan, New York : Radius Press, c1986

• Statistics with Stata (updated for version 9) / Lawrence Hamilton, Thomson Books/Cole, 2006
31
View publication stats PU/DSS/OTR