similar to: questions on ARMA and KPSS

Displaying 20 results from an estimated 5000 matches similar to: "questions on ARMA and KPSS"

2005 Mar 08
2
The null hypothesis in kpss test (kpss.test())
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
2005 Mar 09
1
about kpss.test()
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
2006 Jul 06
2
KPSS 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
2005 May 02
1
Trying to understand kpss.test() in tseries package
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 <-
2009 Mar 06
0
modifying a built in function from the stats package (fixing arima) (CONCLUSIONS)
Thanks a lot to everybody that helped me out with this. Conclusions: (1) In order to edit arima in R: >fix(arima) or alternatively: >arima<-edit(arima) (2) This is not contained in the "Introduction to R" manual. (3) A "productive" fix of arima is attached (arma coefficients printed out and error catched so that it doesn't halt parent loops to search for
2007 Dec 08
2
time series tests
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)
2007 Oct 22
1
Newbie help: Data in an arma fit
I'd like to fit an ARMA(1,1) model to some data (Federal Reserve Bank interest rates) that looks like: ... 30JUN2006, 5.05 03JUL2006, 5.25 04JUL2006, N &lt;---- here! 05JUL2006, 5.25 ... One problem is that holidays have that "N" for their data. As a test, I tried fitting ARMA(1,1) with and without the holidays deleted. In other words, I fit the above data
2007 Nov 02
0
Significance-Problems by using arma/xreg.
Hello. I've got a problem with arma/xreg. I would like to get a better model-fit by implenting some external explanatory variable, so I thought I can implement it by expand the arima-function with an xreg-argument: I have two stationary data vectors y and x of length 201: y <-
2011 Jun 04
0
[R-SIG-Finance] Measure quality of fit for MA(q), ARMA(p, q) and GARCH(p, q)
Thank you so much all for your invaluable inputs. On Sat, Jun 4, 2011 at 3:36 AM, Patrick Burns <patrick at burns-stat.com> wrote: > A common thing to do is the Ljung-Box > test on the residuals. ?For garch it > would be the residuals squared. > > Actually for garch it should be the > rank of the squared residuals -- see >
2009 Oct 13
0
How to specify an ARMA(1, [1,4]) model? Solved
On Tue, Oct 13, 2009 at 5:06 PM, Rolf Turner <r.turner@auckland.ac.nz>wrote: > > Not clear to me what the OP really wants. Perhaps the seasonal > model is what's required; perhaps an arima(1,0,4) model with > theta_2 and theta_3 constrained to be 0. The latter can be > achieved with > > arima(x,order=c(1,0,4),fixed=c(NA,NA,0,0,NA,NA)) > > Or perhaps
2004 Jul 04
1
Re: Seasonal ARMA model
> It might clarify your thinking to note that a seasonal ARIMA model > is just an ``ordinary'' ARIMA model with some coefficients > constrained to be 0 in an efficient way. E.g. a seasonal AR(1) s = > 4 model is the same as an ordinary (nonseasonal) AR(4) model with > coefficients theta_1, theta_2, and theta_3 constrained to be 0. You > can get the same answer as from
2009 Oct 13
1
How to specify an ARMA(1, [1,4]) model?
Hi, I'm trying to model an ARMA(1,[1,4]), i.e. I want only lags 1 and 4 of the Moving Average part. It's the '[1,4]' part that is giving me a problem. I've tried different arma's and arima's in different packages, namely: packages tseries, fArma, FinTS, timeSeries, TSA, Zelig, ds1, forecast For example, with package FinTS: > ( ARIMA(y, order=c(1,0,c(1,4))) )
2004 Feb 12
0
How to predict ARMA models?
Hi all, I am fitting an ARMA(1,(1,4)) model. y(t) = a*y(t-1) + e(t) + b1*e(t-1) + b4*e(t-4) > arma1.14 <- arma(series, lag=list(ar=1, ma=c(1,4)), + include.intercept = F, qr.tol = 1e-07) works fine: Coefficient(s): ar1 ma1 ma4 0.872 -0.445 0.331 I want to forecast 50 periods. I could not find a 'predict' function for ARMA models. I
2004 Oct 13
4
incomplete function output
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)
2008 Jan 10
1
question regarding kpss tests from urca, uroot and tseries packages
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,
2007 Nov 26
3
Time Series Issues, Stationarity ..
Hello, I am very new to R and Time Series. I need some help including R codes about the following issues. I' ll really appreciate any number of answers... # I have a time series data composed of 24 values: myinput = c(n1,n2...,n24); # In order to make a forecasting a, I use the following codes result1 = arima(ts(myinput),order = c(p,d,q),seasonal = list(order=c(P,D,Q))) result2 =
2008 Dec 08
0
ARMA models
Dear ALL: Could you please eamil me how to simulate Mixed Seasonal ARMA (p,q)x(P,Q)12 models [say ARMA(0,1)x(1,0)12 ]from R. With many thanks. Abou ========================== AbouEl-Makarim Aboueissa, Ph.D. Assistant Professor of Statistics Department of Mathematics & Statistics University of Southern Maine 96 Falmouth Street P.O. Box 9300 Portland, ME 04104-9300 Tel: (207)
2008 May 31
0
KPSS test - Lag selection
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
2007 Mar 07
1
good procedure to estimate ARMA(p, q)?
Hi all, I have some residuals from regression, and i suspect they have correlations in them... I am willing to cast the correlation into a ARMA(p, q) framework, what's the best way to identify the most suitable p, and q, and fit ARMA(p, q) model and then correct for the correlations in regression? I know there are functions in R, I have used them before, but I just want to see if I can do
2009 Jul 15
2
storing lm() results and other objects in a list
to clean up some code I would like to make a list of arbitrary length to store?various objects for use in a loop sample code: ############ BEGIN SAMPLE ############## # You can see the need for a loop already linearModel1=lm(modelSource ~ .,mcReg) linearModel2=step(linearModel1) linearModel3=lm(modelSource ~ .-1,mcReg) linearModel4=step(linearModel3) #custom linearModel5=lm(modelSource ~ .