LRurpar.test {partsm}

Likelihood Ratio Test for a Single Unit Root in a PAR(p) Model

Description

This function performs the likelihood ratio test for a single unit root in a PAR(p) model up to order 2.

Usage

    LRurpar.test (wts, detcomp, p)
  

Arguments

wts

a univariate time series object.

detcomp

a vector indicating the deterministic components included in the auxiliary regression. See the corresponding item in fit.ar.par.

p

the order of the PAR model. In this version first and second order models are considered.

Details

In a quarterly time series, the PAR(1) model, y_t = \alpha_{s,1} y_{t-1} + \epsilon_t with \epsilon_t ID(0,1), contains a unit root if g(\alpha) = \Pi_{s=1}^4 \alpha_{s,1} = 1. To test this hypothesis, a likelihood ratio test, LR, is built as the logarithm of the ratio between the residual sum of squares in the unrestricted and the restricted model, weighted by the number of observations.

The unrestricted PAR model is estimated by OLS, whereas the model in which the null hypothesis is imposed, i.e. \Pi_{s=1}^4 \alpha_{s,1} = 1, is estimated by nonlinear least squares.

The critical values are reported in Osterwald-Lenum (1992), table 1.1 (for the case where p-r=1).

In this version, PAR models up to order 2 with seasonal intercepts are considered, since the function fit.piar does not allow for higher orders.

Value

An object of class LRur.partsm-class containing the LR statistics and a one-side test constructed as sign(g(\hat{\alpha}) - 1) * LR^{1/2}, where g(\hat{\alpha}) is the product of the periodic differencing filter parameters estimated under the alternative.

Author(s)

Javier Lopez-de-Lacalle [e-mail].

References

H.P. Boswijk and P.H. Franses (1996), Unit roots in periodic autoregressions. Journal of Time series Analysis, 17, pp. 221-245.

P.H. Franses: Periodicity and Stochastic Trends in Economic Time Series (Oxford University Press, 1996).

Osterwald-Lenum, M. (1992), A Note with Quantiles of the Asymptotic Distribution of the Maximum Likelihood Cointegration Rank Test Statistics: Four Cases. Oxford Bulletin of Economics and Statistics, 54, pp.461-472.

See Also

fit.ar.par, fit.piar, and LRur.partsm-class.

Examples

    ## Test for a single unit root in a PAR(2) model with seasonal intercepts for the
    ## logarithms of the Real GNP in Germany.
    data("gergnp")
    lgergnp <- log(gergnp, base=exp(1))
    detcomp <- list(regular=c(0,0,0), seasonal=c(1,0), regvar=0)
    out <- LRurpar.test(wts=lgergnp, detcomp=detcomp, p=2)