The implementation lives in Pptest, whose auxiliary regression
always includes the linear trend – that is the tseries::pp.test
specification and the one the tabulated critical values belong to.
trend = FALSE is refused rather than silently answered with
trend-case critical values.
References
Phillips and Perron (1988), Biometrika 75:335-346; coded form from tseries::pp.test. See Pptest for the full note.
Examples
D <- data.frame(x = c(1, 2, 3, 4), y = c(2, 4, 5, 9))
Pptrend(D)
#> $statistic
#> [1] -0.8105035
#>
#> $p_value
#> [1] 0.9484554
#>
#> $rho
#> [1] 0.1531532
#>
#> $lags
#> [1] 2
#>
#> $s2
#> [1] 1.777533
#>
#> $lambda2
#> [1] 1.645231
#>
#> $n
#> [1] 7
#>
#> $method
#> [1] "Phillips-Perron unit root test, Z(t_alpha)"
#>
#> attr(,"class")
#> [1] "morie_rich_result" "list"
