similar to: KPSS test - Lag selection

Hi R users! I've come across using kpss tests for time series analysis and i have a question that troubles me since i don't have much experience with time series and the mathematical part underlining it. x<-c(253, 252, 275, 275, 272, 254, 272, 252, 249, 300, 244, 258, 255, 285, 301, 278, 279, 304, 275, 276, 313, 292, 302, 322, 281, 298, 305, 295, 286, 327, 286, 270, 289, 293, 287,

Situation: I had tired a 1000-data generated by random error(i.i.d.), then I sub it into different unit root tests. I got different results among the tests. The following are the test statistics I got: adf.test @ tseries ~ -10.2214 (lag = 9) ur.df @ urca ~ -21.8978 ur.sp @ urca ~ -27.68 pp.test @ tseries ~ -972.3343 (truncation lag =7) ur.pp @ urca ~ -973.2409 ur.kpss @ urca ~ 0.1867 kpss.test @

I'm trying to understand how to use kpss.test() properly. If I have a level stationary series like rnorm() in the help page, shouldn't I get a small p-value with the null hypothesis set to "Trend"? The (condensed) output from kpss.test() for the two possible null hypotheses is given below. I don't see any significant difference between these results. > x <-

Dear R-list member, an update of package 'urca' has been uploaded to CRAN (Mirror: Austria). In the updated release unit root and cointegration tests encountered in applied econometric analysis are implemented. The package is written in 'pure' R and utilises S4 classes. In particular, the Johansen procedure with likelihood ratio tests for the inclusion of a linear trend,

is that 'x' is level or trend stationary. I did this > s<-rnorm(1000) > kpss.test(s) KPSS Test for Level Stationarity data: s KPSS Level = 0.0429, Truncation lag parameter = 7, p-value = 0.1 Warning message: p-value greater than printed p-value in: kpss.test(s) My question is whether p=0.1 is a good number to reject N0? On the other hand, I have a

I'm tring the functions to check the cointegration of a matrix. I'm using **Phillips & Ouliaris Cointegration Test** The function in *tseries* package is **po.test** and **ca.po** in *urca* The results with **URCA** are: > ca.po(prices, demean='none') ######################################## # Phillips and Ouliaris Unit Root Test #

Hi, Am I interpreting the results properly? Are my conclusions correct? > KPSS.test(df) ---- ---- KPSS test ---- ---- Null hypotheses: Level stationarity and stationarity around a linear trend. Alternative hypothesis: Unit root. ---- Statistic for the null hypothesis of level stationarity: 1.089 Critical values: 0.10 0.05 0.025 0.01 0.347 0.463

Hello, I have a question on unit root test with urca toolbox. First, to run a unit root test with lags selected by BIC, I type: > CPILD4UR<-ur.df(x1$CPILD4[5:nr1], type ="drift", lags=12, selectlags ="BIC") > summary(CPILD4UR) The results indicate that the optimal lags selected by BIC is 4. Then I run the same unit root test with drift and 4 lags:

Hi I am using the package urca and I am interested about the KPSS test. That works fine except the method "summary" did not work in the script, only when it is typed direct in the console the results are shown( not a source file). Is there any problem with these method ?

Dear R-Devel subscriber, as a follow up to my yesterday's email: I tested an analogous example with the S4-package "flexclust" by executing the following code: library(flexclust) example(cclust) cl After saving the work space and starting a new R process with the restored work space, the same behaviour (i.e., the methods pertinent to "flexclust" are not used, even after

Dear R-Devel subscriber, I was hinted to the following problem with package 'urca': If one starts R and executes for instance: > library(urca) > example(ur.df) ## output as expected, but omitted here > class(lc.df) [1] "ur.df" attr(,"package") [1] "urca" > class(summary(lc.df)) [1] "sumurca" attr(,"package") [1]

Sorry, to bother the list one more time: but the following worked at least for 'urca': in NAMESPACE I now included explicitly: import(methods) a fix of the 'urca'-package will be uploaded to CRAN on the weekend. Fritz, will this work for ypur package 'flexclust' too? I have in my DESCRIPTION imports: methods and in flexclust it is in depends: methods. However, both

Hi, While doing the ADF test in R using the following command I am getting the error and the result.."> x.ct=ur.df(rev$REVENUE,start=1,end=length(rev$REVENUE),frequency=1) Error in ur.df(rev$REVENUE, start = 1, end = length(rev$REVENUE), frequency = 1) : unused argument(s) (start = 1, end = 4, frequency = 1) >

Hi All, First of all, could you tell me what the "KPSS Level" in the output of the test means? I have a series, x, of periodic data and tried kpss.test() on it to verify its stationarity. The tests gave me the p-value above 0.1. Since the null hypothesis N0 is that the series _is_ stationary, this means that I cannot reject N0. But the series does look periodic! So does all this

Hi, I have been fitting a series of data representing a week of Internet traffic (which is daily seasonal and have a general trend toward lower rate at the weekends). Before I do the ARMA fit (which takes care of seasonality with a lag equal to one day), do I have to make sure the data is stationary? From the results and visually, it seems that this was taken care of. But the residual

Dear all, How to do ADF test ¡¢KPSS¡¢ PP¡¢GLS test in R£¿ Thanks a lot ! [[alternative HTML version deleted]]

Dear All, I am trying to use kpss.test function so as to perform a stationarity test on a data sample. Is it possible to know the associated confidence level for this test? I have not seen any arguments related to it. I had a look at some other tests included in R (adf.test, pp.test, ks.test ...) and I could not find this information for them. Thanks in advanced. Kind regards, Belén

Ukech U. Kidi wrote: > dax<- diff(log(DAX_CAC$DAX[1:1865])) > m1<- garch(dax) > Error: couldn't find function "garch" > m1<- garch(dax[1:1865]) > Error: couldn't find function "garch" > m1<- garch(dax[1:1865]) I am sorry, but I forgot to change the addres to r-help in the reply. Well, I am not sure, wheere do you want to get

Hi all, Can anyone clear my doubts about what conclusions to take with the following what puts of some time series tests: > adf.test(melbmax) Augmented Dickey-Fuller Test data: melbmax Dickey-Fuller = -5.4075, Lag order = 15, p-value = 0.01 alternative hypothesis: stationary Warning message: p-value smaller than printed p-value in: adf.test(melbmax)

Dear Sir, First of all Happy Holidays!,... I am writing to you because I am a bit confused about ARCH estimation. Is there a way to find what garch() exactly does, without the need of reading the source code (because I cannot understand it)? In Eviews (the results at the end) I am getting different results than in R (for those that have the program I do: Quick -> Estimage Equation ->