kpss#

Purpose#

Test for stationarity using a Lagrange Multiplier score statistic.

Format#

{ tstat, crit } = kpss(y, max_lags[, trend, qsk, auto, print_out])#
Parameters:
  • y (Nx1 vector) – data.

  • max_lags (scalar) – Optional input, if max_lags <= 0, maximum lag set using Schwert criterion; if -1, Schwert criterion = 12; if 0, Schwert criterion = 4; else if max_lags > 0, maximum lag = max_lags. Default = 0.

  • trend (scalar) – Optional input, 0 no trend, 1 trend. Default = 0.

  • qsk (scalar) – Optional input, if nonzero, quadratic spectral kernel is used. Default = 0.

  • auto (scalar) – Optional input, if nonzero, automatic maxlags computed. Default = 1.

  • print_out (scalar) – Optional input, if nonzero, intermediate quantities printed to the screen. Default = 1.

Returns:
  • tstat (matrix) – test statistic for each lag.

  • crit (matrix) – Elliot, Rothenberg and Stock (1996) critical values for the GLS detrended unit root test at the 1%, 2.5%, 5%, and 10% significance level.

Example#

new;
cls;
library tsmt;

/********************************************
** LOAD DATA
*********************************************/
npdb= loadd( getGAUSSHome() $+ "pkgs/tsmt/examples/NPDB_orig.csv")';
lrgnp = log(npdb[., 4]);
yt =packr(lrgnp);

// Test using basic KPSS testing: Trend stationary
// Step One: Set-up testing parameters
// Maximum lags to include
max_lags = 5;

// Include trend
trend = 1;

// Use quadratic spectral kernel
qsk = 1;

// Automatic maxlag computation
auto = 1;

// Print results to screen
print_out = 1;

// Step Two: Running KPSS test
{ mat, crit } = kpss(yt, max_lags, trend, qsk, auto, print_out);

print "The tstats for all possible lags:";
mat;

print "Critical values:";
crit;

// Test using basic KPSS testing: Level stationary
// Step One: Set-up testing parameters
// Use defaults for:
// trend = no-trend
// qsk = Bno quadratic spectral kernel used
// printOut = printing output to screen.

// Set maxlags, implies no automatic lag calculation
max_lags = 8;

/ /Running the KPSS test
{ mat, crit } = kpss(yt, max_lags);

print "The tstats for all possible lags:";
mat;

print "Critical values:";
crit;

Library#

tsmt

Source#

kpss.src