pdlm

Purpose

Computes panel data LM tests with 0, 1, or 2 breaks in the level or level and trend.

Format

{ Nlm, Ntb, Np, PD_lm, cvPDLM } = PDLM(y, model[, nbreak, pmax, ic, trimm])
Parameters
  • y (Nx1 matrix) – Time series data to be tested.

  • model (Scalar) –

    Model to be implemented.

    1

    Level break model.

    2

    Level and trend break model.

  • nbreak (Scalar) –

    Optional, number of breaks.

    0

    No break test.

    1

    One break test.

    2

    Two break tests.

    Default = 0.

  • pmax (Scalar) – Optional, maximum number of lags. Default = 8.

  • ic (Scalar) –

    Optional, the information criterion used for choosing lags. Default = 3.

    1

    Akaike.

    2

    Schwarz.

    3

    t-stat significance.

  • trimm (Scalar) – Optional, trimming rate. Default = 0.10.

Returns
  • Nlm (Vector) – Minimum test statistic for each cross-section.

  • Ntb (Vector) – Location of break(s) for each cross-section.

  • Np (Vector) – Number of lags selected by choosen information criterion for each cross-section.

  • PD_lm (Scalar) – Panel LM statistic with N(0,1).

  • pval (Scalar) – P-value of the panel LM statistic.

Examples

new;
cls;
library tspdlib;

// Load date file
y = loadd(getGAUSSHome() $+ "pkgs/tspdlib/examples/pd_full.csv",
                             "Y1 + Y2 + Y3 + Y4 + Y5 + Y6 + Y7 + Y8 +
                              Y9 + Y10 + Y11 + Y12 + Y13 + Y14 + Y15 +
                              Y16 + Y17 + Y18 + Y19 + Y20 + Y21 + Y22 +
                              Y23 + Y24 + date($Date, '%b-%y')");

// Run panel LM test with level shifts
// No breaks
model = 1;

// No breaks
{ Nlm, Ntb, Np, PD_lm, pval } = PDLM(y, model);

Source

pd_iltlevel.src

See also

Functions adf(), lmkpss()