similar to: time series tests

Dear R users, I have a function (below) which encompasses several tests. However, when I run it, only the output of the last test is displayed. How can I ensure that the function root(var) will run and display the output from all tests, and not just the last one? Thank you, b. root <- function(var) { #---Phillips-Perron PP.test(var, lshort = TRUE) PP.test(var, lshort = FALSE)

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

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 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 <-

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,

Hi, How can I access the Values stored in kpssstat-class given by KPSS.test function and store it in a variable. For example: >x <- rnorm(1000) >test <- KPSS.test(ts(x)) >test ---- ---- KPSS test ---- ---- Null hypotheses: Level stationarity and stationarity around a linear trend. Alternative hypothesis: Unit root. ---- Statistic for the null

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 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 ?

Does anyone know of a test for stationarity of a time series, or like all ordination techniques it is a qualitative assessment of a quantitative result. Books, papers, etc. suggestions welcome. thanks Stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are

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

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 @

Dear Users, please help with the following DF test: ===== library(tseries) library(timeSeries) Y=c(3519,3803,4332,4251,4661,4811,4448,4451,4343,4067,4001,3934,3652,3768 ,4082,4101,4628,4898,4476,4728,4458,4004,4095,4056,3641,3966,4417,4367 ,4821,5190,4638,4904,4528,4383,4339,4327,3856,4072,4563,4561,4984,5316 ,4843,5383,4889,4681,4466,4463,4217,4322,4779,4988,5383,5591,5322,5404

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

Short question: Is it possible to use statistical tests, like the Augmented Dickey-Fuller test, in functions with for-loops? If not, are there any alternative ways to scale measures? Detailed explanation: I am working with time-series, and I want to flag curves that are not stationary and which display pulses, trends, or level shifts. >df DATE ID VALUE2012-03-06 1

Hi all, Hope you people do not feel irritated for repeatedly sending mail on Time series. Here I got another problem on the same, and hope I would get some answer from you. I have following dataset: data[,1] [1] 4.96 4.95 4.96 4.96 4.97 4.97 4.97 4.97 4.97 4.98 4.98 4.98 4.98 4.98 4.99 4.99 5.00 5.01 [19] 5.01 5.00 5.01 5.01 5.01 5.01 5.02 5.01 5.02 5.02 5.03 5.03 5.03

Hi All , Could you please tell using which library ,Dickey-Fuller Test can be run? Thanks a lot __________________________________________________ [[alternative HTML version deleted]]

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

Hi all, I was running the adf test in R. CODE 1: adf.test(data$LOSS) Augmented Dickey-Fuller Test data: data$LOSS Dickey-Fuller = -1.9864, Lag order = 2, p-value = 0.5775 alternative hypothesis: stationary CODE 2: adf.test(diff(diff(data$LOSS))) Augmented Dickey-Fuller Test data: diff(diff(data$LOSS)) Dickey-Fuller = -6.9287, Lag order = 2, p-value = 0.01 alternative

Hello everyone! Quite a similar question has been posed here some time ago, but there was no explicit solution offered. So I hope that it is OK to pose it again. I want to perform a KPSS test using the packages "urca" or "tseries". But I neither want to use the predefined lag structures, "short" and "long", nor specify the number of lags arbitrarily by

R has a Dickey-Fuller Test implementation (adf.test) that tests for unit roots in an autoregressive process with a constant and linear trend. Is there a DF implementation that doesn't use the constant or trend? Thanks, Jake. -- View this message in context: http://www.nabble.com/Dickey-Fuller-Tests-with-no-constant-and-no-trend-tp23565210p23565210.html Sent from the R help mailing list