similar to: Significance-Problems by using arma/xreg.

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

Hi all, I want to fit the following model to my data: Y_t= a+bY_(t-1)+cY_(t-2) + Z_t +Z_(t-1) + Z_(t-2) + X_t + M_t i.e. it is an ARMA(2,2) with some additional regressors X and M. [Z_t's are the white noise variables] So, I run the following code: for (i in 1:rep) { index=sample(4,15,replace=T) final<-do.call(rbind,lapply(index,function(i)

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,

Hello everyone, Hope you all are doing great! I have been fitting arima models and performing forecasts pretty straightforwardly in R. However, I wanted to add a couple of regressors to the arima model to see if it could improve the accuracy of the forecasts but have had a hard time trying to do so. I used the following R function: arima(x, order = c(0, 0, 0), seasonal = list(order = c(0, 0,

Hi everyone, I'm trying to include an external regressor in my ARIMA model but am having some problems with the data format in R. I've named my dependent variable of interest "count" and the external regressor "abc". The external regressor is a binary variable. Here are the contents of "abc" and the model I've attemped (along with its error message):

Hi, Using arima on data that are in a data frame, especially when adding xreg, would be much easier if the arima function contained 1) a "data=" argument 2) the possibility to include the covariate(s) in a formula style. Ideally the call could be something like > arima(symptome, order=c(1,0,0), xreg=~trait01*mesure0, data=anxiete) ( or arima(symptome~trait01*mesure0,

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

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

Dear all, I am fitting an arimax (arima with some extra explanatory variables) model to a time series. Say, I have a Y (dependent variable) and an X (explanatory). Y is 100 observations (time series) and X is 100 + 20 (20 to use for the forecast horizon). I can not make xreg work with the forecast function for an arima fit. The "predict" function seems to be working but the

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

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

Hello, The "arma" function in the "tseries" package allows estimation of models with specific "ar" and "ma" lags with its "lag" argument. For example: y[t] = a[0] + a[1]y[t-3] +b[1]e[t-2] + e[t] can be estimated with the following specification : arma(y, lag=list(ar=3,ma=2)). Is this possible with the "arima" function in the

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

When I do ARMA(2,2) using one lag of LCPIH data This is eview result > > *Dependent Variable: DLCPIH > **Method: Least Squares > **Date: 08/12/11 Time: 12:44 > **Sample (adjusted): 1970Q2 2010Q2 > **Included observations: 161 after adjustments > **Convergence achieved after 14 iterations > **MA Backcast: 1969Q4 1970Q1 > ** > **Variable Coefficient Std.

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

Hello,guys: Recently, I am working on a seasonal ARIMA model. And I met some problem in the forecasting. Now I just want to know that How does R perform the predict procedure(the predict formula, the initial setting of errors,etc.)? I run the following commands and get the original code of the "predict" command, but I can't read it. Can anybody explain it to me? Thanks! saji from

I did a regression with ARMA errors using arima0 with ari<-arima0(y,order=c(2,0,2),xreg=reg1,delta=-1) or ari<-arima0(y,order=c(2,0,2),xreg=reg1) where reg1 is the matrix of the regressors and when I see diag(ari$var.coef) I get negative terms. Do you know what this mean ? I try to change transform.pars to 0 or 1 but this crash R on Windows. Is it possible to test the significativity

Hello, Does anybody know if the function arima with an external regressor (xreg) applies the auto correlation on the dependant variable or on the residuals. In comparison with SAS (proc autoreg), it seems that the auto correlation applies on the residuals but i'd like to have the confirmation. I want to estimate: Y[t] = a[1]*X[t] + a[2] + E[t] with E[t]=b[1]*E[t-1] Should I use : arima(Y,

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,