similar to: arima, xreg, and the armax model

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

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 want to model different timeseries with ARMAX models in R because I think that ARMAX models will map best to these data. Besides I don't want to use the order of the AR or MA part but the lag e.g. AR Part =ar1, ar2, ar7; MA Part =ma1, ma3 and I want to use exogenous variables as well. I coudn't find any solutions in the R help and therefore I want to ask all of you. Does anyone

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

When I run predict.Arima in my code, I get warnings like: Warning message: In cbind(intercept = rep(1, n), xreg) : number of rows of result is not a multiple of vector length (arg 1) I think this is because I'm not running predict.Arima in the same environment that I did the fit, so the data object used in the fit is no longer present. Looking at the predict.Arima source,

Pro, I am PhD student , i am usuing Rational Expectations Model, I want to model different timeseries with ARMAX models in R because I think that ARMAX models will map best to these data. I want know the steps to use R for ARMAX models I coudn't find any solutions in the R help and therefore I want to ask all of you. Does anyone know how to solve this problem??? That would be

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,

Hello all, I use arima to fit the model with fit <- arima(y, order = c(1,0,1), xreg = list.indep, include.mean = TRUE) and would like to use predict() to forecast: chn.forecast <- rep(0,times=num.record) chn.forecast[1] <- y[1] for (j in 2:num.record){ indep <- c(aa=chn.forecast[j-1], list.indep[j,2:num.indep]) # this is the newxreg in the

in ?arima (R-1.9.1), the return value component 'convergence' should be 'code'? (it's a pity there is no reliable way to check return value documentation consistency with the code, or is there?) h. ---------------------------------- Hiroyuki Kawakatsu School of Management and Economics 25 University Square Queen's University, Belfast Belfast BT7 1NN Northern Ireland

in ?arima (R-1.9.1), the return value component 'convergence' should be 'code'? (it's a pity there is no reliable way to check return value documentation consistency with the code, or is there?) h. ---------------------------------- Hiroyuki Kawakatsu School of Management and Economics 25 University Square Queen's University, Belfast Belfast BT7 1NN Northern Ireland

I am using arima to develop a time series regression model, I am using arima b/c I have autocorrelated errors. Several of my independent variables are categorical and I have coded them as factors . When I run ARIMA I don't get any warning or error message, but I do not seem to get estimates for all the levels of the factor. Can/how does ARIMA handle factors in xreg? here is some example

Hi, I am trying to learn time series, and I am attending a colleague's course on Econometrics. However, he uses e-views, and I use R. I am trying to reproduce his examples in R, but I am having problems specifying a AR(1) model. Would anyone help me with my code? Thanks in advance! Reproducible code follows: download.file("https://sites.google.com/a/proxima.adm.br/main/ex_32.csv

Hello, Does anyone know how the parameter estimates are calculated for xreg variables when called as part of an arima() command, or know of any literature that provides this info? In particular, I was wondering if there is a quick way to compare different combinations of "xreg" variables in the arima() fit in the same way that you would in multiple regression (using AIC & R^2

hi: you should probably send below to R-Sig-Finance because there are some econometrics people over there who could also possibly give you a good answer and may not see this email ? Also, there's package called mar ( I think that's the name ) that may do what you want ? Finally, I don't know how to do it but I think there are ways of converting a multivariate arima into the

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

Hi I'm trying to measure the effect of a policy intervention (Box and Tiao, 1975). This query has to do with the coding of the model rather than with the particulars of my dataset, so I'm not providing the actual dataset (or a simulated one) in this case, apart from some general description. The time series are of length n=34 (annual observations between 1977 and 2010). The policy

Hello all, I have been using R's time series capabilities to perform analysis for quite some time now and I am having some questions regarding its reliability. In several cases I have had substantial disagreement between R and other packages (such as gretl and the commercial EViews package). I have just encountered another problem and thought I'd post it to the list. In this case,

Hi All, I have a very simple question about a time series object: how to access values for a particular year and quarter (say)? Suppose, following http://www.stat.pitt.edu/stoffer/tsa2/R_time_series_quick_fix.htm I have read in data as a time series; here is how it looks. * Qtr1 Qtr2 Qtr3 Qtr4 1960 0.71 0.63 0.85 0.44 1961 0.61 0.69 0.92 0.55 . . . . .

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

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