similar to: Problems with ARIMA models?

Hi, I have been using this website ( http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm ) to help me to fit ARIMA models to my data. At the moment I have two possible methods to use. Method 1 If I use arima(ts.data, order=c(1,2,0), xreg=1:length(ts.data)) then the wrong value for the intercept/mean is given (checked on SPSS and Minitab) and

Dear all, When I run an arima(1,1,1) on an I(1) variable, y, I get different estimates to when I first difference the variable myself, e.g y2<-diff(y), and then run arima(1,0,1) on y2. Shouldn't these two approaches give the same output? Any help will be much appreciated. george

Hi, My query is related to ARIMA function in stats package. While looking for the time series literature I found following link which highlights discrepancy in "arima" function while dealing with differenced time series. Is there a substitute function similar to "sarima" mentioned in the following website implemened in R? Any pointers would be of great help.

I am trying to understand how to fit an ARMAX model with the arima function from the stats package. I tried the simple data below, where the time series (vector x) is generated by filtering a step function (vector u, the exogenous signal) through a lowpass filter with AR coefficient equal to 0.8. The input gain is 0.3 and there is a 0.01 normal white noise added to the output: x <- u

Hi, I have a suggestion for the fonction arima and arima0. I think you should not call the constant an intercept because it creates confusion. It is not really an intercept but a mean. For an AR(1) the intercept mu should be defined as: X(t)=mu + phi X(t-1) + e(t) What you call intercept mu is rather defined as (X(t)-mu) = phi (X(t-1)-mu)) + e(t) which is not a common way to define an

Hello all, I''m having fun again with the arima function. This time I read in: http://www.stat.pitt.edu/stoffer/tsa2/R_time_series_quick_fix.htm <<It has recently been suggested (by a reliable source) that using xreg in arima() does NOT fit an ARMAX model [insert slap head icon here]. This will be investigated as soon as time permits.>> (by R.H. Shumway & D.S. Stoffer)

Hello,Could somebody tell me what is the difference between  theese 3 calls of functionsarma(x,order=c(1,0)), arima(x,order=c(1,0,0)) ar(x,order=1)?I expected same residuals of theese three models,but unexpectably for the first two R requiredinitial value of something (what?)...Thanks in advance!  [[alternative HTML version deleted]]

Hello everybody! I have a problem when I try to perform a forecast of an ARIMA model produced by an auto.arima function. Here is what I'm doing: c<-auto.arima(fil[[1]],start.p=0,start.q=0,start.P=0,start.Q=0,stepwise=TRUE,stationary=FALSE,trace=TRUE) # fil[[1]] is time series of monthly data ARIMA(0,0,0)(0,1,0)[12] with drift : 1725.272 ARIMA(0,0,0)(0,1,0)[12] with drift

Hi, I am trying to understand exactly what xreg does in arima. The documentation for xreg says:"xreg Optionally, a vector or matrix of external regressors, which must have the same number of rows as x." What does this mean with regard to the action of xreg in arima? Apparently somehow xreg made the following two arima fit equivalent in R: arima(x, order=c(1,1,1), xreg=1:length(x)) is

I have written before, but to no avail. I have found two minor problems with fitting time series models with R. The thing is, they may be solved with MINOR adjustments to the code. I have posted these problems with detailed examples here: http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm Briefly, the problems are (1) When fitting time series models when there is an AR term present, the

Dear R-help-archive.. I am trying to figure out how to make arima prediction when I have a process involving multivariate time series input, and one output time series (output is to be predicted) .. (thus strictly speaking its an ARMAX process). I know that the arima function of R was not designed to handle multivariate analysis (there is dse but it doesnt handle arma multivariate analysis, only

Hi, I want to give a summary or anova for "arima" model in R, as "summary", and "anova" for "lm". As including various intervention factors in arima(xreg = ) part, I want to assess the significancy of thse factors. I can do it using interrupted analysis of time series by linear regression, but want to see whether arima model works for the data first.

David Stoffer describes some challenges with R's output when fitting ARIMA models for different orders (see Issue 2 at http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm). R doesn't fit an intercept in the model if there is any differencing. David describes a workaround using the xreg parameter to force R to calculate an intercept. Assume I have a variable y and 3 explanatory variables a,

Hi, I have been using this site ( http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm) to help me with some ARIMA modelling in R. Unfortunately the methods mentioned do not appear to work with second order differencing; arima(*, 2, *). I have used some dummy data to illustrate my point. When I use the xreg=... method, the estimate of intercept is *way* off. This can be seen by the high s.e but I

Hi all, I have two time series that I would like to correlate but as they are autocorrelated, I am "pre-whitening" them first by fitting ARIMA models, then correlating their residuals....as described in https://onlinecourses.science.psu.edu/stat510/?q=node/75 However, http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm discusses some issues with ARIMA in R. In particular, for issue 2, if

Dear list, Is there any functionality in R that would allow me to parse config files? I have trie ??config and apropos('config') without succes, and also search the R package site. Mvh. Marie [[alternative HTML version deleted]]

Hi all, I would like to know whether there is any function in R were i can find the cross-correlation of two or more multivariate (time series) data. I tried the function ccf() but it seems like to have two univariate datasets. Please let me know. sincerely, sandeep -- Sandeep Joseph PhD Post Doctoral Associate Center for Tropical & Emerging Global Diseases Paul D. Coverdell Center,

Dear list, I have a vector of elements which I want to combined each with each, but none with itself. For example, > v <- c("a", "b", "c") and I need a function 'combine' such that > combine(v) [[1]] [1] "a" "b" [[2]] [1] "a" "b" [[3]] [1] "b" "c" I am not very interested in the orders

Hi, Following convention below: y(t) = Ax(t)+Bu(t)+eps(t) # observation eq x(t) = Cx(t-1)+Du(t)+eta(t) # state eq I modified the following routine (which I copied from: http://www.stat.pitt.edu/stoffer/tsa2/Rcode/Kall.R) to accommodate u(t), an exogenous input to the system. for (i in 2:N){ xp[[i]]=C%*%xf[[i-1]] Pp[[i]]=C%*%Pf[[i-1]]%*%t(C)+Q siginv=A[[i]]%*%Pp[[i]]%*%t(A[[i]])+R

Dear useRs, I have a collection of source file and some of these call others. The files are distribute among a number of directories, and to know how to call some other file they need to know what file is currently executed. As example, I have a file 'search.R' located in directory 'bin' which needs to access the file 'terms.conf' located in the directory 'conf',