Tell me more ×
Cross Validated is a question and answer site for statisticians, data analysts, data miners and data visualization experts. It's 100% free, no registration required.

I am using an ARIMA model to create a model for correlated errors from my regression model. I am using the auto.arima function from the forecast package in R. I am able to get more data at some frequent interval after the regression model is created, so I get more values for the correlated errors.

My question is how do I update the ARIMA model with a gap in time interval between readings.

share|improve this question
When you say update do you mean reestimate the parameters or revise the form of the model? – Michael Chernick Aug 11 '12 at 18:38
Currently I am revising the model and thus the parameters too, as I just call the auto.arima function everytime I get more new information. I do not know what is the best approach either to re-estimate the parameters or revise the model. Also my initial training data set size is around n=30, so it is quite low. I could increase the n value a bit more if required – Mithileash Mohan Aug 11 '12 at 19:02

1 Answer

Suppose x contains the first series of observations and y contains the second series of observations, where there is a gap between them consisting of m time periods. Let z be the complete time series, given by

z <- ts(c(x,rep(NA,m),y), frequency=??, start=??)

I assume you will set the start and frequency to suitable values.

Then you fit the first model using

fit <- auto.arima(x)

After you have obtained the additional observations, you can update the model in several ways. First, you might just compute the residuals on the new observations without revising the parameters or model order:

fit1 <- Arima(z,model=fit)

Alternatively, you might not revise the model order, but you do update the parameter estimates. In that case:

fit2 <- Arima(z,order=fit$arma[c(1,6,2)])

(Here I have assumed it is a non-seasonal ARIMA model. For seasonality, you would need to add a seasonal argument.)

Lastly, you might revise both the model order and the parameter estimates. Then

fit3 <- auto.arima(z)
share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.